Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,547,951 A→G intergenic (+76/+35) gntT → / ← malQ high‑affinity gluconate transporter/4‑alpha‑glucanotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,547,9510AG96.9% 225.0 / ‑2.1 65intergenic (+76/+35)gntT/malQhigh‑affinity gluconate transporter/4‑alpha‑glucanotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (42/21);  total (43/22)
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.76e-01

ACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAA  >  NC_000913/3547819‑3548082
                                                                                                                                    |                                                                                                                                   
aCCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCa                                                                                                                              <  2:976280/140‑1 (MQ=255)
  cATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAg                                                                                                                            >  1:1297627/1‑140 (MQ=255)
    ttATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAgca                                                                                                                          >  2:477394/1‑140 (MQ=255)
     tATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAgcat                                                                                                                         >  2:314644/1‑140 (MQ=255)
         tcGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAgcatcgca                                                                                                                     >  2:1015165/1‑140 (MQ=255)
          cGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAgcatcgcat                                                                                                                    <  2:1185679/140‑1 (MQ=255)
          cGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAgcatcgcat                                                                                                                    <  1:42099/140‑1 (MQ=255)
                 cGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCg                                                                                                                 <  2:532085/136‑1 (MQ=255)
                 cGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCg                                                                                                                 >  1:532085/1‑136 (MQ=255)
                    tCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCaa                                                                                                          <  1:1317111/140‑1 (MQ=255)
                       ggTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCg                                                                                                       >  2:1027894/1‑140 (MQ=255)
                        gTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGc                                                                                                      <  1:1309650/140‑1 (MQ=255)
                            gCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCActc                                                                                                  <  1:887023/140‑1 (MQ=255)
                                tCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTAct                                                                                              >  1:1033180/1‑140 (MQ=255)
                                tCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTAct                                                                                              >  2:219833/1‑140 (MQ=255)
                                tCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTAct                                                                                              >  2:1233745/1‑140 (MQ=255)
                                  tgctgctgAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTActtc                                                                                            >  2:606738/1‑140 (MQ=255)
                                                 tGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGct                                                                             >  2:1187760/1‑140 (MQ=255)
                                                          aaCACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGt                                                                    <  1:1015165/140‑1 (MQ=255)
                                                               tgactgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTc                                                               >  2:587268/1‑140 (MQ=255)
                                                                gactgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTcc                                                              <  2:1297627/140‑1 (MQ=255)
                                                                 actgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                             >  2:98832/1‑140 (MQ=255)
                                                                 actgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                             >  1:1113286/1‑140 (MQ=255)
                                                                  ctgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCaa                                                            >  1:533650/1‑140 (MQ=255)
                                                                         aTGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTc                                                     >  1:64740/1‑140 (MQ=255)
                                                                            cggcgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGc                                                  <  2:533650/140‑1 (MQ=255)
                                                                               cgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAAc                                               >  2:766969/1‑140 (MQ=255)
                                                                                   aCCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGt                                           <  1:471570/140‑1 (MQ=255)
                                                                                    ccGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGtt                                          <  1:1339861/140‑1 (MQ=255)
                                                                                                  cTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCgg                            >  1:1194341/1‑140 (MQ=255)
                                                                                                  cTACGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCgg                            >  1:810945/1‑140 (MQ=255)
                                                                                                    aCGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCa                          >  2:1324327/1‑140 (MQ=255)
                                                                                                    aCGATTCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCa                          >  1:565634/1‑140 (MQ=255)
                                                                                                         tCGGGATGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACAt                     <  1:1324327/140‑1 (MQ=255)
                                                                                                               tGAATTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTc               <  2:598450/140‑1 (MQ=255)
                                                                                                               ccgaatcgtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTc               >  2:697602/8‑140 (MQ=255)
                                                                                                               ccgaatcgtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTc               >  1:67535/8‑140 (MQ=255)
                                                                                                                cgaatcgtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCa              >  2:271563/7‑140 (MQ=255)
                                                                                                                cgaatcgtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCa              >  2:299355/7‑140 (MQ=255)
                                                                                                                  aTTAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAg            >  2:1183109/1‑140 (MQ=255)
                                                                                                                   ttAGTAGGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGc           <  2:384845/140‑1 (MQ=255)
                                                                                                                    tcgtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCg          <  1:271563/138‑1 (MQ=255)
                                                                                                                    tcgtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCg          <  1:697602/138‑1 (MQ=255)
                                                                                                                     cgtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGt         <  2:67535/139‑1 (MQ=255)
                                                                                                                     agtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGt         <  1:314644/140‑1 (MQ=255)
                                                                                                                     agtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGt         <  2:64740/140‑1 (MQ=255)
                                                                                                                     agtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGt         >  1:53290/1‑140 (MQ=255)
                                                                                                                      gtagGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGtt        >  2:876893/1‑140 (MQ=255)
                                                                                                                        agGCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGc      >  1:445006/1‑140 (MQ=255)
                                                                                                                          gCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg    >  1:22047/1‑140 (MQ=255)
                                                                                                                          gCCGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg    >  2:192102/1‑140 (MQ=255)
                                                                                                                           ccGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTc               <  2:149994/128‑1 (MQ=255)
                                                                                                                           ccGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTc               >  1:149994/1‑128 (MQ=255)
                                                                                                                           ccGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGa   >  1:724397/1‑140 (MQ=255)
                                                                                                                           ccGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGa   >  2:1027844/1‑140 (MQ=255)
                                                                                                                           ccGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGa   >  1:3290/1‑140 (MQ=255)
                                                                                                                           ccGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGa   >  1:1119765/1‑140 (MQ=255)
                                                                                                                            cGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGt         >  1:564376/1‑133 (MQ=255)
                                                                                                                            cGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGt         <  2:564376/133‑1 (MQ=255)
                                                                                                                            cGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGaa  >  2:634292/1‑140 (MQ=255)
                                                                                                                            cGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGaa  >  1:763748/1‑140 (MQ=255)
                                                                                                                            cGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGaa  >  1:26888/1‑140 (MQ=255)
                                                                                                                            cGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGaa  >  2:569508/1‑140 (MQ=255)
                                                                                                                            cGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGaa  >  2:220821/1‑140 (MQ=255)
                                                                                                                            cGGATAAGGCGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTATGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGaa  >  1:1192842/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                   
ACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAA  >  NC_000913/3547819‑3548082

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