Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,771,192 G→C 5.3% A63G (GCG→GGG)  yibH ← inner membrane protein YibH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,771,1920GC 5.3% 501.7 / 10.1 148A63G (GCG→GGG) yibHinner membrane protein YibH
Reads supporting (aligned to +/- strand):  ref base G (75/65);  new base C (2/6);  total (77/71)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

CGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCAC  >  NC_000913/3771041‑3771341
                                                                                                                                                        |                                                                                                                                                     
cGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAGCACCT‑‑CCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGc                                                                                                                                                      <  1:462982/150‑1 (MQ=255)
  ccGGCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCg                                                                                                                                                       <  2:941696/150‑1 (MQ=255)
   cGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCNCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCNTCACCTGTGGCGTGATAGGGATCGc                                                                                                                                                      >  2:2051369/1‑150 (MQ=255)
   cGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGc                                                                                                                                                      <  1:234545/150‑1 (MQ=255)
    gTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCt                                                                                                                                                     <  2:1931946/150‑1 (MQ=255)
      catcaGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGAATCGCTAt                                                                                                                                                   >  2:2157772/1‑150 (MQ=255)
       atcaGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATc                                                                                                                                                  <  1:2472098/150‑1 (MQ=255)
           ggTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTg                                                                                                                                              >  1:2033343/1‑150 (MQ=255)
           ggTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGCGATCCCTATCAcgc                                                                                                                                              <  1:1688335/150‑3 (MQ=255)
            gTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGCGATCCCTATCAcgcc                                                                                                                                             <  2:1294427/150‑4 (MQ=255)
            gTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGCGATCCCTATCAcgcc                                                                                                                                             <  2:1433784/150‑4 (MQ=255)
             tCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGTACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCAGTGATAGCGATCCCTATCAcgcc                                                                                                                                             <  2:2146066/150‑4 (MQ=255)
             tCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGcc                                                                                                                                            <  1:1227754/150‑1 (MQ=255)
                gCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCttt                                                                                                                                         >  2:1054642/1‑150 (MQ=255)
                gCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCttt                                                                                                                                         <  2:298186/150‑1 (MQ=255)
                   tGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTgc                                                                                                                                      >  2:425442/1‑150 (MQ=255)
                    gAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTgcg                                                                                                                                     <  2:1729984/150‑1 (MQ=255)
                      aGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGcgcg                                                                                                                                   <  1:79167/150‑1 (MQ=255)
                         cTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTa                                                                                                                                <  2:1222637/150‑1 (MQ=255)
                          tGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTaa                                                                                                                               >  1:815491/1‑150 (MQ=255)
                           gTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTaaa                                                                                                                              >  1:1505140/1‑150 (MQ=255)
                           gTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTaaa                                                                                                                              >  1:1305063/1‑150 (MQ=255)
                            tCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTaaaa                                                                                                                             >  1:376255/1‑150 (MQ=255)
                              aaCTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGTGCGGTAAAAgt                                                                                                                           >  1:2445085/1‑150 (MQ=255)
                               aCTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAgtg                                                                                                                          >  2:1525868/1‑150 (MQ=255)
                                  cgcgCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTaa                                                                                                                       >  2:106889/1‑150 (MQ=255)
                                       cTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGggtg                                                                                                                  <  2:675623/150‑1 (MQ=255)
                                       cTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGggtg                                                                                                                  >  1:1122114/1‑150 (MQ=255)
                                         ggTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCTGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGggtggt                                                                                                                <  1:1301522/150‑1 (MQ=255)
                                             aCGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTa                                                                                                            <  2:77839/150‑1 (MQ=255)
                                              cGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAg                                                                                                           <  2:1935540/150‑1 (MQ=255)
                                              cGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAg                                                                                                           <  2:44749/150‑1 (MQ=255)
                                               gAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGt                                                                                                          >  2:1350426/1‑150 (MQ=255)
                                                 aCCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTtc                                                                                                        <  2:2074183/150‑1 (MQ=255)
                                                   cGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTtcat                                                                                                      >  1:2268951/1‑150 (MQ=255)
                                                    gggTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTtcatc                                                                                                     >  1:1923281/1‑150 (MQ=255)
                                                      gTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAggg                                                                                                                    <  1:1664750/133‑1 (MQ=255)
                                                      gTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAggg                                                                                                                    >  2:1664750/1‑133 (MQ=255)
                                                          aGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACaaa                                                                                               >  1:2206626/1‑150 (MQ=255)
                                                           gCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACaaaa                                                                                              <  2:546622/150‑1 (MQ=255)
                                                            cTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAAt                                                                                             <  1:1613360/150‑1 (MQ=255)
                                                             ttAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATc                                                                                            >  1:2323769/1‑150 (MQ=255)
                                                              tAAAAGGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCCCCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGGGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCa                                                                                           >  1:680239/1‑150 (MQ=255)
                                                              tAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCa                                                                                           >  2:5029/1‑150 (MQ=255)
                                                              tAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCa                                                                                           >  1:198901/1‑150 (MQ=255)
                                                              tAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCa                                                                                           >  2:2389640/1‑150 (MQ=255)
                                                               aaaaaCCTCCTCGCCCTTTTTATTAAGCTTATTATTCTTGTCAGTTATTTTTTTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCCCTATGAGTGCCTTTTGCGAGGTAAATGGGTAAGGGTTCTTGTAGCTCATCATCAAAATCac                                                                                          <  2:1016798/150‑2 (MQ=255)
                                                                 aaaGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAAc                                                                                        <  1:1417488/150‑1 (MQ=255)
                                                                 aaaGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAAc                                                                                        <  2:1292660/150‑1 (MQ=255)
                                                                  aaGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAAcc                                                                                       <  1:916962/150‑1 (MQ=255)
                                                                    gCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCAc                                                                                     <  2:1784648/150‑1 (MQ=255)
                                                                        cTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCAc                                                                                 >  1:1454498/1‑150 (MQ=255)
                                                                         tCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCAcc                                                                                <  2:491521/150‑1 (MQ=255)
                                                                           gCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAg                                                                              <  1:1350426/150‑1 (MQ=255)
                                                                             ccTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGaa                                                                            <  2:1490858/150‑1 (MQ=255)
                                                                                 ttGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAacac                                                                        <  1:1525868/150‑1 (MQ=255)
                                                                                  tGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACg                                                                       <  1:309344/150‑1 (MQ=255)
                                                                                  tGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACg                                                                       >  2:2348785/1‑150 (MQ=255)
                                                                                   gAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGc                                                                      <  1:1040522/150‑1 (MQ=255)
                                                                                      tAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTc                                                                   <  1:2400665/150‑1 (MQ=255)
                                                                                        aGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTccc                                                                 >  1:1362858/1‑150 (MQ=255)
                                                                                        aGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTccc                                                                 <  2:1998740/150‑1 (MQ=255)
                                                                                            gATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAgcg                                                             <  2:1896175/150‑1 (MQ=255)
                                                                                            gATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAgcg                                                             >  1:2305743/1‑150 (MQ=255)
                                                                                            gATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGAGTGATAGGGATAGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAgcg                                                             >  1:1785810/1‑150 (MQ=255)
                                                                                                     tGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCg                                                    <  2:2318039/150‑1 (MQ=255)
                                                                                                         aGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAg                                                >  2:822915/1‑150 (MQ=255)
                                                                                                           tGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCg                                              >  2:2311863/1‑150 (MQ=255)
                                                                                                              cTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTcc                                           >  2:660641/1‑150 (MQ=255)
                                                                                                              cTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTcc                                           >  1:2109114/1‑150 (MQ=255)
                                                                                                               ttCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCa                                          >  1:495940/1‑150 (MQ=255)
                                                                                                                  aGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGCAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTg                                       >  2:764227/1‑150 (MQ=255)
                                                                                                                   gTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGa                                      <  1:1140162/150‑1 (MQ=255)
                                                                                                                     aaCAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGAtt                                    <  1:562439/150‑1 (MQ=255)
                                                                                                                      aCAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGAttt                                   >  2:1689262/1‑150 (MQ=255)
                                                                                                                        aaTTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTAc                                 >  1:125232/1‑150 (MQ=255)
                                                                                                                           tCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCgg                              <  2:1957706/150‑1 (MQ=255)
                                                                                                                             ccGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAaca                                                                         >  2:588628/1‑105 (MQ=255)
                                                                                                                             ccGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAaca                                                                         <  1:588628/105‑1 (MQ=255)
                                                                                                                                tCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGc                         <  1:1054642/150‑1 (MQ=255)
                                                                                                                                 cACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCg                        <  2:872000/150‑1 (MQ=255)
                                                                                                                                  aCCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCCCCAGAAACACGCCTCCCAGCGCCGCCGTCGCCATCGTCAACTCATTTGCCGGAATGCgg                       >  2:2141442/1‑150 (MQ=255)
                                                                                                                                  aCCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCgg                       >  1:87136/1‑150 (MQ=255)
                                                                                                                                  aCCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCgg                       <  2:1122114/150‑1 (MQ=255)
                                                                                                                                   ccTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGa                      >  1:2185776/1‑150 (MQ=255)
                                                                                                                                      gtgGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGa                   >  2:2171938/1‑150 (MQ=255)
                                                                                                                                       aggcagTGATAGCGATCCCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGa                   >  1:892697/6‑150 (MQ=255)
                                                                                                                                        ggCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGAtt                 >  2:724876/1‑150 (MQ=255)
                                                                                                                                           ggTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGAtttt               <  1:224266/149‑1 (MQ=255)
                                                                                                                                             tGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTaa             <  1:660641/150‑1 (MQ=255)
                                                                                                                                             tGATAGGGATCGCAATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTaa             <  2:2100859/150‑1 (MQ=255)
                                                                                                                                              gATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTaaa            >  1:1864767/1‑150 (MQ=255)
                                                                                                                                               aTAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAg           <  1:1999221/150‑1 (MQ=255)
                                                                                                                                               aTAGCGATCCCTATCACTCCCTATTGCGCGCTAAAAGTGAAAGCGTGGTTCTACTTCAACAACAAAAACAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAg           >  2:2439598/1‑150 (MQ=255)
                                                                                                                                                  gggATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGAcc        >  1:920602/1‑150 (MQ=255)
                                                                                                                                                    gATCGCTATCACTGCCTTTTGCGCCGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGc      <  1:476180/150‑1 (MQ=255)
                                                                                                                                                     aTCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGcc     <  2:1725703/150‑1 (MQ=255)
                                                                                                                                                      tCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACATGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGccc    >  1:725289/1‑150 (MQ=255)
                                                                                                                                                      tCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGccc    <  2:74631/150‑1 (MQ=255)
                                                                                                                                                        gCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCAc  >  2:2056886/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
CGCCGTCATCAGGTCAGCCTGAAGTCTGTCAACTCGCGCCTGGTAACGAACCGGGTCGAGCTTAAAAAGCACCTCGCCCTTTTGAATAAGCTGATTATTCTTGTCAGTGACTTCAGTAACAATTCCCGTCACCTGTGGCGTGATAGGGATCGCTATCACTGCCTTTTGCGCGGTAAAAGTGTAAGGGTGGTTGTAGTTCATCAACAAAATCAAACCACTCACCAGAAACACGCCTCCCAGCGCCGCCGTCGCCAGCGTCCACTGATTTACCGGAATGCGGAAGATTTTAAAGACCGCCCAC  >  NC_000913/3771041‑3771341

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