Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 580,302 C→T 6.0% R48R (AGG→AGA ylcI ← DUF3950 family protein, DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,3020CT 6.0% 258.6 / 4.6 84R48R (AGG→AGAylcIDUF3950 family protein, DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base C (40/37);  new base T (3/2);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

AGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAG  >  NC_000913/580155‑580448
                                                                                                                                                   |                                                                                                                                                  
aGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGcc                                                                                                                                                   >  5:136115/1‑149 (MQ=255)
aGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATAGATGTATATGCTCTCTTTTCTGACGTTAGca                                                                                                                                                   >  3:12824/1‑148 (MQ=255)
 ggAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCt                                                                                                                                                  >  8:230034/1‑149 (MQ=255)
  gAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGAAATTAAGTGAATATACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACATTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGACTc                                                                                                                                                 >  7:218051/1‑149 (MQ=255)
          tAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATAGTCCTGTTCGGAGGGAAGAACTCGGGATGTTCATTCGTCATCACTTTTAATTGATGTATATGCTCTCTGTTCTGACGTTAGCCTCCGACgg                                                                                                                                           >  7:173960/1‑147 (MQ=255)
                tGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCtt                                                                                                                                    >  7:208746/1‑148 (MQ=255)
                tGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCtt                                                                                                                                    >  5:188701/1‑148 (MQ=255)
                 gccaCTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGACGTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGCCGTTAGCCTCCGACGGCAGTCTTc                                                                                                                                   >  4:187852/4‑148 (MQ=255)
                        aaTTACTGCAATGCCCTCGTAATTATGTGATTTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGAcc                                                                                                                            <  7:65547/148‑1 (MQ=255)
                        aaTTACTGCAATGCCCTCGTAATTATGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGAcc                                                                                                                            <  7:62033/148‑1 (MQ=255)
                              tGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCt                                                                                                                      >  3:16700/1‑148 (MQ=255)
                              tGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTAGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCt                                                                                                                      >  6:89430/1‑148 (MQ=255)
                                cAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTga                                                                                                                    <  8:32104/148‑1 (MQ=255)
                                cAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTga                                                                                                                    <  6:188701/148‑1 (MQ=255)
                                 aaTGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTgag                                                                                                                   <  2:148220/148‑1 (MQ=255)
                                   tGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGaaa                                                                                                                >  7:9864/1‑149 (MQ=255)
                                     cccTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAAt                                                                                                               <  6:116733/148‑1 (MQ=255)
                                      ccTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTc                                                                                                             >  4:100035/1‑149 (MQ=255)
                                      ccTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTc                                                                                                             >  5:140172/1‑149 (MQ=255)
                                       cTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTGTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTcc                                                                                                            <  4:206798/149‑1 (MQ=255)
                                       cTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACGTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTACGAAGGCAGGCTTCAAAGACCCCGGCTGAGAAATTcc                                                                                                            >  4:139950/1‑149 (MQ=255)
                                          gTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCgg                                                                                                         >  1:123616/1‑149 (MQ=255)
                                            aaTTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAc                                                                                                       >  2:46958/1‑149 (MQ=255)
                                             aTTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAcc                                                                                                      <  8:9864/149‑1 (MQ=255)
                                                     aaTTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTg                                                                                               <  7:121275/148‑1 (MQ=255)
                                                       tttACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCt                                                                                             >  2:174427/1‑148 (MQ=255)
                                                       tttACAATATCGTCCTATTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGTCTGAGAAATTCCCGGACCCTTTTTGCt                                                                                             >  8:11557/1‑148 (MQ=255)
                                                                tcctCCTGTTCGGAGCGAACAACGCCCGATCTTCATTCATCATCACAAATAAAAGATGTATATGCTCCCTTATCTGACCATAGCCTCCGACGCCAGGCTTCAATGACCCAGGCTGAGAAAATCCCGGACCCTTTTTGCTCAAGAGCGAt                                                                                   <  4:156093/146‑1 (MQ=18)
                                                                tCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGa                                                                                    <  3:88654/148‑1 (MQ=255)
                                                                tCGTCCTGTGCTGAGGGAAGCTCGCTGGATGTTCATTCTTCTTCACTTTTAATTGATGTATATGATCTCTTTTCTGACGTTAGACTCAGACGGCAGGCTTCAATGACCAAGGCTGAGAAATTCCCGGAACCTTTTTGCTCAAGAGCgt                                                                                    <  8:218051/148‑2 (MQ=25)
                                                                  gTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATg                                                                                  <  8:225896/148‑1 (MQ=255)
                                                                   tCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGCCCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGt                                                                                 >  4:196043/1‑148 (MQ=255)
                                                                   tCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCgg                                                                                                         <  7:126317/124‑1 (MQ=255)
                                                                   tCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCgg                                                                                                         >  8:126317/1‑124 (MQ=255)
                                                                   tCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGt                                                                                 >  2:222505/1‑148 (MQ=255)
                                                                                gaagaaCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAAt                                                                    >  2:37741/1‑148 (MQ=255)
                                                                                   gaaCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATACCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCAtt                                                                >  7:24215/1‑149 (MQ=255)
                                                                                    atcgctGGATGTTCATTCTTCATCACTTTTAATTTATGTTTATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCGCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGtt                                                                        <  2:121580/138‑1 (MQ=25)
                                                                                    aaCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCATCCGACGGAAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCAttt                                                               <  1:118111/149‑1 (MQ=255)
                                                                                    aaCGCGGGATGTTCATTCTTCATCACTTTTAATTGATCTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTAAAGAGCGATGTTAATTTGTTCagatggg                                                                >  1:121580/1‑142 (MQ=255)
                                                                                     aCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTg                                                              <  7:230034/149‑1 (MQ=255)
                                                                                       gcgGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGt                                                            >  5:157259/1‑149 (MQ=255)
                                                                                          ggATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAg                                                         >  4:115803/1‑149 (MQ=255)
                                                                                            aTGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAgg                                                        >  3:160022/1‑148 (MQ=255)
                                                                                             tGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGa                                                       >  2:242310/1‑148 (MQ=255)
                                                                                               ttCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAg                                                    <  1:174427/149‑1 (MQ=255)
                                                                                               ttCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAg                                                    >  2:59512/1‑149 (MQ=255)
                                                                                                  aTTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTTATAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCgg                                                 <  5:89430/149‑1 (MQ=255)
                                                                                                  aTTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCgg                                                 <  7:85001/149‑1 (MQ=255)
                                                                                                  aTTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCg                                                  >  8:85001/1‑148 (MQ=255)
                                                                                                   ttcttcATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGa                                                <  1:242310/149‑1 (MQ=255)
                                                                                                         atcaCTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGTTTAGGAAAGCGGATGTTGc                                          <  1:37741/149‑1 (MQ=255)
                                                                                                          tcaCTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCg                                         >  8:39755/1‑149 (MQ=255)
                                                                                                           caCTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAg                                                         <  2:63963/132‑1 (MQ=255)
                                                                                                           caCTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAg                                                         >  1:63963/1‑132 (MQ=255)
                                                                                                           caCTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCgg                                        >  8:50058/1‑149 (MQ=255)
                                                                                                              ttttAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGgtt                                     >  8:31095/1‑149 (MQ=255)
                                                                                                              ttttAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGAAGGCTTCAATGCCCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTGAGGAAAGCGGATGTTGCGGgtt                                     >  6:28857/1‑149 (MQ=255)
                                                                                                                   aTTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGgttgttgt                                >  5:186967/1‑149 (MQ=255)
                                                                                                                   aTTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGgttgttg                                 <  1:222505/148‑1 (MQ=255)
                                                                                                                   aTTGATGTATATGCTCTCTTTTCTGACGTTAGCCCCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCAGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGgttgttg                                 <  3:230197/148‑1 (MQ=255)
                                                                                                                      gATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTcg                             >  4:195284/1‑148 (MQ=255)
                                                                                                                      gATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCt                             >  1:13655/1‑149 (MQ=255)
                                                                                                                        tGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCAAGGCTGAGAAATTCCCGGACCCTTTTTACTCAAGAGCGATTTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGc                           <  4:12824/149‑1 (MQ=255)
                                                                                                                          tataTGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCCGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCgg                         >  7:139092/1‑149 (MQ=255)
                                                                                                                          tataTGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCgg                         <  3:115803/149‑1 (MQ=255)
                                                                                                                           atatGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCGGAGAAATTCCCGGACCCTTTTTGCTGAAGAGCGATGTTAATTTGTTCAGTCAGTTGGGTAGGAAAGCGGATGTTGCGGGTGGATGTTCGGCggg                        >  6:229540/1‑149 (MQ=17)
                                                                                                                            tatGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGgt                       >  1:120982/1‑149 (MQ=255)
                                                                                                                            tatGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCGTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGgt                       >  2:224911/1‑149 (MQ=255)
                                                                                                                              tGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGgtt                      <  2:120982/148‑1 (MQ=255)
                                                                                                                                ctctctTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGTCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGgttctg                   <  7:11557/149‑1 (MQ=255)
                                                                                                                                 tctctTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGgttctgt                  >  2:222252/1‑149 (MQ=255)
                                                                                                                                     ttttCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGgttctgtt                 <  4:160022/146‑1 (MQ=255)
                                                                                                                                       ttCTGCCGTTAGCATCCGACGGCAGTCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGGTGTTGTTCTGCGGGTGCTGTTCtt              <  8:173960/147‑1 (MQ=11)
                                                                                                                                       ttCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCg                                                                                     >  3:182029/1‑76 (MQ=255)
                                                                                                                                       ttCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCg                                                                                     <  4:182029/76‑1 (MQ=255)
                                                                                                                                       ttCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTa             <  1:59512/148‑1 (MQ=255)
                                                                                                                                       ttCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTa             <  1:222252/148‑1 (MQ=255)
                                                                                                                                         cTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGt           <  8:208746/148‑1 (MQ=18)
                                                                                                                                         cTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGt           <  4:16700/148‑1 (MQ=255)
                                                                                                                                           gACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGa        >  2:229706/1‑149 (MQ=255)
                                                                                                                                            aCGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGAc       <  1:46958/149‑1 (MQ=14)
                                                                                                                                              gTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACAt     <  3:46265/149‑1 (MQ=39)
                                                                                                                                              gTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACAt     <  1:229706/149‑1 (MQ=39)
                                                                                                                                               ttAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGAAAt     >  8:127614/1‑148 (MQ=38)
                                                                                                                                                  gCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAg  <  8:196286/148‑1 (MQ=35)
                                                                                                                                                   |                                                                                                                                                  
AGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAG  >  NC_000913/580155‑580448

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