Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,916,386 C→A 5.6% intergenic (+88/‑27) PP_1715 → / → PP_1716 hypothetical protein/Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,916,3860CA 5.6% 227.4 / 3.4 71intergenic (+88/‑27)PP_1715/PP_1716hypothetical protein/Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base A (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

GGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAGG  >  NC_002947/1916244‑1916533
                                                                                                                                              |                                                                                                                                                   
ggTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACGTCAACTGCACAGACGTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAAt                                                                                                                                               >  1:39868/1‑149 (MQ=255)
    gAGCGCCTTCCTGCTACCCCGGAAGCGCTGCGGCTGTAGCGCAACAAGTAATGAAACCTGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACAGCATGTAATccgt                                                                                                                                           <  6:156382/149‑3 (MQ=255)
    gAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCtt                                                                                                                                           >  5:298718/1‑149 (MQ=255)
    gAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACAGCATGTAATccgt                                                                                                                                           >  7:175866/1‑147 (MQ=255)
     aGCGCCTGCTGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCttt                                                                                                                                          >  5:91567/1‑149 (MQ=255)
     aGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTTGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCttt                                                                                                                                          <  5:259367/149‑1 (MQ=255)
     aGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACAGCATGTAATccgta                                                                                                                                          <  8:175866/149‑4 (MQ=255)
      gcgcCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCtttt                                                                                                                                         <  2:25934/149‑1 (MQ=255)
              cGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAgggg                                                                                                                                 <  6:239767/149‑1 (MQ=255)
                    ccccGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTg                                                                                                                             <  1:974/147‑1 (MQ=255)
                    ccccGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGt                                                                                                                            >  2:233643/1‑148 (MQ=255)
                     cccGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTaa                                                                                                                          >  6:254210/1‑149 (MQ=255)
                     cccGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTaa                                                                                                                          >  2:95470/1‑149 (MQ=255)
                     cccGGAAGCGCTGCGGCTGGAGCGCAACAACTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTaa                                                                                                                          >  1:30078/1‑149 (MQ=255)
                      ccGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAAt                                                                                                                         >  3:61635/1‑149 (MQ=255)
                        ggAAGCGCGGCGGCTGGAGCCCAACAAGTAATGAAACCGCCCTGAATCTCACCTCAACTGCACACACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATgg                                                                                                                       <  4:108888/149‑1 (MQ=255)
                           aGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTa                                                                                                                    >  4:161945/1‑149 (MQ=255)
                              gcTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTc                                                                                                                 <  1:233643/149‑1 (MQ=255)
                              gcTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTc                                                                                                                 <  6:91567/149‑1 (MQ=255)
                                  cGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAg                                                                                                               <  6:298718/147‑1 (MQ=255)
                                  cGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGt                                                                                                             <  4:180845/149‑1 (MQ=255)
                                        gAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAg                                                                                                       <  3:161945/149‑1 (MQ=255)
                                        gAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAg                                                                                                       <  3:174779/149‑1 (MQ=255)
                                               acaaGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGAc                                                                                                >  5:218932/1‑149 (MQ=255)
                                                  aGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTTGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGt                                                                                             <  5:360936/149‑1 (MQ=255)
                                                      aTGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAAc                                                                                         >  1:61764/1‑149 (MQ=255)
                                                          aaCCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTTGCATAAGGGTTATGCCTACTGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGa                                                                                     <  4:173759/149‑1 (MQ=255)
                                                          aaCCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGa                                                                                     <  1:95470/149‑1 (MQ=255)
                                                           aCCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAg                                                                                    <  6:218932/149‑1 (MQ=255)
                                                                ccTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGc                                                                               <  2:30078/149‑1 (MQ=255)
                                                                 cTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGcc                                                                              >  6:366032/1‑149 (MQ=255)
                                                                    aaTCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCg                                                                             <  5:366032/147‑1 (MQ=255)
                                                                        tcACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTc                                                                       >  3:324689/1‑149 (MQ=255)
                                                                        tcACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTc                                                                       >  1:148070/1‑149 (MQ=255)
                                                                         cACCTCAACTGCACAGACCTCGGATTTTGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTc                                                                       <  4:258667/148‑1 (MQ=255)
                                                                            cTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGc                                                                               >  8:272328/1‑137 (MQ=255)
                                                                            cTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGc                                                                               <  7:272328/137‑1 (MQ=255)
                                                                             tCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTgg                                                                  >  3:291414/1‑149 (MQ=255)
                                                                              cAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTgg                                                                  >  3:91312/1‑148 (MQ=255)
                                                                               aaCTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGa                                                                <  4:61635/149‑1 (MQ=255)
                                                                               aaCTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGa                                                                <  4:324689/149‑1 (MQ=255)
                                                                                      caGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTTTTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATAc                                                         <  5:348700/149‑1 (MQ=255)
                                                                                         aCCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTgg                                                                  <  8:2131/137‑1 (MQ=255)
                                                                                         aCCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTgg                                                                  >  7:2131/1‑137 (MQ=255)
                                                                                           cTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAAACCGAGTACGCCGTATTTCCCTGGGAGCAATACCAGGc                                                    >  7:163637/1‑149 (MQ=255)
                                                                                               gAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGtt                                                 >  4:238715/1‑148 (MQ=255)
                                                                                               gAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTg                                                <  5:254210/149‑1 (MQ=255)
                                                                                                  ttGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTc                                             <  4:91312/149‑1 (MQ=255)
                                                                                                         tataATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTg                                      >  6:85447/1‑149 (MQ=255)
                                                                                                          ataATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTg                                      >  5:380686/1‑148 (MQ=255)
                                                                                                           taATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTg                                      >  4:299860/1‑147 (MQ=255)
                                                                                                             aTTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTgcagg                                  >  7:156440/1‑149 (MQ=255)
                                                                                                              ttGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACACGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTgcaggg                                 >  1:37806/1‑149 (MQ=255)
                                                                                                                 gCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCa                                            <  1:163199/135‑1 (MQ=255)
                                                                                                                 gCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCa                                            >  2:163199/1‑135 (MQ=255)
                                                                                                                          ttATGCCTACGGATTACATGCTGTTATCCTTTTTAAAGGGGTTTGTAATGGGTATTAAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCACTGGTAGCAATACCAGGCGTTGCTCAGCGATTCAGTGCTGGTACCAGcc                     >  7:168259/1‑149 (MQ=255)
                                                                                                                             tGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAg                  >  8:10236/1‑149 (MQ=255)
                                                                                                                               ccTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGc                <  3:299860/149‑1 (MQ=255)
                                                                                                                                cTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCg               >  7:129775/1‑149 (MQ=255)
                                                                                                                                aaaaggATTACAGCATGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCg               >  6:190155/5‑149 (MQ=255)
                                                                                                                                   ccgaatACATGCTGTAATACTTTGTAAAGGGGTTTTTAATGGGTATTCAGGTCCTCAGCCGCGACGGTCCACCCGAGTACGCCGTAGTTCCCTGGGATCAATACAAGGCGTTGCTCAGCGCTGCAGGGCAGGTACCAGCCGAGGCGGtt            <  8:168259/144‑1 (MQ=255)
                                                                                                                                     gATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTc          <  8:156440/149‑1 (MQ=255)
                                                                                                                                     gATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTc          <  5:85447/149‑1 (MQ=255)
                                                                                                                                          cATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGt      >  3:242087/1‑148 (MQ=255)
                                                                                                                                           aTGCTGTAATCCTATTTACAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGa    <  2:39868/149‑1 (MQ=255)
                                                                                                                                            tGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAg   <  5:65428/149‑1 (MQ=255)
                                                                                                                                             gCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACTGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAg   <  6:212559/148‑1 (MQ=255)
                                                                                                                                             gCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAg       <  2:177425/144‑1 (MQ=255)
                                                                                                                                             gCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAg       >  1:177425/1‑144 (MQ=255)
                                                                                                                                             gCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAgg  >  3:300414/1‑149 (MQ=255)
                                                                                                                                             gCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAg   >  5:212559/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAGG  >  NC_002947/1916244‑1916533

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: