Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 25.6% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC25.6% 227.4 / 51.2 193intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (60/81);  new base C (7/45);  total (67/126)
Fisher's exact test for biased strand distribution p-value = 1.43e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                       |                                                                                                                                                     
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                       <  2:5410923/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:1687869/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:5308621/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:4217213/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:6870721/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:7019471/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                               >  2:3990884/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCACAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgcg                                                                                                                                               >  2:6522288/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGAGTTATCCCAAACA‑TGTCTGGACGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                             >  1:6026056/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACg                                                                                                                                           >  2:7460188/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                           >  1:6993587/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                          >  1:7354329/1‑150 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgc                                                                                                                                        >  1:7467687/1‑148 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                                        >  1:3441430/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcgccc                                                                                                                                    >  1:2979721/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                                    >  1:6112502/1‑150 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccc                                                                                                                                   >  1:570281/1‑143 (MQ=255)
                      aatgggCTTTATCCAAAACA‑TGTCGGGCCGGACGCCCCGCAAAA‑TACTCCCCCGCCGCGACGGACTAAAATTTTTTCTTATGTTTGTTTTTCACCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                 <  2:2178275/147‑1 (MQ=255)
                        aGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtggcccacacc                                                                                                                               >  2:352982/1‑150 (MQ=255)
                        aGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGggacccgcacc                                                                                                                               >  1:6848742/1‑150 (MQ=255)
                         gggCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                              <  2:5311254/150‑1 (MQ=255)
                          ggCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTggacccgcccccc                                                                                                                             >  2:6061707/1‑149 (MQ=38)
                              ttATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:5566231/150‑1 (MQ=255)
                              ttATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCAAAcc                                                                                                                         >  1:4537206/1‑150 (MQ=38)
                                aTCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgggcccccaccccccaccc                                                                                                                       >  2:3661375/1‑143 (MQ=38)
                                aTCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAacca                                                                                                                       >  2:4692926/1‑149 (MQ=255)
                                 acccAAAAATTGTCGGGCCGGACGCCCCGCAAAA‑AACGGCCCCGCCGCGCCGGACTAAAATTTATTCTTAGGGTGGTTTTTCCCCCTCCCCCGGGGGGAGGGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                       <  1:6030172/149‑1 (MQ=255)
                                   ccAAACT‑TGTCGGGCCGGACGACCCGTCAAA‑TCCTGCACCGACGCGACGGACGAAAATGTTTTGGTATGGTTGGTTTTCTCCCTCTCCCGGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:6090762/150‑1 (MQ=37)
                                   caaaaCT‑TGTCGGGCCGTACGACCCGGAAAA‑TACTCCACCGACGCGACGTACTAAAATTTTTTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:5909500/148‑1 (MQ=25)
                                   caaaaCA‑TGTGGGGCGGGCCGACCGGTCAAT‑TACTGCCCCGACGCGCCGGACTAAAATGTTTTGTTATGGTTGGTTTACCCCCTCCCCGTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:4310688/148‑1 (MQ=255)
                                     aaaCA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGTCCCCCACCCCAAaccacaacc                                                                                                                  >  2:245186/1‑144 (MQ=37)
                                     aaaCA‑TGTCTGGCCGGACGACCCGCCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGACGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                  <  2:4585300/150‑1 (MQ=255)
                                      aaccctGTCGGGCCGGACGCCCCGGCAAA‑TACTGCACCGACGCGACGGACGAAATATTATCGGAATGGTGGGTTTGCTCCCTCCCCCGGGGGTTGGGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                  <  2:5041571/145‑1 (MQ=255)
                                           tGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACACACACCCAAACCCACTCCCACAAg                                                                                                             >  1:6300272/1‑150 (MQ=38)
                                           ggtggggCCGGACGACCGGGCACA‑TACGGCCCCGCCGCGCCGGCCGAAAAAGTATTGGTATGGTGGGTTTTCTCCCTCGCCCGGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:293317/145‑1 (MQ=255)
                                            gTCTGGCCGGACGACCCGGCAAA‑AACTGCACCGACGCGACGTACCAAAAATTATTCTTAAGGTTGGTTTTCCCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                            <  2:1789872/150‑1 (MQ=38)
                                             tCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGCACCCGCACCCAAACCCACACCCGCCAggg                                                                                                           >  2:8171921/1‑150 (MQ=37)
                                              cTGGCCGGCCGCCCCGCCAAT‑TACTCCACCGACGCGCCGGCCTAAAAAGTTTTCGTATGGTTGGTTTACTCCCTCACCCGGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                          <  1:6873427/150‑1 (MQ=255)
                                               tGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCACCGCCCCCCAACCCACGCCCACAAGGgag                                                                                                         >  2:4702660/1‑150 (MQ=37)
                                                ggCCGGACGACCCGGCAAA‑TACGGCCCCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTGCGCCCGCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                        <  2:4637577/150‑1 (MQ=38)
                                                ggCCGCCCGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATGCGTACGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                        <  1:6103928/150‑1 (MQ=38)
                                                 gCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCAAACCCTCAAGGGagag                                                                                                       >  2:7287283/1‑150 (MQ=255)
                                                 gCCGGACGACCCGGAAAC‑TCCTGCACCGACGCTACGCACGAAAATTTATTCGTATGGTTGTTTTACCCCCTCCCCCTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                       <  1:4873490/150‑1 (MQ=255)
                                                 cccGGACGACCCGGCCCC‑TACTGCACCCACGCGCCGGACGACAACGTATGCGTATGGTTGGTTTTCCCCCCCCCCCGGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  2:795209/149‑1 (MQ=255)
                                                   cGGACGCCCGGTAAAA‑TACTGCACCGCCGCGACGTACTAAAATTTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                     <  1:2853805/150‑1 (MQ=255)
                                                   cGGACGACCCGGCAAA‑TACTGCACCGAGGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCACCCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  2:1931074/150‑1 (MQ=255)
                                                   cGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCGCACCCACACCCACACCCACAACGGGGAggg                                                                                                     >  1:5093514/1‑150 (MQ=37)
                                                    ggACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGAACCACACCCAAACCCACACCCCCCACGGGGAgggg                                                                                                    >  1:1365179/1‑150 (MQ=37)
                                                        gacCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCGCTCCCACCAGGGAGAGGGGCCCg                                                                                                >  2:4881150/1‑150 (MQ=38)
                                                          cccGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCACCCAAACCCCCACCCCCCAGGGAGGGGGGGCCgcg                                                                                              >  1:6356193/1‑148 (MQ=25)
                                                           ccGGCAAA‑TACTGCCCCGACGCGACGGACTAAAAATTATTCGTATGGTTGGTTTCCCCCCTCCCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                             <  2:8148659/150‑1 (MQ=255)
                                                           ccGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGGTCGCGCCCCCCCCAACCCACGCCCGCAAGGGCggggggggggcgc                                                                                             >  1:2367247/1‑143 (MQ=25)
                                                            ccgCAAATTACTGCCCCGAGGCTCCGGACTAAATTTTTTTGTTATGGTTGGTTTTCCCCCTCCCCGTGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                             <  1:6570261/148‑1 (MQ=255)
                                                            cGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCTAACCCACACCCACCAGGGGGAGGGGGCCCATCg                                                                                            >  1:6259963/1‑150 (MQ=37)
                                                              gaaaa‑TCCTGCACTGCCGCGACGTACTAAAACTTATTGTTCTGTTTGGTTTTCCCCCCCTCCTTGTGGGTTAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  1:7019471/148‑1 (MQ=255)
                                                              gaaaa‑CACTGCACCGCCGCGACGGACGAAATAGTATTCTTCCGGTTGTTTTCCCCCCACACCCTGTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  1:1327397/148‑1 (MQ=255)
                                                               cAAA‑TCCTGCCCCGACGCGCCGGCCTAAAAAGTATTCGTAAGGTTGGTTTTCGCCCTCCCCCTGGGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                         <  2:6810620/150‑1 (MQ=255)
                                                                aaa‑TACTGCACCGACGCGACGGACGAAAATGTTTTGGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                        <  2:4930660/150‑1 (MQ=255)
                                                                      cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCGAACCCGAAACCGCAAGGGTGTGGGGTGCGATAGAGGACAAAt                                                                                   >  1:7499453/1‑150 (MQ=25)
                                                                      cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAAACACAACCACAAGGGAGATGGGACCCATCGAGAACAAAt                                                                                   >  2:7734226/1‑150 (MQ=25)
                                                                       tGCACCGACGCTACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                  <  2:4762457/150‑1 (MQ=255)
                                                                       ggccccgccgCGACGCACTAAAATTTATTCTTTTGGTTGTTTTTCCCCCCCCCCGGGTGGGCGCGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                  <  1:797503/146‑1 (MQ=255)
                                                                          cccGACGCGAGGGACCAAAATTTATTGGTAAGGTTGTTTTTCTCCCGCACCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                               <  1:2605190/149‑1 (MQ=255)
                                                                           ccgccgCGACGTACTAAAATTTATTGGTATGGTTGTTTTTCTCCCCCGCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:8199528/150‑1 (MQ=255)
                                                                           ccgccgCGACGGACGAAAATGAATTCTTACGGTTGGTTTTCGCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:6683029/150‑1 (MQ=255)
                                                                           ccgccgCGACGGACCAAAATTTATTGTTAGGGTGGTTTTTCGCCCTCGCCCTGGGGGCGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:2438516/150‑1 (MQ=255)
                                                                             gacgcgacgGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  2:8320199/150‑1 (MQ=255)
                                                                                     ggacgAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                    <  1:7698925/150‑1 (MQ=255)
                                                                                     ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCACCCAAAACCACACCCCCAACGGGGGGGGGGCCGAACGCGCCCAAAATTTTTGAGATAtatt                                                                    >  1:974418/1‑150 (MQ=18)
                                                                                     ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCGCCCCAAACCACAACCACAAAGGAGAGGGGGCCGCGAGAGCACCAAATTTgtgaggtttaat                                                                    >  1:3780794/1‑143 (MQ=16)
                                                                                     ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTCCACCCGCACCCAAAAACACAACCAAAAAGGGGGGGGGGCCCATCGAGCCCAAAATTTTATAGATAtatt                                                                    >  1:4267577/1‑150 (MQ=16)
                                                                                         gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCCCCCGCCCCCAAACCCACACCCACAAGGGCGCGGGGACCCCTCGCGCCCAAATTTTGAGAGATAtcttcttc                                                                >  2:2084662/1‑150 (MQ=25)
                                                                                           aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCGCACCCACAACGGCGAGGGGGCCCCTACAGCACAAATTATGAGAGAAATCTTATTATg                                                              >  2:4527098/1‑150 (MQ=21)
                                                                                               gTATTCGTATGGTTGGTTTTCGCCCTCGCCCTGTGGGAGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:5025023/150‑1 (MQ=255)
                                                                                                tagccGTAAGGTTGGTTTTCGCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                         <  2:4490259/146‑1 (MQ=255)
                                                                                                tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                         <  1:8369731/150‑1 (MQ=255)
                                                                                                   tCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  2:6305760/150‑1 (MQ=255)
                                                                                                     gTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  2:7210112/150‑1 (MQ=255)
                                                                                                     gTAAGGTTGGTTTTCTCCCTCTCCCCGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  2:4934525/150‑1 (MQ=255)
                                                                                                        tggttgtttTTCTCCCTCGCCGTGTGGGCTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  2:3199411/149‑1 (MQ=255)
                                                                                                        agtttggttTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  1:4441059/147‑1 (MQ=255)
                                                                                                           ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCGCACCCAAAACCACACCCCCCAGGGGGTGGGGCCCGAGCGAGCCCAAAAGTTGTGAGATAAGTTATTGTGTGTTGTTAGAgaagaa                                              >  2:850550/1‑150 (MQ=12)
                                                                                                           ttggttTTCGCCCCCTCCCTGTGGGATAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                              <  1:1878006/150‑1 (MQ=255)
                                                                                                            tggttTCCTCCCCCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  1:5874214/150‑1 (MQ=255)
                                                                                                                    tccctctccctGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  2:5490067/150‑1 (MQ=255)
                                                                                                                     ccctctccctGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:8343241/150‑1 (MQ=255)
                                                                                                                      cctctccctGTGGTAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  2:115883/150‑1 (MQ=255)
                                                                                                                      cctcgccatGTGGTAGTGGGCCGGGGTGTGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  2:7354329/150‑1 (MQ=255)
                                                                                                                                ggggcgtgggCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  1:7287283/143‑1 (MQ=255)
                                                                                                                                  ggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:3990884/150‑1 (MQ=255)
                                                                                                                                    agagGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCAAAACCGCACCCACAAGGGGGGGGGGTCCCGTAGAGCACAAAATTTGTGGGATATAATATTATGTCTTGTAACACAAAAACAAAAAAAAGGGCGGTTGGGCCAGa                     >  2:1372014/1‑150 (MQ=255)
                                                                                                                                        ggCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCCCCCCCAACCCACACCCACAACGGACAGGGGCCCGAACCAGCCCCAAATTTGAGATATAAAATAATATGTCTTGTAAAACAACAAAACAAAAACCGGCGTTTTCGccagtaaat                 >  2:1282281/1‑143 (MQ=255)
                                                                                                                                         gcgggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  2:570281/147‑1 (MQ=255)
                                                                                                                                         gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  1:7915150/150‑1 (MQ=255)
                                                                                                                                          ccGGGGTGCGGGCATCAGCGCGCACAGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  1:1853170/150‑1 (MQ=255)
                                                                                                                                           cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  2:947537/150‑1 (MQ=255)
                                                                                                                                            ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:4537206/150‑1 (MQ=255)
                                                                                                                                              ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:5246802/150‑1 (MQ=255)
                                                                                                                                                    ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:2068452/150‑1 (MQ=255)
                                                                                                                                                      cATCAGCGCGCACGTTCACACTCACCCTATCCCTCTGCCTAAAGGGAGAGGGGACCGATCGAGCACTAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGTAATAGc   <  2:6026056/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAA‑TACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑