Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 267,148 Δ1 bp 5.1% intergenic (‑181/+36) ykfK ← / ← yafZ CP4‑6 prophage; putative DNA‑binding transcriptional regulator;Phage or Prophage Related/CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913267,1480G. 5.1% 256.2 / 6.9 78intergenic (‑181/+36)ykfK/yafZCP4‑6 prophage; putative DNA‑binding transcriptional regulator;Phage or Prophage Related/CP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base . (2/2);  total (39/39)
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.20e-01

TTACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGC  >  NC_000913/267010‑267294
                                                                                                                                          |                                                                                                                                                  
ttACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGCGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGata                                                                                                                                          >  3:353573/1‑149 (MQ=255)
  aCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGatat                                                                                                                                         <  3:82753/148‑1 (MQ=255)
   cTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATAtt                                                                                                                                        <  1:29769/148‑1 (MQ=255)
    tttAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCCGCGGGGGTTGCCGGCCTTCCGTCGGTGAGTGAGGAGCGTCAGCACGGTGCGCTGCGGACAGGTGGTGTCCGATATTat                                                                                                                                      >  2:101752/1‑148 (MQ=37)
    tttAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTaa                                                                                                                                      >  2:106418/1‑149 (MQ=255)
    tttAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTaa                                                                                                                                      >  2:100719/1‑149 (MQ=255)
    tttAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTaa                                                                                                                                      >  4:155028/1‑149 (MQ=255)
    tttAAAGAGGGGGTCGGACAGGGTATGGACACAACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTaa                                                                                                                                      >  5:69848/1‑149 (MQ=255)
        aaGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACaa                                                                                                                                   <  1:107766/148‑1 (MQ=255)
         agagGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAag                                                                                                                                 >  8:156549/1‑149 (MQ=255)
          gagGGGGTCGGACAGGGTATGGACACCACAGAAACTTTTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCGTTCGCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGGGAGGTGGTGTGCGATATTTACAAag                                                                                                                                 >  2:197951/1‑148 (MQ=255)
               ggTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCaa                                                                                                                           >  3:388054/1‑149 (MQ=255)
                   ggACAGGGTATGCACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGTACCCTCAGCACTGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACAtt                                                                                                                       <  1:197951/149‑1 (MQ=255)
                      cAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTaaa                                                                                                                    <  5:172383/149‑1 (MQ=255)
                      cAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTaaa                                                                                                                    <  7:156549/149‑1 (MQ=255)
                         ggTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGat                                                                                                                                           <  4:340582/123‑1 (MQ=255)
                         ggTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGat                                                                                                                                           >  3:340582/1‑123 (MQ=255)
                             tGGACACCACAGTAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCa                                                                                                              <  6:287121/148‑1 (MQ=255)
                             tGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAggtggt                                                                                                                                                  >  4:317819/1‑112 (MQ=255)
                             tGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAggtggt                                                                                                                                                  <  3:317819/112‑1 (MQ=255)
                             tGGACACCACAGAAAATATTTTTCATTTCTGCAAAAGCCAGATGGAGCTGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTTGTGTCCGATATTAACAAATCAGCAACATTAAACTTTCAc                                                                                                             <  8:268892/149‑1 (MQ=255)
                              ggACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGt                                                                                                                >  8:164240/1‑145 (MQ=255)
                              ggACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGt                                                                                                                <  7:164240/145‑1 (MQ=255)
                              ggACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACt                                                                                                            >  1:158114/1‑149 (MQ=255)
                               gagagcacaGATACTATTTTTCTTTTCTGCAAAAGGAAGAGTCAGCGGGGGTTACAGGCGTTCCGTCGGTGAGTGAGGACCGTGAGCACGGTGCGATGCGGACAGGTGGTGTCGGATATTAACAAatcag                                                                                                                              <  8:240614/125‑1 (MQ=37)
                               gACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGGGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTctct                                                                                                            <  8:327715/148‑1 (MQ=255)
                               gACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAagcag                                                                                                                              >  7:240614/1‑130 (MQ=255)
                               gACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGAGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTg                                                                                                           >  2:153575/1‑149 (MQ=255)
                                   ccacAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGt                                                                                                       >  7:8166/1‑149 (MQ=255)
                                     acaGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTg                                                                                                     >  1:249570/1‑149 (MQ=255)
                                      caGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGc                                                                                                    >  4:379480/1‑149 (MQ=255)
                                       aGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTcc                                                                                                                                              <  2:304854/106‑1 (MQ=255)
                                       aGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTcc                                                                                                                                              >  1:304854/1‑106 (MQ=255)
                                        gAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGt                                                                                                  <  4:386878/149‑1 (MQ=255)
                                         aaaCTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGt                                                                                                  >  6:154576/1‑148 (MQ=255)
                                          aaCTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGa                                                                                                <  7:344835/149‑1 (MQ=255)
                                             tATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAg                                                                                               >  5:26020/1‑147 (MQ=255)
                                                      tttCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGAAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCg                                                                                    >  5:404544/1‑149 (MQ=255)
                                                        tCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGc                                                                                   <  1:153575/148‑1 (MQ=255)
                                                          tGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGCGCAGGGCTTCCGCCAt                                                                                >  8:191936/1‑149 (MQ=255)
                                                          tGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCAt                                                                                >  5:347310/1‑149 (MQ=255)
                                                          tGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCAt                                                                                >  8:421551/1‑149 (MQ=255)
                                                          tGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCAt                                                                                >  8:415907/1‑149 (MQ=255)
                                                            cAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCa                                                                              <  1:100719/149‑1 (MQ=255)
                                                            cAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGGGTCCGATATTAACAAAGCAGCAACATTAAACGGTCCCTGCAGTTGCGTGAGCAGGGCTTCCGCCATCa                                                                              <  1:106418/149‑1 (MQ=255)
                                                              aaaGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCAcc                                                                            <  7:31972/149‑1 (MQ=255)
                                                                   cAGAGGCAGCGGGGGCTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAg                                                                       >  5:327879/1‑149 (MQ=255)
                                                                      agGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGc                                                                    <  6:69848/149‑1 (MQ=255)
                                                                            cGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGtt                                                              <  6:26020/149‑1 (MQ=255)
                                                                            cGGGGGCTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGtt                                                              <  6:327879/149‑1 (MQ=255)
                                                                                  ttACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAg                                                           <  4:337635/146‑1 (MQ=255)
                                                                                  ttACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAg                                                           >  3:337635/1‑146 (MQ=255)
                                                                                  ttACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGttt                                                        >  6:351713/1‑149 (MQ=255)
                                                                                  ttACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGttt                                                        >  6:388550/1‑149 (MQ=255)
                                                                                    aaaGGCCTTCCCTCGGTGAGTGAGGACCGTCTGCTAGGTGCGATGCGGTCAGGTGGTGTCCGATATTAACAAATGAGCAAGATTTAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCTGCATCACCCACAGTGCGCTGTTAAGTTTc                                                       <  6:404544/146‑1 (MQ=25)
                                                                                     cAGGCCTTCCCTCGGGGAGGGAGGACCGTCAGCACGGTGCGATGCGGACAGGCGGTGTCCGATATTAACAAAGCCGCAACATTAACCGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCAc                                                     <  2:295912/149‑1 (MQ=25)
                                                                                           ttCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTcccc                                               >  4:118361/1‑149 (MQ=255)
                                                                                            tCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCg                                              >  7:7399/1‑149 (MQ=255)
                                                                                            tCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGATAAGTTTCACGTCCCCg                                              >  6:85239/1‑149 (MQ=255)
                                                                                              ccTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTc                                            >  1:29818/1‑149 (MQ=255)
                                                                                                  ggTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGa                                          <  5:388550/147‑1 (MQ=255)
                                                                                                  ggTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGa                                          <  2:249570/147‑1 (MQ=255)
                                                                                                  ggTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATa                                        >  7:268279/1‑149 (MQ=255)
                                                                                                  ggTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATa                                        <  3:118361/149‑1 (MQ=255)
                                                                                                           ggACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGGCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAAc                                <  3:222641/148‑1 (MQ=255)
                                                                                                           ggACCGTCAGCACGGTGCGATGCGGACAGGTGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACgg                              >  2:150198/1‑149 (MQ=255)
                                                                                                           ggACCGTCAGCACGGTGCGATGCGGACAGGTGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACgg                              >  6:281978/1‑149 (MQ=255)
                                                                                                              ccGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCa                            <  8:7399/149‑1 (MQ=255)
                                                                                                                gTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACg                          >  5:83763/1‑149 (MQ=255)
                                                                                                                    gCACGGTGCGATGCGGACAGGTGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTag                     <  5:281978/149‑2 (MQ=255)
                                                                                                                    gCACGGTGCGATGCGGACAGGTGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTAt                     <  1:150198/149‑1 (MQ=255)
                                                                                                                       cGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTAtgt                   <  5:159721/149‑1 (MQ=255)
                                                                                                                        ggTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTAtgtg                  >  7:237356/1‑149 (MQ=255)
                                                                                                                        ggTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTAtgtg                  <  8:268279/149‑1 (MQ=255)
                                                                                                                         gTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGa                 <  2:158114/149‑1 (MQ=255)
                                                                                                                            cGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGAcc               <  3:155028/148‑1 (MQ=255)
                                                                                                                              aTGCGTACAGGTGGGGTCTGATATTAACAAAGCAGCATCATTAAACGGTCACTGCAGTTGCGTGCGCAGGGCTTCGGCGATCACCCACTGTGCGAGGTGAAGTTTGACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGAACGtc            <  4:112551/149‑1 (MQ=255)
                                                                                                                              aTGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGtc            <  6:83084/149‑1 (MQ=255)
                                                                                                                                  ggACAGGTAGTGTCCGATATTAACAAAGCAGCAACCTTAAACGGTCTCTGCAGTTGCGTGCGCAGGGCGTCCGCACTCACCCACAGTGACCGGTTAAGTGTCACGTCCCCGGCGCTACCGCGAACGGCACGCGTATGTGCCCGtcctcc        >  1:133216/1‑149 (MQ=255)
                                                                                                                                     cAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGGCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCttt     <  5:85239/149‑1 (MQ=255)
                                                                                                                                        gtggtgTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTgg   >  3:68643/1‑148 (MQ=255)
                                                                                                                                        gtggtgTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGc  <  5:154576/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
TTACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAGCAACATTAAACGGTCACTGCAGTTGCGTGAGCAGGGCTTCCGCCATCACCCACAGTGCCCGGTTAAGTTTCACGTCCCCGTCGATACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGC  >  NC_000913/267010‑267294

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: