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

GCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGC  >  NC_000913/4364561‑4365032
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
gcgGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTg                                                                                                                                                                                                                               <  1:437865/251‑1 (MQ=255)
   gTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCAc                                                                                                                                                                                                                            <  1:407371/251‑1 (MQ=255)
      ggCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTc                                                                                                                                                                                                                         <  2:120856/251‑1 (MQ=255)
          gctCTTGCTGCGGAACAGACACAGGCTGTGGTGCTGCGTTGTGGGCGACCACTACGCTTAACGGAACGGTTTTGGTGTCTGGCGGATAACAGAAACCGGCATCAGTTCAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCCCTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGc                                                                                                                                                                                                                      <  2:352021/250‑1 (MQ=255)
           ctctTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGCGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGcc                                                                                                                                                                                                                     >  1:223436/1‑249 (MQ=255)
            tctTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCCACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACCGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGCCGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAAt                                                                                                                                                                                                                   >  1:283392/1‑251 (MQ=255)
                   gcGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTcgcg                                                                                                                                                                                                            <  2:255630/251‑1 (MQ=255)
                    cGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCTCgt                                                                                                                                                                                                           >  2:64483/1‑251 (MQ=255)
                          aGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCg                                                                                                                                                                                                     <  2:218482/251‑1 (MQ=255)
                               caGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTa                                                                                                                                                                                                <  1:179775/251‑1 (MQ=255)
                                       ggCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGAt                                                                                                                                                                                        >  1:368582/1‑251 (MQ=255)
                                          gcTGCGTTGTTGGCGACCACTTCGCTTAAGGGCACGGTTTTGGTTTCTGGCGGATCACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTg                                                                                                                                                                                     >  2:31702/1‑251 (MQ=255)
                                          gcTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTg                                                                                                                                                                                     <  1:404971/251‑1 (MQ=255)
                                           cTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGt                                                                                                                                                                                    >  1:348276/1‑251 (MQ=255)
                                             gCGTTGTTGGCGCCCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGAGCTGtt                                                                                                                                                                                   <  2:225565/250‑1 (MQ=255)
                                                  gttgGCGACCACTTCGCTTAACGGAACGTTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCCGCTGCACGTAGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACggg                                                                                                                                                                             >  2:19040/1‑250 (MQ=255)
                                                       cGACCACTTAGCTTACCGGAACGGTTTTGGTTTCTGGCGGATAACAGAACCCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAACTCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAgg                                                                                                                                                                        >  2:155918/1‑251 (MQ=255)
                                                            aCTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGTCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGgtagta                                                                                                                                                                   <  1:284797/251‑1 (MQ=255)
                                                              ttCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGAGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGCTCTGTTTACGGTAGAGGTAGTAAc                                                                                                                                                                 >  1:406713/1‑251 (MQ=255)
                                                               tctcTTAACGGAACGGTTTGGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGGTAACGTCGCTCCCGCACTCGCCTGGTTGATGGGGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGGCTTGCGGCAGCTGCACGTCGGCAATGTTCGCGTGTTCCGGCGTAAGGCGTATCTGTTTACGGTAGAGGTAGTAAc                                                                                                                                                                 <  2:272899/247‑1 (MQ=255)
                                                                cgctttaccgacgggTTTTGGTTGCTGGCGGATAACAGAACCCGGCAGGAGCACAGCCCTGCTAGGTGCACGTTAACGTCGCTCCCGGACTCGCCTGGTTGCTGGTTACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCGCGTCGGCAATTTTCGCGTGTTCCGGGGTAATGCTGATCTGTTGCCGGTAGAGgttgtcacc                                                                                                                                                                <  2:125528/237‑8 (MQ=255)
                                                                      acggaacggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTg                                                                                                                                                         <  2:213330/251‑1 (MQ=255)
                                                                       cggaacgtTGTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTTACAGTTAACGTCTCTCCCGCTCTCGCCTGGTTGATGGTGACGTGAACCGTCAGCCGATCGCTGTAATTCTCGCTTTTGCCGTAAAACTCATCTTCAGGCCAGCCGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTTACCGTCTTTGa                                                                                                                                                        <  2:2152/251‑1 (MQ=255)
                                                                        ggaacggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGAt                                                                                                                                                       <  1:289050/251‑1 (MQ=255)
                                                                            cggTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCCGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTCCGGTAGCGGTAGTAACCGGCTTTGATCTGc                                                                                                                                                   >  2:237260/1‑251 (MQ=255)
                                                                             ggTTTTGGATTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGGGACAGTTCACGGCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGCTCGCGGTAACTCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGGCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGCTCTGTTTACGGTAGAGGAAGTAAACGTCTTTGTTCTGcc                                                                                                                                                  >  2:225477/1‑251 (MQ=255)
                                                                                tgtGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGGCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAg                                                                                                                                                <  2:368582/248‑1 (MQ=255)
                                                                                    gTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTc                                                                                                                                             >  1:39942/1‑249 (MQ=255)
                                                                                    gTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAg                                                                                                                                           <  1:16094/251‑1 (MQ=255)
                                                                                      ttCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGTTCAGa                                                                                                                                          <  1:328663/250‑1 (MQ=255)
                                                                                          ggCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGAtt                                                                                                                                        <  1:122126/248‑1 (MQ=255)
                                                                                                    aGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAAGCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCGTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGAGCTGCCAGGTCAGATTAAGGTCATGTTgg                                                                                                                           <  2:380738/251‑1 (MQ=255)
                                                                                                     gAAACCGGCAGCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACGCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCGTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCTTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGt                                                                                                                          <  1:225623/251‑1 (MQ=255)
                                                                                                      aacccGGCATCAGCACAGCCCTGGTAGGTTACAGTTAACGTCGCCCCCGCACTCGCCTGGTGGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCATACGTCTTGCGGCAGCTGCACGGCGGCAATTTTCGCGTGTTCCGTCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTTCCAGGTCAGAGTAAGGTCATGTTGGtt                                                                                                                         <  2:335045/248‑1 (MQ=255)
                                                                                                        cccGGCCTCAGCGCAGGCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTTATGGTGACGGGAAGCGTCAGCCTATCTCGGTAAATCGCGCTTTTGCCCTAAAACTCATCCTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCTTGTTCCGGCGTAATGCGTATCTTTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGGTGGtttt                                                                                                                       <  2:200381/250‑1 (MQ=255)
                                                                                                        aCCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGtttt                                                                                                                       <  2:240200/251‑1 (MQ=255)
                                                                                                            gCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTg                                                                                                                   <  1:130241/251‑1 (MQ=255)
                                                                                                                    acaGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCCTGTTGGTTTTGCTGAAAATCa                                                                                                            >  2:88480/1‑250 (MQ=255)
                                                                                                                     caGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATTCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCaaa                                                                                                          <  1:422152/251‑1 (MQ=255)
                                                                                                                       gtacTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGGGACGGTAAGCGTCATCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGTCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTCAGGTCAGGTTGGTTTTGCCGAAAATCAAAAg                                                                                                        <  2:379689/248‑1 (MQ=255)
                                                                                                                          cTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCa                                                                                                      >  1:446871/1‑250 (MQ=255)
                                                                                                                             gTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAgg                                                                                                  <  1:64483/251‑1 (MQ=255)
                                                                                                                              tAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGttt                                                                                                                        >  2:165477/1‑228 (MQ=255)
                                                                                                                              tAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGTCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGttt                                                                                                                        <  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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
GCGGTGGGCTGCTCTTGCTGCGGAACAGACACAGGCTGTGGCGCTGCGTTGTTGGCGACCACTTCGCTTAACGGAACGGTTTTGGTTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTGATCCGCGGGGACAAATTGTGAACGTCCCGGCGCGTCGAATAATCCGGCAAAAACGGAAGTGCTGCAAAGTAGCAGGATCAGCGTAAAGATGC  >  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: