Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,916,384 T→G 6.2% intergenic (+86/‑29) 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,3840TG 6.2% 199.3 / 4.4 65intergenic (+86/‑29)PP_1715/PP_1716hypothetical protein/Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base G (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

GGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAG  >  NC_002947/1916244‑1916532
                                                                                                                                            |                                                                                                                                                    
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                                                                                                                         >  2:95470/1‑149 (MQ=255)
                     cccGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTaa                                                                                                                         >  6:254210/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                                                                                                                <  6:91567/149‑1 (MQ=255)
                              gcTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTc                                                                                                                <  1:233643/149‑1 (MQ=255)
                                  cGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAg                                                                                                              <  6:298718/147‑1 (MQ=255)
                                  cGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGt                                                                                                            <  4:180845/149‑1 (MQ=255)
                                        gAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAg                                                                                                      <  3:174779/149‑1 (MQ=255)
                                        gAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAg                                                                                                      <  3:161945/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         <  5:85447/149‑1 (MQ=255)
                                                                                                                                     gATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTc         <  8:156440/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)
                                                                                                                                            |                                                                                                                                                    
GGTCGAGCGCCTGCCGGCTACCCCGGAAGCGCTGCGGCTGGAGCGCAACAAGTAATGAAACCGGCCTGAATCTCACCTCAACTGCACAGACCTCGGAGTTGGTAATATAATTGGCATAAGGGTTATGCCTACGGATTACATGCTGTAATCCTTTTTAAAGGGGTTTGTAATGGGTATTCAGGTCATCAGCCGGGACGGTCAACCCGAGTACGCCGTAGTTCCCTGGGAGCAATACCAGGCGTTGCTCAGCGCTGCAGGGCAGGCACCAGCCGAGGCGGTTTCCAGTGAG  >  NC_002947/1916244‑1916532

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