Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,610,105 C→G G154A (GGC→GCC)  PP_4929 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,610,1050CG98.7% 237.0 / ‑4.5 74G154A (GGC→GCC) PP_4929LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (35/38);  total (35/39)
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.20e-01

TGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGC  >  NC_002947/5609961‑5610250
                                                                                                                                                |                                                                                                                                                 
tGCCGGTGGCGGGCCAATTGTTGGCGGGAGCCCCGGCCTTCGCCGGCCAGCGGGGGGCCGACCGCGCACACCGTGACCCTCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGCAGGAAGAACAGCGCCAgg                                                                                                                                              >  1:1111965/1‑149 (MQ=255)
    ggTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCgg                                                                                                                                          <  1:652817/150‑1 (MQ=255)
         cGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGgggcga                                                                                                                                     <  2:412000/150‑2 (MQ=255)
             ccAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCt                                                                                                                                 >  2:485685/1‑150 (MQ=255)
               aaTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTcc                                                                                                                               <  1:642269/150‑1 (MQ=255)
                    ttgGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAg                                                                                                                          >  2:542130/1‑150 (MQ=255)
                       gCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGACGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGAtt                                                                                                                       >  2:395065/1‑150 (MQ=255)
                         gggTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTgc                                                                                                                     >  1:241512/1‑150 (MQ=255)
                          ggTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTgcg                                                                                                                    >  1:368728/1‑150 (MQ=255)
                           gTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTgcgc                                                                                                                   >  1:830678/1‑150 (MQ=255)
                             gACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTgcgcgc                                                                                                                 <  2:830678/150‑1 (MQ=255)
                              aCCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGcc                                                                                                                >  1:1148560/1‑150 (MQ=255)
                               cccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCca                                                                                                               >  1:712375/1‑150 (MQ=255)
                               cccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTCTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGccc                                                                                                               >  2:412068/1‑149 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTccac                                                                                                                             <  2:608012/135‑1 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTccac                                                                                                                             >  1:608012/1‑135 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTcc                                                                                                                               >  1:974643/1‑133 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTcc                                                                                                                               <  2:974643/133‑1 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCt                                                                                                                        >  1:1115960/1‑140 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCt                                                                                                                        <  2:1115960/140‑1 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTgcgcgc                                                                                                                 >  2:1191699/1‑147 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTgcgcgc                                                                                                                 <  1:1191699/147‑1 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCcac                                                                                                              <  2:956865/150‑1 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCcac                                                                                                              >  2:602496/1‑150 (MQ=255)
                                ccGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCCGCGCACCGCCATCCTCGCAAATGTCCAAAGAAGCCGGTAGGCTTTCGTGGCGGAAGCAAAGGGCGAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCcac                                                                                                              <  1:602496/150‑1 (MQ=255)
                                 cGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGNAGAACAGGGNCAGATCGGNGCGCCGCTCCACCAGCTTGCGCGCcaca                                                                                                             <  1:395065/150‑1 (MQ=255)
                                      ttCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAACCCGGGATGCTTTCGTGGAGCAAGAACAGGGCCAGACAGGCGTTCCGCTCCACCTTCGTGCGCGCCACATCGcc                                                                                                        >  2:457211/1‑150 (MQ=255)
                                          cccgccAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTgg                                                                                                    <  2:1022103/150‑1 (MQ=255)
                                                aGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCg                                                                                              >  2:854857/1‑150 (MQ=255)
                                                          cGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGtt                                                                                    <  1:458335/150‑1 (MQ=255)
                                                                      ccGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGa                                                                        >  2:760574/1‑150 (MQ=255)
                                                                        gTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAAc                                                                      <  1:542130/150‑1 (MQ=255)
                                                                        gTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAAc                                                                      <  2:712375/150‑1 (MQ=255)
                                                                        gTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAAc                                                                      >  1:668949/1‑150 (MQ=255)
                                                                           aCCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACggg                                                                   >  2:392748/1‑150 (MQ=255)
                                                                           aCCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACggg                                                                   >  2:725270/1‑150 (MQ=255)
                                                                            ccATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGcc                                                                                                        >  2:15483/1‑112 (MQ=255)
                                                                            ccATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGcc                                                                                                        <  1:15483/112‑1 (MQ=255)
                                                                            ccATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGt                                                                  <  2:857664/150‑1 (MQ=255)
                                                                                  ccACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGc                                                            <  1:642661/150‑1 (MQ=255)
                                                                                   cACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGc                                                            <  1:301613/149‑1 (MQ=255)
                                                                                   cACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGc                                                            >  2:301613/1‑149 (MQ=255)
                                                                                   cACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCt                                                           >  2:1211488/1‑150 (MQ=255)
                                                                                   cACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCt                                                           >  1:511822/1‑150 (MQ=255)
                                                                                            cAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTtc                                                  >  1:1170287/1‑150 (MQ=255)
                                                                                            cAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTtc                                                  >  2:233229/1‑150 (MQ=255)
                                                                                               cGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTtcgtc                                               <  1:392748/150‑1 (MQ=255)
                                                                                                  aCGGCGCTCCAAGTAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCGCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAg                                            <  1:615152/150‑1 (MQ=255)
                                                                                                  aCGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAg                                            <  2:253893/150‑1 (MQ=255)
                                                                                                  aCGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAg                                            <  2:368728/150‑1 (MQ=255)
                                                                                                  aCGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAg                                            <  1:167477/150‑1 (MQ=255)
                                                                                                   cGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGCCCAGGGACATTTTGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAgg                                           <  1:85038/150‑1 (MQ=255)
                                                                                                   cGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGGGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAgg                                           <  1:85037/150‑1 (MQ=255)
                                                                                                   cGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAgg                                           <  2:241512/150‑1 (MQ=255)
                                                                                                        cTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTc                                      <  2:1148560/150‑1 (MQ=255)
                                                                                                         tCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGAGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCg                                     >  2:647470/1‑150 (MQ=255)
                                                                                                           cAAGGAAGCCTGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGa                                    >  2:1158589/1‑149 (MQ=255)
                                                                                                           cAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACggg                                                                   >  2:988323/1‑118 (MQ=255)
                                                                                                           cAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACggg                                                                   <  1:988323/118‑1 (MQ=255)
                                                                                                           cAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGa                                    <  1:1158589/149‑1 (MQ=255)
                                                                                                            aaGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATc                                  <  1:1178699/150‑1 (MQ=255)
                                                                                                               gAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCa                                                                           <  2:1241312/106‑1 (MQ=255)
                                                                                                               gAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCa                                                                           >  1:1241312/1‑106 (MQ=255)
                                                                                                                   ccGGGATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCt                           <  1:725270/150‑1 (MQ=255)
                                                                                                                     gggATGCTTTCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGa                         <  1:819664/150‑1 (MQ=255)
                                                                                                                            tttCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATa                       >  2:68007/1‑145 (MQ=255)
                                                                                                                            tttCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATa                       <  1:68007/145‑1 (MQ=255)
                                                                                                                            tttCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCa                  >  2:1188834/1‑150 (MQ=255)
                                                                                                                            tttCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCa                  >  1:460248/1‑150 (MQ=255)
                                                                                                                              tCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATg                <  2:1170287/150‑1 (MQ=255)
                                                                                                                              tCGTGGAGGAAGAACAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATg                <  1:1211488/150‑1 (MQ=255)
                                                                                                                                   gaggaAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCtta           <  2:813247/150‑2 (MQ=255)
                                                                                                                                      gaagaaCAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTc        <  2:511822/150‑1 (MQ=255)
                                                                                                                                            cAGGGCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTgcgc  >  2:593687/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTCGCCAGTTCGTCGAGGCTGTCGATCACCGGCTGATACGGCATGGCTTCGTCCTGCGC  >  NC_002947/5609961‑5610250

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: ‑