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.

AGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGA  >  NC_000913/164447‑164743
                                                                                                                                                       |                                                                                                                                                    
agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtc                                                                                                                                                     >  2:2787680/1‑150 (MQ=255)
  cgatgatTGGGCGTTAACCCAAACATGTCGGGCGGGACGACCCGGA‑AAATACTGCACCGCCGCGAGGCACTAAAATGTTTTCTTATGGTTGG‑TTTTCTCCCGCTCCGTGGGGT‑TGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcac                                                                                                                                                   <  2:1883129/149‑1 (MQ=255)
   ttttaatgggCTTTATCCCAAAAATGTCGGGCCTTACGACCGGCA‑AAATATGCCACCGACGCGAGGGACGAAAATTTATTCGTTTGTTTGT‑TTTTCTCCCTCTCCGTGTGGG‑TGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                                  <  2:114620/143‑1 (MQ=255)
   tatcaaaGGGCTTTAACCCAAACATTACTGGCCGCACCACCCGCA‑AAAAACTCCACCCACGCGACGCACGAAAAATTATTCCTAAGTTTGT‑TTTTCTCCCACACCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                                  <  2:1994936/146‑1 (MQ=19)
        ttgggCGTTACCCAAAACTTGTGGGGCCGCACG‑CCCGGCAAAATACGCCACCGACGCTCCGCACCAAAATTTTTTCTTTTGGTTGT‑TTTTCCCCCTCTCCCTGTGGG‑TGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                             <  2:4460452/148‑1 (MQ=12)
        acgggCTTTAGCCCAAACATGTGTGGCCGGACGACCCGGA‑AATTAGTGCACCGACGCGACGGACGAAAATGTATTCTAACGGTTGT‑TTTTCTCCCTCCCCCTGTGGG‑TGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTtcaccctcac                                                                                                                                             <  2:295530/148‑1 (MQ=255)
           ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccccaaccc                                                                                                                                          >  2:87242/1‑149 (MQ=255)
             cGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCCaa                                                                                                                                        >  2:646505/1‑150 (MQ=255)
                tATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCAAAacc                                                                                                                                     >  1:771054/1‑150 (MQ=38)
                 aTCCCAAACATGTCTGGCCGGACGACCCGGA‑AAATACTGCACGGACGCGACGGACGAAAATGTATTCGTATGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                                    <  2:3898376/150‑1 (MQ=38)
                  tCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                   <  2:3568958/150‑1 (MQ=255)
                   cccAACCCTGTCTGGCCGGACGACCCGCA‑AAATAGTGCCCCGACGCGAGGGACTAAAATTTTTTGGTAGGTTTGT‑TTTTCTCCGTCTCCGGGTGGT‑CGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                                  <  2:3705476/150‑1 (MQ=255)
                    ccAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAaccacaa                                                                                                                                 >  1:3346910/1‑146 (MQ=38)
                    ccAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCAAAcccacac                                                                                                                                 >  1:2308812/1‑146 (MQ=38)
                           tGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:3176065/150‑1 (MQ=255)
                           ggccggcccGGACGCCCCGCC‑AACTCCGGCACCGCCCCGACGGACCACCAAGTATTCGCACGGTGGG‑TTGGCGCCCGCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:2739614/143‑1 (MQ=255)
                           ggTCGGGCCGGACCCCCCGCC‑AAATACGGCACCGAGGCGCCGCACGAAAATTTATTGGTATGGTTGT‑TTTTCGCCCACCCCGGGGGGT‑TGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                          <  2:1729422/149‑1 (MQ=255)
                             tCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAACCCACACCCACCAggg                                                                                                                        >  2:4001971/1‑150 (MQ=38)
                             tCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCGAAACCACTcccgcccggg                                                                                                                        >  1:277500/1‑143 (MQ=37)
                                 gccggccgCCCCGGA‑AAATCCTCCACCGACGCGACGGTCTAAAATTTTTTGTTATGTTTGG‑TTTTCTCCCTCGGCGTGTGGG‑TGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                    <  1:2634351/150‑1 (MQ=255)
                                 gCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTGCCCCCGCACCCCAACCCACACCCCCAAggggggg                                                                                                                    >  2:2954720/1‑146 (MQ=37)
                                  ccGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGACCCACACCCAAAACCACACCCACAAcggcgggg                                                                                                                   >  2:4101536/1‑145 (MQ=25)
                                  ccGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAACCCACACCCACAACGGGGAgg                                                                                                                   >  2:4101359/1‑150 (MQ=37)
                                   cGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCACCCCCACCCAAAACCACACCCACCAGGGAGgggg                                                                                                                  >  1:1480559/1‑150 (MQ=37)
                                   cGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                                  <  2:4723510/150‑1 (MQ=255)
                                   cGGACG‑CCCGGC‑AAAATCCGCCCCGCCGGGGCGGCCCTAAAAATTTTGCTTAGGTTGGTTTTTCACCCACCCCGTGGGGT‑CGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                                  <  2:1272845/150‑1 (MQ=255)
                                    ggACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCGCCACACCCAAAACCACAACCACAAGGGCGAgggg                                                                                                                 >  1:1426749/1‑150 (MQ=25)
                                        gacCCGTA‑AAAAACTGCCCCGCCGCGACGTACTAAATTTTATTCTTACGGTGGG‑TTTTCTCCCTCGCCGTGTGGG‑AGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                                             <  1:3500626/150‑6 (MQ=255)
                                        gacCCGGC‑AACTCCTCCACCGACGCGACGGCCGAAAAAGTATTCGTAAGGTTGG‑TTTACTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttagg                                                                                                             <  2:2386457/150‑6 (MQ=37)
                                          cccGGC‑AAAAACGCCCCCGCGGCGACGGACTAAAACTCTTTCGTCAGGTTGG‑TTTTCCCCCGCCCCCGGGGGG‑CGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                           <  2:2183810/150‑1 (MQ=255)
                                           ccGCC‑AAATCCGGCCCCGCCGCGACGGACGAAAACGTTTTCGTAAGGTTGG‑TTTTCTCCCCCCCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                          <  2:2521740/150‑1 (MQ=37)
                                              gcaaaaTACTGCCCCGACGGGCGGGACTAAAATGTTTTCGTATGTTTGG‑TTTTCGCCGTGTCCCTGTGGG‑TGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                                        <  2:4058975/147‑1 (MQ=255)
                                              gC‑AAATACTGCCCCGCCGCGCCGGACGAAAACTTATTCGTATGGTTGG‑TTTTCACCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                       <  2:2512286/150‑1 (MQ=255)
                                              gC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCCCACCCGAAACCACGCCCTCAAGGGGGAGGGGGCCGGTCGAg                                                                                                       >  1:746039/1‑150 (MQ=25)
                                              gC‑AAATACGGCACCGCCGCGCCGGACGAAACAGTATTCGTATGGTTGG‑TTTTCGCCCTCTCCCTGTGGG‑AGTGGGCGGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                       <  2:3895376/150‑1 (MQ=38)
                                              cc‑AAAGACTGCCCCGACGGGCCGGACCAAAAAGTATTCGTACGGTGGG‑TTTTCGCCCTCTCCCTGGGGG‑AGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                       <  1:2755647/149‑1 (MQ=255)
                                               c‑AAATCCTCCACCGAGGCGCCGGACCAAAAAGTATTCGAACGTTTGT‑TTTACCCCCACTCACTGGGGG‑AGCGGGCGGGGGTTTGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                      <  1:342522/150‑1 (MQ=255)
                                                 aaaTCCGGCACCGCGGCGCCGGTTGAAAATGTTTTGGTCTGGTTGG‑TTTTCGCCCTCGCCGGGGGGG‑AGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                     <  1:3458360/150‑1 (MQ=255)
                                                         cACCGCCGCTACGCACAAAAAATTATTCTTATGGTTTT‑TTTTCCCCCTCTCCGTGTGTT‑AGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                             <  1:3262724/150‑1 (MQ=255)
                                                           ccgccgCGACGGCCCAAAATTTTTTCTTATGGTTGG‑TTTTCCCCCCCCCCGTGGGGG‑TGTGGGGGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  1:1469459/150‑1 (MQ=255)
                                                           ccgccgCGACGGACGAAAACGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  2:3588470/150‑1 (MQ=255)
                                                           ccGACGCGCCGGACCAAAATTTATTCGTACGTTTGG‑TTTTCTCCCACTCCCTGTGGT‑TGCGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                           <  2:3834234/150‑1 (MQ=255)
                                                            cgacgcgacgTCCTATATATTATTCTTAGGCTTGG‑TTTTCGCCCGCGCCCTTTGGT‑TGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                          <  2:1473763/150‑1 (MQ=255)
                                                             gacgcgacgGACTAAAATGTATTCGTATGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                                         <  2:777288/150‑1 (MQ=255)
                                                              acgcgacgGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCCAAACCCCACCCACACCGGGGAGGGGGCCGCGACAGCACCAAATTTgagtg                                                                                        >  2:3052020/1‑148 (MQ=21)
                                                                     gtaggAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCGGTGGG‑AGTGGGCGGGGGTTAGGGCCTCAGCGCGCACGTTCACCCTCACGCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGTGATAtctt                                                                                 <  2:2737265/146‑1 (MQ=255)
                                                                         taaaaaGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                                             <  2:1153873/149‑1 (MQ=255)
                                                                         gAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCCCAAGGGAGAGGGGACCGTTCGAGCACAAATTTTGAGAGATAtcttcttc                                                                             <  2:3912934/150‑1 (MQ=255)
                                                                         gAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACACGAAAACACTAACGCCAGGGAGAGGGGGCCCATCGAGCACAAAATTTGTGAGAAAtattctta                                                                             >  2:747288/1‑149 (MQ=21)
                                                                          aaaaaGTATTCGTACGGTTGG‑TTTACTCCCTCACCATGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                                            <  2:1075304/150‑1 (MQ=255)
                                                                            aaaGTATTCGTACGGTTGG‑TTTTCTCCCTCTCCCGGTGGG‑AGAGGGCCGGGGTGAGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                          <  2:2731739/150‑1 (MQ=255)
                                                                             aaGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCGCACCCAAAACCACTACCACAACGGGGAGGGGGCCCATCGAGCACAAAATTTGTGAGATATATTATTCTGtg                                                                         >  2:4284166/1‑149 (MQ=21)
                                                                               ttATTCGTAAGGTTGG‑TCTGCCCCCCCCCCCGGTGGG‑TGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                       <  1:428416/149‑1 (MQ=255)
                                                                               gTATTCGTATGGTTGT‑TTTTCTCCATCTCCCTGTGGG‑ATTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACTAATTTTGAGAGATATCTTCTTCTGTCtt                                                                       <  2:380627/150‑1 (MQ=255)
                                                                                ttttCGTTTGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑TGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  2:2797492/148‑1 (MQ=255)
                                                                                ttttCGTATGTTTGG‑TTTTCTCCCTCTCCCGGGGGG‑AGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  1:923438/148‑1 (MQ=255)
                                                                                tttgcGTAAGGTTTG‑TTTTCTCCCTCTCCCTTTGGG‑TGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  2:2292383/146‑1 (MQ=39)
                                                                                tATTCGTTTGGTTGT‑TTTTCTCCCGCTCCGTGGGGG‑AGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                      <  1:923419/150‑1 (MQ=255)
                                                                                   tggTATGGTTGG‑TTTCCGCCCCCGCCCTGGGGG‑TGCGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  2:1816158/148‑1 (MQ=255)
                                                                                   tCGTATGGTTGG‑TTTTCTCCCTCGCCCGGTGGG‑TGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  2:4296628/150‑1 (MQ=39)
                                                                                   tCGTATGGTTGG‑TTTTCCCCCTCTCCCTGTGGT‑TGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                   <  2:1816316/150‑1 (MQ=255)
                                                                                     gTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                 <  2:3158360/150‑1 (MQ=255)
                                                                                       taGGTTGT‑TTTTCTCCCTCCCCGTGTGGG‑TGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                               <  1:2637701/149‑1 (MQ=255)
                                                                                       atgtttgg‑ttTTCGCCCTCCCCCTGTGGT‑TGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                               <  1:2637504/146‑1 (MQ=255)
                                                                                            tgg‑ttTTCTCCCTCTCCCTGTGGG‑AGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                          <  1:756927/150‑1 (MQ=255)
                                                                                            tgg‑ttTACTCCCTCTCCCTGGGGG‑CGTGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                          <  1:2297629/150‑1 (MQ=255)
                                                                                              g‑ttTTCCCCCCCCCCCGGGGGGTAGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGACGAACa                                                         <  1:1281904/150‑1 (MQ=255)
                                                                                                ttttCTCCCTCTCCCTGTGGG‑AGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                       <  1:1174195/150‑1 (MQ=255)
                                                                                                ttttCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                                       <  1:1174443/150‑1 (MQ=255)
                                                                                                     tccctctccctGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                  <  2:3125609/150‑1 (MQ=255)
                                                                                                     tccctctcccgGTGGG‑AGAGGGCGGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCACTCAAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCTGTATTGTAACAGAAGAACAGAATATc                                                  <  2:2394215/150‑1 (MQ=255)
                                                                                                     gccctctccctGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                  <  2:3125751/149‑1 (MQ=255)
                                                                                                                 tgGG‑TGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                      <  2:4174898/150‑1 (MQ=255)
                                                                                                                 gggg‑AGTGGGCGGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTTCCCCTCTCCCTCAAGGGAGAGGGGACCGATCGCGCACAAATTTTGAGAGATATCTTCTTCTGCCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                      <  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)
                                                                                                                                                       |                                                                                                                                                    
AGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGC‑AAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGA  >  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: ‑