Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,304,808 C→T 6.7% intergenic (+18/+14) 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,8080CT 6.7% 176.3 / 2.9 60intergenic (+18/+14)cysS/PP_2906cysteine‑‑tRNA ligase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (24/32);  new base T (2/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCGGCA  >  NC_002947/3304661‑3304957
                                                                                                                                                    |                                                                                                                                                     
ccAGGCGCGGCTGAAGGGGCGGGCGGACAAGAACTGGACTGAGTCTGACCGCATTCGCAATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGc                                                                                                                                                     >  2:516249/1‑150 (MQ=255)
 cAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCa                                                                                                                                                    <  1:772638/150‑1 (MQ=255)
  aGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAg                                                                                                                                                   >  1:598692/1‑150 (MQ=255)
     cgcgcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGccc                                                                                                                                                <  1:165528/150‑1 (MQ=255)
        gcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCGGTGCAGCCCCAg                                                                                                                                             <  2:128454/150‑1 (MQ=255)
             aaGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCtttt                                                                                                                                        <  2:242310/150‑1 (MQ=255)
                   cgTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGc                                                                                                                                  <  2:404978/150‑1 (MQ=255)
                     tGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCtt                                                                                                                                >  1:1258917/1‑150 (MQ=255)
                      gCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCCTa                                                                                                                               <  1:1281397/150‑1 (MQ=255)
                          aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTa                                                                                                                               <  1:594352/146‑1 (MQ=255)
                          aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTa                                                                                                                               >  2:312553/1‑146 (MQ=255)
                          aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTa                                                                                                                               <  1:312553/146‑1 (MQ=255)
                          aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTa                                                                                                                               >  2:594352/1‑146 (MQ=255)
                          aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAg                                                                                                                              <  1:1133806/147‑1 (MQ=255)
                          aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAg                                                                                                                              >  2:1133806/1‑147 (MQ=255)
                          aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAcc                                                                                                                           >  1:298388/1‑150 (MQ=255)
                               aaCTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGc                                                                                                                                  >  1:981522/1‑138 (MQ=255)
                               aaCTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGc                                                                                                                                  <  2:981522/138‑1 (MQ=255)
                                    ggCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCggg                                                                                                                 <  1:490219/150‑1 (MQ=255)
                                       tGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAg                                                                                                              <  1:740778/150‑1 (MQ=255)
                                         aGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTg                                                                                                            <  2:298388/150‑1 (MQ=255)
                                            cTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCaa                                                                                                         <  1:653705/150‑1 (MQ=255)
                                             tGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAAt                                                                                                        >  2:1147813/1‑150 (MQ=255)
                                               aCCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGCGTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATAc                                                                                                      >  1:922326/1‑150 (MQ=255)
                                               aCCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATAc                                                                                                      <  1:87838/150‑1 (MQ=255)
                                                     ttCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGAACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGTCAGTGCAATACCAGCCt                                                                                                <  2:158092/150‑1 (MQ=255)
                                                         cgATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGcc                                                                                            <  1:1147813/150‑1 (MQ=255)
                                                           aTCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGcccc                                                                                          <  1:378058/150‑1 (MQ=255)
                                                                cTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCGGTGCAGCCCCAg                                                                                                                                             >  1:1281629/1‑94 (MQ=255)
                                                                cTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCGGTGCAGCCCCAg                                                                                                                                             <  2:1281629/94‑1 (MQ=255)
                                                                       ccATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCt                                                                              >  2:459948/1‑150 (MQ=255)
                                                                       ccATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCt                                                                              <  2:618182/150‑1 (MQ=255)
                                                                             gtgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCaa                                                                                    <  2:976806/138‑1 (MQ=255)
                                                                             gtgtAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCaa                                                                                    >  1:976806/1‑138 (MQ=255)
                                                                                   tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGTC‑GCAGTAAGCGGGGGCGGCACCGCAGCCCAACCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtctc                                                                  <  1:287865/150‑1 (MQ=255)
                                                                                     cTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGTGGACACCCAACTCGCTTCCCCCAATCTCCa                                                                >  1:324907/1‑150 (MQ=255)
                                                                                           gaCAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGa                                                          >  1:1178035/1‑150 (MQ=255)
                                                                                             cAGCAAGGGCACGACTACCTGGCGTCTCACTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGaaa                                                        <  1:643994/150‑1 (MQ=255)
                                                                                              aGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGc                                                                                                  >  1:17639/1‑107 (MQ=255)
                                                                                              aGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGc                                                                                                  <  2:17639/107‑1 (MQ=255)
                                                                                                cAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAg                                                               >  2:1051788/1‑140 (MQ=255)
                                                                                                cAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAg                                                               <  1:1051788/140‑1 (MQ=255)
                                                                                                cAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCa                                                     >  1:890473/1‑150 (MQ=255)
                                                                                                   gggCACGACTACCTGGCGTCTGGGTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCACCCTGGCCCCACCCAAATCGGTTCCACAAGTCTCCAACCGGTAACCATGc                                                  <  2:472254/150‑1 (MQ=255)
                                                                                                          aCTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTcga                                           <  2:598692/150‑1 (MQ=255)
                                                                                                                     tCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTg                                                                                               >  2:860721/1‑87 (MQ=255)
                                                                                                                     tCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTg                                                                                               <  1:860721/87‑1 (MQ=255)
                                                                                                                       tGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGt                                              >  2:1255065/1‑134 (MQ=255)
                                                                                                                       tGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGt                                              <  1:1255065/134‑1 (MQ=255)
                                                                                                                         gCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                            <  2:890473/150‑1 (MQ=255)
                                                                                                                         gCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                            <  2:665497/150‑1 (MQ=255)
                                                                                                                          ctgactgaAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGAt                           <  2:1258917/150‑1 (MQ=255)
                                                                                                                             actgaAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGa                        >  2:1235517/1‑150 (MQ=255)
                                                                                                                                gaAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACaa                     >  2:762530/1‑150 (MQ=255)
                                                                                                                                 aaGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAc                    <  2:1178035/150‑1 (MQ=255)
                                                                                                                                  taaaagcTGGGGCTGCGGTGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAc                    >  1:698957/7‑150 (MQ=255)
                                                                                                                                             gctccgccGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGaa         <  1:193816/144‑1 (MQ=255)
                                                                                                                                             gCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGaa         >  2:101758/1‑150 (MQ=255)
                                                                                                                                             gCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGaa         >  1:1111354/1‑150 (MQ=255)
                                                                                                                                               tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGc       <  1:459948/150‑1 (MQ=255)
                                                                                                                                                 cACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCg     >  2:133484/1‑150 (MQ=255)
                                                                                                                                                    cGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCGGcc  >  2:95451/1‑149 (MQ=255)
                                                                                                                                                    cGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCGGCa  >  1:363619/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAA‑CTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCGGCA  >  NC_002947/3304661‑3304957

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