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.

GAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGC  >  NC_000913/164420‑164718
                                                                                                                                                      |                                                                                                                                                      
gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                        >  1:147290/1‑150 (MQ=255)
 aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                                       >  2:3057059/1‑150 (MQ=255)
    ggAGCGTACCGTGAGGTGCAAAACGAGATGAACGGGCGTTATCCCACCCAGGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                                                                                                                                                    <  2:9633433/150‑1 (MQ=255)
     gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                                   >  1:9451287/1‑150 (MQ=255)
     gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                                   >  2:10744086/1‑150 (MQ=255)
       gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                                 <  2:412684/150‑1 (MQ=255)
        cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAACCATGTCTGGCCGGACGATCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                                <  2:454241/150‑1 (MQ=255)
        cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                                >  2:8003841/1‑150 (MQ=255)
         gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                               >  1:10040554/1‑150 (MQ=255)
          tACCGTGCGGTGCAAATCGAGATGAACGGGCGTTATCCCAACCAGGTCTGGCCGGACGCCCCGGCAAATACTGCACCGACGCGCCGGACGAAAACGTTTTCGTACGGTTGGTTTTCGCCCGCTCCCGGTGGGAGAGGGCCGGGGTGAggg                                                                                                                                              <  2:966417/150‑1 (MQ=255)
            ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                            <  2:51811/150‑1 (MQ=255)
             cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGcct                                                                                                                                           >  1:5489630/1‑148 (MQ=255)
             cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGATTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                           >  2:3240074/1‑150 (MQ=255)
                  ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                      >  1:1018042/1‑150 (MQ=255)
                   gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                     >  1:7764284/1‑150 (MQ=255)
                    tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                    >  2:9434793/1‑150 (MQ=255)
                     gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgcgcgc                                                                                                                                   >  1:5380533/1‑150 (MQ=255)
                       aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTAGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                 <  2:2397362/150‑1 (MQ=255)
                         aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                               >  1:2964139/1‑150 (MQ=255)
                          aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGGGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                              >  2:743108/1‑149 (MQ=255)
                              tattaatgggCTTTAACCAAAACCTGTCGGGCGGTACGACCCGGAAATTACTGCCCCGACGCGAGTTAATAAAATTTTTTCTTATGGTTGTTTTTCCCCCTCCCCCTGTGGTTTGGGGCGGGGTTGCGGGCATCAGCGCGCACGTtcacc                                                                                                                          <  2:3142937/143‑1 (MQ=255)
                                 gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                     <  2:2599767/120‑1 (MQ=255)
                                 gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                     >  1:2599767/1‑120 (MQ=255)
                                   atgggCTTTAACCCAAACTGGTCGGGCCGGACGACCCGGAAAATACGGCACCGACGCTACGGACGAAAATTTATTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                     <  2:709183/148‑1 (MQ=25)
                                   aaGGGCGTTATCCCAACCTGGTCTGGCCGGACGACCCGTCAAATCCTGCACCGACGCGACGGACGAAAATTTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                     <  2:617059/150‑1 (MQ=37)
                                    tgggCGTTAGCCCAAACATGGCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAATGTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                    <  2:4574145/149‑1 (MQ=37)
                                          ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGACCCACACCCCAAcc                                                                                                              >  1:6369974/1‑150 (MQ=38)
                                          ttATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACGGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              <  2:5380533/150‑1 (MQ=38)
                                          ttATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACGGCCCCGAGGCTAGTGACTAAAATTTTTTCTTATGTTTGTTTTTCTCCCTACCCCGGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              <  2:293729/150‑1 (MQ=255)
                                          ttATCCAAACCTTCTCGGGCCGCACGACCCGGCAAATACTGCACCGCCGCGCCGGACTAAAAAGTTTTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              <  2:5179103/150‑1 (MQ=25)
                                          taaccccAAACAGGTGGGGCCGGCCGACCCGGCAAATACTGCACCGCGGCGACGGACTAAAAAGTTTTGGTATGGTTGGTTTTCGCCCTCCCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              <  2:5616464/146‑1 (MQ=255)
                                           atattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                              >  2:1353786‑M1/98‑150 (MQ=255)
                                              cccAACCTTGTGGGGCCGGACGCCCCGGCATTTCCTGCCCCGCCGCGCGGGACTAAAATTTTTTCGTATGTTTGGTTTTCTCCCTCTCCGGGGGGGTGTGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                          <  2:6678617/150‑1 (MQ=255)
                                              cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACCCCCAAAccaaca                                                                                                          >  1:9552340/1‑146 (MQ=38)
                                               ccAATCTTGGCTGGCCGGACCACCGGGAAATTCCGGCCCCTCCGCGCCGGACGAAAATGTATTCTTTGGGTTGGTTTTCGCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                         <  1:9849500/150‑1 (MQ=255)
                                               ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAAcccacac                                                                                                         >  1:8093122/1‑146 (MQ=255)
                                               ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCCAAccccccc                                                                                                         >  1:2089773/1‑146 (MQ=38)
                                                cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCACCCCCAAAcccacaac                                                                                                        >  1:6112602/1‑145 (MQ=38)
                                                     ttGTCTGGCGGCACGACCCGCAAAATACTCCACCGACGCGACGCACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                   <  2:5030470/149‑1 (MQ=25)
                                                     aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                   <  2:2684807/150‑1 (MQ=255)
                                                       gTCTGTCCGGCCGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                 <  2:7296340/150‑1 (MQ=38)
                                                       gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAACCTCTCCCTCAAgg                                                                                                 >  2:1744447/1‑150 (MQ=255)
                                                          tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                                <  2:10289832/100‑1 (MQ=255)
                                                          tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCACAAGGgag                                                                                              >  1:9825867/1‑150 (MQ=38)
                                                           ggCCGGACCACCCGGCAAATACTGCACCTACGCGACGGACTAAAATGTTTTCGTATGGTTGTTTTTCTCCAACGCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                             <  2:573528/150‑1 (MQ=255)
                                                            gCCGGACGCCCCGCAAAATCCGCCCCCTACGCTCCGGACTAAATATTATGCGTTTGGTTGTTTTTCTCCCGCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                            <  1:2387064/150‑1 (MQ=255)
                                                             ccGGACTCCCCGCAAATTACTCCCCCGCCGCGCCGTACGAAAATTTATAGTTATGTTGGTTGTCCCCCCTCCCCCTGGGGGTGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                           <  1:6808108/150‑1 (MQ=255)
                                                              cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCACCCACACCCACACCCCCCAGGGAGAggg                                                                                          >  1:4366184/1‑150 (MQ=37)
                                                                gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCACTCCCACAAGGGGGAggggg                                                                                        >  2:3662965/1‑149 (MQ=255)
                                                                gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCACAACCACAAGGGAggggggg                                                                                        >  1:1841835/1‑149 (MQ=37)
                                                                 acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCGCACCCCAACCCACACCCACAAGGGGGAgggggg                                                                                       >  1:10818156/1‑148 (MQ=37)
                                                                   gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCGCACCCAAAACCACAACCCCAAGCGGGAGGGGCCCg                                                                                     >  1:8476966/1‑150 (MQ=25)
                                                                   gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCCCCGAAACCCCAACCCCAAGGGGGTGGGGGCCg                                                                                     >  1:1865563/1‑150 (MQ=25)
                                                                      ccGGCAAATACGGCACCGACGCGACGGACTAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                  <  2:8660648/150‑1 (MQ=38)
                                                                         gCAAATCCGGCCCCGACGCGCCGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCCCCCCCGGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                               <  2:5811129/150‑1 (MQ=255)
                                                                         gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGCACCCACGCCCGCCCGGGAGAGGGGGCCGATCGAg                                                                               >  2:9814813/1‑150 (MQ=37)
                                                                          cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                              <  2:10042544/150‑1 (MQ=255)
                                                                                 tGCCCCGACGCGACGGACGAAAAATTTTTCTTATGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                      <  2:1539184/150‑1 (MQ=38)
                                                                                  gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCACACCCACACCCGCTCCCGCCCGGGAGAGGGGGCCGAGCGCGCCCAAAatt                                                                     >  2:10851448/1‑150 (MQ=25)
                                                                                   cacccccgcgacgGACGAAAATTTATCCTCATGGTTTTTTAACCCCCTAGCCCTGTGTTGGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                    <  1:4045932/144‑1 (MQ=255)
                                                                                     ccgccgCGCGGGACTAAAATTTTTTCGTATGGTGGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                  <  1:7252915/150‑1 (MQ=255)
                                                                                     ccgccgCGACGGACTAAAATGTTTTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                  <  2:8476966/150‑1 (MQ=38)
                                                                                     cccacgcgacgCACCAAAAAGTATTCGAAAGGTTGGTTTACACCCACCCCCTGTGGCAGAGGGCCGGGGTGAGGCAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                  <  2:8796251/147‑1 (MQ=255)
                                                                                     ccGACGCGACGTACGAAAAATTTTTCTTATGTTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                  <  1:2683272/150‑1 (MQ=38)
                                                                                       gcggcgacgTACTAAAAAGTTTTCTTATGGTTGTTTTCCGCCCTCCCCCGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                <  1:6852933/147‑1 (MQ=255)
                                                                                        acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACAACCACAAGGGAGAGTGGACCCAACGAGCACAAAATTTtagag                                                               >  2:2789576/1‑150 (MQ=25)
                                                                                            gagggaccAAAATTTATTCTTACGTTTGGTTTTCCCCCACCCCGGGTGGGTGGGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                           <  2:1972041/147‑1 (MQ=255)
                                                                                                 acgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAAAACACACCCACAAGGGGGCGGGGCCCGATCGAGCCCAAATTTTTAGATATAtcatat                                                      >  1:9941741/1‑146 (MQ=21)
                                                                                                  cgAAAAAGTATTCTTACGGTTGGTTTCCCCCCTCTCCCTGTGGCAGAGGCCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcta                                                     <  1:5734634/150‑2 (MQ=255)
                                                                                                    aaaaaTTATTCTAATGGTTGGTTATCACCCTCCCCCTGTGGGAGTGGCCCGGGGTGAGCCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                   <  1:5063565/150‑1 (MQ=255)
                                                                                                     aaaaGTAATCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCACCACCCCCAAAACCACTCCCACAAGGGCGACGGGCCCCATCCAGCAACAAATTTTCGAGATAtcttcttcgg                                                  >  2:6043175/1‑148 (MQ=16)
                                                                                                       aaGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                <  1:7958573/150‑1 (MQ=255)
                                                                                                         gTATTCGTAAGGTTGTTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                              <  2:3618661/150‑1 (MQ=255)
                                                                                                          tATTCTTATGGTTGTTTTTCTCCCTCTCCCGGGGGGAGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGAGATCTTCTTCTGTCTTg                                             <  2:4908414/150‑1 (MQ=255)
                                                                                                           aTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                            <  1:7361234/150‑1 (MQ=255)
                                                                                                            ttCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGTAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                           <  1:5805344/150‑1 (MQ=255)
                                                                                                             tcttatggttggttTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                          <  2:3501642/144‑1 (MQ=39)
                                                                                                             tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCAACCGAACCCACTCCCACAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGTGATATCTTATTCTGTGTTGTaa                                          >  1:1001509/1‑150 (MQ=37)
                                                                                                               gtttggttgtttTTCTCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                        <  2:683384/146‑1 (MQ=39)
                                                                                                               gTATGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                        <  1:10740866/150‑1 (MQ=255)
                                                                                                                  cggttggttTACTCCCTCCCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGTTATCTTCTTCTGTCTTGTAACagaa                                     <  1:921224/149‑1 (MQ=255)
                                                                                                                    gttggtcTTCCCCCCCCCCCTGGGGGTGGGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                   <  1:3201849/150‑1 (MQ=255)
                                                                                                                         ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                              <  2:5898068/150‑1 (MQ=255)
                                                                                                                          ttgtcgccctctccctGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                              <  1:1213530/144‑1 (MQ=255)
                                                                                                                                ccccccccctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGAGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                        <  1:1031213/144‑1 (MQ=255)
                                                                                                                                 cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                       <  1:3667211/150‑1 (MQ=255)
                                                                                                                                  ctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                      <  2:9049874/150‑1 (MQ=255)
                                                                                                                                     tccctGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCGCACCCAAACCCCCGCCCACCAGGGCGGGGGGGCCCCGCGCGCCCAAAGTGTTGGTGGTATAAAGATCTGGGGTGTAACAGAAGAACACGAAATCgggggg                   >  2:7910296/1‑147 (MQ=255)
                                                                                                                                        ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAACCCACACCCTCAAGGGAGAGGGGACCCAACGAGCACAAATTTTGAGAGATATATTCTTCTGTCGTTTAACAGAAGAAAAGAAAATCGGGCGTTTg                >  2:9101082/1‑150 (MQ=255)
                                                                                                                                         tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGGGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc               >  2:7051134/1‑150 (MQ=255)
                                                                                                                                           ggggtgtgggCGGGGTTGTGGGCATCAGCGCGCACGTTCACACTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGCTATCTTCTTCTGTCTTGTAACAGAAGAACAGAATATCGGGCTTTTgcgc             <  2:10099539/143‑1 (MQ=255)
                                                                                                                                           ggggcgcgggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc             <  1:4665968/143‑1 (MQ=255)
                                                                                                                                           gGGGAGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc             >  2:2304225‑M1/2‑150 (MQ=255)
                                                                                                                                             gtcgtGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt           >  1:10563457‑M1/6‑150 (MQ=255)
                                                                                                                                             ggtgagGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt           <  1:6459175/147‑1 (MQ=255)
                                                                                                                                              GAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg          >  1:10392241‑M1/1‑150 (MQ=255)
                                                                                                                                                GAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCCAACCCACTCCCACAAGGGGGAGGGGACCCATCGAGCACAAAATTTGAGTGATAAATTATTATGTCTTGTAACAGAAGAACAGAAAATAGGGCTTTTGCGCCGGaa        <  2:8203458‑M1/150‑1 (MQ=255)
                                                                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCAAAACCACACCCGCAATGGGGAGGGGGCCGATCGAGCCCAAAATTTGAGATATATATTATTATGTCTTTTAAAAGAAGAAAAGAAAATACGGCTTTTTCGCCAGAaaat     >  2:6401717/1‑150 (MQ=255)
                                                                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTACCACCAGGGAGAGGGGACCGATCGAGCACAAATTTTTAGAGTTATAATATTAAGTAATTTAACAAAAGACACCAAAAACGGGCTTTTGAGCCGGAatat     >  2:8221209/1‑150 (MQ=255)
                                                                                                                                                    GCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                         <  2:2441762‑M1/50‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGC  >  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: