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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC23.2% 345.5 / 49.3 253intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (79/111);  new base C (8/55);  total (87/166)
Fisher's exact test for biased strand distribution p-value = 1.48e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164433‑164727
                                                                                                                                                     |                                                                                                                                                  
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                   >  2:8881506/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                             >  1:7645810/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                             >  1:10607746/1‑150 (MQ=255)
            aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                       >  2:9828039/1‑150 (MQ=255)
              agagATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAATTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                                     <  2:1200912/150‑1 (MQ=39)
                 ttttaaggggCTTTACCCCAAACTGGGCGGGCCGGACGCCCCGCAAAATAGTGCACCGACGCGACGGACTAAAATTTTTTGTTATGTTTTTTTTTCCCCCTCCCCGGGGGGGTGTGGGCGGGGGTTTGGGCATCAGCGCGCACGTtcacc                                                                                                                                  <  2:9390445/143‑1 (MQ=255)
                   tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggccca                                                                                                                                >  2:5771373/1‑143 (MQ=255)
                    gAATGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGAGGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                               <  2:6747533/150‑1 (MQ=255)
                     aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                              <  2:10063433/150‑1 (MQ=255)
                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATCCTGCCCCGACGCGACGGACGAAAAAGTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                            <  2:2561919/150‑1 (MQ=255)
                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                            <  2:6887150/150‑1 (MQ=255)
                          gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccacccccga                                                                                                                         >  2:1676694/1‑148 (MQ=255)
                            tttATCCCAAACTGGGGGGGCGGGACGACCGGGAAAATACTGCACCGACGCGAGGGACTAAAAAGTATTGTTATGGTTGGTTTTCTCCCTCCCCCTGTGGTTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                       <  2:9790605/149‑1 (MQ=255)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGTGACcccaaccccaac                                                                                                                       >  2:373265/1‑150 (MQ=38)
                             ttATCCCAACCAGGTCTGGCCGGCCGACCGGGCAAATACTGCACCGCCGCGCCGGACTAAAAATTTTTCGTACGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:7899557/150‑1 (MQ=25)
                             ttATCCCAAACTTGTCGGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTATTCTTACGGTTGGTTTTCGCCCACCCCCGGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:6976897/150‑1 (MQ=255)
                             ttAGCCCAAACTGGTCTGGCCGGACGACCGGGCAAATACTGCACCGCCGCGAGGGACGAAAATGTTTTCGTACGGTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                      <  2:5787398/150‑1 (MQ=37)
                                  ccAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                 <  2:11009585/150‑1 (MQ=255)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCGaagacacta                                                                                                                 >  2:5469871/1‑143 (MQ=255)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCACACCCGAAcccacac                                                                                                                 >  2:3186467/1‑146 (MQ=38)
                                    aaaCATGTCTGGCAGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACAACCAAAaccacacac                                                                                                               >  1:6623919/1‑144 (MQ=37)
                                       ctggTCGGGCCGGACCCCCCGGAAATCCCGGCCCCGACGCGCGGGCCGAAAAAGTTTTGGTTGGGTGGGTTTTCTCCCCCCCCCGGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                            <  2:2028699/147‑1 (MQ=255)
                                        aTGTCGGGCCGGACGACCCGGAAAATACTCCACCGACGCCACGCACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                           <  2:7989721/150‑1 (MQ=38)
                                         tggcgggCCGGACGCCCCGGAAAATACTGCACCGACGCGACGGACTAAAATGTATTCGTATGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                          <  2:8711731/145‑1 (MQ=255)
                                           tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAACCCCCACCCACACggg                                                                                                        >  2:4793137/1‑150 (MQ=37)
                                              ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAATTTTTTGGTATGGTTGTTTTTCTCCCTCGCGCTGTGGTCGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  1:9157925/150‑1 (MQ=255)
                                              ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCACCCTCTCCCGGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                     <  2:821401/150‑1 (MQ=255)
                                               gcgggccgacCCGGCAAATCCTGCACCGCCGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgagat                                                                                                    <  2:5423386/144‑2 (MQ=255)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCCAACCCACACCCACaaccgggag                                                                                                    >  2:10358322/1‑143 (MQ=37)
                                                cgggccgacCCGCAAAAACCGCCACCGAGGCGCCGGACGAAAATGTTTTCGTATGTTGGGTTTTCGCCCCCCCCGGGGGGGTTTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                   <  1:4417494/145‑1 (MQ=255)
                                                ccGGACGCCCCGCAAAATCCTCCACCGCCGCGCCGCACTAAAATTTTTGCTTTTGGTTGGTTTCCCCCCTCTCCGTGTGGGAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                   <  1:6398284/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAACCACACCCCCCAgggggggg                                                                                                   >  1:1202921/1‑145 (MQ=25)
                                                  ggtcgccccGCCAAATACGGCCCCGCCGCGCCGCACTAAAATGTTTTCTTATGGTGGGTTTACTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggt                                                                                                 <  1:10416188/144‑2 (MQ=255)
                                                     cgacCCGGCAAATACGGCCCCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCGCCCCCTCCCGGTGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                              <  2:10717633/150‑1 (MQ=255)
                                                       tcccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCACCCCCAAACCCACGCCCACCCGGGCGGGGGGGCCgc                                                                                            >  1:1802127/2‑149 (MQ=25)
                                                       acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCCCACCCGCACCCACACGGGGGCGGGGGcccg                                                                                            >  1:8600067/1‑148 (MQ=25)
                                                         ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCACAAGGGGGCGGGGGCCGCTc                                                                                          >  1:10193711/1‑150 (MQ=37)
                                                         ccGGCAAAACCACCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  2:9692227/150‑1 (MQ=38)
                                                         ccGGCAAAAACTGCACCGACGCGAGGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  2:6451638/150‑1 (MQ=255)
                                                          cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACTCCCTCAAGGGAGAGGGGGCCGATCg                                                                                         >  2:7363665/1‑150 (MQ=38)
                                                          cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAACCACAACCCCAACGGCGAGGGGGcccaaag                                                                                         >  2:4322513/1‑145 (MQ=25)
                                                            gaaaaTCCTGCCCCGACGCGCCGGACTAAAAATTATTCGTAAGGTTGGTTTACTCCCTCGCCCTGTGGTAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                       <  1:6783149/148‑1 (MQ=255)
                                                            gCAAATACTGCACCGACGCGACGGACGAAAAAGTGTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCACAACCCCACCCCCACCGGGGAGGGGGCCgggaggg                                                                                       >  1:2340931/1‑144 (MQ=21)
                                                                  aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACCCCCCAAACCGCACCCACAAGGGGGAGGGGGCCGATAGAGCCCaaa                                                                                 >  1:5416018/1‑150 (MQ=25)
                                                                  aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAACCCACACCCCCACGGGGGAGGGGCCCGCTCGTGCACaaa                                                                                 >  1:7103696/1‑150 (MQ=25)
                                                                       cACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCCCCCGCACCCGAACCCACACCCCCCCGGGGGGGGGGCCCCGTCGGGCCCCAAtttt                                                                             >  2:7298613/1‑150 (MQ=21)
                                                                         ccgccgCGCCGTACGAAAAAGTATTCTTCCGGTGGGTTGTCCCCCGCGCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:4262030/150‑1 (MQ=255)
                                                                         ccgccgCGAGGTACTAAAAAGTATTGGTAGGGTGGTTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:4298478/150‑1 (MQ=37)
                                                                         ccgccgCGACGGCCGAAAAATTATTGGTAGGGTTGGTTTTCTCCCTCCCCGGGTGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:5231854/150‑1 (MQ=255)
                                                                         ccgccgCGACGGACTAAAATGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:10866627/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACTAAAAAGTTTTCGTCCGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:9889806/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACTAAAAAGTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:3525055/150‑1 (MQ=255)
                                                                         ccGACGCGACGTACTAAAATTTATTCTTATGGTGGTTTTCCGCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:6938040/150‑1 (MQ=37)
                                                                         ccGACGCGACGGACTAAAACTTTTTCGTAAGTTTGGTTTTCTCCCTCCCCCTGGGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:7379325/150‑1 (MQ=38)
                                                                         ccGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:11093008/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCAAACCCAAAACCACACACACAAGGGGGAGGGGGCCGAAAGAGCCCAAAatttta                                                                           >  1:9098699/1‑148 (MQ=21)
                                                                              gcgacgGACTAAAAAGTATTCGTAAGGTTGGATTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  1:9973545/150‑1 (MQ=255)
                                                                              gcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  2:10056570/150‑1 (MQ=255)
                                                                               cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAACAGCGCGCACGTTCACCCGCACCCGAACCCACTACCGCAAGGGCGGGGGGCCCGATCGAGCACAAAATTTgagagaaa                                                                     >  2:9164962/1‑148 (MQ=37)
                                                                                   ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCCCACCCAAAACCGCAACCACCAGGGAGGGGGGGCCCGGCGAGCACAAAATTTTAGTGGTAtatt                                                                 >  1:8162743/1‑150 (MQ=16)
                                                                                    gaccAAAATGTTTTCTTACGGTGGGTTCGCCCCCTCGCCGGGGGGGTGCGGGCGGGGGGTCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                <  1:3016227/150‑1 (MQ=255)
                                                                                       gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAACCCACACCCACACCGGAGAGGGGACCCAGCGACCACCAAAATTTAGAGATAtcttattc                                                             >  2:4516856/1‑150 (MQ=21)
                                                                                          aatttATTCGTTTGGTTGTTTTTCGCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                          <  1:1637537/146‑1 (MQ=255)
                                                                                               tttCTTACGGTGGGTTTACTCCCTCTCCCGGTGGGTGTGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                     <  2:9405387/149‑1 (MQ=255)
                                                                                                  cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                  <  2:7111709/150‑1 (MQ=255)
                                                                                                   ttttggttgtttTTCTCCCTCGCCGGGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  2:8384/146‑1 (MQ=255)
                                                                                                   gtttggtgggttTGCTCCCTCGCCGGGGGGGTGGGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  1:9501670/146‑1 (MQ=255)
                                                                                                      aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  2:5941244/150‑1 (MQ=255)
                                                                                                         tgggttTTCTCCCTCTCCCTGTGGGCGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:3423140/148‑1 (MQ=255)
                                                                                                          tggttTTCTCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:6135264/150‑1 (MQ=255)
                                                                                                          gggttGTCTCCCTCCCCCTGTGGGTGCGGGCCGGGGTGAGCGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:8967225/149‑1 (MQ=255)
                                                                                                                tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                    <  2:1355471/150‑1 (MQ=255)
                                                                                                                  gccctctcggtGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                  <  2:7645810/149‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCGAACCCGATCCCGCCAGGGAGGGGGGGCCGGTCGGGCCCCAATTTTGTGAGGTATCTTCTTCAGTCTTGTTACAGGAGCACAGAAAATCg                                 >  2:2696207/1‑150 (MQ=21)
                                                                                                                    cctcgccctGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:10038001/150‑1 (MQ=255)
                                                                                                                         ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAACCACTCCCGCAAGGGGGAGGGGGCCGGTCGGGCACAAAATTTGTGGGATATATTATTATGTCTTTTAACAGAAGAAAAGAAAATCGGGCgtt                           >  2:7262620/1‑150 (MQ=19)
                                                                                                                          cctGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                          <  1:10720002/150‑1 (MQ=255)
                                                                                                                              tgggtgggGGCGGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  2:6835720/150‑1 (MQ=255)
                                                                                                                                ggcgggggCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:4932802/145‑1 (MQ=255)
                                                                                                                                 gagtGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                   <  2:284490/150‑1 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAACCCACTCCCACAAGGGAGAGGGGGCCGAACGAGCACAAATTTTGTGAGATATCATCTTCTGTCGTGTAACAGAAGAACCGAAAAAAGGGCTTTTGCGCCGGa                  >  2:10987088/1‑150 (MQ=255)
                                                                                                                                   gagGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTTCGCCTGaa                 <  2:2641772/150‑1 (MQ=255)
                                                                                                                                     gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata               <  1:11131818/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACACCCACAAGGGGGAGGGGCCCCACAGAGCACAAAAGTTTGGAGAAATAATATTCTTGGGTGTAAAAGAAGTACAGAAAATCGGGCGTTTTCGCCAGtataa              >  2:2578007/1‑149 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCCCACCCGAACCCGCTCCCGCCAGGGCGAGCGGCCCCCTCCCGCCCCAATGGTGAGAGATAAATTCTTCTGGCGGGGAACAGCAGAACCGCACAGCGGGCGGTTGCGCCGGacgag              >  2:10122746/1‑146 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATGTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  2:8997364/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:7653905/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:7535004/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc           <  1:10113578/150‑1 (MQ=255)
                                                                                                                                          ggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGGAGAACTGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:5110124/150‑1 (MQ=255)
                                                                                                                                          ggggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:11005171/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACGCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:7557684/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:5923829/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:8331718/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  2:3329714/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:4676592/150‑1 (MQ=255)
                                                                                                                                            ggggcgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:1070607/145‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:8106545/150‑1 (MQ=255)
                                                                                                                                              tGAGGGCCTCAGCGCGCACGTTCACCCACACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga      <  1:5645625/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCACACCCAAAACCACAACCACAAGGGAGAGGGGACCCATCGAGCACAAAATTTGAGAGAAATAATATTAATTATTGTAACAGAAGAAAACAAAAACGGGCTTTTTCGCCAGAATATTTCGGAGaaaaa  >  2:8279411/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164433‑164727

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: