Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,648 T→A 12.8% intergenic (‑262/+65) sufA ← / ← rydB Fe‑S cluster assembly protein/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,6480TA12.8% 239.8 / 5.3 74intergenic (‑262/+65)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (43/21);  new base A (2/8);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCGTATAACAGCATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGG  >  NC_000913/1764410‑1764883
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
tCAGCGTATAACAGCATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGt                                                                                                                                                                                                                                              >  1:467131/1‑239 (MQ=255)
 cAGCGTATAACAGCATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTa                                                                                                                                                                                                                                             <  1:257965/239‑1 (MQ=255)
              cATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCg                                                                                                                                                                                                                               >  1:1209306/1‑239 (MQ=255)
              cATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCg                                                                                                                                                                                                                               >  2:1178236/1‑239 (MQ=255)
                tGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGtt                                                                                                                                                                                                                             <  1:972545/239‑1 (MQ=255)
                tGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGt                                                                                                                                                                                                                              >  1:1199533/1‑238 (MQ=255)
                tGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATAAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGtt                                                                                                                                                                                                                             <  1:920266/239‑1 (MQ=255)
                   tAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAg                                                                                                                                                                                                                          >  2:593038/1‑239 (MQ=255)
                    aGTGATAATGATTATCAGTTCAACCCAGCAACCGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGa                                                                                                                                                                                                                         >  1:772928/1‑239 (MQ=255)
                        aTAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGa                                                                                                                                                                                                                         <  1:714797/235‑1 (MQ=255)
                        aTAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGa                                                                                                                                                                                                                         >  2:714797/1‑235 (MQ=255)
                             gagtATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTACAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGaaa                                                                                                                                                                                                                <  2:218656/236‑1 (MQ=255)
                             gATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGaaa                                                                                                                                                                                                                <  2:772928/239‑1 (MQ=255)
                                      ttCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGa                                                                                                                                                                                                                         <  2:1208056/221‑1 (MQ=255)
                                      ttCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGa                                                                                                                                                                                                                         >  1:1208056/1‑221 (MQ=255)
                                           cccAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCtgt                                                                                                                                                                                                  <  2:467131/239‑1 (MQ=255)
                                           cccAGCAAACGCAGGAGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGAAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCATTAAGTCtgt                                                                                                                                                                                                  >  2:1180907/1‑239 (MQ=255)
                                              aGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCa                                                                                                                                                                                                                                   <  1:353174/204‑1 (MQ=255)
                                              aGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCa                                                                                                                                                                                                                                   >  2:353174/1‑204 (MQ=255)
                                                cAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGac                                                                                                                                                                                             >  1:1119622/1‑239 (MQ=255)
                                                               tACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATaa                                                                                                                                                                                                                    >  2:715684/1‑201 (MQ=255)
                                                               tACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATaa                                                                                                                                                                                                                    <  1:715684/201‑1 (MQ=255)
                                                                   gTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCt                                                                                                                                                                           >  1:261981/1‑238 (MQ=255)
                                                                            gCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCg                                                                                                                                                                 <  1:115967/239‑1 (MQ=255)
                                                                                      tGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAtct                                                                                                                                                       <  1:593038/239‑1 (MQ=255)
                                                                                       gCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAtct                                                                                                                                                       >  1:592219/1‑238 (MQ=255)
                                                                                        cTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTcc                                                                                                                                                     >  1:657368/1‑239 (MQ=255)
                                                                                                             gATGTGTGGGTTAACAGGCCCGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGAAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGAGCTTATCTCCGGCACTCTCAGTGGCTTAGc                                                                                                                                 >  1:151267/1‑238 (MQ=255)
                                                                                                                  gtgGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGa                                                                                                                           >  1:707015/1‑239 (MQ=255)
                                                                                                                  gtgGGTTAACAGGCACGCTAAGAATTAGGGATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGa                                                                                                                           >  1:3249/1‑239 (MQ=255)
                                                                                                                    gggTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                                                                                                                         >  1:841726/1‑239 (MQ=255)
                                                                                                                           cAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGAt                                                                                                                  >  2:374688/1‑239 (MQ=255)
                                                                                                                                    tAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa                                                                                                          <  2:592219/238‑1 (MQ=255)
                                                                                                                                         aTTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCa                                                                                                     <  1:530194/238‑1 (MQ=255)
                                                                                                                                                  ctcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGATGGGGCGATAAGAATAATCTCATAAAGCTaa                                                                                            >  2:861087/1‑238 (MQ=255)
                                                                                                                                                  ctcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTaa                                                                                            >  1:1012491/1‑238 (MQ=255)
                                                                                                                                                       tGTTAGATAATGGTTATTATCTAAGAAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                                                                      >  1:86467/1‑239 (MQ=255)
                                                                                                                                                           aGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGttt                                                                                  >  1:223792/1‑239 (MQ=255)
                                                                                                                                                                 tGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAacac                                                                            >  2:388411/1‑239 (MQ=255)
                                                                                                                                                                        tatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata                                                                                                           >  2:194810/1‑201 (MQ=255)
                                                                                                                                                                        tatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata                                                                                                           <  1:194810/201‑1 (MQ=255)
                                                                                                                                                                             aaGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGAt                                                                                                                  >  2:928864/1‑189 (MQ=255)
                                                                                                                                                                             aaGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGAt                                                                                                                  <  1:928864/189‑1 (MQ=255)
                                                                                                                                                                             aaGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTa                                                                >  2:534299/1‑239 (MQ=255)
                                                                                                                                                                             aaGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACACCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTa                                                                >  1:236/1‑239 (MQ=255)
                                                                                                                                                                               ggAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAt                                                                                                                                                         >  2:652431/1‑148 (MQ=255)
                                                                                                                                                                               ggAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAt                                                                                                                                                         <  1:652431/148‑1 (MQ=255)
                                                                                                                                                                                   cATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtatta                                                          >  1:306587/1‑239 (MQ=255)
                                                                                                                                                                                    aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat                                                         >  1:173793/1‑239 (MQ=255)
                                                                                                                                                                                    aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat                                                         >  1:554824/1‑239 (MQ=255)
                                                                                                                                                                                    aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat                                                         >  1:567609/1‑239 (MQ=255)
                                                                                                                                                                                      aaaGGCTGTGATGAAAAAGAAAAAGTTGTATTAAAAATAACTATTAAAGAGATAGGTAAAAAAGTGAAGCGTTCAGAAAAAAGAAACCCCTTAAGTCTGTGCGACACAGGCTTAAGGTTTCCTACCCCATCCGGCGTTTATCTCCGGCCCTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt                                                       <  1:1178236/238‑1 (MQ=255)
                                                                                                                                                                                           atgtaaaaaaaaagaaaaagTTGTATAGAAAATGAATATTAAAGAGAGAGGAAAAAAAGGACAGCTCAAAAAAAAAAGAAACCCCTAAAGTCGGGTCAAACCCGCCTTAGGTTTTTCCCCCCCCCCGGGCTCTTTTCCCCGGCACTCTCAGTGCCTTACCTTTAGAGGGGGTAATAAAAATAATCAAAAAAACCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                    <  2:554824/229‑1 (MQ=255)
                                                                                                                                                                                            tgaaaaaaaaaagaaaaagTGTAATGAAAAATAATAATTAAAGAGAAAGAAAAAAAAGGGCGGCGTTAAAAAAAAAAAAACCCCTAAAGTCTGGGCGACACAGGCTTAAGGTTTTCTCCCCCACCCGCCGCTAACCCCCGGCACTCTCAGTGGTTTAGCTTTGGAAGGGGCGATAAGAATAATCTAATAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                  <  2:223792/231‑1 (MQ=255)
                                                                                                                                                                                            ggtaataaaaaaaaaaaagTTGTTTAAAAAATGATTATTAAAAAGATAGATAAAAAATGCCCCTATTAAAAAAAAAAAAACCCCTTAAGTCTGGGCGCAACAGGCTTAAGGTTTTCCCCCCCCGCCGGCGTTTCTCTCCGCCCCTCTGAGTGCTTTGTCTTTTGAAGGGGCAATAAAAATAATCTCAAAAAGTCACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                  <  2:173793/231‑1 (MQ=255)
                                                                                                                                                                                             gtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATACAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGCATAAtttttacaca                                                >  2:901599/1‑234 (MQ=255)
                                                                                                                                                                                             gtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa                                                >  1:283618/1‑239 (MQ=255)
                                                                                                                                                                                              tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                         >  1:1021342/1‑181 (MQ=255)
                                                                                                                                                                                              tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                         <  2:1021342/181‑1 (MQ=255)
                                                                                                                                                                                              tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACaa                                               >  2:782613/1‑239 (MQ=255)
                                                                                                                                                                                              tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACaa                                               >  1:173976/1‑239 (MQ=255)
                                                                                                                                                                                                     aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAg                                         >  1:767235/1‑238 (MQ=255)
                                                                                                                                                                                                      aagaaaaatTTGGATTAAAAAGGACTATTTAAGAGAAAGCAAAAAAAGGACAGCTTTAAGAAAAAAAAACCCCCTTAAGTCTGTGCACCACAGCCTTAAGGGTCTCTGCCCCATCCGGGCCTAATCGCCGACCCTCTCGGGGGCTGAGATCTTGAGGGGGTAATAAGAAGATTTCAAAAAAACTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                       <  2:1209306/239‑1 (MQ=255)
                                                                                                                                                                                                      aaaaaaaagTTGGATTAAAAATGACTATTTAAGAGATACAAAAAAAAGTCCATGGTACAGAAAAAAGAACACCCTTAAGCCTGCGCGACCCAGTCTTAAGGGTTTCCCCCCCACCCCGCGCTAACCTCCGGCCCTCTCGGTGGCTTAGCTCTTGAAGGGGCAAAAAGAATAATCTAATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                       <  2:657368/236‑1 (MQ=255)
                                                                                                                                                                                                      aaaaaaaagTTATATGAAAAAAGATAATTAAAGACATAGATAAAAAAGTGGGGCCTTAAAAAAAAAGAAACCCCTTAGTTCTGCGGCACACGGTCTAAAGGTTTCCACCCCCCTCCGGCGCTTAGCCCCGCCCCTCTCTGTGGAGCAGCCTTTAAAGGGGATAAAAGAATAATCTAAAAAAGAAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                       <  2:3249/236‑1 (MQ=255)
                                                                                                                                                                                                        aaaaaagTGGTAATAAAAATAATTATGAAAGAAAGAGGAAAAAAAGTGCAGGGTTAAAAAAAAAAAACACCTTTAAGTCTGTGCGACAAAGTCTAAGGGGTTTCTACCCCACCCGGCGATTATCCCCGCCCTCCTAGGTGGTTTAGCTCTAGAAGGGGGAAAAAGAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCgg                                     <  2:385526/238‑1 (MQ=255)
                                                                                                                                                                                                                tGTATTGAAAATGATTATTTAAGAGATAGGAAAAAAAGTGCAGCGTTAAGAAATAAGAAACCCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTCCCCCATCCGGCGCTTATCTCCGGCACTCTAAGTGGCTTAGCTCTTGAAGGGGCGATAAGAAAAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAg                             <  1:193351/239‑1 (MQ=255)
                                                                                                                                                                                                                  tATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCa                           >  1:568225/1‑239 (MQ=255)
                                                                                                                                                                                                                                   taaaGAGACAGAAAAAAAGGGGAGGCATTAAGAAAAAAAAAACCCTTTATGTCGGGCCGACCCGGGCTAAAGGTTTTCACCCCCCCCCGACTTTTATCCCCGCCACTCTCGGTGCTTTACTTTTTGAAGGGACAATAAAAATAATCACATAAATCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCAt           <  2:1315861/237‑1 (MQ=255)
                                                                                                                                                                                                                                     aaGAGAAAGAAAAAAAAGGGCAGCTTTAAAAAATAAGAAACCCTAGAAGTCTGTGGGACACAGGCTTAAGGGTTTCTACCCCACCCGGCTTTTATCTCCGGCACTCGAGGGGGTTTAGCTTTTGAAGGGGAAAAAAGAATAATCTAATAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCg         <  2:1057336/239‑1 (MQ=255)
                                                                                                                                                                                                                                          aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGaa       >  2:159109/1‑236 (MQ=255)
                                                                                                                                                                                                                                          aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGaa       >  1:1216728/1‑236 (MQ=255)
                                                                                                                                                                                                                                            aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTgg  >  1:316585/1‑239 (MQ=255)
                                                                                                                                                                                                                                            aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTgg  >  1:1189700/1‑239 (MQ=255)
                                                                                                                                                                                                                                             ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTgg  >  1:1057336/1‑238 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
TCAGCGTATAACAGCATGTTAGTGATAATGATTATCAGTTCAACCCAGCAAACGCAGGGGCTTTACCGTAAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGG  >  NC_000913/1764410‑1764883

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