Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,304,807 C→G 7.9% intergenic (+17/+15) 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,8070CG 7.9% 234.7 / 10.2 77intergenic (+17/+15)cysS/PP_2906cysteine‑‑tRNA ligase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base G (2/4);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

ATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGC  >  NC_002947/3304659‑3304952
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                    <  1:1545030/150‑1 (MQ=255)
                                              cTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCaa                                                                                                    <  2:2179825/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                                                                                                          >  2:933541/1‑110 (MQ=255)
                                                                                tgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCa                                                                                                          <  1:933541/110‑1 (MQ=255)
                                                                                    gTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTccc                                                                     >  1:235829/1‑143 (MQ=255)
                                                                                    gTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTccc                                                                     <  2:235829/143‑1 (MQ=255)
                                                                                     tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             >  1:706723/1‑150 (MQ=255)
                                                                                     tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             >  2:2392199/1‑150 (MQ=255)
                                                                                     tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             <  2:446429/150‑1 (MQ=255)
                                                                                     tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             >  1:1988050/1‑150 (MQ=255)
                                                                                     tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                             <  2:1563490/150‑1 (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                                                                                                                                       <  1:20745/60‑1 (MQ=255)
                                                                                                     gggCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGc                                                                                                                                       >  2:20745/1‑60 (MQ=255)
                                                                                                     gggCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGCGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGc                                             <  1:21461/150‑1 (MQ=255)
                                                                                                       gCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAg                                           <  2:2307984/150‑1 (MQ=255)
                                                                                                       gCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAg                                           >  1:2034564/1‑150 (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                       >  2:352204/1‑150 (MQ=255)
                                                                                                                           gCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                       >  1:1328293/1‑150 (MQ=255)
                                                                                                                           gCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                       >  1:1018318/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                >  2:860809/1‑137 (MQ=255)
                                                                                                                                               cTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACaa                <  1:860809/137‑1 (MQ=255)
                                                                                                                                                tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGc  >  1:447679/1‑150 (MQ=255)
                                                                                                                                                tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGc  >  1:79534/1‑150 (MQ=255)
                                                                                                                                                tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATACAGGTCGACGATCGCCGCCACGATCGCCAACCCCAGCCCGAAGc  >  2:2009855/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGC  >  NC_002947/3304659‑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: ‑