Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,203,427 A→G 5.6% S60G (AGC→GGC)  PP_2810 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,203,4270AG 5.6% 207.4 / 2.2 71S60G (AGC→GGC) PP_2810hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base G (2/2);  total (36/35)
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.53e-01

GTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTGGCACAGCATGCGAGCAA  >  NC_002947/3203283‑3203572
                                                                                                                                                |                                                                                                                                                 
gTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAgcgg                                                                                                                                               >  8:338269/1‑149 (MQ=255)
 tATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAgcggc                                                                                                                                              >  2:169309/1‑149 (MQ=255)
      gcccgccGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                                                                                         <  7:187272/149‑1 (MQ=255)
      gcccgccGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                                                                                         >  1:244071/1‑149 (MQ=255)
      gcccgccGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                                                                                         >  3:300598/1‑149 (MQ=255)
      acccgccGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                                                                                         <  2:349366/148‑1 (MQ=255)
           ccgTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGCCCGCGATCGGGTCGGTCCAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTGCGAAGGTGACGCCGCTGGGGGGCGAGCGCGCCGCCAGGTCTAGCCGCGATAGTCCAg                                                                                                                                    >  6:283525/1‑149 (MQ=255)
                ttttCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTgg                                                                                                                               <  5:364176/149‑1 (MQ=255)
                   tCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTggg                                                                                                                              >  7:241248/1‑147 (MQ=255)
                   tCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTggg                                                                                                                              <  8:241248/147‑1 (MQ=255)
                     aGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCTAGAACTGgcggcgc                                                                                                                                            >  4:62659/1‑130 (MQ=255)
                     aGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCTAGAACTGgcggcgc                                                                                                                                            <  3:62659/131‑2 (MQ=255)
                      ggCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCATGGGCACAg                                                                                                                         >  2:139879/1‑149 (MQ=255)
                       gCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGc                                                                                                                        >  3:249367/1‑149 (MQ=255)
                             gctgGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                                                                                         >  6:14531/1‑126 (MQ=255)
                             gctgGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTGCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGatat                                                                                                                                         <  5:14531/126‑1 (MQ=255)
                               tgGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGGTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGgcag                                                                                                                >  8:289895/1‑149 (MQ=255)
                                 gTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGgcag                                                                                                                >  8:287188/1‑147 (MQ=255)
                                           cGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAgg                                                                                                     <  4:300598/148‑1 (MQ=255)
                                              ccTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATg                                                                                                  >  5:15876/1‑148 (MQ=255)
                                                  gACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAg                                                                                                                     >  1:250612/1‑125 (MQ=255)
                                                  gACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAg                                                                                                                     <  2:250612/125‑1 (MQ=255)
                                                       cgATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGc                                                                                        <  1:252671/149‑1 (MQ=255)
                                                          tCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATg                                                                                                  <  2:269158/136‑1 (MQ=255)
                                                          tCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATg                                                                                                  >  1:269158/1‑136 (MQ=255)
                                                              ttCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAggg                                                                                                           >  5:186644/1‑123 (MQ=255)
                                                              ttCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAggg                                                                                                           <  6:186644/123‑1 (MQ=255)
                                                                cGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGcc                                                                                       >  8:242270/1‑141 (MQ=255)
                                                                cGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGcc                                                                                       <  7:242270/141‑1 (MQ=255)
                                                                cGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTAc                                                                                 >  8:289593/1‑147 (MQ=255)
                                                                 ggTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAg                                                                              <  7:287188/149‑1 (MQ=255)
                                                                  gTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGa                                                                                                <  2:304777/130‑1 (MQ=255)
                                                                  gTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGa                                                                                                >  1:304777/1‑130 (MQ=255)
                                                                  gTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAgg                                                                             >  4:320728/1‑149 (MQ=255)
                                                                  gTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAgg                                                                             >  1:210034/1‑149 (MQ=255)
                                                                   tACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAggg                                                                            <  7:289895/149‑1 (MQ=255)
                                                                   tACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAggg                                                                            <  3:69953/149‑1 (MQ=255)
                                                                      aGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAg                                                                         >  3:259191/1‑149 (MQ=255)
                                                                        gCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCa                                                                       >  1:289933/1‑149 (MQ=255)
                                                                         ctccgcCCCCTGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCAc                                                                      >  1:37009/4‑149 (MQ=255)
                                                                          tGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACt                                                                     >  7:209776/1‑149 (MQ=255)
                                                                             gcCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGa                                                                  <  3:320728/149‑1 (MQ=255)
                                                                              ccccAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTgg                                                                   >  5:60860/1‑147 (MQ=255)
                                                                              ccccAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGaa                                                                 <  7:289593/149‑1 (MQ=255)
                                                                               cccAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAg                                                                >  4:42183/1‑149 (MQ=255)
                                                                               cccAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAg                                                                >  3:56515/1‑149 (MQ=255)
                                                                               cccAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAg                                                                <  1:169309/149‑1 (MQ=255)
                                                                               cccAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCCAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAg                                                                >  2:205908/1‑149 (MQ=255)
                                                                                ccAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGt                                                               <  6:15876/149‑1 (MQ=255)
                                                                                  aGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTTGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTc                                                             <  1:139879/149‑1 (MQ=255)
                                                                                  aGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGtt                                                              >  7:7539/1‑148 (MQ=255)
                                                                                     gggCGGAGGCGCTGGGTTCGAAGTTGCCGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAATTCAAg                                                          <  2:244071/149‑1 (MQ=255)
                                                                                     gggCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAg                                                          <  4:249367/149‑1 (MQ=255)
                                                                                       gCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAggg                                                        <  2:210034/149‑1 (MQ=255)
                                                                                           aGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGCTGAt                                                    <  8:7539/149‑1 (MQ=255)
                                                                                               gcTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGc                                                >  1:143348/1‑149 (MQ=255)
                                                                                                     ttCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCctct                                          <  8:209776/149‑1 (MQ=255)
                                                                                                        gAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTac                                        <  2:37009/148‑1 (MQ=255)
                                                                                                        gAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTac                                        <  2:289933/148‑1 (MQ=255)
                                                                                                               aCGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGa                                >  2:197841/1‑149 (MQ=255)
                                                                                                                cgccgcTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGAc                               <  4:259191/149‑1 (MQ=255)
                                                                                                                 gccgcTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACt                              >  3:25411/1‑149 (MQ=255)
                                                                                                                    gcTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGa                           <  2:307251/149‑1 (MQ=255)
                                                                                                                             cGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTgg                  <  6:60860/149‑1 (MQ=255)
                                                                                                                                  tcgccgcTAGAACTGGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGa                           <  6:303562/134‑1 (MQ=255)
                                                                                                                                  tcgccgcTAGAACTGGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGa                           >  5:303562/2‑135 (MQ=255)
                                                                                                                                   cgccgcCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAAt                      >  2:54031/1‑139 (MQ=255)
                                                                                                                                   cgccgcCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAAt                      <  1:54031/139‑1 (MQ=255)
                                                                                                                                    gccgccAGTTCCAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTGGCACAGCa           <  1:205908/149‑1 (MQ=255)
                                                                                                                                               tAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTGGCACAGCATGCGAGCaa  >  7:305514/1‑147 (MQ=255)
                                                                                                                                               tAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTGGCACAGCATGCGAGCaa  <  8:305514/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTATCCGCCCGCCGTATTTTCAGGCAGCTGCTGGTTACCCTTGCGGCCTCGACCGCGATCTGTTCGGTACAGGCTGCGCCCCAGGGGGCGGAGGCGCTGGGTTCGAAGTTGACGCCGCTGGGGGGCGAGAGCGCCGCCAGTTCTAGCGGCGATATTCCAGCCTGGGCACAGCCAGGGCAGCAGGGGGCAGGATGGAGTTACGGCCAGGTACGAGGGGAGCACTGGAAGTTCAAGGGTGATAAGCCCCTCTACAACATCGACTCGAACAATGTGGCACAGCATGCGAGCAA  >  NC_002947/3203283‑3203572

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