Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,569 C→G 21.9% intergenic (+35/‑161) 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,5690CG21.9% 638.3 / 63.0 268intergenic (+35/‑161)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (110/94);  new base G (4/58);  total (115/153)
Fisher's exact test for biased strand distribution p-value = 1.69e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

GAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGC  >  NC_000913/164420‑164718
                                                                                                                                                      |                                                                                                                                                      
gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                        >  1:147290/1‑150 (MQ=255)
 aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                                       >  2:3057059/1‑150 (MQ=255)
    ggAGCGTACCGTGAGGTGCAAAACGAGATGAACGGGCGTTATCCCACCCAGGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                                                                                                                                                    <  2:9633433/150‑1 (MQ=255)
     gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                                   >  1:9451287/1‑150 (MQ=255)
     gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                                   >  2:10744086/1‑150 (MQ=255)
       gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                                 <  2:412684/150‑1 (MQ=255)
        cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAACCATGTCTGGCCGGACGATCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                                <  2:454241/150‑1 (MQ=255)
        cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                                >  2:8003841/1‑150 (MQ=255)
         gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                               >  1:10040554/1‑150 (MQ=255)
          tACCGTGCGGTGCAAATCGAGATGAACGGGCGTTATCCCAACCAGGTCTGGCCGGACGCCCCGGCAAATACTGCACCGACGCGCCGGACGAAAACGTTTTCGTACGGTTGGTT‑TTCGCCCGCTCCCGGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                              <  2:966417/150‑1 (MQ=255)
            ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                            <  2:51811/150‑1 (MQ=255)
             cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGcct                                                                                                                                           >  1:5489630/1‑148 (MQ=255)
             cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGATT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                           >  2:3240074/1‑150 (MQ=255)
                  ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                      >  1:1018042/1‑150 (MQ=255)
                   gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                     >  1:7764284/1‑150 (MQ=255)
                    tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                    >  2:9434793/1‑150 (MQ=255)
                     gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgcgcgc                                                                                                                                   >  1:5380533/1‑150 (MQ=255)
                       aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTAGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                 <  2:2397362/150‑1 (MQ=255)
                         aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                               >  1:2964139/1‑150 (MQ=255)
                          aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGGGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                              >  2:743108/1‑149 (MQ=255)
                              tattaatgggCTTTAACCAAAACCTGTCGGGCGGTACGACCCGGAAATTACTGCCCCGACGCGAGTTAATAAAATTTTTTCTTATGGTTGTTT‑TTCCCCCTCCCCCTGTGGTTTGGGGCGGGGTTGCGGGCATCAGCGCGCACGTtcacc                                                                                                                          <  2:3142937/143‑1 (MQ=255)
                                 gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                     <  2:2599767/120‑1 (MQ=255)
                                 gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                     >  1:2599767/1‑120 (MQ=255)
                                   atgggCTTTAACCCAAACTGGTCGGGCCGGACGACCCGGAAAATACGGCACCGACGCTACGGACGAAAATTTATTCGTATGGTTGGTT‑TTCTCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                     <  2:709183/148‑1 (MQ=25)
                                   aaGGGCGTTATCCCAACCTGGTCTGGCCGGACGACCCGTCAAATCCTGCACCGACGCGACGGACGAAAATTTATTCGTATGGTTGGTT‑TTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                     <  2:617059/150‑1 (MQ=37)
                                    tgggCGTTAGCCCAAACATGGCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAATGTTTTCTTATGGTTGGTT‑TTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                    <  2:4574145/149‑1 (MQ=37)
                                          ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGACCCACACCCCAAcc                                                                                                              >  1:6369974/1‑150 (MQ=38)
                                          ttATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACGGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              <  2:5380533/150‑1 (MQ=38)
                                          ttATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACGGCCCCGAGGCTAGTGACTAAAATTTTTTCTTATGTTTGTTT‑TTCTCCCTACCCCGGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              <  2:293729/150‑1 (MQ=255)
                                          ttATCCAAACCTTCTCGGGCCGCACGACCCGGCAAATACTGCACCGCCGCGCCGGACTAAAAAGTTTTCGTATGGTTGGTT‑TTCCCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              <  2:5179103/150‑1 (MQ=25)
                                          taaccccAAACAGGTGGGGCCGGCCGACCCGGCAAATACTGCACCGCGGCGACGGACTAAAAAGTTTTGGTATGGTTGGTT‑TTCGCCCTCCCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              <  2:5616464/146‑1 (MQ=255)
                                           atattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              >  2:1353786‑M1/98‑150 (MQ=255)
                                              cccAACCTTGTGGGGCCGGACGCCCCGGCATTTCCTGCCCCGCCGCGCGGGACTAAAATTTTTTCGTATGTTTGGTT‑TTCTCCCTCTCCGGGGGGGTGTGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                          <  2:6678617/150‑1 (MQ=255)
                                              cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACCCCCAAAccaaca                                                                                                          >  1:9552340/1‑146 (MQ=38)
                                               ccAATCTTGGCTGGCCGGACCACCGGGAAATTCCGGCCCCTCCGCGCCGGACGAAAATGTATTCTTTGGGTTGGTT‑TTCGCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                         <  1:9849500/150‑1 (MQ=255)
                                               ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAAcccacac                                                                                                         >  1:8093122/1‑146 (MQ=255)
                                               ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCCAAccccccc                                                                                                         >  1:2089773/1‑146 (MQ=38)
                                                cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCACCCCCAAAcccacaac                                                                                                        >  1:6112602/1‑145 (MQ=38)
                                                     ttGTCTGGCGGCACGACCCGCAAAATACTCCACCGACGCGACGCACTAAAATTTATTCGTATGTTTGTTT‑TTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                   <  2:5030470/149‑1 (MQ=25)
                                                     aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                   <  2:2684807/150‑1 (MQ=255)
                                                       gTCTGTCCGGCCGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                 <  2:7296340/150‑1 (MQ=38)
                                                       gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAACCTCTCCCTCAAgg                                                                                                 >  2:1744447/1‑150 (MQ=255)
                                                          tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                                <  2:10289832/100‑1 (MQ=255)
                                                          tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCACAAGGgag                                                                                              >  1:9825867/1‑150 (MQ=38)
                                                           ggCCGGACCACCCGGCAAATACTGCACCTACGCGACGGACTAAAATGTTTTCGTATGGTTGTTT‑TTCTCCAACGCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                             <  2:573528/150‑1 (MQ=255)
                                                            gCCGGACGCCCCGCAAAATCCGCCCCCTACGCTCCGGACTAAATATTATGCGTTTGGTTGTTT‑TTCTCCCGCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                            <  1:2387064/150‑1 (MQ=255)
                                                             ccGGACTCCCCGCAAATTACTCCCCCGCCGCGCCGTACGAAAATTTA‑TAGTTATGTTGGTTGTCCCCCCTCCCCCTGGGGGTGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                           <  1:6808108/150‑1 (MQ=255)
                                                              cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCACCCACACCCACACCCCCCAGGGAGAggg                                                                                          >  1:4366184/1‑150 (MQ=37)
                                                                gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCACTCCCACAAGGGGGAggggg                                                                                        >  2:3662965/1‑149 (MQ=255)
                                                                gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCACAACCACAAGGGAggggggg                                                                                        >  1:1841835/1‑149 (MQ=37)
                                                                 acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCGCACCCCAACCCACACCCACAAGGGGGAgggggg                                                                                       >  1:10818156/1‑148 (MQ=37)
                                                                   gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCACCCAAAACCACAACCCCAAGCGGGAGGGGCCCg                                                                                     >  1:8476966/1‑150 (MQ=25)
                                                                   gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCCCCGAAACCCCAACCCCAAGGGGGTGGGGGCCg                                                                                     >  1:1865563/1‑150 (MQ=25)
                                                                      ccGGCAAATACGGCACCGACGCGACGGACTAAAATGTATTCGTAAGGTTGGTT‑TTCTCCCTCGCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                  <  2:8660648/150‑1 (MQ=38)
                                                                         gCAAATCCGGCCCCGACGCGCCGGACGAAAAAGTATTCGTACGGTTGGTT‑TTCTCCCCCCCCCGGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                               <  2:5811129/150‑1 (MQ=255)
                                                                         gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGCACCCACGCCCGCCCGGGAGAGGGGGCCGATCGAg                                                                               >  2:9814813/1‑150 (MQ=37)
                                                                          cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                              <  2:10042544/150‑1 (MQ=255)
                                                                                 tGCCCCGACGCGACGGACGAAAAATTTTTCTTATGGTTGGTT‑TTCTCCCTCCCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAAtt                                                                      <  2:1539184/150‑1 (MQ=38)
                                                                                  gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCACACCCACACCCGCTCCCGCCCGGGAGAGGGGGCCGAGCGCGC‑CCAAAatt                                                                     >  2:10851448/1‑150 (MQ=25)
                                                                                   cacccccgcgacgGACGAAAATTTATCCTCATGGTTTTTT‑AACCCCCTAGCCCTGTGTTGGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAAtttt                                                                    <  1:4045932/144‑1 (MQ=255)
                                                                                     ccgccgCGCGGGACTAAAATTTTTTCGTATGGTGGGTT‑TTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTga                                                                  <  1:7252915/150‑1 (MQ=255)
                                                                                     ccgccgCGACGGACTAAAATGTTTTCGTATGTTTGTTT‑TTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTga                                                                  <  2:8476966/150‑1 (MQ=38)
                                                                                     cccacgcgacgCACCAAAAAGTATTCGAAAGGTTGGTT‑TACACCCACCCCCTGTGGCAGAGGGCCGGGGTGAGGCAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTga                                                                  <  2:8796251/147‑1 (MQ=255)
                                                                                     ccGACGCGACGTACGAAAAATTTTTCTTATGTTTGTTT‑TTCTCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTga                                                                  <  1:2683272/150‑1 (MQ=38)
                                                                                       gcggcgacgTACTAAAAAGTTTTCTTATGGTTGTTT‑TCCGCCCTCCCCCGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTgaga                                                                <  1:6852933/147‑1 (MQ=255)
                                                                                        acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACAACCACAAGGGAGAGTGGACCCAACGAGC‑ACAAAATTTtagag                                                               >  2:2789576/1‑150 (MQ=25)
                                                                                            gagggaccAAAATTTATTCTTACGTTTGGTT‑TTCCCCCACCCCGGGTGGGTGGGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGatat                                                           <  2:1972041/147‑1 (MQ=255)
                                                                                                 acgAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAAAACACACCCACAAGGGGGCGGGGCCCGATCGAGC‑CCAAATTTTTAGATATAtcatat                                                      >  1:9941741/1‑146 (MQ=21)
                                                                                                  cgAAAAAGTATTCTTACGGTTGGTT‑TCCCCCCTCTCCCTGTGGCAGAGGCCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATAtcttcta                                                     <  1:5734634/150‑2 (MQ=255)
                                                                                                    aaaaaTTATTCTAATGGTTGGTT‑ATCACCCTCCCCCTGTGGGAGTGGCCCGGGGTGAGCCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATAtcttcttct                                                   <  1:5063565/150‑1 (MQ=255)
                                                                                                     aaaaGTAATCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCACCACCCCCAAAACCACTCCCACAAGGGCGACGGGCCCCATCCAGCAACAAATTTT‑CGAGATAtcttcttcgg                                                  >  2:6043175/1‑148 (MQ=16)
                                                                                                       aaGTATTCGTATGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTc                                                <  1:7958573/150‑1 (MQ=255)
                                                                                                         gTATTCGTAAGGTTGTTT‑TTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCtt                                              <  2:3618661/150‑1 (MQ=255)
                                                                                                          tATTCTTATGGTTGTTT‑TTCTCCCTCTCCCGGGGGGAGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAAGTTTGAGAGAGATCTTCTTCTGTCTTg                                             <  2:4908414/150‑1 (MQ=255)
                                                                                                           aTTCGTATGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                            <  1:7361234/150‑1 (MQ=255)
                                                                                                            ttCGTATGGTTGGTT‑TTCTCCCTCTCCCTGTGGTAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                           <  1:5805344/150‑1 (MQ=255)
                                                                                                             tcttatggttggtt‑TTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                          <  2:3501642/144‑1 (MQ=39)
                                                                                                             tCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCAACCGAACCCACTCCCACAAGGGAGAGGGGACCGATCGAGC‑ACAAAATTTGAGTGATATCTTATTCTGTGTTGTaa                                          >  1:1001509/1‑150 (MQ=37)
                                                                                                               gtttggttgttt‑TTCTCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                        <  2:683384/146‑1 (MQ=39)
                                                                                                               gTATGGTTGGTT‑TTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                        <  1:10740866/150‑1 (MQ=255)
                                                                                                                  cggttggtt‑TACTCCCTCCCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGTTATCTTCTTCTGTCTTGTAACagaa                                     <  1:921224/149‑1 (MQ=255)
                                                                                                                    gttggtc‑TTCCCCCCCCCCCTGGGGGTGGGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                   <  1:3201849/150‑1 (MQ=255)
                                                                                                                         tt‑ttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                              <  2:5898068/150‑1 (MQ=255)
                                                                                                                          ttgtcgccctctccctGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                              <  1:1213530/144‑1 (MQ=255)
                                                                                                                                ccccccccctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGAGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                        <  1:1031213/144‑1 (MQ=255)
                                                                                                                                 cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                       <  1:3667211/150‑1 (MQ=255)
                                                                                                                                  ctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                      <  2:9049874/150‑1 (MQ=255)
                                                                                                                                     tccctGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCGCACCCAAACCCCCGCCCACCAGGGCGGGGGGGCCCCGCGCGC‑CCAAAGTGTTGGTGGTATAAAGATCTGGGGTGTAACAGAAGAACACGAAATCgggggg                   >  2:7910296/1‑147 (MQ=255)
                                                                                                                                        ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAACCCACACCCTCAAGGGAGAGGGGACCCAACGAGC‑ACAAATTTTGAGAGATATATTCTTCTGTCGTTTAACAGAAGAAAAGAAAATCGGGCGTTTg                >  2:9101082/1‑150 (MQ=255)
                                                                                                                                         tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGGGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc               >  2:7051134/1‑150 (MQ=255)
                                                                                                                                           ggggtgtgggCGGGGTTGTGGGCATCAGCGCGCACGTTCACACTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGCTATCTTCTTCTGTCTTGTAACAGAAGAACAGAATATCGGGCTTTTgcgc             <  2:10099539/143‑1 (MQ=255)
                                                                                                                                           ggggcgcgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc             <  1:4665968/143‑1 (MQ=255)
                                                                                                                                           gGGGAGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc             >  2:2304225‑M1/2‑150 (MQ=255)
                                                                                                                                             gtcgtGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt           >  1:10563457‑M1/6‑150 (MQ=255)
                                                                                                                                             ggtgagGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt           <  1:6459175/147‑1 (MQ=255)
                                                                                                                                              GAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg          >  1:10392241‑M1/1‑150 (MQ=255)
                                                                                                                                                GAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCCAACCCACTCCCACAAGGGGGAGGGGACCCATCGAGC‑ACAAAATTTGAGTGATAAATTATTATGTCTTGTAACAGAAGAACAGAAAATAGGGCTTTTGCGCCGGaa        <  2:8203458‑M1/150‑1 (MQ=255)
                                                                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAAACCACACCCGCAATGGGGAGGGGGCCGATCGAGC‑CCAAAATTTGAGATATATATTATTATGTCTTTTAAAAGAAGAAAAGAAAATACGGCTTTTTCGCCAGAaaat     >  2:6401717/1‑150 (MQ=255)
                                                                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTACCACCAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTTAGAGTTATAATATTAAGTAATTTAACAAAAGACACCAAAAACGGGCTTTTGAGCCGGAatat     >  2:8221209/1‑150 (MQ=255)
                                                                                                                                                    GCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                         <  2:2441762‑M1/50‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTT‑TTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGC‑ACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGC  >  NC_000913/164420‑164718

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: ‑