Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,772,788 C→G 7.3% intergenic (+71/‑30) PP_1580 → / → PP_1581 hypothetical protein/BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,772,7880CG 7.3% 207.8 / 6.6 69intergenic (+71/‑30)PP_1580/PP_1581hypothetical protein/BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base G (2/3);  total (33/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.28e-01

GGCATTTCCCAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGTCGATC  >  NC_002947/1772641‑1772929
                                                                                                                                                   |                                                                                                                                              
ggCATTTCCCAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCg                                                                                                                                               >  4:152783/1‑149 (MQ=255)
ggCATTTCCCAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCg                                                                                                                                               >  2:335360/1‑149 (MQ=255)
  cATTTCCCAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTg                                                                                                                                             <  6:285357/149‑1 (MQ=255)
   aTTTCCCAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGc                                                                                                                                            <  2:380692/149‑1 (MQ=255)
      tCCCAGAGATCGTAGTCGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCgg                                                                                                                                          <  2:354840/148‑1 (MQ=255)
      tCCCAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCggg                                                                                                                                         >  3:329778/1‑149 (MQ=255)
        ccAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCACGgggcgggc                                                                                                                                        <  2:62449/148‑1 (MQ=255)
         cAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCtt                                                                                                                                      >  5:381483/1‑149 (MQ=255)
         cAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCACGgggcgggcag                                                                                                                                      <  3:354985/149‑3 (MQ=255)
         cAGAGATCGTAGTAGGTGCGCCAATCATTGACTGGGCGATCAACAAGCAATCGCCGGAACCCACATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCtt                                                                                                                                      <  5:387979/149‑1 (MQ=255)
          agagATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCACGgggcgggcagg                                                                                                                                     <  8:123061/149‑4 (MQ=255)
                   agtagGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCt                                                                                                                            >  4:48173/1‑149 (MQ=255)
                    gtagGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCt                                                                                                                            <  1:335360/148‑1 (MQ=255)
                    gtagGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTg                                                                                                                           <  4:140420/149‑1 (MQ=255)
                               aaTCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATg                                                                                                                >  2:183661/1‑149 (MQ=255)
                                 tCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATgcg                                                                                                              <  6:30379/149‑1 (MQ=255)
                                     tgattgGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTa                                                                                                           <  4:329778/148‑1 (MQ=255)
                                         tgGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTacaaca                                                                                                      >  8:313031/1‑149 (MQ=255)
                                            gCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGTTCCCATGCGCTACAACACGa                                                                                                   <  5:290287/149‑1 (MQ=255)
                                            gCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGCTCCCATGCGCTACAACACg                                                                                                    >  7:105005/1‑148 (MQ=255)
                                                     aaGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAAccccg                                                                                          >  3:275167/1‑149 (MQ=255)
                                                     aaGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAAccccg                                                                                          >  5:54256/1‑149 (MQ=255)
                                                         aaTGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAAccccgcccc                                                                                      >  3:283742/1‑149 (MQ=255)
                                                             gCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGag                                                                                   >  3:14604/1‑148 (MQ=255)
                                                                ggAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGagcagc                                                                               <  6:216149/149‑1 (MQ=255)
                                                                     ccAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATCCGCTAAAACACGAATAACCCCGCACCGAGCAGCGAccc                                                                          >  7:25814/1‑149 (MQ=255)
                                                                       aaaTAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGAccc                                                                          >  4:193694/1‑147 (MQ=255)
                                                                       aaaTAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCga                                                                        <  4:275167/149‑1 (MQ=255)
                                                                                  ccGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCaa                                                             <  1:183661/148‑1 (MQ=255)
                                                                                  ccGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAAc                                                            >  8:142204/1‑149 (MQ=255)
                                                                                  ccGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAAc                                                            <  5:439936/149‑1 (MQ=255)
                                                                                  ccGCAGAGCCGAACCAGCAACGAGATTCAGCCGGTCAGTCCGGAGCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCAGGAGATCCCCTGCGCTACAACAGGAATAGCCGCGCCCCGAGCGGCGACCCGAGAGACTTCAAc                                                            >  4:421921/1‑149 (MQ=255)
                                                                                       tAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATaa                                                       <  4:283742/149‑1 (MQ=255)
                                                                                          ccGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAActc                                                    >  7:192930/1‑149 (MQ=255)
                                                                                           cGAAGCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAActct                                                   <  3:421921/149‑1 (MQ=255)
                                                                                            gAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAActctc                                                  <  4:14604/149‑1 (MQ=255)
                                                                                                cAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAAc                                                                                              >  6:236786/1‑102 (MQ=255)
                                                                                                cAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAAc                                                                                              <  5:236786/102‑1 (MQ=255)
                                                                                                 aGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAgcg                                             >  5:86319/1‑149 (MQ=255)
                                                                                                    aaCGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCCACGCTACCTCTCTAGCGCtt                                          >  2:11253/1‑149 (MQ=255)
                                                                                                     aCGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTGCAACGATAAGTCTCTAGCGCtt                                          >  2:95964/1‑148 (MQ=255)
                                                                                                       gagaTTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTtgat                                       <  2:52663/149‑1 (MQ=255)
                                                                                                         gaTTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTtgatga                                     <  8:281494/149‑1 (MQ=255)
                                                                                                              aGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGc                                 <  7:397949/148‑1 (MQ=255)
                                                                                                               gCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCAt                               <  7:142204/149‑1 (MQ=255)
                                                                                                                 cGGTCAATCCGGATCCTATCGCGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATCCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGCCTTCAACGATAACTCTCTAGACCTAGCAGAGGGCATgg                             >  4:304273/1‑149 (MQ=255)
                                                                                                                  ggTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGc                            <  7:313031/149‑1 (MQ=255)
                                                                                                                        tCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCt                      <  6:54256/149‑1 (MQ=255)
                                                                                                                          cGGATCCTATCGAGCCTGCCCGCCCCGTGCTGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTgg                    <  6:86319/149‑1 (MQ=255)
                                                                                                                           ggATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACCTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTgg                    <  1:95964/149‑1 (MQ=255)
                                                                                                                            gATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGaa                  <  3:152783/149‑1 (MQ=255)
                                                                                                                               ccTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGaa                  >  1:209341/1‑146 (MQ=255)
                                                                                                                               ccTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGaa                  <  2:209341/146‑1 (MQ=255)
                                                                                                                                 tATCGAGGCTGCCGGCCCCGTGCGTGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGacac             <  8:192930/149‑1 (MQ=255)
                                                                                                                                 tATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGacac             >  3:280487/1‑149 (MQ=255)
                                                                                                                                 tATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGacac             >  5:78282/1‑149 (MQ=255)
                                                                                                                                  aTCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACt            >  1:67623/1‑149 (MQ=255)
                                                                                                                                     aaaaaagcccgcacGGGGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACtttt         >  5:167798/7‑149 (MQ=255)
                                                                                                                                      aGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTg        >  8:31469/1‑149 (MQ=255)
                                                                                                                                       gCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACtt           <  3:48173/145‑1 (MQ=255)
                                                                                                                                       gCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGt       <  7:391606/149‑1 (MQ=255)
                                                                                                                                        ccTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGTc      >  7:389379/1‑149 (MQ=255)
                                                                                                                                         cTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGTc      >  5:172360/1‑148 (MQ=255)
                                                                                                                                         aagcccgcacGGGGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGTc      >  2:49222/3‑148 (MQ=255)
                                                                                                                                          tGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGt       <  3:188652/146‑1 (MQ=255)
                                                                                                                                          tGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGt       >  4:188652/1‑146 (MQ=255)
                                                                                                                                            cccgcccCGTGCGTGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCTCCCCGAGCAGCGACCCGAGAGACTTCAAGGATAACTCTCTAGCGGTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGTCGATc  <  8:25814/149‑1 (MQ=255)
                                                                                                                                               gcccCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAActct                                                   >  2:89985/1‑97 (MQ=255)
                                                                                                                                               gcccCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAActct                                                   <  1:89985/97‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GGCATTTCCCAGAGATCGTAGTAGGTGCGCCAATCATTGATTGGGCGATCAACAAGGAATGGCCGGAACCCAAATAGAACTGCCGCATAGCCGAACCAGCAACGAGATTCAGCCGGTCAATCCGGATCCTATCGAGCCTGCCCGCCCCGTGCGGGCTTTTTTGTGCCTGGAGATCCCATGCGCTACAACACGAATAACCCCGCCCCGAGCAGCGACCCGAGAGACTTCAACGATAACTCTCTAGCGCTTGATGAGGGCATGGCCAGTCTGGAAGACACTTTTGTCGATC  >  NC_002947/1772641‑1772929

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: