Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,595 T→G 17.8% 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,5950TG17.8% 397.6 / 48.2 168intergenic (+61/‑135)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (16/120);  new base G (30/2);  total (46/122)
Fisher's exact test for biased strand distribution p-value = 2.81e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

AGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGA  >  NC_000913/164447‑164743
                                                                                                                                                    |                                                                                                                                                    
agttaTGAAAGGGCTTTATCCCAAACATGTCTGGCCGGACGACCGGTAAAATACTGCACCGACGCGAGGTACTAAAATTTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTtc                                                                                                                                                     <  2:480800/146‑1 (MQ=255)
  ataTGAAAGGGCTTTATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                                   <  2:2887253/148‑1 (MQ=38)
    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                                                 >  2:5700416/1‑150 (MQ=255)
    aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                                                 >  2:6265976/1‑150 (MQ=255)
     tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTgcccccc                                                                                                                                                >  1:5274309/1‑149 (MQ=255)
       aaaGGGCTTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                              <  2:5860616/150‑1 (MQ=255)
        aaGGGCGTTATCCCAAACATGTCTGGCGGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTtcaccctcac                                                                                                                                             <  2:2964667/150‑1 (MQ=255)
         tgggCTTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGAGGGACTAAAAAGTTTTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                            <  2:591227/149‑1 (MQ=37)
         cgggCTTTACCCCAAAATGGTGTGGCCGGACGACCCGGAAAATACTGCACCGACGCGCGGGAGGAAAATGTTGTCGTATGTTGGTTTTTCTCCCTCGCCGTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                                            <  2:6237700/149‑1 (MQ=255)
             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCaaa                                                                                                                                        >  1:270325/1‑150 (MQ=255)
             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCACCCCaaa                                                                                                                                       >  1:270348/1‑149 (MQ=255)
               ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTTCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCGaaac                                                                                                                                      >  1:3343894/1‑148 (MQ=38)
                  tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTGGCCCCCCCCCCCAAcccac                                                                                                                                   >  1:5689461/1‑148 (MQ=38)
                   cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCGCCCGCCCCCGAAccccca                                                                                                                                  >  1:1996976/1‑147 (MQ=38)
                    ccAAACTTGTCGGGCGGGCCGACCGGGAAAATACTGCCCCGACGCGACGTACGAAAACGTATTGGTATGTTTGGTTTCCCCCCTCTCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                                 <  2:5116368/150‑1 (MQ=255)
                    ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGCCCCCCACCCAAAcccacac                                                                                                                                 >  1:6792010/1‑146 (MQ=37)
                    caaaaaATTTCTGGCCGGACGACCGGCAAAATACTCCACCGACGCGACGTACTAAAATTTTTTCGTACGGTTGGTTTCCCCCCTCTCCCGGGGGTTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                                 <  1:5833073/148‑1 (MQ=255)
                     caaccctGTCTGGCCGGACGACCCGCCAAATCCGGCCCCGACGCGACGGACGAAACAGTTTTCGTTTGTTGGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                                <  2:6931764/144‑1 (MQ=255)
                     cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                                <  2:7189205/150‑1 (MQ=255)
                     cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAAccccaacc                                                                                                                                >  2:1395437/1‑145 (MQ=255)
                           tGTCTGGCCGGCCGACCCGCAAAATACCCCACCGCCGCGCCGCCCGAAAAAGTTTTCGAAGGGTGGGTTTCCTCCCTCTCCCTGTGGGTTTGGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:6132617/150‑1 (MQ=255)
                           tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:2308009/150‑1 (MQ=255)
                           tGTCTGGCCGGACGACCCGCCAAATACTGCACGGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:5939162/150‑1 (MQ=255)
                           tGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGCAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:2356929/150‑1 (MQ=255)
                           ggccgggccggccGCCCCGCAAAATAGTGCACCGACGCGCGGGACTAAAAAGTATTGGTATGGTTGTTTTTCCCCCGCCCCGGGTGGGAGCGGGCCGGGGTGGGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:2562804/145‑1 (MQ=255)
                           ggTCGGGCGGCACGACCCGGAAATTACTGCACCGACGCGACGGACGAAAAAGTTTTCGTACGTTTGTTTTTCACCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:5182982/149‑1 (MQ=37)
                             tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCCCGGGGCCCCGCACCCCAACCCACAcccccccggg                                                                                                                        >  2:2169217/1‑143 (MQ=25)
                             tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCAGCAGCGCGCACGGGGGCCCGCCCCCCAACCCACTCCCCCaccgg                                                                                                                        >  1:4188499/1‑146 (MQ=37)
                                gcccGGACCACCCGGCAAAGACTCCCCCGACGCGCGGGCCGAAAATGTTTTCTTATGGTTGGTTTTCTCCCCCTCCCTGGGGGTGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                     <  2:6771636/148‑1 (MQ=255)
                                 gccggccgACCCGGCAAATACTGCACCGACGCGCCGGACGAAAATGTTTTCGTATGGTTGGTTTTCTCCCCCCCCCTGTGGGAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                    <  2:5154517/150‑1 (MQ=255)
                                 gccgcccgacCCGCCAATCACTCACCCGTCGCGACGGACCAAAACTTATTTTAATGGTTGTTTTACCCCCTCCCCCGGTGGGAGCGGGCGGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                    <  1:1636483/144‑1 (MQ=255)
                                   cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCGAACACCCACCCCCACGgggggggg                                                                                                                  >  1:4962556/1‑150 (MQ=25)
                                     gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCAAAACCACACCCCCAACGGCGAGGGGa                                                                                                                >  2:1699595/1‑150 (MQ=25)
                                      acgacCGGGAAAATACTGCACCGACGCGAGGGACGAAAAATTATTCTTACGGTTGGTTTTCGCCCTCTCCCTGGGGGATAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                               <  2:5608188/150‑1 (MQ=255)
                                       ctacCCGGAAAATACTGCCCCGACGGTAGGGACTAAAATTTTTTGTTATGGTTGTTTTGCTCCCTCTGCGTGTGGGTTTGGGCGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                              <  1:4165252/148‑1 (MQ=255)
                                        gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCGCACCCAAAACCCAGCCCCCCAGGGAGAGGGGACCg                                                                                                             >  1:2960841/1‑150 (MQ=37)
                                         acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCCAACCCACACCCACAACGGCGCGGGGGccca                                                                                                            >  2:5252051/1‑148 (MQ=37)
                                           ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCACACCCGAACCCACACCCGCAAGGGAGGGGGGGcccagc                                                                                                          >  2:2215407/1‑146 (MQ=25)
                                           ccGGCAAAGACTGCACCGACGCGACGTACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                          <  2:6577497/150‑1 (MQ=255)
                                              gCAAATACTGCACCGACGCGACGTACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                       <  2:7142547/150‑1 (MQ=255)
                                                       gCACGGACGCGACGGACGAAAATGTTTTCTTTCGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                              <  1:5409374/150‑1 (MQ=37)
                                                        cACCGCCGCGACGCACGAAAACGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                             <  2:6848919/150‑1 (MQ=38)
                                                        cACCGACGCGACGGACGAAAAAGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                             <  2:2443529/150‑1 (MQ=255)
                                                         aCCTACGCTACGGACGAAAATGTATTCTTATGGTTGGTTTTCTCCCGCTCCCTTTGGTTGAGGGCCGGGGTTTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGTACCGATCGAGCACAAATTTTg                                                                                            <  1:3943243/150‑1 (MQ=255)
                                                          cccccgcgccgCACTAAAAAGTATGCTTAAGGTTTTTTTACGCCCTCCCGTGGGGGGTGAGGGGGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  1:1302730/146‑1 (MQ=255)
                                                          ccGACGCTACGGACTAAAATGTTTTGGTATGGTTGGTTTTCTCCCACTCCGTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  1:4847951/150‑1 (MQ=255)
                                                          ccGACGCTACGCACTAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  1:6438290/150‑1 (MQ=255)
                                                          ccGACGCGCCGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  2:3979770/150‑1 (MQ=255)
                                                          ccGACGCGACGGACGAAAATGTATTCGTTAGGTTGGTTTTCTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  1:7072439/150‑1 (MQ=255)
                                                          ccGACGCGACGGACGAAAAAGTATTCGGACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  2:4193762/150‑1 (MQ=255)
                                                                cgacgTACTAAAAATTATTCGTATGGTTGTTTTTCCCCCCCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                                                     <  1:155692/150‑1 (MQ=255)
                                                                 gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGGCCCCCACCCAAACCCCCTACCCCACAGGCGAGGGGGCCGATCGAGCCCAAAATTTTAGGGaaat                                                                                    >  2:1595694/1‑150 (MQ=21)
                                                                     gactAAAATTTATTGGTTTGGTTGTTTTTCTCCCTCCCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                                <  1:1157855/150‑1 (MQ=37)
                                                                       cgAAAATTTATTCGTTTGGTTGTTTTTCGCCCTCCCCGGGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                              <  2:2806739/150‑1 (MQ=255)
                                                                         aaaatGTATTCTTAAGTTTGGTTTACTCCAACTCCCTGTGGGAGTGGGCCGGGGTGCGGGAATCAGCGCGAACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                                            <  2:713498/150‑1 (MQ=255)
                                                                            atgttttCGTTTGGTTGGTTTTCTCCCACGCCCGGGGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                         <  1:5159772/145‑1 (MQ=255)
                                                                            aagaactCGTACGGTGGTATTACTCCCTCGCCGTGTGGGTGAGGGCCGTGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCATGGGTGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                         <  1:5159745/144‑1 (MQ=255)
                                                                            aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                         <  2:7149690/150‑1 (MQ=255)
                                                                              tttttCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                       <  2:4989802/147‑1 (MQ=255)
                                                                              ttATTCGTATGTTTGGTTTCCACCCTCTCCCTGTGGGATGGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                       <  1:5877418/149‑1 (MQ=255)
                                                                              gttttCGTTGGGTGGGTTTTCTCCCTCGCCCTGGGGGCGGGGGGGGGGGTGTGGGCATCCGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                       <  1:925509/147‑1 (MQ=255)
                                                                              gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCACACCCAACAGGGAGCGGGGGCCCAGCGAGCCCCAAAATTGAGAGATATAATATTCTGTGtt                                                                       >  2:3472963/1‑150 (MQ=16)
                                                                               ttttCTTATGGTTGGTTTACGCCCTCTCGCTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  1:2599987/148‑1 (MQ=255)
                                                                               tattgttAAGGGTGTTTTACTCCCTCCCCCTGTGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  1:6831484/144‑1 (MQ=255)
                                                                               tATTCTTATGGTTGTTTGTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  1:6831894/150‑1 (MQ=39)
                                                                                 ttCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                                    <  1:6237093/150‑1 (MQ=255)
                                                                                  tcttatggttggttTGCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  1:6372632/144‑1 (MQ=255)
                                                                                  tCGTATGTTTGGTTTCCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  1:5962616/150‑1 (MQ=39)
                                                                                    gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGACCCCCCCCCAAAACCACACCCCACAGGGAGAGGGGGCCGAACGAGCACAAAAATTTTTAGATATATTATTATGTATTTTAACa                                                                 >  1:6561091/1‑150 (MQ=12)
                                                                                      ttggttggttTTCTCCCTCTCCGTGTGGGCGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                               <  1:797472/148‑1 (MQ=255)
                                                                                       tggttggttTTCTCCCTCTCCGTGGGGGAGGGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                              <  1:2381360/149‑1 (MQ=255)
                                                                                          ttggttTTCTCCCTCTCCGGGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                                           <  2:1487028/150‑1 (MQ=255)
                                                                                              ttttCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                       <  2:1094025/150‑1 (MQ=255)
                                                                                              ttttCCCCCCCTCCGTGTGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                       <  1:3175407/150‑1 (MQ=255)
                                                                                                ttCCCCCCCGCCGTGTGGTTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                                     <  2:319096/150‑1 (MQ=255)
                                                                                                    ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGCACCCCCCCCCGAAACCTCGCCCGCAAGGGCGGGGGGGCCGCTCGCGCCCACAGTTTGAGAGATATCTTCGTCTGGCGTGTAACAGCAGCACAGAAAAACg                                                 >  2:5346661/1‑150 (MQ=16)
                                                                                                    ccctctccctGGGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                 <  1:7031549/150‑1 (MQ=255)
                                                                                                    cccactccctGGGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGGTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                 <  1:1154335/150‑1 (MQ=255)
                                                                                                     cccctcccgGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                <  1:5843845/147‑1 (MQ=255)
                                                                                                             tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCAAAACCACAACCACAAGGGCGAGGGGACCCAAAGAGCACCAAATTTTTGAGATATAATATTATTTGTTGTTAAAGAAGAACACAAAAAAGGGCgtttta                                        >  2:829920/1‑148 (MQ=255)
                                                                                                                gggcgcgggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                                     <  1:764362/144‑1 (MQ=255)
                                                                                                                gggCGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCATGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTAGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                                     <  1:2148265/150‑1 (MQ=255)
                                                                                                                          gggggggggggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                           <  1:3740611/149‑1 (MQ=255)
                                                                                                                          cGGGGTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                           <  2:3343894/150‑1 (MQ=255)
                                                                                                                            gggTGAGGGCATCCGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg                         <  2:4093177/150‑1 (MQ=255)
                                                                                                                            gggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCACAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGAACTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg                         <  1:1803394/150‑1 (MQ=255)
                                                                                                                             ggTGGGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                        <  2:2728911/150‑1 (MQ=255)
                                                                                                                             ggTGCGGGCCTCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCGAGCGAGCACAAATTTTGAGAGATCTCTTCTTCTGTCTTGTAACAGAAGAACAGAAACTCGGGCTTTTGCGCCTGAATATTGCGga                        <  2:5830948/150‑1 (MQ=255)
                                                                                                                             ggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                        <  1:1699595/150‑1 (MQ=255)
                                                                                                                                gAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa                     <  2:5082732/150‑1 (MQ=255)
                                                                                                                                  gggCATCAGCGCGCACGTTCACCCACACACAAACCCTAACCCACCACGGCCAGGGGACCGATCCAGCACCAAATTTTAGAGATATAATCATATGTCTTTTAACAAAAGAAAAGAAAAAAGGGCGTTTGCGCCGGAATATTTCGGAGaaaa                   >  2:3461702/1‑150 (MQ=255)
                                                                                                                                   ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                  <  1:5070034/150‑1 (MQ=255)
                                                                                                                                      aTCAGCGCGCACGTTCACCCAAACCCAAAACCAAAACCTCAAGGGAGAGGGGACCCATAGATCACAAAATTTGATATATATATTATTAAGTCTTTTAACAGAAGAACAGCAAAACGGGCTTTTTCGCCGGAAAAATTCGGGGAAAAAGCa               >  2:4682484/1‑150 (MQ=255)
                                                                                                                                            gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg         <  1:6634779/150‑1 (MQ=255)
                                                                                                                                             cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg        <  2:5956506/150‑1 (MQ=255)
                                                                                                                                             cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCCTCTGTCATGTAACAGAAGAACAGAAAAACGGGCTTTTGCGCCAGAATATTGCGGAGAAAAAGCATGGCCgg        <  1:4103472/150‑1 (MQ=255)
                                                                                                                                              gcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCggg       <  1:366673/150‑1 (MQ=255)
                                                                                                                                               cACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa      <  2:1216054/150‑1 (MQ=255)
                                                                                                                                                 cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt    <  2:4567010/150‑1 (MQ=255)
                                                                                                                                                  gTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTACTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAGCATGGCCGGGAATGa  >  2:4471546/1‑150 (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: