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.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                       |                                                                                                                                                     
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                       <  2:5410923/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:1687869/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:5308621/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:4217213/150‑1 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:6870721/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:7019471/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                               >  2:3990884/1‑150 (MQ=255)
        tGCAAAAAGAGATGAAAGGGCGTTATCACAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgcg                                                                                                                                               >  2:6522288/1‑150 (MQ=255)
          cAAAAAGAGATGAAAGGGAGTTATCCCAAACATGTCTGGACGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                             >  1:6026056/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACg                                                                                                                                           >  2:7460188/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                           >  1:6993587/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                          >  1:7354329/1‑150 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgc                                                                                                                                        >  1:7467687/1‑148 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                                        >  1:3441430/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcgccc                                                                                                                                    >  1:2979721/1‑150 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                                    >  1:6112502/1‑150 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccc                                                                                                                                   >  1:570281/1‑143 (MQ=255)
                      aatgggCTTTATCCAAAACATGTCGGGCCGGACGCCCCGCAAAATACTCCCCCGCCGCGACGGACTAAAATTTTTTCTTATGTTTGTTTTTCACCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                 <  2:2178275/147‑1 (MQ=255)
                        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtggcccacacc                                                                                                                               >  2:352982/1‑150 (MQ=255)
                        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGggacccgcacc                                                                                                                               >  1:6848742/1‑150 (MQ=255)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                              <  2:5311254/150‑1 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTggacccgcccccc                                                                                                                             >  2:6061707/1‑149 (MQ=38)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:5566231/150‑1 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCAAAcc                                                                                                                         >  1:4537206/1‑150 (MQ=38)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgggcccccaccccccaccc                                                                                                                       >  2:3661375/1‑143 (MQ=38)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAacca                                                                                                                       >  2:4692926/1‑149 (MQ=255)
                                 acccAAAAATTGTCGGGCCGGACGCCCCGCAAAAAACGGCCCCGCCGCGCCGGACTAAAATTTATTCTTAGGGTGGTTTTTCCCCCTCCCCCGGGGGGAGGGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                       <  1:6030172/149‑1 (MQ=255)
                                   ccAAACTTGTCGGGCCGGACGACCCGTCAAATCCTGCACCGACGCGACGGACGAAAATGTTTTGGTATGGTTGGTTTTCTCCCTCTCCCGGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:6090762/150‑1 (MQ=37)
                                   caaaaCTTGTCGGGCCGTACGACCCGGAAAATACTCCACCGACGCGACGTACTAAAATTTTTTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:5909500/148‑1 (MQ=25)
                                   caaaaCATGTGGGGCGGGCCGACCGGTCAATTACTGCCCCGACGCGCCGGACTAAAATGTTTTGTTATGGTTGGTTTACCCCCTCCCCGTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:4310688/148‑1 (MQ=255)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGTCCCCCACCCCAAaccacaacc                                                                                                                  >  2:245186/1‑144 (MQ=37)
                                     aaaCATGTCTGGCCGGACGACCCGCCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGACGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                  <  2:4585300/150‑1 (MQ=255)
                                      aaccctGTCGGGCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAATATTATCGGAATGGTGGGTTTGCTCCCTCCCCCGGGGGTTGGGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                  <  2:5041571/145‑1 (MQ=255)
                                           tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACACACACCCAAACCCACTCCCACAAg                                                                                                             >  1:6300272/1‑150 (MQ=38)
                                           ggtggggCCGGACGACCGGGCACATACGGCCCCGCCGCGCCGGCCGAAAAAGTATTGGTATGGTGGGTTTTCTCCCTCGCCCGGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:293317/145‑1 (MQ=255)
                                            gTCTGGCCGGACGACCCGGCAAAAACTGCACCGACGCGACGTACCAAAAATTATTCTTAAGGTTGGTTTTCCCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                            <  2:1789872/150‑1 (MQ=38)
                                             tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGCACCCGCACCCAAACCCACACCCGCCAggg                                                                                                           >  2:8171921/1‑150 (MQ=37)
                                              cTGGCCGGCCGCCCCGCCAATTACTCCACCGACGCGCCGGCCTAAAAAGTTTTCGTATGGTTGGTTTACTCCCTCACCCGGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                          <  1:6873427/150‑1 (MQ=255)
                                               tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCACCGCCCCCCAACCCACGCCCACAAGGgag                                                                                                         >  2:4702660/1‑150 (MQ=37)
                                                ggCCGGACGACCCGGCAAATACGGCCCCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTGCGCCCGCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                        <  2:4637577/150‑1 (MQ=38)
                                                ggCCGCCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTACGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                        <  1:6103928/150‑1 (MQ=38)
                                                 gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCAAACCCTCAAGGGagag                                                                                                       >  2:7287283/1‑150 (MQ=255)
                                                 gCCGGACGACCCGGAAACTCCTGCACCGACGCTACGCACGAAAATTTATTCGTATGGTTGTTTTACCCCCTCCCCCTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                       <  1:4873490/150‑1 (MQ=255)
                                                 cccGGACGACCCGGCCCCTACTGCACCCACGCGCCGGACGACAACGTATGCGTATGGTTGGTTTTCCCCCCCCCCCGGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  2:795209/149‑1 (MQ=255)
                                                   cGGACGCCCGGTAAAATACTGCACCGCCGCGACGTACTAAAATTTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                     <  1:2853805/150‑1 (MQ=255)
                                                   cGGACGACCCGGCAAATACTGCACCGAGGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCACCCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  2:1931074/150‑1 (MQ=255)
                                                   cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCGCACCCACACCCACACCCACAACGGGGAggg                                                                                                     >  1:5093514/1‑150 (MQ=37)
                                                    ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGAACCACACCCAAACCCACACCCCCCACGGGGAgggg                                                                                                    >  1:1365179/1‑150 (MQ=37)
                                                        gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCGCTCCCACCAGGGAGAGGGGCCCg                                                                                                >  2:4881150/1‑150 (MQ=38)
                                                          cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCACCCAAACCCCCACCCCCCAGGGAGGGGGGGCCgcg                                                                                              >  1:6356193/1‑148 (MQ=25)
                                                           ccGGCAAATACTGCCCCGACGCGACGGACTAAAAATTATTCGTATGGTTGGTTTCCCCCCTCCCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                             <  2:8148659/150‑1 (MQ=255)
                                                           ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGGTCGCGCCCCCCCCAACCCACGCCCGCAAGGGCggggggggggcgc                                                                                             >  1:2367247/1‑143 (MQ=25)
                                                            ccgCAAATTACTGCCCCGAGGCTCCGGACTAAATTTTTTTGTTATGGTTGGTTTTCCCCCTCCCCGTGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                             <  1:6570261/148‑1 (MQ=255)
                                                            cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCTAACCCACACCCACCAGGGGGAGGGGGCCCATCg                                                                                            >  1:6259963/1‑150 (MQ=37)
                                                              gaaaaTCCTGCACTGCCGCGACGTACTAAAACTTATTGTTCTGTTTGGTTTTCCCCCCCTCCTTGTGGGTTAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  1:7019471/148‑1 (MQ=255)
                                                              gaaaaCACTGCACCGCCGCGACGGACGAAATAGTATTCTTCCGGTTGTTTTCCCCCCACACCCTGTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  1:1327397/148‑1 (MQ=255)
                                                               cAAATCCTGCCCCGACGCGCCGGCCTAAAAAGTATTCGTAAGGTTGGTTTTCGCCCTCCCCCTGGGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                         <  2:6810620/150‑1 (MQ=255)
                                                                aaaTACTGCACCGACGCGACGGACGAAAATGTTTTGGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                        <  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)
                                                                                                                                                       |                                                                                                                                                     
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  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: