Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,417,011 C→T 5.5% E126K (GAA→AAA)  recE ← Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,417,0110CT 5.5% 225.7 / 2.3 73E126K (GAA→AAA) recERac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base T (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

AGGCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGA  >  NC_000913/1416866‑1417153
                                                                                                                                                 |                                                                                                                                              
aGGCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTACTGCCTGGAGTCCGTCGAGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTTAGCAAGCCAACGAATGGGCAGTTCCTg                                                                                                                                             <  3:368620/149‑1 (MQ=255)
aGGCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTg                                                                                                                                             <  7:42272/149‑1 (MQ=255)
aGGCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTg                                                                                                                                             <  7:157480/149‑1 (MQ=255)
 ggCTGGTTTTGTCGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTgg                                                                                                                                            >  4:463124/1‑149 (MQ=255)
  gCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGc                                                                                                                                           >  7:272346/1‑149 (MQ=255)
     ggTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTg                                                                                                                                             <  8:153809/144‑1 (MQ=255)
     ggTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTg                                                                                                                                             >  7:153809/1‑144 (MQ=255)
              ggAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTgg                                                                                                                                 >  6:385602/1‑147 (MQ=255)
              ggAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTgg                                                                                                                                 <  5:385602/147‑1 (MQ=255)
              ggAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGa                                                                                                                               <  3:426845/149‑1 (MQ=255)
                 aaCAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTa                                                                                                                            <  1:444821/149‑1 (MQ=255)
                               tAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGCTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTctcctc                                                                                                               >  3:65651/1‑148 (MQ=255)
                                    gCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAAtctc                                                                                                         <  5:177047/149‑1 (MQ=255)
                                     ccGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCa                                                                                                        <  7:441723/149‑1 (MQ=255)
                                       ggTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTgg                                                                                                                                 >  8:258616/1‑122 (MQ=255)
                                       ggTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTgg                                                                                                                                 <  7:258616/122‑1 (MQ=255)
                                        gTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGt                                                                                                      >  2:28145/1‑148 (MQ=255)
                                        gTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTc                                                                                                     <  2:292327/149‑1 (MQ=255)
                                        gTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTCCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGt                                                                                                      >  3:29973/1‑148 (MQ=255)
                                             tcttcAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCAt                                                                                                   <  7:234568/146‑1 (MQ=255)
                                             tcttcAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCAt                                                                                                   >  8:234568/1‑146 (MQ=255)
                                             tcttcAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTc                                                                                                >  1:205948/1‑149 (MQ=255)
                                               ttcAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCACGCCAACGAATGGGCAGTTACTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCtt                                                                                              >  8:481654/1‑149 (MQ=255)
                                                   gCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGcc                                                                                          >  2:440637/1‑149 (MQ=255)
                                                     ccGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGt                                                                                        <  4:65651/149‑1 (MQ=255)
                                                      cGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGtt                                                                                       >  3:136086/1‑149 (MQ=255)
                                                      cGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGtt                                                                                       >  4:34750/1‑149 (MQ=255)
                                                            aTGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTACTCAATCTCAGTCATGTCTTCGCCGTTGACGtt                                                                                 >  7:504091/1‑149 (MQ=255)
                                                               tgtAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTgg                                                                               >  1:262019/1‑148 (MQ=255)
                                                                 tAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTa                                                                             <  8:272346/148‑1 (MQ=255)
                                                                   aTGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTg                                                                          <  7:13858/149‑1 (MQ=255)
                                                                     gCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGcc                                                                        >  6:113869/1‑149 (MQ=255)
                                                                             tCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCTATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGt                                                                <  3:463124/149‑1 (MQ=255)
                                                                               cGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGcc                                                                        <  1:47907/139‑1 (MQ=255)
                                                                               cGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGcc                                                                        >  2:47907/1‑139 (MQ=255)
                                                                               cGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGt                                                                >  1:363166/1‑147 (MQ=255)
                                                                                      tGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGt                                                       <  3:34750/149‑1 (MQ=255)
                                                                                            gTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGtg                                                 <  2:205948/149‑1 (MQ=255)
                                                                                             tGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCGTGGCCAGAAATTGGGAGTAGCATTTTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGtgc                                                <  1:141974/149‑1 (MQ=255)
                                                                                              gCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGtgct                                               >  2:142270/1‑149 (MQ=255)
                                                                                              gCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGtgct                                               >  6:187159/1‑149 (MQ=255)
                                                                                                gTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGTAGTAGCATATTCTCCTCAATCTCATTCATTTCTTCTCCGTTGACGTTGGTATTGACTTGATATTGAGCGTTGTCTGGtgatgc                                             >  4:37767/1‑149 (MQ=255)
                                                                                                gTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGt                                                       >  6:76662/1‑139 (MQ=255)
                                                                                                gTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGt                                                       <  5:76662/139‑1 (MQ=255)
                                                                                                gTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGtgctgc                                             <  2:363166/149‑1 (MQ=255)
                                                                                                   aCCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTTCTGGCCAGAAATTGGGa                                                                                                                               <  3:115072/64‑1 (MQ=255)
                                                                                                   aCCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTTCTGGCCAGAAATTGGGa                                                                                                                               <  3:108840/64‑1 (MQ=255)
                                                                                                   aCCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTTCTGGCCAGAAATTGGGa                                                                                                                               >  4:115072/1‑64 (MQ=255)
                                                                                                   aCCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTTCTGGCCAGAAATTGGGa                                                                                                                               >  4:108840/1‑64 (MQ=255)
                                                                                                    ccGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTcc                                          <  1:266394/148‑1 (MQ=255)
                                                                                                      gggttttCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCATGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCgg                                       >  5:141337/4‑149 (MQ=255)
                                                                                                            tCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTa                                  <  4:99564/148‑1 (MQ=255)
                                                                                                                 gCCGTGTTGAGCATGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAACATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTcc                             >  7:5970/1‑148 (MQ=255)
                                                                                                                  ccGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTccc                            <  4:84407/148‑1 (MQ=255)
                                                                                                                      gtTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGt                        >  1:277209/1‑148 (MQ=255)
                                                                                                                      gtTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGGCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTGTAGTTCCCATggc                       <  3:37767/149‑3 (MQ=255)
                                                                                                                          aGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATg                    <  2:277209/148‑1 (MQ=255)
                                                                                                                           gCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATg                    >  7:374961/1‑147 (MQ=255)
                                                                                                                           gCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGa                  <  1:28145/149‑1 (MQ=255)
                                                                                                                                ccAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCtt             <  4:136086/149‑1 (MQ=255)
                                                                                                                                ccAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCtt             >  8:454591/1‑149 (MQ=255)
                                                                                                                                  aaCGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTg           >  5:410879/1‑149 (MQ=255)
                                                                                                                                  aaCGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATCGAGTCTTTg           >  4:316629/1‑149 (MQ=255)
                                                                                                                                   aCGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGc          <  1:440637/149‑1 (MQ=255)
                                                                                                                                     gAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTGATGGAGTCTTTGCTg        >  5:199609/1‑149 (MQ=255)
                                                                                                                                     gAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCt         >  2:86236/1‑148 (MQ=255)
                                                                                                                                       aTGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAg      <  8:374961/149‑1 (MQ=255)
                                                                                                                                          ggCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGtt    <  5:113869/148‑1 (MQ=255)
                                                                                                                                           gCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGa                  <  7:22542/133‑1 (MQ=255)
                                                                                                                                           gCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGa                  >  8:22542/1‑133 (MQ=255)
                                                                                                                                           gCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGa  >  2:457200/1‑149 (MQ=255)
                                                                                                                                              gTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGtt    >  7:477553/1‑144 (MQ=255)
                                                                                                                                              gTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGtt    <  8:477553/144‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AGGCTGGTTTTGTGGGAAACAGCCAGGGCAGTAACAGCCGGTAGTTCTTCAGCCCGAGCAATGTGTAATGCCTGGAGTCCGTCGCGTGAAACGTGCGTTACCGGTTTTTCGCTGCCGTGTTGAGCAAGCCAACGAATGGGCAGTTCCTGGCCAGAAATTGGGAGTAGCATATTCTCCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGA  >  NC_000913/1416866‑1417153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: