Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,501,098 C→T 6.6% P138P (CCG→CCA PP_2194 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,501,0980CT 6.6% 173.4 / 3.1 61P138P (CCG→CCAPP_2194LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (36/21);  new base T (2/2);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

CCGTACAGCACATGCCCCCTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAACGCCGTTGGC  >  NC_002947/2500952‑2501247
                                                                                                                                                  |                                                                                                                                                     
ccGTACAGCACATGCCCCCTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGt                                                                                                                                                    <  1:1364654/150‑1 (MQ=255)
       gCACATGCCCCCTGCCGCAATACAGCCGCTGTCGTTCGGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTgc                                                                                                                                             >  1:1325023/1‑149 (MQ=255)
               ccccTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCCGCTCCGCCGGCGGTGCACTGACAAGCTcc                                                                                                                                     >  2:122797/1‑150 (MQ=255)
               ccccTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTcc                                                                                                                                     <  2:481133/150‑1 (MQ=255)
               ccccTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTcc                                                                                                                                     >  1:178637/1‑150 (MQ=255)
                cccTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCa                                                                                                                                    <  1:571192/150‑1 (MQ=255)
                cccTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCa                                                                                                                                    >  2:1174486/1‑150 (MQ=255)
                 ccTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTg                                                                                                                                                   <  2:1617314/134‑1 (MQ=255)
                 ccTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTg                                                                                                                                                   >  1:1617314/1‑134 (MQ=255)
                   tGCCGCAATACAGCCGCCGTCGTGCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTGCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGtt                                                                                                                                 <  2:692233/150‑1 (MQ=255)
                      cGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGcc                                                                                                                              >  2:655178/1‑149 (MQ=255)
                      cGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCa                                                                                                                              <  1:1154078/150‑1 (MQ=255)
                              aGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCCTGCCCCATCCGCTGTTCCAGCAACAGGCGTTGCATCTTCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTcc                                                                                                                      >  1:774061/1‑150 (MQ=255)
                              aGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTcc                                                                                                                      >  2:490850/1‑150 (MQ=255)
                                        gTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGc                                                                                                            >  2:762645/1‑150 (MQ=255)
                                                gaCAAGTGCCAGTAATCAAAGGGCGCCTGGTTGACGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGcc                                                                                                    <  2:418372/149‑1 (MQ=255)
                                                aaCAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCAAATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGcc                                                                                                    >  1:1116605/1‑150 (MQ=255)
                                                 aCAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTACACCGATTCCTCTCGCCTGCGAAATGCCg                                                                                                   >  2:462688/1‑150 (MQ=255)
                                                         tGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGgcc                                                                                           >  2:1318369/1‑150 (MQ=255)
                                                               tCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCGAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGAc                                                                                                                                            >  2:1061054/1‑95 (MQ=255)
                                                               tCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGAc                                                                                                                                            <  1:1061054/95‑1 (MQ=255)
                                                                cAAAGGCCGCCTCGTTGCCGCAACAATAGCTGATCACATAGTCCCGCCGCTGTTCCAACACCAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCCCCGATTCCTCTCGCCTGCGCAATCCCGCCATGGCCGCGGCCa                                                                                    >  2:439880/1‑150 (MQ=255)
                                                                 aaaGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAg                                                                                   >  1:899926/1‑150 (MQ=255)
                                                                         ccTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGc                                                                           >  2:1462888/1‑150 (MQ=255)
                                                                          cTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGTGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCt                                                                          >  1:803895/1‑150 (MQ=255)
                                                                          cTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGTGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCt                                                                          >  1:804064/1‑150 (MQ=255)
                                                                             gTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCAc                                                                       >  2:1076285/1‑150 (MQ=255)
                                                                                  cGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGa                                                                                                                          <  2:65085/94‑1 (MQ=255)
                                                                                  cGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGa                                                                                                                          >  1:65085/1‑94 (MQ=255)
                                                                                  cGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGATCACGCACCGCCGCCACTGCACTGACAAGCTGCATGTCCACCGAGTCTTAGCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCt                                                                  >  2:184364/1‑150 (MQ=255)
                                                                                  cGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAAAAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCGATCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCt                                                                  >  2:184402/1‑150 (MQ=255)
                                                                                    gAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAgg                                                                            >  2:458566/1‑138 (MQ=255)
                                                                                    gAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAgg                                                                            <  1:458566/138‑1 (MQ=255)
                                                                                      aaaaTAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGTGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTc                                                              <  2:803895/150‑1 (MQ=255)
                                                                                      aaaaTAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGTGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTc                                                              <  2:804064/150‑1 (MQ=255)
                                                                                       aaaTAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTACACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTcc                                                             <  1:462688/150‑1 (MQ=255)
                                                                                          tAGCTGATGGCATAGTCCCAGAGGTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACACGCTCCAGGTGCACCGATTTCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGa                                                          <  1:1318369/150‑1 (MQ=255)
                                                                                          tAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGa                                                          <  1:122797/150‑1 (MQ=255)
                                                                                                 tGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCg                                                   >  1:248699/1‑150 (MQ=255)
                                                                                                    cATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCGGCCa                                                >  2:600624/1‑150 (MQ=255)
                                                                                                       aGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCAGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGc                                             <  1:1076285/150‑1 (MQ=255)
                                                                                                           cAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTAt                                         <  2:159318/150‑1 (MQ=255)
                                                                                                           cAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTAt                                         >  1:159318/1‑150 (MQ=255)
                                                                                                           cAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTAt                                         >  1:212754/1‑150 (MQ=255)
                                                                                                             gCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGa                                                                                                          >  1:1166599/1‑83 (MQ=255)
                                                                                                             gCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGa                                                                                                          <  2:1166599/83‑1 (MQ=255)
                                                                                                              ccGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGtt                                                                                                                                 >  2:980287/1‑59 (MQ=255)
                                                                                                              ccGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGtt                                                                                                                                 <  1:980287/59‑1 (MQ=255)
                                                                                                                gCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCTTGATTTTCCGCCATGCCTATGCGca                                    >  1:415588/1‑150 (MQ=255)
                                                                                                                 cTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGcac                                   >  1:889579/1‑150 (MQ=255)
                                                                                                                  tGTTCCAGCAACAGTCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGcaca                                  <  2:706980/150‑1 (MQ=255)
                                                                                                                    ttCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGt                                <  1:1174486/150‑1 (MQ=255)
                                                                                                                        aGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTcccc                            >  1:1078877/1‑150 (MQ=255)
                                                                                                                           aaCAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCcagc                         >  1:1466409/1‑150 (MQ=255)
                                                                                                                            aCAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCcagca                        >  1:982204/1‑150 (MQ=255)
                                                                                                                               gCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAgcc                     >  2:1535063/1‑150 (MQ=255)
                                                                                                                                    cccAGCTCCGCCCCCTGAGCACCGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCGATGGCCGCCGCCAGTTTCAGACTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCg                <  2:1116605/149‑1 (MQ=255)
                                                                                                                                        gCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCaa            >  2:288365/1‑150 (MQ=255)
                                                                                                                                         cTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAAc           >  1:285229/1‑150 (MQ=255)
                                                                                                                                            cgccgGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAACGcc        >  1:1260541/1‑150 (MQ=255)
                                                                                                                                              ccgGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAACGCCGt      <  1:1462888/150‑1 (MQ=255)
                                                                                                                                                ggcggTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAACGCCGTTg    <  2:415588/150‑1 (MQ=255)
                                                                                                                                                ggcggTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAACGCCGTTg    >  2:879792/1‑150 (MQ=255)
                                                                                                                                                 gcggTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAACGCCGTTgg   <  2:1078877/150‑1 (MQ=255)
                                                                                                                                                  cggTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAACGCCGTTGGc  <  2:774061/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCGTACAGCACATGCCCCCTGCCGCAATACAGCCGCTGTCGTTCCGCGAACAAGGGCTGGTAATCAAAGGCCGCCTGGTTGCCGGAAAAATAGCTGATGGCATAGTCCAGCCGCTGTTCCAGCAACAGCCGTTCCAGCTCCGCCGGCGGTGCACTGACAAGCTCCAGTTGCACCGATTCCTCTCGCCTGCGAAATGCCGCAATGGCCGCCGCCAGTTTCAGGCTCACCGCCTTGTCCTGATTTTCCGCCATGCCTATGCGCACAGTCCCCAGCAGCCGCCCGGCAACGCCGTTGGC  >  NC_002947/2500952‑2501247

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: