Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,279,857 T→G 8.8% L338R (CTG→CGG)  PP_3750 → GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,279,8570TG 8.8% 186.5 / 9.8 68L338R (CTG→CGG) PP_3750GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base G (3/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

TACGACTGTGAGTTCCGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCACGCCGCGAGCGCATCC  >  NC_002947/4279709‑4280006
                                                                                                                                                    |                                                                                                                                                     
tacgactGTGAGTTCCGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTg                                                                                                                                                      >  1:287699/1‑150 (MQ=255)
tacgactGTGAGTTCCGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTg                                                                                                                                                      >  1:1983038/1‑150 (MQ=255)
tacgactGTGAGTTCCGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTg                                                                                                                                                      >  1:1317159/1‑150 (MQ=255)
 acgagtGGGGGTTCCGCTACCACGGGCGTGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGc                                                                                                                                                     >  2:32856/1‑150 (MQ=255)
     ctGTGAGTTCCGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGc                                                                                                                                                 >  1:1017489/1‑150 (MQ=255)
     ctGTGAGTTCCGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGc                                                                                                                                                 <  1:904855/150‑1 (MQ=255)
      tgtgAGTTCCGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGcc                                                                                                                                                >  1:322703/1‑150 (MQ=255)
       gtgAGTTCCGCTACCACGGGACAGCATCGGAAGCGGTGCAGCCTCTGGAGCGACATGGGCTGGTTGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCAGTACGGCAGGcgg                                                                                                                                               <  2:1129774/150‑3 (MQ=255)
              ccGCTACCACGGGCCCGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGTGCTGGTCGCCTACGTCGGGACCTTCCCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGt                                                                                                                                        <  2:542488/150‑1 (MQ=255)
              ccGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGt                                                                                                                                        <  2:2354707/150‑1 (MQ=255)
                   aCCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCa                                                                                                                                   <  2:1168173/150‑1 (MQ=255)
                    ccACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCaa                                                                                                                                  <  2:1965034/150‑1 (MQ=255)
                      aCGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGc                                                                                                                                >  2:2391806/1‑150 (MQ=255)
                         ggCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCAcc                                                                                                                             <  2:287699/150‑1 (MQ=255)
                              gCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGc                                                                                                                        >  1:381343/1‑150 (MQ=255)
                                cTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGc                                                                                                                                          >  1:205403/1‑130 (MQ=255)
                                cTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGc                                                                                                                                          <  2:205403/130‑1 (MQ=255)
                                cTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTGCTCCAAGACCTTGCTGCCTTTGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGa                                                                                                                      <  2:1771361/150‑1 (MQ=255)
                                 tCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGAc                                                                                                                     >  2:523071/1‑150 (MQ=255)
                                          cTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGc                                                                                                            <  2:1017489/150‑1 (MQ=255)
                                          cTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGc                                                                                                            <  1:523071/150‑1 (MQ=255)
                                               gCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTAc                                                                                                       <  1:2391806/150‑1 (MQ=255)
                                                     ggACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTCAGCGACTGGCACAGCCCTACCctgct                                                                                                 <  1:1596570/150‑1 (MQ=255)
                                                           gCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTg                                                                                           >  2:204128/1‑150 (MQ=255)
                                                                  cTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTg                                                                                    >  2:2399222/1‑150 (MQ=255)
                                                                   tGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTgg                                                                                   <  2:1983038/150‑1 (MQ=255)
                                                                     gTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGgcg                                                                                 >  1:1628565/1‑150 (MQ=255)
                                                                        gCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCAc                                                                                                                              <  2:1020662/102‑1 (MQ=255)
                                                                        gCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCAc                                                                                                                              >  1:1020662/1‑102 (MQ=255)
                                                                          cTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCtt                                                                            >  1:307396/1‑150 (MQ=255)
                                                                          cTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCtt                                                                            >  2:1026571/1‑150 (MQ=255)
                                                                          cTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCAGTACGGCAGGcgg                                                                                                                                               >  2:298991/1‑81 (MQ=255)
                                                                          cTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCAGTACGGCAGGcgg                                                                                                                                               <  1:298991/83‑3 (MQ=255)
                                                                           tACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCa                                                                                                                                                    >  1:1309318/1‑77 (MQ=255)
                                                                           tACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCa                                                                                                                                                    <  2:1309318/77‑1 (MQ=255)
                                                                              gTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCtt                                                                        <  2:1628565/150‑1 (MQ=255)
                                                                                cGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTgg                                                                      <  1:133044/150‑1 (MQ=255)
                                                                                  ggACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTCTGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGcc                                                                    >  2:590337/1‑149 (MQ=255)
                                                                                  ggACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTgg                                                                                          >  2:2347585/1‑128 (MQ=255)
                                                                                  ggACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTgg                                                                                          <  1:2347585/128‑1 (MQ=255)
                                                                                     ccTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAg                                                                 <  1:2016952/150‑1 (MQ=255)
                                                                                       ttCCCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGc                                                               <  1:410396/150‑1 (MQ=255)
                                                                                          tcCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCAc                                                            <  1:2399222/150‑1 (MQ=255)
                                                                                           ccAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCAcc                                                           >  2:739935/1‑150 (MQ=255)
                                                                                            cAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCt                                                          >  2:534907/1‑150 (MQ=255)
                                                                                            cAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCATTGGCGAAgagatcgg                                                          >  1:404252/1‑143 (MQ=255)
                                                                                                  cTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACaa                                                    <  2:472592/150‑1 (MQ=255)
                                                                                                   ttGCTGCCCGAGCTGCGCCTGGGCTATGGCCGATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAg                                                   >  1:674440/1‑150 (MQ=255)
                                                                                                   ttGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAg                                                   <  1:1026571/150‑1 (MQ=255)
                                                                                                    tgctgcCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCt                                                                                                                            >  2:1469189/1‑76 (MQ=255)
                                                                                                    tgctgcCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCt                                                                                                                            <  1:1469189/76‑1 (MQ=255)
                                                                                                    tgctgcCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGTCGAAGGCCACCTGAACAAGc                                                  <  2:1294218/150‑1 (MQ=255)
                                                                                                           cGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCGTCCTTGGCGAAGGCCACCTGCACAAGCACATccg                                           >  2:180539/1‑150 (MQ=255)
                                                                                                                     cTGGGCTATGCCGTATTACGGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCAc                                 <  2:1584895/150‑1 (MQ=255)
                                                                                                                      tGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACg                                >  2:855391/1‑150 (MQ=255)
                                                                                                                       gggCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGa                               <  2:2252497/150‑1 (MQ=255)
                                                                                                                           tATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTc                           >  1:1079832/1‑150 (MQ=255)
                                                                                                                             tGCCGTATTACCGCCTGCCCTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTa                         <  1:534907/150‑1 (MQ=255)
                                                                                                                               ccGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCCTAGCTAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACa                       <  1:1833435/150‑1 (MQ=255)
                                                                                                                                       ccGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTc                           >  1:1806424/1‑138 (MQ=255)
                                                                                                                                       ccGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTc                           <  2:1806424/138‑1 (MQ=255)
                                                                                                                                        ggCCTGCAGTACGGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGaa                                                                  >  2:371441/2‑98 (MQ=255)
                                                                                                                                        ggCCTGCAGTACGGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGaa                                                                  <  1:371441/97‑1 (MQ=255)
                                                                                                                                        ggCCTGCAGTACGGCAGGCCGCTTGCGTGGCAACGGACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCAcgcc              >  1:2371181/2‑150 (MQ=255)
                                                                                                                                            tGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCa                  >  1:1036999/1‑142 (MQ=255)
                                                                                                                                            tGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCa                  <  2:1036999/142‑1 (MQ=255)
                                                                                                                                              ccGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCGCGCCGCGAgc        >  1:1306840/1‑150 (MQ=255)
                                                                                                                                                   cTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCGCAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTACTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCACGCCGCGAGCGCATc   >  1:2168748/1‑150 (MQ=255)
                                                                                                                                                    tGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCACGCCGCGATCGCATcc  <  2:1079832/150‑1 (MQ=255)
                                                                                                                                                    tGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCACGCCGCGAGCGCATcc  >  2:2108822/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACGACTGTGAGTTCCGCTACCACGGGCCAGCCTCGGACGCGCTGCAGCGTCTGGACCGGCATGGGCTGGTCGCCTACGTCGGGACCTTCTCCAAGACCTTGCTGCCCGAGCTGCGCCTGGGCTATGCCGTATTACCGCCTGCCGTACTGCAGGCCGCTTGCGTGGCCAAGCACCTTAGCGACTGGCACAGCCCTACCCTGCTGCAGTGGGCGCTGGCGCGCTTCCTTGGCGAAGGCCACCTGAACAAGCACATCCGCCGCTGCCACGAGGTCTACAGCGCACGCCGCGAGCGCATCC  >  NC_002947/4279709‑4280006

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: