Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,595 T→G 11.3% intergenic (+61/‑135) 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,5950TG11.3% 356.9 / 5.4 123intergenic (+61/‑135)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (7/99);  new base G (15/2);  total (22/101)
Fisher's exact test for biased strand distribution p-value = 2.92e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGA  >  NC_000913/164447‑164743
                                                                                                                                                       |                                                                                                                                                    
agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtc                                                                                                                                                     >  2:2787680/1‑150 (MQ=255)
  cgatgatTGGGCGTTAACCCAAACATGTCGGGCGGGACGACCCGGAAAATACTGCACCGCCGCGAGGCACTAAAATGTTTTCTTATGGTTGGTTTTCTCCCGCTCCGTGGGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcac                                                                                                                                                   <  2:1883129/149‑1 (MQ=255)
   ttttaatgggCTTTATCCCAAAAATGTCGGGCCTTACGACCGGCAAAATATGCCACCGACGCGAGGGACGAAAATTTATTCGTTTGTTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                                  <  2:114620/143‑1 (MQ=255)
   tatcaaaGGGCTTTAACCCAAACATTACTGGCCGCACCACCCGCAAAAAACTCCACCCACGCGACGCACGAAAAATTATTCCTAAGTTTGTTTTTCTCCCACACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                                  <  2:1994936/146‑1 (MQ=19)
        ttgggCGTTACCCAAAACTTGTGGGGCCGCACGCCCGGCAAAATACGCCACCGACGCTCCGCACCAAAATTTTTTCTTTTGGTTGTTTTTCCCCCTCTCCCTGTGGGTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                             <  2:4460452/148‑1 (MQ=12)
        acgggCTTTAGCCCAAACATGTGTGGCCGGACGACCCGGAAATTAGTGCACCGACGCGACGGACGAAAATGTATTCTAACGGTTGTTTTTCTCCCTCCCCCTGTGGGTGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTtcaccctcac                                                                                                                                             <  2:295530/148‑1 (MQ=255)
           ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccccaaccc                                                                                                                                          >  2:87242/1‑149 (MQ=255)
             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCCaa                                                                                                                                        >  2:646505/1‑150 (MQ=255)
                tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCAAAacc                                                                                                                                     >  1:771054/1‑150 (MQ=38)
                 aTCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACGGACGCGACGGACGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                                    <  2:3898376/150‑1 (MQ=38)
                  tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                   <  2:3568958/150‑1 (MQ=255)
                   cccAACCCTGTCTGGCCGGACGACCCGCAAAATAGTGCCCCGACGCGAGGGACTAAAATTTTTTGGTAGGTTTGTTTTTCTCCGTCTCCGGGTGGTCGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                                  <  2:3705476/150‑1 (MQ=255)
                    ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAaccacaa                                                                                                                                 >  1:3346910/1‑146 (MQ=38)
                    ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCAAAcccacac                                                                                                                                 >  1:2308812/1‑146 (MQ=38)
                           tGTCTGGCCGGACGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:3176065/150‑1 (MQ=255)
                           ggccggcccGGACGCCCCGCCAACTCCGGCACCGCCCCGACGGACCACCAAGTATTCGCACGGTGGGTTGGCGCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:2739614/143‑1 (MQ=255)
                           ggTCGGGCCGGACCCCCCGCCAAATACGGCACCGAGGCGCCGCACGAAAATTTATTGGTATGGTTGTTTTTCGCCCACCCCGGGGGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:1729422/149‑1 (MQ=255)
                             tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAACCCACACCCACCAggg                                                                                                                        >  2:4001971/1‑150 (MQ=38)
                             tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCGAAACCACTcccgcccggg                                                                                                                        >  1:277500/1‑143 (MQ=37)
                                 gccggccgCCCCGGAAAATCCTCCACCGACGCGACGGTCTAAAATTTTTTGTTATGTTTGGTTTTCTCCCTCGGCGTGTGGGTGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                    <  1:2634351/150‑1 (MQ=255)
                                 gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTGCCCCCGCACCCCAACCCACACCCCCAAggggggg                                                                                                                    >  2:2954720/1‑146 (MQ=37)
                                  ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGACCCACACCCAAAACCACACCCACAAcggcgggg                                                                                                                   >  2:4101536/1‑145 (MQ=25)
                                  ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAACCCACACCCACAACGGGGAgg                                                                                                                   >  2:4101359/1‑150 (MQ=37)
                                   cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCACCCCCACCCAAAACCACACCCACCAGGGAGgggg                                                                                                                  >  1:1480559/1‑150 (MQ=37)
                                   cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                                  <  2:4723510/150‑1 (MQ=255)
                                   cGGACGCCCGGCAAAATCCGCCCCGCCGGGGCGGCCCTAAAAATTTTGCTTAGGTTGGTTTTTCACCCACCCCGTGGGGTCGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                                  <  2:1272845/150‑1 (MQ=255)
                                    ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCGCCACACCCAAAACCACAACCACAAGGGCGAgggg                                                                                                                 >  1:1426749/1‑150 (MQ=25)
                                        gacCCGTAAAAAACTGCCCCGCCGCGACGTACTAAATTTTATTCTTACGGTGGGTTTTCTCCCTCGCCGTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                                             <  1:3500626/150‑6 (MQ=255)
                                        gacCCGGCAACTCCTCCACCGACGCGACGGCCGAAAAAGTATTCGTAAGGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                                             <  2:2386457/150‑6 (MQ=37)
                                          cccGGCAAAAACGCCCCCGCGGCGACGGACTAAAACTCTTTCGTCAGGTTGGTTTTCCCCCGCCCCCGGGGGGCGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                           <  2:2183810/150‑1 (MQ=255)
                                           ccGCCAAATCCGGCCCCGCCGCGACGGACGAAAACGTTTTCGTAAGGTTGGTTTTCTCCCCCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                          <  2:2521740/150‑1 (MQ=37)
                                              gcaaaaTACTGCCCCGACGGGCGGGACTAAAATGTTTTCGTATGTTTGGTTTTCGCCGTGTCCCTGTGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                                        <  2:4058975/147‑1 (MQ=255)
                                              gCAAATACTGCCCCGCCGCGCCGGACGAAAACTTATTCGTATGGTTGGTTTTCACCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                       <  2:2512286/150‑1 (MQ=255)
                                              gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCCCACCCGAAACCACGCCCTCAAGGGGGAGGGGGCCGGTCGAg                                                                                                       >  1:746039/1‑150 (MQ=25)
                                              gCAAATACGGCACCGCCGCGCCGGACGAAACAGTATTCGTATGGTTGGTTTTCGCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                       <  2:3895376/150‑1 (MQ=38)
                                              ccAAAGACTGCCCCGACGGGCCGGACCAAAAAGTATTCGTACGGTGGGTTTTCGCCCTCTCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                       <  1:2755647/149‑1 (MQ=255)
                                               cAAATCCTCCACCGAGGCGCCGGACCAAAAAGTATTCGAACGTTTGTTTTACCCCCACTCACTGGGGGAGCGGGCGGGGGTTTGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                      <  1:342522/150‑1 (MQ=255)
                                                 aaaTCCGGCACCGCGGCGCCGGTTGAAAATGTTTTGGTCTGGTTGGTTTTCGCCCTCGCCGGGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                     <  1:3458360/150‑1 (MQ=255)
                                                         cACCGCCGCTACGCACAAAAAATTATTCTTATGGTTTTTTTTCCCCCTCTCCGTGTGTTAGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                             <  1:3262724/150‑1 (MQ=255)
                                                           ccgccgCGACGGCCCAAAATTTTTTCTTATGGTTGGTTTTCCCCCCCCCCGTGGGGGTGTGGGGGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  1:1469459/150‑1 (MQ=255)
                                                           ccgccgCGACGGACGAAAACGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  2:3588470/150‑1 (MQ=255)
                                                           ccGACGCGCCGGACCAAAATTTATTCGTACGTTTGGTTTTCTCCCACTCCCTGTGGTTGCGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  2:3834234/150‑1 (MQ=255)
                                                            cgacgcgacgTCCTATATATTATTCTTAGGCTTGGTTTTCGCCCGCGCCCTTTGGTTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                          <  2:1473763/150‑1 (MQ=255)
                                                             gacgcgacgGACTAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                                         <  2:777288/150‑1 (MQ=255)
                                                              acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCCAAACCCCACCCACACCGGGGAGGGGGCCGCGACAGCACCAAATTTgagtg                                                                                        >  2:3052020/1‑148 (MQ=21)
                                                                     gtaggAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCGGGGGTTAGGGCCTCAGCGCGCACGTTCACCCTCACGCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGTGATAtctt                                                                                 <  2:2737265/146‑1 (MQ=255)
                                                                         taaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                                             <  2:1153873/149‑1 (MQ=255)
                                                                         gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCCCAAGGGAGAGGGGACCGTTCGAGCACAAATTTTGAGAGATAtcttcttc                                                                             <  2:3912934/150‑1 (MQ=255)
                                                                         gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACACGAAAACACTAACGCCAGGGAGAGGGGGCCCATCGAGCACAAAATTTGTGAGAAAtattctta                                                                             >  2:747288/1‑149 (MQ=21)
                                                                          aaaaaGTATTCGTACGGTTGGTTTACTCCCTCACCATGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                                            <  2:1075304/150‑1 (MQ=255)
                                                                            aaaGTATTCGTACGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                          <  2:2731739/150‑1 (MQ=255)
                                                                             aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCGCACCCAAAACCACTACCACAACGGGGAGGGGGCCCATCGAGCACAAAATTTGTGAGATATATTATTCTGtg                                                                         >  2:4284166/1‑149 (MQ=21)
                                                                               ttATTCGTAAGGTTGGTCTGCCCCCCCCCCCGGTGGGTGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                       <  1:428416/149‑1 (MQ=255)
                                                                               gTATTCGTATGGTTGTTTTTCTCCATCTCCCTGTGGGATTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACTAATTTTGAGAGATATCTTCTTCTGTCtt                                                                       <  2:380627/150‑1 (MQ=255)
                                                                                ttttCGTTTGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  2:2797492/148‑1 (MQ=255)
                                                                                ttttCGTATGTTTGGTTTTCTCCCTCTCCCGGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  1:923438/148‑1 (MQ=255)
                                                                                tttgcGTAAGGTTTGTTTTCTCCCTCTCCCTTTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  2:2292383/146‑1 (MQ=39)
                                                                                tATTCGTTTGGTTGTTTTTCTCCCGCTCCGTGGGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  1:923419/150‑1 (MQ=255)
                                                                                   tggTATGGTTGGTTTCCGCCCCCGCCCTGGGGGTGCGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  2:1816158/148‑1 (MQ=255)
                                                                                   tCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  2:4296628/150‑1 (MQ=39)
                                                                                   tCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGTTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  2:1816316/150‑1 (MQ=255)
                                                                                     gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                 <  2:3158360/150‑1 (MQ=255)
                                                                                       taGGTTGTTTTTCTCCCTCCCCGTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                               <  1:2637701/149‑1 (MQ=255)
                                                                                       atgtttggttTTCGCCCTCCCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                               <  1:2637504/146‑1 (MQ=255)
                                                                                            tggttTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                          <  1:756927/150‑1 (MQ=255)
                                                                                            tggttTACTCCCTCTCCCTGGGGGCGTGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                          <  1:2297629/150‑1 (MQ=255)
                                                                                              gttTTCCCCCCCCCCCGGGGGGTAGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGACGAACa                                                         <  1:1281904/150‑1 (MQ=255)
                                                                                                ttttCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                       <  1:1174195/150‑1 (MQ=255)
                                                                                                ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                       <  1:1174443/150‑1 (MQ=255)
                                                                                                     tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                  <  2:3125609/150‑1 (MQ=255)
                                                                                                     tccctctcccgGTGGGAGAGGGCGGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCACTCAAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCTGTATTGTAACAGAAGAACAGAATATc                                                  <  2:2394215/150‑1 (MQ=255)
                                                                                                     gccctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                  <  2:3125751/149‑1 (MQ=255)
                                                                                                                 tgGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                      <  2:4174898/150‑1 (MQ=255)
                                                                                                                 ggggAGTGGGCGGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTTCCCCTCTCCCTCAAGGGAGAGGGGACCGATCGCGCACAAATTTTGAGAGATATCTTCTTCTGCCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                      <  2:1879098/149‑1 (MQ=255)
                                                                                                                            ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg                            <  2:3140408/150‑1 (MQ=255)
                                                                                                                              ggggggAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  2:812043/150‑1 (MQ=255)
                                                                                                                              ggggTGCGGGCCTCCGCGCGCCCGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:4548320/150‑1 (MQ=255)
                                                                                                                              ggggTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:4548262/150‑1 (MQ=255)
                                                                                                                              ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:4296447/150‑1 (MQ=255)
                                                                                                                              ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAGCCCAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAACTCGGGCTTTTGCGCCTGAATATTGCg                          <  2:2246729/150‑1 (MQ=255)
                                                                                                                              ggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  2:2928212/150‑1 (MQ=255)
                                                                                                                              ggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:1799984/150‑1 (MQ=255)
                                                                                                                              ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:4599114/150‑1 (MQ=255)
                                                                                                                              ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:2787680/150‑1 (MQ=255)
                                                                                                                              ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:1521026/150‑1 (MQ=255)
                                                                                                                              ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  2:3571517/150‑1 (MQ=255)
                                                                                                                              ggggTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTACTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                          <  1:3052020/150‑1 (MQ=255)
                                                                                                                                   gAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa                     <  1:4075084/150‑1 (MQ=255)
                                                                                                                                      ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                  <  2:679432/150‑1 (MQ=255)
                                                                                                                                        cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                <  2:4219898/150‑1 (MQ=255)
                                                                                                                                        cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                <  2:3701742/150‑1 (MQ=255)
                                                                                                                                            aGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCTAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTATTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATgg            <  2:4529217/150‑1 (MQ=255)
                                                                                                                                                cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg        <  2:3945529/150‑1 (MQ=255)
                                                                                                                                                cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg        <  2:277500/150‑1 (MQ=255)
                                                                                                                                                cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg        <  2:1226662/150‑1 (MQ=255)
                                                                                                                                                    cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt    <  2:745934/150‑1 (MQ=255)
                                                                                                                                                    cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt    <  2:746039/150‑1 (MQ=255)
                                                                                                                                                    cGTTCACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt    >  2:1106870/1‑150 (MQ=255)
                                                                                                                                                     gTTCACCCACACCCAAACCCACACCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGAGATATATTCATATGTCTTGTAAAAGAAGAACAGAAAAACGGGCTTTTGCGCCGGAATATTTCGGAGAAAAAGCATGGCCGGGAATg   >  2:598412/1‑150 (MQ=255)
                                                                                                                                                      gtcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa  <  1:4001971/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
AGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGA  >  NC_000913/164447‑164743

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: