Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,564 A→T 19.1% intergenic (+30/‑166) 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,5640AT19.1% 650.6 / 33.4 294intergenic (+30/‑166)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (124/92);  new base T (2/55);  total (127/167)
Fisher's exact test for biased strand distribution p-value = 4.19e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 319 total aligned reads displayed.

TTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAA  >  NC_000913/164415‑164712
                                                                                                                                                     |                                                                                                                                                    
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                      >  1:8361053/1‑150 (MQ=255)
ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                      <  1:10430320/150‑1 (MQ=255)
gcgtccaattggctcgcggtcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGgaga                                                                                                                                                      >  1:6046642‑M1/145‑150 (MQ=255)
gcgtccaattggctcgcggtcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGgaga                                                                                                                                                      >  1:6564157‑M1/145‑150 (MQ=255)
gcgtccaattggctcgcggtcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGgaga                                                                                                                                                      >  2:8201144‑M1/145‑150 (MQ=255)
gcgtccaattggctcgcggtcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGgaga                                                                                                                                                      >  1:1230997‑M1/145‑150 (MQ=255)
 tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGagag                                                                                                                                                     <  1:6788889/150‑1 (MQ=255)
 tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGagag                                                                                                                                                     <  2:9548341/150‑1 (MQ=255)
 cgtccaattggctcgcggtcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGagag                                                                                                                                                     <  2:9297903‑M1/7‑1 (MQ=255)
 cgtccaattggctcgcggtcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGagag                                                                                                                                                     >  1:6116233‑M1/144‑150 (MQ=255)
  gtccaagtggctcgcgggcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAgg                                                                                                                                                    >  1:4289189‑M1/143‑150 (MQ=255)
  cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                                    >  1:3534802/1‑150 (MQ=255)
  cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                                    <  2:10744167/150‑1 (MQ=255)
   tccaattggctcgcggtcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctccctcaaGGGAGAggg                                                                                                                                                   >  2:687780‑M1/142‑150 (MQ=255)
   tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                                   >  1:4639352/1‑150 (MQ=255)
   tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                                   >  1:4783731/1‑150 (MQ=255)
   tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                                   <  1:6458424/150‑1 (MQ=255)
    ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                                  >  1:3860258/1‑150 (MQ=255)
    ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAggga                                                                                                                                                  <  1:9399495/150‑2 (MQ=255)
     gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAGGGATTCGTAAGGTTGGTGTTCTCCCTCTCCCTGTGGGAGAGGGcc                                                                                                                                                 >  2:5716001/1‑150 (MQ=255)
     caattggctcgcggtcattcccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctccctcaaGGGAgagggga                                                                                                                                                 >  2:482142‑M1/140‑148 (MQ=255)
      aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                                <  2:5394194/150‑1 (MQ=255)
      aaaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                                >  2:3057059/1‑150 (MQ=255)
       aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                               <  2:7223008/150‑1 (MQ=255)
            gCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                          <  2:8470439/150‑1 (MQ=255)
             cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                         <  2:4449235/150‑1 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        <  2:3233338/150‑1 (MQ=255)
              gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                        <  2:9674471/150‑1 (MQ=255)
               tACCGTGAGGTGCAAAAAGAGATGAACGGGCGTTATCCCAAACAGGTCTGGCCGGACGACCGGGCAAATCCTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGCGGGCCGGGGTGAggg                                                                                                                                       <  2:966188/150‑1 (MQ=255)
                 ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGACAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                     >  1:5463539/1‑150 (MQ=255)
                 ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                     <  2:51811/150‑1 (MQ=255)
                 ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                     >  2:2197027/1‑150 (MQ=255)
                 ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                     >  2:8515143/1‑150 (MQ=255)
                  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCag                                                                                                                                    >  1:5489914/1‑149 (MQ=255)
                  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGcct                                                                                                                                    >  1:1451571/1‑148 (MQ=255)
                    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                  >  2:3112173/1‑150 (MQ=255)
                        gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                              >  1:7764147/1‑150 (MQ=255)
                         ccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGagag                                                                                                                                                     <  2:8704012‑M1/7‑1 (MQ=255)
                         ccggccatgctttttctccgcaatattcaggcgcaaaagcccgattttctgttcttctgttacaagacagaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGagag                                                                                                                                                     >  1:8704012‑M1/120‑126 (MQ=255)
                            aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTAGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                          <  2:10361790/150‑1 (MQ=255)
                            aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTAGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                          <  2:2397362/150‑1 (MQ=255)
                                  tgttGAATGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTACGGTTGGTTTTCTCCCCCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtcac                                                                                                                    <  2:6712214/147‑1 (MQ=255)
                                    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTgcgccc                                                                                                                  >  1:7958154/1‑150 (MQ=255)
                                      gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTtgaccccc                                                                                                                >  1:261030/1‑148 (MQ=255)
                                      gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgc                                                                                                                >  1:5284801/1‑148 (MQ=255)
                                      gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGgcacccgc                                                                                                                >  2:6443313/1‑148 (MQ=255)
                                          gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggtcccccaccc                                                                                                            >  1:1017066/1‑150 (MQ=255)
                                               ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCCCACCCCaaac                                                                                                       >  2:7819694/1‑148 (MQ=255)
                                               ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGACCCACACCCCAAcc                                                                                                       >  1:6369974/1‑150 (MQ=38)
                                               ttATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACGGCCCCGAGGCTAGTGACTAAAATTTTTTCTTATGTTTGTTTTTCTCCCTACCCCGGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                       <  2:293729/150‑1 (MQ=255)
                                               taaccccAAACAGGTGGGGCCGGCCGACCCGGCAAATACTGCACCGCGGCGACGGACTAAAAAGTTTTGGTATGGTTGGTTTTCGCCCTCCCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                       <  2:5616464/146‑1 (MQ=255)
                                                 aTCCCAACCATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTACTCCCTCTCCCTGTGGGATAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                     <  2:6337069/150‑1 (MQ=255)
                                                 aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCGCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                     <  2:9004316/150‑1 (MQ=255)
                                                   cccAAACATGTCTGGCCGGCCGACCCGGCAAATACTGCACCGCCGCGACGGACTAAAAATTATTCTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                   <  2:2572175/150‑1 (MQ=38)
                                                   cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACCCCCAAAccaaca                                                                                                   >  1:9552340/1‑146 (MQ=38)
                                                    ccAATCTTGGCTGGCCGGACCACCGGGAAATTCCGGCCCCTCCGCGCCGGACGAAAATGTATTCTTTGGGTTGGTTTTCGCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                  <  1:9849500/150‑1 (MQ=255)
                                                    ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCCAAccccccc                                                                                                  >  1:2089773/1‑146 (MQ=38)
                                                     caaccttGTCGGGCGGGACGCCCGGGCAAATCCGGCCCCGCCGCGTGGGCCGAAAATGTTTTCGTTGGGTGGGTTTTCTCCCTCGCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCtctcn                                                                                                 <  1:10373179/144‑2 (MQ=255)
                                                     cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAAacccacac                                                                                                 >  2:7402179/1‑145 (MQ=37)
                                                          ctgcgtGGCCGGACGACCCGGCAAATACGGCCCCGACGCGACGGACTAAAAAGTTTTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGGGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                            <  2:5029656/145‑1 (MQ=37)
                                                              cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                        <  2:2665951/150‑1 (MQ=255)
                                                               tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                         >  1:10289832/1‑100 (MQ=255)
                                                               tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCACAAGGgag                                                                                       >  1:9825867/1‑150 (MQ=38)
                                                               gggccggccGCCCCGGCAAATACGGCCCCGCCGCGCCGGCCGAACAAGTATTCGTACGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgag                                                                                       <  2:10037604/149‑1 (MQ=25)
                                                                 gCCGGACGCCCCGCAAAATCCGCCCCCTACGCTCCGGACTAAATATTATGCGTTTGGTTGTTTTTCTCCCGCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                     <  1:2387064/150‑1 (MQ=255)
                                                                 gCCGGACGACCCGGAAAATACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                     <  2:7060138/150‑1 (MQ=255)
                                                                      acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCGCACCCCAACCCACACCCACAAGGGGGAgggggg                                                                                >  1:10818156/1‑148 (MQ=37)
                                                                      acgacCCGGAAAATACTGCCCCGCGGCGACGGACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                <  2:5997348/150‑1 (MQ=255)
                                                                        gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCACCCCAAAACCACAACCACCAGGGGGAGGGGGCCg                                                                              >  1:1865546/1‑150 (MQ=25)
                                                                           ccgtaaaaTCCTGCACCGCCGCGCGGGACCAAAACTTTTTGGTAGGGTTGGTTTGCCCCCCCCCCGTGGGGGTGCGGGCGGGGGGGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                           <  2:3860258/145‑1 (MQ=255)
                                                                              gCAAATCCGGCCCCGACGCGCCGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCCCCCCCGGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                        <  2:5811129/150‑1 (MQ=255)
                                                                                       gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCACACCCACACCCGCTCCCGCCCGGGAGAGGGGGCCGAGCGCGCCCAAAatt                                                               >  2:10851448/1‑150 (MQ=25)
                                                                                        cacccccgcgacgGACGAAAATTTATCCTCATGGTTTTTTAACCCCCTAGCCCTGTGTTGGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                              <  1:4045932/144‑1 (MQ=255)
                                                                                          ccgccgCGCGGGACTAAAATTTTTTCGTATGGTGGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                            <  1:7252915/150‑1 (MQ=255)
                                                                                          cccacgcgacgCACCAAAAAGTATTCGAAAGGTTGGTTTACACCCACCCCCTGTGGCAGAGGGCCGGGGTGAGGCAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                            <  2:8796251/147‑1 (MQ=255)
                                                                                            gcggcgacgTACTAAAAAGTTTTCTTATGGTTGTTTTCCGCCCTCCCCCGGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                          <  1:6852933/147‑1 (MQ=255)
                                                                                             agaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                                                                                          >  2:1535925‑M1/52‑101 (MQ=255)
                                                                                             agaagaagatatctctcaaaatttgtgctcgatcggtcccctctcccttgaGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                                                                                          <  1:1535925‑M1/50‑1 (MQ=255)
                                                                                             acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACAACCACAAGGGAGAGTGGACCCAACGAGCACAAAATTTtagag                                                         >  2:2789576/1‑150 (MQ=25)
                                                                                                     gactAAAATGTTTTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                 <  2:800747/150‑1 (MQ=38)
                                                                                                      acgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAAAACACACCCACAAGGGGGCGGGGCCCGATCGAGCCCAAATTTTTAGATATAtcatat                                                >  1:9941741/1‑146 (MQ=21)
                                                                                                       cgAAAAAGTATTCTTACGGTTGGTTTCCCCCCTCTCCCTGTGGCAGAGGCCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcta                                               <  1:5734634/150‑2 (MQ=255)
                                                                                                         aaaaaGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                             <  2:10046251/150‑1 (MQ=255)
                                                                                                            atgTATTCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                          <  2:7958154/148‑1 (MQ=255)
                                                                                                              gTATTCGTAAGGTTGTTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                        <  2:3618661/150‑1 (MQ=255)
                                                                                                                  tcttatggttggttTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                    <  2:3501642/144‑1 (MQ=39)
                                                                                                                  tCGTATGGTTGGTTTCCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCACGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGATCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                    <  1:3112173/150‑1 (MQ=255)
                                                                                                                  tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCAACCGAACCCACTCCCACAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGTGATATCTTATTCTGTGTTGTaa                                    >  1:1001509/1‑150 (MQ=37)
                                                                                                                    gTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                  <  1:10740861/150‑1 (MQ=255)
                                                                                                                       tggttggttTTCTCCATCCCCGTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                               <  1:921273/149‑1 (MQ=255)
                                                                                                                       cggttggttTTCCCCCTCTCCCTGTTGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                               <  1:2910042/149‑1 (MQ=255)
                                                                                                                       aGGTTGTTTTTCGCCCTCTCCCTGTGGGTTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                               <  2:1281411/150‑1 (MQ=255)
                                                                                                                        gtttggttTTCGCCCGCTCCGTGGGGGTGTGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                              <  2:187998/148‑1 (MQ=255)
                                                                                                                         gttgtttTTCCCCCGCGCCCGGTGGGCGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                             <  2:5155676/150‑1 (MQ=255)
                                                                                                                         gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCGAACCCACCCCCTCAAGGGAGAGGGGACCGATCGCGCCCAAAGTTTGTGAGATATCTTCTTCTGTCTTGTGACagaaga                             >  1:8167551/1‑150 (MQ=37)
                                                                                                                              ttttCTCCCTCTCCGTGGGGGTGTGGGGGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                        <  1:7402179/150‑1 (MQ=255)
                                                                                                                              ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                        <  2:5898068/150‑1 (MQ=255)
                                                                                                                                  ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCGCAAGGGGGAGGGGCCCCCTAGCGCACAAATTTTGGGAGATATAATATTCAGTCTTGGTACAGAAGAACAGtaaaa                    >  2:4560660/1‑149 (MQ=19)
                                                                                                                                            ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAACCCACACCCTCAAGGGAGAGGGGACCCAACGAGCACAAATTTTGAGAGATATATTCTTCTGTCGTTTAACAGAAGAAAAGAAAATCGGGCGTTTg          >  2:9101082/1‑150 (MQ=255)
                                                                                                                                               gGGGAGTGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc       >  2:2304225‑M1/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAA  >  NC_000913/164415‑164712

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: