Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,230 T→A 5.8% intergenic (‑120/+193) sapA ← / ← ymjA antimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,2300TA 5.8% 198.1 / 2.8 78intergenic (‑120/+193)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base T (40/33);  new base A (3/2);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

CCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACG  >  NC_000913/1357083‑1357377
                                                                                                                                                   |                                                                                                                                                   
ccAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTc                                                                                                                                                    >  2:941399/1‑149 (MQ=255)
ccAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTc                                                                                                                                                    >  2:946056/1‑149 (MQ=255)
ccAAAAGAGACGATAATACCTGGCGCATAATATTAACGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTc                                                                                                                                                    >  1:275896/1‑149 (MQ=255)
 cAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                                                                                   <  2:451913/149‑1 (MQ=255)
      gagaCGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTATTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACTTCGCCATc                                                                                                                                              >  1:1059438/1‑149 (MQ=255)
       agaCGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTGAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCt                                                                                                                                             >  1:690155/1‑149 (MQ=255)
           gATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTAGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTc                                                                                                                                         >  1:950821/1‑149 (MQ=255)
              aataCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAc                                                                                                                                      >  1:990928/1‑149 (MQ=255)
              aataCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTAGCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAc                                                                                                                                      <  2:1165039/149‑1 (MQ=255)
               ataCCTGGCGCATACTATTAACGGATTTTATGTATAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                                                                     <  1:259311/149‑2 (MQ=255)
                  ccTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTc                                                                                                                                  <  2:923957/149‑1 (MQ=255)
                   cTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                                                                                 >  2:108474/1‑149 (MQ=255)
                      gcgcATAATATTAAGGGATTTTATGTAAAGTGGCTATATTACTAATATTTAATGACAATTGACATTACCGTTTTTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTAGCCATCTGTTCAACATTCGTAc                                                                                                                              <  2:782960/149‑1 (MQ=255)
                      gcgcATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTATTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAc                                                                                                                              <  1:881110/149‑1 (MQ=255)
                         cATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTGCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGACGCCATCTGTTCAACATTCGTACCCg                                                                                                                           >  2:759185/1‑149 (MQ=255)
                           taataTTAAGGGATTTTATCTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTATGTTCAGGGGTCGTAATGAGGCCCGTGTGGTTGGTCTTTGAACAGGTGACGTCGCCATCTGTTAAACATTCTTACACGtt                                                                                                                         >  2:441868/1‑148 (MQ=255)
                            aataTTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGa                                                                                                                        <  1:590101/149‑1 (MQ=255)
                             atatTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAt                                                                                                                       >  1:200727/1‑149 (MQ=255)
                              tatTAAGGGATTTTATTTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGCAATGAGGGCCGTGTGGTTTGTCTTTGAACAGGTGACGTCGACATCTGTTAAACATTCGTACCCGTTAtt                                                                                                                      >  2:848514/1‑149 (MQ=255)
                                 tAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACAGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtct                                                                                                                   >  1:393415/1‑149 (MQ=255)
                                   tcgtaTTTTATGTAAAGAGTCGATCTTTCTTATATTTAATCACATTTGCCATTACCGTTTGTGTTCCGGGTGCGTATTGAGGGCGGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                                                                                 <  2:852359/145‑1 (MQ=255)
                                       aTTTTATGTAAAGAGGCTATCTTACTAATCTTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataa                                                                                                             <  2:198728/149‑1 (MQ=255)
                                       aTTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataa                                                                                                             >  2:536577/1‑149 (MQ=255)
                                        ttttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCTCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                                            >  2:550073/1‑149 (MQ=255)
                                        ttttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                                            >  2:594222/1‑149 (MQ=255)
                                        ttttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataaa                                                                                                            >  2:269057/1‑148 (MQ=255)
                                         tttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCAGTTCAACATTCGTACACGTGATTTCTctataact                                                                                                           >  2:1164994/1‑149 (MQ=255)
                                          ttATGTACAGAGGCTATCATACTAAAATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGGCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                                                                          <  1:207151/149‑1 (MQ=255)
                                                 acgagGCTATCTTACTAATATTTAATGCCATTTGCCATTACCGATTGTGGTCCGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACCGGGGACGTCGCCATCTGTTCAACATTCGTACCCGGGATTTCTCTATAACTATAACTca                                                                                                   <  2:198957/147‑1 (MQ=255)
                                                 aaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTAGCGTTTGTGTTCAGGTGTCGTAATGAGGGCCGTGTTGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTAAAACTca                                                                                                   <  2:458943/149‑1 (MQ=255)
                                                 aaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTca                                                                                                   >  1:979745/1‑149 (MQ=255)
                                                  agatGCTATATTACTAATATTTAATGACATTTGCCATTACGGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACTTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcac                                                                                                  <  1:686422/149‑1 (MQ=255)
                                                   gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcaca                                                                                                 >  2:979448/1‑149 (MQ=255)
                                                      ggtATCTTACTAATATTTAATGAAATTTGCCATTAACGTTTGTGTTGAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                              <  1:441868/147‑1 (MQ=255)
                                                      gcgaTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                              <  1:361316/146‑1 (MQ=255)
                                                      gCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                              >  2:729113/1‑149 (MQ=255)
                                                         aTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaa                                                                                           >  1:533941/1‑149 (MQ=255)
                                                           cTTACTAATATTTAATTACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACt                                                                                         >  2:1123525/1‑149 (MQ=255)
                                                              aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                                                                       >  1:508501/1‑148 (MQ=255)
                                                              aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                                                                       <  2:508501/148‑1 (MQ=255)
                                                              aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                                                                      >  1:572308/1‑149 (MQ=255)
                                                              aCTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                                                                      >  2:120093/1‑149 (MQ=255)
                                                                 aataatTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGGCGTCATGAGGGCCGTGTGGTTGGTCGTTCAACACGTGTCGTCGCCATCTGTTCAACAGTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCga                                                                                   <  2:903555/149‑1 (MQ=255)
                                                                 aaTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCga                                                                                   >  2:467369/1‑149 (MQ=255)
                                                                         aTGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCaga                                                                           >  2:312119/1‑149 (MQ=255)
                                                                           gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                                                                         >  2:406926/1‑149 (MQ=255)
                                                                            aCATTTGCCATTACCGTTTTTGTGCAGGGGTCTTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGGGACGTCGCCTTCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                                                                        <  1:1154335/149‑1 (MQ=255)
                                                                            aCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                                                                        >  1:180043/1‑149 (MQ=255)
                                                                            aCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                                                                        <  1:754205/149‑1 (MQ=255)
                                                                               tttGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTATTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTa                                                                     >  2:601059/1‑149 (MQ=255)
                                                                                      ttACCGTTTGTGTTCAGGGGTCGTAATGTGGGCGGTGTGGTTGGTCGTTGAACTGGTGACGTCGCCATGTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGATCAACTTAGCGAGGAGCAGAGCGGTACGATCGt                                                              <  2:954279/149‑1 (MQ=255)
                                                                                      ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt                                                              <  2:327590/149‑1 (MQ=255)
                                                                                      ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTCAACAGGTGACGTCGCCATCTGTTCAACATACGTACCCGTGATTTCTCTATAACTATATCTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt                                                              >  2:921344/1‑149 (MQ=255)
                                                                                           gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                                 >  1:1103302/1‑125 (MQ=255)
                                                                                           gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                                 <  2:1103302/125‑1 (MQ=255)
                                                                                            tttGTGTTGAGGGGTCTTAATGATGGGGGTGTGTTATGGAGTTGAACAGTTGAAGTCTCCATCTGTTCAACATTCGTTCCCGTGATTTCTCTCTAACTAAAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaaa                                                        <  1:820946/149‑1 (MQ=255)
                                                                                                  ttCAGGGGTCGTAATGAGGGACGTGTTGTTGGTCGTTGAACAGGTGACGTCGACATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGAGGTACGATCGTGAGAAAAACAtt                                                  >  2:1084781/1‑149 (MQ=255)
                                                                                                     tggggTCTTTATGATGGCCGTGTGGGTGGGCGTTGTACAGGTGTCGTCGCCATCTGTTCAACATTCGTACCGGTGATTTCTCTATTTCTATAACTCAGAGAACAACTTTGCGAGTTGGTGAGCGGTTCGATCGTGAGAAACGGATTaaa                                               <  2:1102608/148‑1 (MQ=255)
                                                                                                     aGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaaa                                               <  1:941399/149‑1 (MQ=255)
                                                                                                     aGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaaa                                               <  1:946056/149‑1 (MQ=255)
                                                                                                        ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCa                                            >  1:436659/1‑149 (MQ=255)
                                                                                                          tCGTAATGAGCGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACg                                          <  2:275896/149‑1 (MQ=255)
                                                                                                                tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAGCCCATTAAAGCACGCTGCAg                                    >  1:34702/1‑149 (MQ=255)
                                                                                                                tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAg                                    <  1:729113/149‑1 (MQ=255)
                                                                                                                    ggCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCAATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCg                                >  1:93643/1‑149 (MQ=255)
                                                                                                                           tggttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAGTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCAtt                         >  2:994224/1‑149 (MQ=255)
                                                                                                                             gttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGt                       <  1:536577/149‑1 (MQ=255)
                                                                                                                                 gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGc                   <  2:341889/149‑1 (MQ=255)
                                                                                                                                  tCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGcc                  >  2:687251/1‑149 (MQ=255)
                                                                                                                                  tCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGcc                  <  1:108474/149‑1 (MQ=255)
                                                                                                                                  gcGTTTAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTACAGCACGCTGCAGTGCGTGTCATTGTTAGcc                  <  1:1114520/148‑1 (MQ=255)
                                                                                                                                    gTTGAACAGATGGCGACGTCACCTGTTCAACAGTGGTCCCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAg                >  1:695675/1‑149 (MQ=255)
                                                                                                                                     ttGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATACCTCCCAGAACAACTTAGCGAGGAGCAGAGCGGGACGCTCGTGAGAAACCCATTAACGCACGCTGCAGTGCGTGTCATTGTTAGCCAGa               <  2:871893/149‑1 (MQ=255)
                                                                                                                                         aCAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATgcg           >  1:467508/1‑149 (MQ=255)
                                                                                                                                          cAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATgcgg          >  1:141132/1‑149 (MQ=255)
                                                                                                                                             gTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGt       <  2:572308/149‑1 (MQ=255)
                                                                                                                                             gTGACGTCGCCATCTGGTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACCGAACTACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGt       <  2:1168557/149‑1 (MQ=255)
                                                                                                                                                  gTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACg  >  2:1030214/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACG  >  NC_000913/1357083‑1357377

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: