Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,557,086 C→A 12.4% 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,0860CA12.4% 180.8 / 11.4 65L145L (CTG→CTT) gntZ6‑phosphogluconate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base A (5/3);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTTG  >  NC_002947/4556941‑4557227
                                                                                                                                                 |                                                                                                                                             
cGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGTAGCGGCCTGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTAGGCGGGATGGCGCCCACGCCCGGCGCCAGCGCCTGGAACAGCGGTTCCAGgc                                                                                                                                           >  1:416719/1‑150 (MQ=255)
cGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGgc                                                                                                                                           >  1:850578/1‑150 (MQ=255)
 gTTATCCACCCTTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCGCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGgcg                                                                                                                                          <  2:416719/150‑1 (MQ=255)
     ttcACCATTTTCACGAAGTGCCAGGAGCACGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGGCGCtgg                                                                                                                                      <  2:589009/148‑3 (MQ=255)
       cACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGaa                                                                                                                                    >  2:360388/1‑150 (MQ=255)
       cACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGaa                                                                                                                                    <  2:618612/150‑1 (MQ=255)
       cACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTCGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGGctctggcg                                                                                                                                    >  1:589009/1‑146 (MQ=255)
        aCCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAac                                                                                                                                   >  1:1131225/1‑150 (MQ=255)
        aCCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAac                                                                                                                                   >  1:70863/1‑150 (MQ=255)
        aCCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAac                                                                                                                                   <  1:697347/150‑1 (MQ=255)
          cATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGTGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAacac                                                                                                                                 <  1:299786/150‑1 (MQ=255)
                   gTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTc                                                                                                                        <  2:70863/150‑1 (MQ=255)
                                 ccGGTGGACCGGCATGGATGTAGCCGTGCGCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGAtcatca                                                                                                          >  1:212092/1‑150 (MQ=255)
                                   ggTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATg                                                                                                        >  2:155562/1‑150 (MQ=255)
                                    gTGGACCGGCATGGATGTAGCCGTACTAGGAATGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGc                                                                                                       <  2:1125569/150‑1 (MQ=255)
                                          cGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGACGATCAGCTTACGGCTATACCGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGc                                                                                                 >  1:348846/1‑150 (MQ=255)
                                            gCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCcg                                                                                               >  2:1099087/1‑150 (MQ=255)
                                            gCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCcg                                                                                               >  2:178640/1‑150 (MQ=255)
                                             cATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCcgc                                                                                              >  1:1052360/1‑150 (MQ=255)
                                                      aGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGc                                                                                     <  2:783254/150‑1 (MQ=255)
                                                        ccGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGTCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGccc                                                                                   <  1:807119/150‑1 (MQ=255)
                                                             cTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCacac                                                                              >  2:231391/1‑150 (MQ=255)
                                                                    cGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGccgccg                                                                       <  1:193472/150‑1 (MQ=255)
                                                                     gCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGGCGCtg                                                                                                                                       <  1:677923/85‑2 (MQ=37)
                                                                     gCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGGCGCtg                                                                                                                                       >  2:677923/1‑84 (MQ=37)
                                                                      cTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAg                                                                                                                                             <  1:1094837/78‑1 (MQ=255)
                                                                      cTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAg                                                                                                                                             >  2:1094837/1‑78 (MQ=255)
                                                                       tGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGaa                                                                    <  1:202403/150‑1 (MQ=255)
                                                                            cTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTccc                                                               <  2:348846/150‑1 (MQ=255)
                                                                            cTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTccc                                                               <  1:547044/150‑1 (MQ=255)
                                                                             tCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTACAGGCCCCACACGCCGCCGGAAGTCCCg                                                              >  1:493188/1‑150 (MQ=255)
                                                                                        gCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCacac                                                                              <  1:539057/123‑1 (MQ=255)
                                                                                        gCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCacac                                                                              >  2:539057/1‑123 (MQ=255)
                                                                                         ccATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCGAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATaa                                                  >  2:831592/1‑150 (MQ=255)
                                                                                         ccATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATaa                                                  >  1:691822/1‑150 (MQ=255)
                                                                                             gggTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCa                                              <  2:850578/150‑1 (MQ=255)
                                                                                              ggTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAg                                             <  1:1021875/150‑1 (MQ=255)
                                                                                               gTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGAGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAgg                                            >  2:1194562/1‑150 (MQ=255)
                                                                                               gTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGGCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAgg                                            >  2:353333/1‑150 (MQ=255)
                                                                                               gTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGgcgc                                                                                                                                         <  2:767370/57‑1 (MQ=18)
                                                                                               gTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTGGAACCGCTGTTCAAGgcgc                                                                                                                                         >  1:767370/1‑57 (MQ=18)
                                                                                                tACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGc                                           <  2:1052360/150‑1 (MQ=255)
                                                                                                   gcgGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACCCGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGAt                                                    <  1:671298/138‑1 (MQ=255)
                                                                                                   gcgGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGAt                                                    >  2:671298/1‑138 (MQ=255)
                                                                                                   gcgGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCcgc                                        >  2:935958/1‑150 (MQ=255)
                                                                                                        aTGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTg                                   <  2:314516/150‑1 (MQ=255)
                                                                                                          gTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCCTGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGc                                 <  2:1041482/150‑1 (MQ=255)
                                                                                                             gCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTg                                   <  2:945734/145‑1 (MQ=255)
                                                                                                             gCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTg                                   >  1:945734/1‑145 (MQ=255)
                                                                                                                cACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTc                           >  2:769904/1‑150 (MQ=255)
                                                                                                                 aCGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTccc                                                               <  1:662952/113‑1 (MQ=255)
                                                                                                                 aCGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTccc                                                               >  2:662952/1‑113 (MQ=255)
                                                                                                                 aCGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCt                          <  1:360388/150‑1 (MQ=255)
                                                                                                                     ccGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAg                      >  1:620098/1‑150 (MQ=255)
                                                                                                                      cGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAgc                     <  1:1099087/150‑1 (MQ=255)
                                                                                                                           cgaGCGCCTGGAACCGCTGTTCAAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCgcggc                >  2:353370/3‑150 (MQ=255)
                                                                                                                           ccAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCgcggc                >  1:219222/1‑150 (MQ=255)
                                                                                                                            cAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATaa                                                  <  1:300359/115‑1 (MQ=255)
                                                                                                                            cAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATaa                                                  >  2:300359/1‑115 (MQ=255)
                                                                                                                              gcgcCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCa             >  2:479209/1‑150 (MQ=255)
                                                                                                                                gcCTGGAACCGCTGTTCAAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATg           >  1:317633/1‑150 (MQ=255)
                                                                                                                                  cTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATgtc         >  1:453695/1‑150 (MQ=255)
                                                                                                                                     gAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATgtcgtc      <  2:620098/150‑1 (MQ=255)
                                                                                                                                        cAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCtt   <  2:219222/150‑1 (MQ=255)
                                                                                                                                         aGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTTg  <  1:975334/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CGTTATGCACCATTTTCACGTAGTGCCCGGCGCCCGGTGGACCGGCATGGATGTAGCCGTGCTCGGCACGCTGGTGCTCGCCGCTGCGGCCATGGGTACGCGGGATGTCGCCCACGCCCGGCGCCAGCGCCTTGAACAGCGGTTCCAGGCGCTCGAACACGTCCTTTTCGCCGCCGATCATCATGCAGTAGCCGCGGTCCAGGCCCCACACGCCGCCGGAAGTCCCGACGTCCAGATAATGCAGGCCGCGCTTGGCCAGTTCTGCAGCGCGGCGCATGTCGTCCTTG  >  NC_002947/4556941‑4557227

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