Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,304,806 A→G 7.7% intergenic (+16/+16) cysS → / ← PP_2906 cysteine‑‑tRNA ligase/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,304,8060AG 7.7% 224.1 / 6.5 78intergenic (+16/+16)cysS/PP_2906cysteine‑‑tRNA ligase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base G (2/4);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

TGATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGC  >  NC_002947/3304657‑3304952
                                                                                                                                                     |                                                                                                                                                  
tGATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCa                                                                                                                                                    <  1:1700956/150‑1 (MQ=255)
  aTCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCAcc                                                                                                                                                  <  1:844688/150‑1 (MQ=255)
       ggCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGc                                                                                                                                             <  1:1934541/150‑1 (MQ=255)
        gcgcgcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGcc                                                                                                                                            <  1:2290573/150‑1 (MQ=255)
          gcgcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGccac                                                                                                                                          >  1:1157196/1‑148 (MQ=255)
            gcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAg                                                                                                                                        <  2:1091472/150‑1 (MQ=255)
            gcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAg                                                                                                                                        <  1:594773/150‑1 (MQ=255)
              cTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCt                                                                                                                                      >  2:326005/1‑150 (MQ=255)
                  gggCGCGTGCGCAATAGAACTAGGCATAGTCTTACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTa                                                                                                                                  <  2:1718492/149‑1 (MQ=255)
                  aGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTa                                                                                                                                  <  1:2133897/150‑1 (MQ=255)
                    gcccgTGCGGACAAGAACTGGGCTGAGTCTGACCGACTGCGCGATCAGCTTACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTTGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTAtt                                                                                                                                <  1:254765/147‑1 (MQ=255)
                              aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCAACGACACACCAGCTTTTATTTGCTTAGAcc                                                                                                                      >  1:1698738/1‑150 (MQ=255)
                                aaGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCg                                                                                                                    >  2:2082816/1‑150 (MQ=255)
                                          cTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCa                                                                                                          <  1:326005/150‑1 (MQ=255)
                                           tGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAg                                                                                                         <  2:117810/150‑1 (MQ=255)
                                             aGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTg                                                                                                       <  2:489393/150‑1 (MQ=255)
                                              gTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGc                                                                                                      <  1:1251784/150‑1 (MQ=255)
                                              gTCTGACCGCATTCGCGATCAGCTGACCGCCAAGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGc                                                                                                      >  2:1349696/1‑150 (MQ=255)
                                                cTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCaa                                                                                                    <  2:2179825/150‑1 (MQ=255)
                                                cTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCaa                                                                                                    <  1:1545030/150‑1 (MQ=255)
                                                   aCCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCa                                                                                                                                               <  1:2033287/104‑1 (MQ=255)
                                                   aCCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCa                                                                                                                                               >  2:2033287/1‑104 (MQ=255)
                                                       cATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGc                                                                                             <  1:2156382/150‑1 (MQ=255)
                                                         ttCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCt                                                                                           >  2:2194968/1‑150 (MQ=255)
                                                         ttCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCt                                                                                           >  1:446429/1‑150 (MQ=255)
                                                            gcgATCAGCTGACCGCCATGGGTGTAGTGCTGGCGGGCAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGc                                                                                        <  2:2409704/150‑1 (MQ=255)
                                                              gATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGccc                                                                                      >  1:1517572/1‑150 (MQ=255)
                                                                tCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCccca                                                                                    <  1:286244/150‑1 (MQ=255)
                                                                  aGCTGACCACCATGGGTGTAGTGCTGGCGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCcccacc                                                                                  <  2:1309202/150‑1 (MQ=255)
                                                                    cTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCcccaccca                                                                                >  1:1019919/1‑150 (MQ=255)
                                                                       cccGCCAGGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTTCCTGGCGTCTGGATGACTGAAGTAACTGGGGCTGCACCGCGGCCACGTCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACt                                                                             <  2:1157196/149‑1 (MQ=255)
                                                                          gCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGc                                                                          <  2:2422447/150‑1 (MQ=255)
                                                                              tGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAg                                                                                                                                        >  1:2413689/1‑84 (MQ=255)
                                                                              tGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAg                                                                                                                                        <  2:2413689/84‑1 (MQ=255)
                                                                                 gtgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCgg                                                                                                             <  1:585678/108‑1 (MQ=255)
                                                                                 gtgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCgg                                                                                                             >  2:585678/1‑108 (MQ=255)
                                                                                  tgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTAtt                                                                                                                                >  2:1785790/1‑88 (MQ=255)
                                                                                  tgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTAtt                                                                                                                                <  1:1785790/88‑1 (MQ=255)
                                                                                  tgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCa                                                                                                          <  1:933541/110‑1 (MQ=255)
                                                                                  tgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCa                                                                                                          >  2:933541/1‑110 (MQ=255)
                                                                                      gTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTccc                                                                     >  1:235829/1‑143 (MQ=255)
                                                                                      gTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTccc                                                                     <  2:235829/143‑1 (MQ=255)
                                                                                       tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             <  2:446429/150‑1 (MQ=255)
                                                                                       tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             >  1:706723/1‑150 (MQ=255)
                                                                                       tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             >  2:2392199/1‑150 (MQ=255)
                                                                                       tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             <  2:1563490/150‑1 (MQ=255)
                                                                                       tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             >  1:1988050/1‑150 (MQ=255)
                                                                                         cTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCa                                                           >  2:616995/1‑150 (MQ=255)
                                                                                          tGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTNCAg                                                          <  2:1517572/150‑1 (MQ=255)
                                                                                           ggaggaCAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGc                                                         >  2:307208/1‑150 (MQ=255)
                                                                                            gaggaCAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGcc                                                        >  1:1830488/1‑150 (MQ=255)
                                                                                                 cAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGaaa                                                   <  1:2392199/150‑1 (MQ=255)
                                                                                                   gCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAAcc                                                 >  2:2061404/1‑150 (MQ=255)
                                                                                                      aGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATg                                              <  2:1988050/150‑1 (MQ=255)
                                                                                                       gggCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGc                                                                                                                                       >  2:20745/1‑60 (MQ=255)
                                                                                                       gggCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGc                                                                                                                                       <  1:20745/60‑1 (MQ=255)
                                                                                                       gggCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGCGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGc                                             <  1:21461/150‑1 (MQ=255)
                                                                                                         gCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAg                                           >  1:2034564/1‑150 (MQ=255)
                                                                                                         gCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAg                                           <  2:2307984/150‑1 (MQ=255)
                                                                                                           aCGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCa                                                                                                                                         >  1:2122643/1‑54 (MQ=25)
                                                                                                           aCGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCa                                                                                                                                         <  2:2122643/54‑1 (MQ=25)
                                                                                                           aCGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTCTTTGCTTAGACCGGCAAGCCTGCAGTTGAATACCATCCTTGCCCCACCCATCCCGCTTCCCCCAAGCTCCAGCCCGAAACCATGCAGGt                                         >  2:1705949/1‑150 (MQ=255)
                                                                                                            cGACTACCTGGCGTCTGGCTGACTGAAGTAACTGAGGCTGAACCGCGGCGCCATCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTc                                        >  1:2254402/1‑150 (MQ=255)
                                                                                                                    tGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATcg                                <  1:1349696/150‑1 (MQ=255)
                                                                                                                             gCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                       >  1:1018318/1‑150 (MQ=255)
                                                                                                                             gCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                       >  2:352204/1‑150 (MQ=255)
                                                                                                                             gCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                       >  1:1328293/1‑150 (MQ=255)
                                                                                                                               tgactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                     <  1:2061404/150‑1 (MQ=255)
                                                                                                                               tgactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                     <  1:628223/150‑1 (MQ=255)
                                                                                                                                   tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCATAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCg                                                                           >  2:1583301/1‑92 (MQ=255)
                                                                                                                                   tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCATAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCg                                                                           <  1:1583301/92‑1 (MQ=255)
                                                                                                                                        tAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAGCCGCCACCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAcccc            >  2:1242825/1‑150 (MQ=255)
                                                                                                                                           cTGGGGCTGCGGTGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGc         <  2:647744/150‑1 (MQ=255)
                                                                                                                                           cTGGGGCTGCGGTGCAGCCCCAGCTTTTATTTGCTAAGCAGCGACAACGCGCGACACAAGACCAGCCTGGCGCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGc         >  1:647744/1‑150 (MQ=255)
                                                                                                                                                 cTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACaa                <  1:860809/137‑1 (MQ=255)
                                                                                                                                                 cTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACaa                >  2:860809/1‑137 (MQ=255)
                                                                                                                                                  tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGc  >  1:79534/1‑150 (MQ=255)
                                                                                                                                                  tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGc  >  1:447679/1‑150 (MQ=255)
                                                                                                                                                  tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATACAGGTCGACGATCGCCGCCACGATCGCCAACCCCAGCCCGAAGc  >  2:2009855/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGC  >  NC_002947/3304657‑3304952

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