Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,364,798 C→T G82D (GGC→GAC)  dsbD ← thiol‑disulfide exchange protein DsbD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,364,7980CT94.6% 165.6 / NA 57G82D (GGC→GAC) dsbDthiol‑disulfide exchange protein DsbD
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (18/35);  minor base G (1/2);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGC  >  NC_000913/4364561‑4365032
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
gcgGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTg                                                                                                                                                                                                                               <  1:437865/251‑1 (MQ=255)
   gTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCAc                                                                                                                                                                                                                            <  1:407371/251‑1 (MQ=255)
      ggCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTc                                                                                                                                                                                                                         <  2:120856/251‑1 (MQ=255)
          gctCTTGCTGCGGAACAGACACAGGCTGTGGTGCTGCGTTGTGGGCGACCACTACGCTTAACGGAACGGTTTTGGTGTCTGGCGGATAACAGAAACCGGCATCAGTTCAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCCCTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGc                                                                                                                                                                                                                      <  2:352021/250‑1 (MQ=255)
           ctctTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGCGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGcc                                                                                                                                                                                                                     >  1:223436/1‑249 (MQ=255)
            tctTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCCACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACC‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGCCGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAAt                                                                                                                                                                                                                   >  1:283392/1‑251 (MQ=255)
                   gcGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTcgcg                                                                                                                                                                                                            <  2:255630/251‑1 (MQ=255)
                    cGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCTCgt                                                                                                                                                                                                           >  2:64483/1‑251 (MQ=255)
                          aGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCg                                                                                                                                                                                                     <  2:218482/251‑1 (MQ=255)
                               caGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTa                                                                                                                                                                                                <  1:179775/251‑1 (MQ=255)
                                       ggCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGAt                                                                                                                                                                                        >  1:368582/1‑251 (MQ=255)
                                          gcTGCGTTGTTGGCGACCACTTCGCTTAAGGGCACGGTTTTGGTTTCTGGCGGATCACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTg                                                                                                                                                                                     >  2:31702/1‑251 (MQ=255)
                                          gcTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTg                                                                                                                                                                                     <  1:404971/251‑1 (MQ=255)
                                           cTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGt                                                                                                                                                                                    >  1:348276/1‑251 (MQ=255)
                                             gCGTTGTTGGCGCCCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGAGCTGtt                                                                                                                                                                                   <  2:225565/250‑1 (MQ=255)
                                                  gttgGCGACCACTTCGCTTAACGGAACGTTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCCGCTGCACGTAGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACggg                                                                                                                                                                             >  2:19040/1‑250 (MQ=255)
                                                       cGACCACTTAGCTTACCGGAACGGTTTTGGTTTCTGGCGGATAACAGAACCCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAACTCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAgg                                                                                                                                                                        >  2:155918/1‑251 (MQ=255)
                                                            aCTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGTCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGgtagta                                                                                                                                                                   <  1:284797/251‑1 (MQ=255)
                                                              ttCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGAGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGCTCTGTTTACGGTAGAGGTAGTAAc                                                                                                                                                                 >  1:406713/1‑251 (MQ=255)
                                                               tctcTTAACGGAACGGTTTGGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GGTAACGTCGCTCCCGCACTCGCCTGGTTGATGGGGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGGCTTGCGGCAGCTGCACGTCGGCAATGTTCGCGTGTTCCGGCGTAAGGCGTATCTGTTTACGGTAGAGGTAGTAAc                                                                                                                                                                 <  2:272899/247‑1 (MQ=255)
                                                                cgctttaccgacgggTTTTGGTTGCTGGCGGATAACAGAACCCGGCAGGAGCACAGCCCTGCTAGGTG‑CACGTTAACGTCGCTCCCGGACTCGCCTGGTTGCTGGTTACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCGCGTCGGCAATTTTCGCGTGTTCCGGGGTAATGCTGATCTGTTGCCGGTAGAGgttgtcacc                                                                                                                                                                <  2:125528/237‑8 (MQ=255)
                                                                      acggaacggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTg                                                                                                                                                         <  2:213330/251‑1 (MQ=255)
                                                                       cggaacgtTGTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTTACA‑GTTAACGTCTCTCCCGCTCTCGCCTGGTTGATGGTGACGTGAACCGTCAGCCGATCGCTGTAATTCTCGCTTTTGCCGTAAAACTCATCTTCAGGCCAGCCGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTTACCGTCTTTGa                                                                                                                                                        <  2:2152/251‑1 (MQ=255)
                                                                        ggaacggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGAt                                                                                                                                                       <  1:289050/251‑1 (MQ=255)
                                                                            cggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCCGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTCCGGTAGCGGTAGTAACCGGCTTTGATCTGc                                                                                                                                                   >  2:237260/1‑251 (MQ=255)
                                                                             ggTTTTGGATTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGGGACA‑GTTCACGGCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGCTCGCGGTAACTCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGGCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGCTCTGTTTACGGTAGAGGAAGTAAACGTCTTTGTTCTGcc                                                                                                                                                  >  2:225477/1‑251 (MQ=255)
                                                                                tgtGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGGCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAg                                                                                                                                                <  2:368582/248‑1 (MQ=255)
                                                                                    gTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTc                                                                                                                                             >  1:39942/1‑249 (MQ=255)
                                                                                    gTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAg                                                                                                                                           <  1:16094/251‑1 (MQ=255)
                                                                                      ttCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGTTCAGa                                                                                                                                          <  1:328663/250‑1 (MQ=255)
                                                                                          ggCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGAtt                                                                                                                                        <  1:122126/248‑1 (MQ=255)
                                                                                                    aGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAAGCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCGTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGAGCTGCCAGGTCAGATTAAGGTCATGTTgg                                                                                                                           <  2:380738/251‑1 (MQ=255)
                                                                                                     gAAACCGGCAGCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACGCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCGTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCTTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGt                                                                                                                          <  1:225623/251‑1 (MQ=255)
                                                                                                      aacccGGCATCAGCACAGCCCTGGTAGGTTACA‑GTTAACGTCGCCCCCGCACTCGCCTGGTGGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCATACGTCTTGCGGCAGCTGCACGGCGGCAATTTTCGCGTGTTCCGTCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTTCCAGGTCAGAGTAAGGTCATGTTGGtt                                                                                                                         <  2:335045/248‑1 (MQ=255)
                                                                                                        cccGGCCTCAGCGCAGGCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTTATGGTGACGGGAAGCGTCAGCCTATCTCGGTAAATCGCGCTTTTGCCCTAAAACTCATCCTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCTTGTTCCGGCGTAATGCGTATCTTTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGGTGGtttt                                                                                                                       <  2:200381/250‑1 (MQ=255)
                                                                                                        aCCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGtttt                                                                                                                       <  2:240200/251‑1 (MQ=255)
                                                                                                            gCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTg                                                                                                                   <  1:130241/251‑1 (MQ=255)
                                                                                                                    acaGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCCTGTTGGTTTTGCTGAAAATCa                                                                                                            >  2:88480/1‑250 (MQ=255)
                                                                                                                     caGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATTCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCaaa                                                                                                          <  1:422152/251‑1 (MQ=255)
                                                                                                                       gtacTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGGGACGGTAAGCGTCATCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGTCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTCAGGTCAGGTTGGTTTTGCCGAAAATCAAAAg                                                                                                        <  2:379689/248‑1 (MQ=255)
                                                                                                                          cTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCa                                                                                                      >  1:446871/1‑250 (MQ=255)
                                                                                                                             gTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAgg                                                                                                  <  1:64483/251‑1 (MQ=255)
                                                                                                                              tAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGttt                                                                                                                        >  2:165477/1‑228 (MQ=255)
                                                                                                                              tAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGttt                                                                                                                        <  1:165477/228‑1 (MQ=255)
                                                                                                                                                              ccTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTccc                                                                  <  1:286171/251‑1 (MQ=255)
                                                                                                                                                                  gTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGgcg                                                              <  2:388300/251‑1 (MQ=255)
                                                                                                                                                                                gAAGCGTCAGCCGATCGCGGTAAAGCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCTGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGGCGAATAATCCg                                                <  1:349346/251‑1 (MQ=255)
                                                                                                                                                                                                           cGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCGGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATTCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGAGCTGCAAAGTAGc                     >  1:18074/1‑251 (MQ=255)
                                                                                                                                                                                                           cGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAg                         <  1:228300/247‑1 (MQ=255)
                                                                                                                                                                                                           cGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTGTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGTTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGACGAATAATCCGGCAAAAACGGAAGTTCTGCAAAg                         <  2:89525/247‑1 (MQ=255)
                                                                                                                                                                                                             cTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGGAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCGGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCTGGGACAAATTGTGAACTTCCCGGCGCGTCGACTAATCCGGCAAAAACGGAAGAGCTGCAAAGtcgc                     >  1:184011/1‑246 (MQ=255)
                                                                                                                                                                                                                 tGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGAt                <  1:244348/250‑1 (MQ=255)
                                                                                                                                                                                                                  gCCGTAAAAATAATCTTCATGCCAGACGTCTTGCGGCAGCTGCCCGACGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGACGCGTCGAATACTCCGCCAAAAACGGAAGTGCTGCAAAGTAACAGGATc               >  2:71155/1‑250 (MQ=255)
                                                                                                                                                                                                                   ccGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAg             <  1:387567/251‑1 (MQ=255)
                                                                                                                                                                                                                          actccaattCATGCCAGCCGGCTGGCGCCAGCTGCACGGCAGCAATTTTCGCGTGTACCGGCGTAATGCGGATCTTTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCCGGTCAGCCTAAGGTCATGTTGGTTTTGCTGAAAATAACAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCACGGCGCGTCGAATAATCCGGCAAAAACGGAAGAGCTGCAAAGTAGCAGGATCAGCGTaa        >  2:231358/8‑249 (MQ=255)
                                                                                                                                                                                                                          aCTCATCTTCATGCCAGACGTCCTGCGGCCGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCgaaa        >  1:364995/1‑246 (MQ=255)
                                                                                                                                                                                                                              atcgtCATGCCAGACGTCGTGCGGCAGCTGCACGTAGGCAATGTTCGCTTGTTCCGGCGTAATGCGGATCTTTTTACGGTAGAGGTAGTAACCGTCTTTGATCGGACAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGAGGGGACAAATTGTGAACGTACCGGAGAGTCGAATAATCCGGCAAAAATGGAAGTGCTGCAAAGTAGCAGGATAAGCGTAAAGaggc  <  2:203199/247‑4 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
GCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACA‑GTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGC  >  NC_000913/4364561‑4365032

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