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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5810AC20.9% 180.2 / 26.6 130intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (27/74);  new base C (6/23);  total (33/97)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729
                                                                                                                                                       |                                                                                                                                                    
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  2:3683047/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGAAt                                                                                                                                                     >  1:889766/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  2:1373834/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  1:3158360/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACg                                                                                                                                          >  1:4549101/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                         >  2:1957046/1‑150 (MQ=255)
              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                        >  1:3834234/1‑149 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtc                                                                                                                                       >  2:2787680/1‑150 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                       >  2:1799984/1‑148 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgc                                                                                                                                       >  1:3571517/1‑148 (MQ=255)
                 cgatgatTGGGCGTTAACCCAAACATGTCGGGCGGGACGACCCGGAAAATACTGCACCGCCGCGAGGCACTAAAATGTTTTCTTATGGTTGG‑TTTTCTCCCGCTCCGTGGGGT‑TGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcac                                                                                                                                     <  2:1883129/149‑1 (MQ=255)
                  ttttaatgggCTTTATCCCAAAAATGTCGGGCCTTACGACCGGCAAAATATGCCACCGACGCGAGGGACGAAAATTTATTCGTTTGTTTGT‑TTTTCTCCCTCTCCGTGTGGG‑TGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                    <  2:114620/143‑1 (MQ=255)
                  tatcaaaGGGCTTTAACCCAAACATTACTGGCCGCACCACCCGCAAAAAACTCCACCCACGCGACGCACGAAAAATTATTCCTAAGTTTGT‑TTTTCTCCCACACCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                    <  2:1994936/146‑1 (MQ=19)
                  gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGtgcccc                                                                                                                                    >  2:4102342/1‑145 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccccaaccc                                                                                                                            >  2:87242/1‑149 (MQ=255)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCCaa                                                                                                                          >  2:646505/1‑150 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTATGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:720442/150‑1 (MQ=255)
                               tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCAAAacc                                                                                                                       >  1:771054/1‑150 (MQ=38)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACGGACGCGACGGACGAAAATGTATTCGTATGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                      <  2:3898376/150‑1 (MQ=38)
                                 tCCCAACCATGTCTGGCAGGAGTACGCGTAATTTAGTGCACAGTCGCGACTGACGAAATTGTATTCGTATGGTTGT‑TTTTCTCGCTATCACTGTGGT‑AGTGGGCAGGGGTGTGTGAATCAGCGCGCACGTTCACCCTCACCCTaaaactc                                                                                                                     <  2:3568887/150‑6 (MQ=255)
                                 tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                     <  2:3568958/150‑1 (MQ=255)
                                  cccAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTACGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                    <  2:564081/150‑1 (MQ=255)
                                  cccAAACTTGTGGGGCCGGACGACCCGCAAAACACGGCCCCGACGCGACGGACGAAAATGTATTCGTATGGTTGG‑TTTTCCCCCCCTCCCTGTGGG‑AGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                    <  1:3987728/150‑1 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCAAAcccacac                                                                                                                   >  1:2308812/1‑146 (MQ=38)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCCCCCCCCCCCCCAcccacacc                                                                                                                  >  1:3362530/1‑145 (MQ=37)
                                          ggccggcccGGACGCCCCGCCAACTCCGGCACCGCCCCGACGGACCACCAAGTATTCGCACGGTGGG‑TTGGCGCCCGCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:2739614/143‑1 (MQ=255)
                                          ggTCGGGCCGGACCCCCCGCCAAATACGGCACCGAGGCGCCGCACGAAAATTTATTGGTATGGTTGT‑TTTTCGCCCACCCCGGGGGGT‑TGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:1729422/149‑1 (MQ=255)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCGAAACCACTcccgcccggg                                                                                                          >  1:277500/1‑143 (MQ=37)
                                                gccggccgCCCCGGAAAATCCTCCACCGACGCGACGGTCTAAAATTTTTTGTTATGTTTGG‑TTTTCTCCCTCGGCGTGTGGG‑TGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  1:2634351/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCACCCCCAAACCCCCTCCCACCAggggggg                                                                                                      >  2:3926941/1‑146 (MQ=37)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTGCCCCCGCACCCCAACCCACACCCCCAAggggggg                                                                                                      >  2:2954720/1‑146 (MQ=37)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAACCCACACCCACAACGGGGAgg                                                                                                     >  2:4101359/1‑150 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCACCCCCACCCAAAACCACACCCACCAGGGAGgggg                                                                                                    >  1:1480559/1‑150 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                    <  2:4723510/150‑1 (MQ=255)
                                                  cGGACG‑CCCGGCAAAATCCGCCCCGCCGGGGCGGCCCTAAAAATTTTGCTTAGGTTGGTTTTTCACCCACCCCGTGGGGT‑CGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                    <  2:1272845/150‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCGCCACACCCAAAACCACAACCACAAGGGCGAgggg                                                                                                   >  1:1426749/1‑150 (MQ=25)
                                                          ccGCCAAATCCGGCCCCGCCGCGACGGACGAAAACGTTTTCGTAAGGTTGG‑TTTTCTCCCCCCCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:2521740/150‑1 (MQ=37)
                                                            ggCATAAACGGCACCGTCGCGCCTGACGAAATTGTATTCGTACGTTGTT‑TTTTCACCTTCTCGAGGGGTG‑AGTGGGCGGGGTTGAGGGCATCAGCGCGCAGGTTCACCCTCACCCTAACCCTCTCCCTCAATGGAGAGGGGACCGTTCGa                                                                                          <  2:446001/150‑1 (MQ=255)
                                                            gccAAATCCTGCCCCGACGCGCCGGACCAAAATGTATCCTCACGTTTGG‑TTGACCCCCTCCCCGTGGGGG‑TGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:4059325/148‑1 (MQ=255)
                                                            gcaaaaTACTGCCCCGACGGGCGGGACTAAAATGTTTTCGTATGTTTGG‑TTTTCGCCGTGTCCCTGTGGG‑TGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:4058975/147‑1 (MQ=255)
                                                             gCAAATACTGCCCCGCCGCGCCGGACGAAAACTTATTCGTATGGTTGG‑TTTTCACCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  2:2512286/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCCCACCCGAAACCACGCCCTCAAGGGGGAGGGGGCCGGTCGAg                                                                                         >  1:746039/1‑150 (MQ=25)
                                                             ccAAAGACTGCCCCGACGGGCCGGACCAAAAAGTATTCGTACGGTGGG‑TTTTCGCCCTCTCCCTGGGGG‑AGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  1:2755647/149‑1 (MQ=255)
                                                              cAAATCCTCCACCGAGGCGCCGGACCAAAAAGTATTCGAACGTTTGT‑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)
                                                                         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)
                                                                                       gAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCCCAAGGGAGAGGGGACCGTTCGAGCACAAATTTTGAGAGATAtcttcttc                                                               <  2:3912934/150‑1 (MQ=255)
                                                                                       gAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACACGAAAACACTAACGCCAGGGAGAGGGGGCCCATCGAGCACAAAATTTGTGAGAAAtattctta                                                               >  2:747288/1‑149 (MQ=21)
                                                                                        aaataGTATCCGTTCTGTTGT‑TTTTCCCCCTCACCCTTTGGG‑TGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTAACCCTAACCCTATCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTGAGAGATAtcttcttat                                                              <  2:1077360/150‑3 (MQ=255)
                                                                                        aaaaaGTATTCGTACGGTTGG‑TTTACTCCCTCACCATGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                              <  2:1075304/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)
                                                                                             gTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:2819206/150‑1 (MQ=255)
                                                                                              ttttggTATGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  2:2370489/145‑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)
                                                                                                   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)
                                                                                                     aaGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  2:363033/150‑1 (MQ=255)
                                                                                                         ttgg‑ttTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                             <  2:3049672/150‑1 (MQ=255)
                                                                                                          tgg‑ttTTCTCCCTCTCCCTGTGGG‑AGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                            <  1:756927/150‑1 (MQ=255)
                                                                                                          tgg‑ttTACTCCCTCTCCCTGGGGG‑CGTGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                            <  1:2297629/150‑1 (MQ=255)
                                                                                                           tgtttttCTCCCTCTCCCTGGGGG‑AGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                            <  2:1867597/147‑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)
                                                                                                                  ctccctctccctGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                     <  2:2028665/150‑1 (MQ=255)
                                                                                                                   tccctctccctGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                    <  2:3125609/150‑1 (MQ=255)
                                                                                                                   tccctctcccgGTGGG‑AGAGGGCGGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCACTCAAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCTGTATTGTAACAGAAGAACAGAATATc                                    <  2:2394215/150‑1 (MQ=255)
                                                                                                                               tgGG‑TGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  2:4174898/150‑1 (MQ=255)
                                                                                                                               gggg‑AGTGGGCGGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTTCCCCTCTCCCTCAAGGGAGAGGGGACCGATCGCGCACAAATTTTGAGAGATATCTTCTTCTGCCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  2:1879098/149‑1 (MQ=255)
                                                                                                                               gggg‑AGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  2:4174916/149‑1 (MQ=255)
                                                                                                                                          ccGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg              <  2:3140162/150‑1 (MQ=255)
                                                                                                                                          ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg              <  2:3140408/150‑1 (MQ=255)
                                                                                                                                            ggggggAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:812043/150‑1 (MQ=255)
                                                                                                                                            ggggTGCGGGCCTCCGCGCGCCCGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:4548320/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)
                                                                                                                                            ggggTGAGGGCCTCAGCGCGCCCGTGCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:4296388/150‑1 (MQ=255)
                                                                                                                                            ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:2787680/150‑1 (MQ=255)
                                                                                                                                            ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:4599114/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       <  2:3827117/150‑1 (MQ=255)
                                                                                                                                                    ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:679432/150‑1 (MQ=255)
                                                                                                                                                      cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc  <  2:3701742/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGG‑TTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729

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: ‑