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.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729
                                                                                                                                                       |                                                                                                                                                    
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  2:3683047/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGAAt                                                                                                                                                     >  1:889766/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  2:1373834/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  1:3158360/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACg                                                                                                                                          >  1:4549101/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                         >  2:1957046/1‑150 (MQ=255)
              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                        >  1:3834234/1‑149 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtc                                                                                                                                       >  2:2787680/1‑150 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                       >  2:1799984/1‑148 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgc                                                                                                                                       >  1:3571517/1‑148 (MQ=255)
                 cgatgatTGGGCGTTAACCCAAACATGTCGGGCGGGACGACCCGGAAAATACTGCACCGCCGCGAGGCACTAAAATGTTTTCTTATGGTTGGTTTTCTCCCGCTCCGTGGGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcac                                                                                                                                     <  2:1883129/149‑1 (MQ=255)
                  ttttaatgggCTTTATCCCAAAAATGTCGGGCCTTACGACCGGCAAAATATGCCACCGACGCGAGGGACGAAAATTTATTCGTTTGTTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                    <  2:114620/143‑1 (MQ=255)
                  tatcaaaGGGCTTTAACCCAAACATTACTGGCCGCACCACCCGCAAAAAACTCCACCCACGCGACGCACGAAAAATTATTCCTAAGTTTGTTTTTCTCCCACACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                    <  2:1994936/146‑1 (MQ=19)
                  gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGtgcccc                                                                                                                                    >  2:4102342/1‑145 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccccaaccc                                                                                                                            >  2:87242/1‑149 (MQ=255)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCCaa                                                                                                                          >  2:646505/1‑150 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:720442/150‑1 (MQ=255)
                               tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCAAAacc                                                                                                                       >  1:771054/1‑150 (MQ=38)
                                aTCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACGGACGCGACGGACGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                      <  2:3898376/150‑1 (MQ=38)
                                 tCCCAACCATGTCTGGCAGGAGTACGCGTAATTTAGTGCACAGTCGCGACTGACGAAATTGTATTCGTATGGTTGTTTTTCTCGCTATCACTGTGGTAGTGGGCAGGGGTGTGTGAATCAGCGCGCACGTTCACCCTCACCCTaaaactc                                                                                                                     <  2:3568887/150‑6 (MQ=255)
                                 tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                     <  2:3568958/150‑1 (MQ=255)
                                  cccAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                    <  2:564081/150‑1 (MQ=255)
                                  cccAAACTTGTGGGGCCGGACGACCCGCAAAACACGGCCCCGACGCGACGGACGAAAATGTATTCGTATGGTTGGTTTTCCCCCCCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                    <  1:3987728/150‑1 (MQ=37)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCAAAcccacac                                                                                                                   >  1:2308812/1‑146 (MQ=38)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCCCCCCCCCCCCCAcccacacc                                                                                                                  >  1:3362530/1‑145 (MQ=37)
                                          ggccggcccGGACGCCCCGCCAACTCCGGCACCGCCCCGACGGACCACCAAGTATTCGCACGGTGGGTTGGCGCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:2739614/143‑1 (MQ=255)
                                          ggTCGGGCCGGACCCCCCGCCAAATACGGCACCGAGGCGCCGCACGAAAATTTATTGGTATGGTTGTTTTTCGCCCACCCCGGGGGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  2:1729422/149‑1 (MQ=255)
                                            tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCCCCCGAAACCACTcccgcccggg                                                                                                          >  1:277500/1‑143 (MQ=37)
                                                gccggccgCCCCGGAAAATCCTCCACCGACGCGACGGTCTAAAATTTTTTGTTATGTTTGGTTTTCTCCCTCGGCGTGTGGGTGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  1:2634351/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCACCCCCAAACCCCCTCCCACCAggggggg                                                                                                      >  2:3926941/1‑146 (MQ=37)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTGCCCCCGCACCCCAACCCACACCCCCAAggggggg                                                                                                      >  2:2954720/1‑146 (MQ=37)
                                                 ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAACCCACACCCACAACGGGGAgg                                                                                                     >  2:4101359/1‑150 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCACCCCCACCCAAAACCACACCCACCAGGGAGgggg                                                                                                    >  1:1480559/1‑150 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                    <  2:4723510/150‑1 (MQ=255)
                                                  cGGACGCCCGGCAAAATCCGCCCCGCCGGGGCGGCCCTAAAAATTTTGCTTAGGTTGGTTTTTCACCCACCCCGTGGGGTCGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                    <  2:1272845/150‑1 (MQ=255)
                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCGCCACACCCAAAACCACAACCACAAGGGCGAgggg                                                                                                   >  1:1426749/1‑150 (MQ=25)
                                                          ccGCCAAATCCGGCCCCGCCGCGACGGACGAAAACGTTTTCGTAAGGTTGGTTTTCTCCCCCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:2521740/150‑1 (MQ=37)
                                                            ggCATAAACGGCACCGTCGCGCCTGACGAAATTGTATTCGTACGTTGTTTTTTCACCTTCTCGAGGGGTGAGTGGGCGGGGTTGAGGGCATCAGCGCGCAGGTTCACCCTCACCCTAACCCTCTCCCTCAATGGAGAGGGGACCGTTCGa                                                                                          <  2:446001/150‑1 (MQ=255)
                                                            gccAAATCCTGCCCCGACGCGCCGGACCAAAATGTATCCTCACGTTTGGTTGACCCCCTCCCCGTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:4059325/148‑1 (MQ=255)
                                                            gcaaaaTACTGCCCCGACGGGCGGGACTAAAATGTTTTCGTATGTTTGGTTTTCGCCGTGTCCCTGTGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:4058975/147‑1 (MQ=255)
                                                             gCAAATACTGCCCCGCCGCGCCGGACGAAAACTTATTCGTATGGTTGGTTTTCACCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  2:2512286/150‑1 (MQ=255)
                                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCCCACCCGAAACCACGCCCTCAAGGGGGAGGGGGCCGGTCGAg                                                                                         >  1:746039/1‑150 (MQ=25)
                                                             ccAAAGACTGCCCCGACGGGCCGGACCAAAAAGTATTCGTACGGTGGGTTTTCGCCCTCTCCCTGGGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                         <  1:2755647/149‑1 (MQ=255)
                                                              cAAATCCTCCACCGAGGCGCCGGACCAAAAAGTATTCGAACGTTTGTTTTACCCCCACTCACTGGGGGAGCGGGCGGGGGTTTGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                        <  1:342522/150‑1 (MQ=255)
                                                               aaaTCCGGCACCGCGGCGCCGGTTGAAAATGTTTTGGTCTGGTTGGTTTTCGCCCTCGCCGGGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                       <  1:3458360/150‑1 (MQ=255)
                                                                       cACCGCCGCTACGCACAAAAAATTATTCTTATGGTTTTTTTTCCCCCTCTCCGTGTGTTAGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  1:3262724/150‑1 (MQ=255)
                                                                         ccgccgCGACGGCCCAAAATTTTTTCTTATGGTTGGTTTTCCCCCCCCCCGTGGGGGTGTGGGGGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:1469459/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGACCAAAATTTATTCGTACGTTTGGTTTTCTCCCACTCCCTGTGGTTGCGGGCGGGGTTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:3834234/150‑1 (MQ=255)
                                                                          cgacgcgacgTCCTATATATTATTCTTAGGCTTGGTTTTCGCCCGCGCCCTTTGGTTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                            <  2:1473763/150‑1 (MQ=255)
                                                                           gacgcgacgGACTAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                           <  2:777288/150‑1 (MQ=255)
                                                                            acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCCAAACCCCACCCACACCGGGGAGGGGGCCGCGACAGCACCAAATTTgagtg                                                                          >  2:3052020/1‑148 (MQ=21)
                                                                                   gtaggAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCGGGGGTTAGGGCCTCAGCGCGCACGTTCACCCTCACGCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGTGATAtctt                                                                   <  2:2737265/146‑1 (MQ=255)
                                                                                       gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCCCAAGGGAGAGGGGACCGTTCGAGCACAAATTTTGAGAGATAtcttcttc                                                               <  2:3912934/150‑1 (MQ=255)
                                                                                       gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACACGAAAACACTAACGCCAGGGAGAGGGGGCCCATCGAGCACAAAATTTGTGAGAAAtattctta                                                               >  2:747288/1‑149 (MQ=21)
                                                                                        aaataGTATCCGTTCTGTTGTTTTTCCCCCTCACCCTTTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTAACCCTAACCCTATCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTGAGAGATAtcttcttat                                                              <  2:1077360/150‑3 (MQ=255)
                                                                                        aaaaaGTATTCGTACGGTTGGTTTACTCCCTCACCATGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                              <  2:1075304/150‑1 (MQ=255)
                                                                                           aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCGCACCCAAAACCACTACCACAACGGGGAGGGGGCCCATCGAGCACAAAATTTGTGAGATATATTATTCTGtg                                                           >  2:4284166/1‑149 (MQ=21)
                                                                                             ttATTCGTAAGGTTGGTCTGCCCCCCCCCCCGGTGGGTGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  1:428416/149‑1 (MQ=255)
                                                                                             gTATTCGTATGGTTGTTTTTCTCCATCTCCCTGTGGGATTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACTAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:380627/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  2:2819206/150‑1 (MQ=255)
                                                                                              ttttggTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  2:2370489/145‑1 (MQ=255)
                                                                                              ttttCGTATGTTTGGTTTTCTCCCTCTCCCGGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  1:923438/148‑1 (MQ=255)
                                                                                              tttgcGTAAGGTTTGTTTTCTCCCTCTCCCTTTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  2:2292383/146‑1 (MQ=39)
                                                                                              tATTCGTTTGGTTGTTTTTCTCCCGCTCCGTGGGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  1:923419/150‑1 (MQ=255)
                                                                                                 tggTATGGTTGGTTTCCGCCCCCGCCCTGGGGGTGCGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                     <  2:1816158/148‑1 (MQ=255)
                                                                                                 tCGTATGGTTGGTTTTCTCCCTCGCCCGGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                     <  2:4296628/150‑1 (MQ=39)
                                                                                                   gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                   <  2:3158360/150‑1 (MQ=255)
                                                                                                     taGGTTGTTTTTCTCCCTCCCCGTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  1:2637701/149‑1 (MQ=255)
                                                                                                     atgtttggttTTCGCCCTCCCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  1:2637504/146‑1 (MQ=255)
                                                                                                     aaGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  2:363033/150‑1 (MQ=255)
                                                                                                         ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                             <  2:3049672/150‑1 (MQ=255)
                                                                                                          tggttTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                            <  1:756927/150‑1 (MQ=255)
                                                                                                          tggttTACTCCCTCTCCCTGGGGGCGTGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                            <  1:2297629/150‑1 (MQ=255)
                                                                                                           tgtttttCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                            <  2:1867597/147‑1 (MQ=255)
                                                                                                            gttTTCCCCCCCCCCCGGGGGGTAGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGACGAACa                                           <  1:1281904/150‑1 (MQ=255)
                                                                                                              ttttCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                         <  1:1174195/150‑1 (MQ=255)
                                                                                                              ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                         <  1:1174443/150‑1 (MQ=255)
                                                                                                                  ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                     <  2:2028665/150‑1 (MQ=255)
                                                                                                                   tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                    <  2:3125609/150‑1 (MQ=255)
                                                                                                                   tccctctcccgGTGGGAGAGGGCGGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCACTCAAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATATCTTCTTCTGTATTGTAACAGAAGAACAGAATATc                                    <  2:2394215/150‑1 (MQ=255)
                                                                                                                               tgGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  2:4174898/150‑1 (MQ=255)
                                                                                                                               ggggAGTGGGCGGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTTCCCCTCTCCCTCAAGGGAGAGGGGACCGATCGCGCACAAATTTTGAGAGATATCTTCTTCTGCCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  2:1879098/149‑1 (MQ=255)
                                                                                                                               ggggAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  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)
                                                                                                                                                       |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  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: