Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,557,086 C→A 10.3% L145L (CTG→CTT gntZ ← 6‑phosphogluconate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,557,0860CA10.3% 196.9 / 9.4 68L145L (CTG→CTTgntZ6‑phosphogluconate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (3/4);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

TGCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTTGTAGAAGGT  >  NC_002947/4556938‑4557235
                                                                                                                                                    |                                                                                                                                                     
tGCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCa                                                                                                                                                      <  1:405887/150‑1 (MQ=255)
 gCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAg                                                                                                                                                     <  1:214096/150‑1 (MQ=255)
   cGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGgc                                                                                                                                                   <  2:1046879/150‑1 (MQ=255)
    gTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCGGgcg                                                                                                                                                  >  2:469532/1‑150 (MQ=255)
     ttATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCATGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGgcgc                                                                                                                                                 >  1:855155/1‑150 (MQ=255)
     ttATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGgcgc                                                                                                                                                 <  2:285333/150‑1 (MQ=255)
      tATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGGCGCt                                                                                                                                                <  2:1158167/150‑1 (MQ=255)
          cACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTCCGCGGCATGTCTCCCACGCACGCCGCCAGGTCCTTGAACAGCGGTTCCAGGCGCTTGaa                                                                                                                                            >  2:584971/1‑150 (MQ=255)
          cACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGaa                                                                                                                                            <  2:57304/150‑1 (MQ=255)
            ccATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAaca                                                                                                                                          >  1:602096/1‑150 (MQ=255)
             cATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAacac                                                                                                                                         <  2:855155/150‑1 (MQ=255)
               ttttCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGt                                                                                                                                       >  1:984297/1‑150 (MQ=255)
                tttCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTc                                                                                                                                      <  2:984297/150‑1 (MQ=255)
                 ttCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTcc                                                                                                                                     <  2:790925/150‑1 (MQ=255)
                    aCGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCttt                                                                                                                                  <  1:880230/150‑1 (MQ=255)
                            cccGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCagcgcc                                                                                                                          >  1:170744/1‑150 (MQ=255)
                             ccGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGgc                                                                                                                                                   <  2:119096/124‑1 (MQ=255)
                             ccGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGgc                                                                                                                                                   >  1:119096/1‑124 (MQ=255)
                             ccGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCACCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGccgccg                                                                                                                         >  1:840957/1‑150 (MQ=255)
                             ccGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCCGGCTGTTGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGccgccg                                                                                                                         >  1:510192/1‑150 (MQ=255)
                                  gcCCGGTGGACCGGCATGGATGTACATCTGGACGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGatcat                                                                                                                    >  1:937441/1‑150 (MQ=255)
                                      ggTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCGGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATg                                                                                                                <  1:469532/150‑1 (MQ=255)
                                       gTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGc                                                                                                               >  1:1085007/1‑150 (MQ=255)
                                        tGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCACCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCa                                                                                                              <  2:840957/150‑1 (MQ=255)
                                            ccGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATTTCTCCCACGCCCGGCACCAGCGCCTTGATCAGCGGTTCCAGGCGCTCGAAAACCGCCTGGTCGCCGCCGATCATCATGCAGTAg                                                                                                          <  1:634580/150‑1 (MQ=255)
                                              ggCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGAGGTCCGCCGCACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGcc                                                                                                        <  2:937441/150‑1 (MQ=255)
                                               gCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAg                                                                                                             >  2:1049021/1‑144 (MQ=255)
                                               gCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAg                                                                                                             <  1:1049021/144‑1 (MQ=255)
                                                    gATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTc                                                                                                  >  2:826154/1‑150 (MQ=255)
                                                         aGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGc                                                                                             <  2:829128/150‑1 (MQ=255)
                                                            cGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCttt                                                                                                                                  >  2:378313/1‑110 (MQ=255)
                                                            cGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCttt                                                                                                                                  <  1:378313/110‑1 (MQ=255)
                                                             gTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCca                                                                                         <  1:1126710/150‑1 (MQ=255)
                                                                cTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCacac                                                                                      <  2:602096/150‑1 (MQ=255)
                                                                 tCGGCACGCTGGTGCTCGCCGCTGCTTCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACAcg                                                                                     >  1:766323/1‑150 (MQ=255)
                                                                 tCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACAcg                                                                                     >  2:821473/1‑150 (MQ=255)
                                                                        gCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGAAGCGGTCCAGGCCCCACACgccgccag                                                                              >  1:622650/1‑148 (MQ=255)
                                                                            gTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATACAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGt                                                                          >  2:276901/1‑150 (MQ=255)
                                                                              gCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGgcgc                                                                                                                                                 <  2:482276/77‑1 (MQ=37)
                                                                              gCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGgcgc                                                                                                                                                 >  1:482276/1‑77 (MQ=37)
                                                                                tCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCg                                                                      >  2:540529/1‑150 (MQ=255)
                                                                                tCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCg                                                                      >  1:526285/1‑150 (MQ=255)
                                                                                     gcTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTc                                                                 <  2:105592/150‑1 (MQ=255)
                                                                                      cTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTcc                                                                >  2:224412/1‑150 (MQ=255)
                                                                                         cGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGa                                                             <  1:826154/150‑1 (MQ=255)
                                                                                          ggCCATGGGTACGAGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGAt                                                            >  1:24071/1‑150 (MQ=255)
                                                                                            ccATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATaa                                                          <  2:766323/150‑1 (MQ=255)
                                                                                                   tACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGTCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGc                                                   <  2:643699/150‑1 (MQ=255)
                                                                                                                gCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCCCCGCCGATCATGGGGTGGTAGCGGCGGTCCAGGCCCCACACGGCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCGTGGCCAg                                      >  2:346791/1‑150 (MQ=255)
                                                                                                                        ccGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAg                              >  1:199451/1‑150 (MQ=255)
                                                                                                                              ccAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTTCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGAAAGTCCCGACACCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCgcggc                        <  2:170744/150‑1 (MQ=255)
                                                                                                                              ccAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCgcggc                        >  2:332215/1‑150 (MQ=255)
                                                                                                                               cAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAg                                                     >  2:764185/1‑120 (MQ=255)
                                                                                                                               cAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAg                                                     <  1:764185/120‑1 (MQ=255)
                                                                                                                                   gcCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATg                   <  2:24071/150‑1 (MQ=255)
                                                                                                                                   gcCTGGAACCGCTGTTCAAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGt                                                                                                   >  1:1070393/1‑70 (MQ=25)
                                                                                                                                   gcCTGGAACCGCTGTTCAAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGt                                                                                                   <  2:1070393/70‑1 (MQ=25)
                                                                                                                                   gcCTGGAACCGCTGTTCAAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTcc                                                                                                 <  2:230721/72‑1 (MQ=37)
                                                                                                                                   gcCTGGAACCGCTGTTCAAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTcc                                                                                                 >  1:230721/1‑72 (MQ=37)
                                                                                                                                    ccTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGaa                                                                            <  2:1093853/92‑1 (MQ=255)
                                                                                                                                    ccTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGaa                                                                            >  1:1093853/1‑92 (MQ=255)
                                                                                                                                     cTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCa                               >  2:257405/1‑136 (MQ=255)
                                                                                                                                     cTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCa                               <  1:257405/136‑1 (MQ=255)
                                                                                                                                         aaCAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTcc             <  2:510192/150‑1 (MQ=255)
                                                                                                                                         aaCAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTcc             <  1:821473/150‑1 (MQ=255)
                                                                                                                                          aCAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCt            >  1:334670/1‑150 (MQ=255)
                                                                                                                                            aGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTTg          >  2:992274/1‑150 (MQ=255)
                                                                                                                                             gCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTTGt         <  2:1085007/150‑1 (MQ=255)
                                                                                                                                                    cAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTTGTAGAAGgt  >  1:607879/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCCGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTTGTAGAAGGT  >  NC_002947/4556938‑4557235

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