Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,157,442 G→T 6.5% G274C (GGC→TGC)  dkgA → 2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,157,4420GT 6.5% 220.1 / 2.3 61G274C (GGC→TGC) dkgA2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base G (22/35);  new base T (2/2);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATG  >  NC_000913/3157304‑3157578
                                                                                                                                          |                                                                                                                                        
gtgATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGAGGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcgg                                                                                                                                         >  2:561501/1‑140 (MQ=255)
   aTCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcggcgg                                                                                                                                      >  1:1017575/1‑140 (MQ=255)
   aTCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcggcgg                                                                                                                                      >  2:385381/1‑140 (MQ=255)
   aTCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcggcgg                                                                                                                                      >  1:537498/1‑140 (MQ=255)
   aTCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcggcgg                                                                                                                                      >  1:312457/1‑140 (MQ=255)
   aTCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcggcgg                                                                                                                                      >  2:594042/1‑140 (MQ=255)
   aTCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcggcgg                                                                                                                                      >  2:642459/1‑140 (MQ=255)
   aTCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcggcgg                                                                                                                                      >  1:1506499/1‑140 (MQ=255)
   aTCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTcggcgg                                                                                                                                      >  1:1017508/1‑140 (MQ=255)
     cccGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCt                                                                                                                                    <  1:1326530/140‑1 (MQ=255)
       cGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTaa                                                                                                                                  <  1:935681/140‑1 (MQ=255)
       cGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTaa                                                                                                                                  >  2:779806/1‑140 (MQ=255)
               gTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCaaa                                                                                                                          >  1:703373/1‑140 (MQ=255)
                        tCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTTACCAGTTCTGCGGCTAACATGCAAATTCTCCCgg                                                                                                                 >  2:510754/1‑140 (MQ=255)
                        tCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTTACCAGTTCTGCGGCTAACATGCAAATTCTCCCgg                                                                                                                 >  2:511112/1‑140 (MQ=255)
                           cGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTgg                                                                                                              <  1:1649571/140‑1 (MQ=255)
                                 gCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCTAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTaa                                                                                                        >  2:1605836/1‑140 (MQ=255)
                                 gCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTaa                                                                                                        >  2:897097/1‑140 (MQ=255)
                                 gCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTaa                                                                                                        >  2:1073845/1‑140 (MQ=255)
                                 gCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTaa                                                                                                        >  2:761269/1‑140 (MQ=255)
                                      aaaCTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTc                                                                                                   <  1:1394985/140‑1 (MQ=255)
                                      aaaCTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTc                                                                                                   <  1:799314/140‑1 (MQ=255)
                                            tGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTAc                                                                                             >  2:348572/1‑140 (MQ=255)
                                              ctGTCTGGGATTTCCGTATCGACAAAGTAGAACTCGGCGAAATTGCAAAACTCGATCTGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGGTCCGCTACCg                                                                                           <  1:1638109/139‑1 (MQ=255)
                                              aTGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCg                                                                                           <  1:1623329/140‑1 (MQ=255)
                                                        tttCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAg                                                                                 <  1:149123/140‑1 (MQ=255)
                                                          tCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGaa                                                                               <  1:187632/140‑1 (MQ=255)
                                                          tCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGaa                                                                               <  1:1502140/140‑1 (MQ=255)
                                                          tCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGaa                                                                               >  1:1075012/1‑140 (MQ=255)
                                                          tCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGaa                                                                               <  1:1616609/140‑1 (MQ=255)
                                                          tCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGaa                                                                               <  2:587989/140‑1 (MQ=255)
                                                          tCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGaa                                                                               <  2:590040/140‑1 (MQ=255)
                                                          tCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGaa                                                                               <  2:869291/140‑1 (MQ=255)
                                                             gTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATc                                                                            <  2:504615/140‑1 (MQ=255)
                                                                tcGACAAAGTCGATCTCTGCGAAAGTGCAAAACTCGTTCAGGGCAAGCGTCTCGGTCCCGATCCTGATCTGTTCGGCGGCTAACATGCATTTTCTCCCGGTGGAGGTAATGTTCCGCTACCTGACTTTTCAGAATTCAtt                                                                         <  2:861216/140‑1 (MQ=255)
                                                                tcGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTACCTTGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGATATCAtt                                                                         <  2:949482/140‑1 (MQ=255)
                                                                tcGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCAtt                                                                         <  2:1448547/140‑1 (MQ=255)
                                                                       aGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTccc                                                                  <  1:836860/140‑1 (MQ=255)
                                                                                       aTTGCAAACCTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCgt                                                  <  1:1094846/140‑1 (MQ=255)
                                                                                                   gATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTcc                                      <  1:48257/140‑1 (MQ=255)
                                                                                                          gCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTGGTTACTCTTCCTTGTTCa                               >  2:1478418/1‑140 (MQ=255)
                                                                                                             aGCGTCTCGGTCCCGATCCTTACCAGTTCTGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGa                            <  1:510754/140‑1 (MQ=255)
                                                                                                             aGCGTCTCGGTCCCGATCCTTACCAGTTCTGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGa                            <  1:511112/140‑1 (MQ=255)
                                                                                                             aGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGa                            <  1:1194753/140‑1 (MQ=255)
                                                                                                             aGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGa                            <  1:1677898/140‑1 (MQ=255)
                                                                                                             aGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGa                            <  1:1678073/140‑1 (MQ=255)
                                                                                                                gTCTCGGTCCCGTTCCTGACCAGTTCGGCGGCTAACATGCAATTTCTCCCGGTGGCGGTATTGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGGGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATg                         <  2:1082534/140‑1 (MQ=255)
                                                                                                                gTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGTCTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATg                         <  2:555692/140‑1 (MQ=255)
                                                                                                                   tcGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCa                      >  1:566607/1‑140 (MQ=255)
                                                                                                                   tcGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCa                      >  1:567077/1‑140 (MQ=255)
                                                                                                                   tcGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCa                      >  1:333168/1‑140 (MQ=255)
                                                                                                                   tcGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCa                      >  1:1309234/1‑140 (MQ=255)
                                                                                                                   tcGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCa                      >  1:1327018/1‑140 (MQ=255)
                                                                                                                   tcGGTCCAGATCCTGAGGAGTTCGGCGGCTAAGTTTCAAATTCTCCCGGTGGCGGTAATGTTGCGCTACCGGACTTTTCAGAAATCTTTTATTCCCCTCGCGTCCCGCGCGTTGTTACTCTTCATTGTTCAGGAATGCCa                      <  2:419864/140‑1 (MQ=255)
                                                                                                                      gTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAat                   >  1:823915/1‑140 (MQ=255)
                                                                                                                        cccGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAatat                 <  2:312457/140‑1 (MQ=255)
                                                                                                                        cccGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAatat                 <  2:409914/140‑1 (MQ=255)
                                                                                                                        cccGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAatat                 <  2:858036/140‑1 (MQ=255)
                                                                                                                        cccGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAatat                 <  2:857891/140‑1 (MQ=255)
                                                                                                                         ccGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTGCTGGACTGCCAAAtata                >  1:1604819/1‑140 (MQ=255)
                                                                                                                             tCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGa            <  1:837164/140‑1 (MQ=255)
                                                                                                                             tCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGa            <  1:836486/140‑1 (MQ=255)
                                                                                                                             tCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGa            <  1:1270876/140‑1 (MQ=255)
                                                                                                                             tCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGa            <  2:686700/140‑1 (MQ=255)
                                                                                                                                 gACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCTGAAATCATTTATTCCCCTCGCGTCCCGCCCGTGGTTACTCTTCCTTGTTCAGGAATGCCAATTATAAGGAcatc        >  2:9035/1‑140 (MQ=255)
                                                                                                                                 gACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAATTATAAGGAcatc        >  2:8894/1‑140 (MQ=255)
                                                                                                                                 gACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcatc        >  2:1171371/1‑140 (MQ=255)
                                                                                                                                      gTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcatcatcat   >  1:407130/1‑140 (MQ=255)
                                                                                                                                      gTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGAcatcatcat   >  1:575220/1‑140 (MQ=255)
                                                                                                                                       ttCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATg  >  1:298833/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                        
GTGATCCCGAAATCGGTCACACCTTCACGTATTGCCGAAAACTTTGATGTCTGGGATTTCCGTCTCGACAAAGACGAACTCGGCGAAATTGCAAAACTCGATCAGGGCAAGCGTCTCGGTCCCGATCCTGACCAGTTCGGCGGCTAACATGCAAATTCTCCCGGTGGCGGTAATGTTCCGCTACCGGACTTTTCAGAAATCATTTATTCCCCTCGCGTCCCGCCCGTTGTTACTCTTCCTTGTTCAGGAATGCCAAATATAAGGACATCATCATG  >  NC_000913/3157304‑3157578

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