Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,159,523 T→C 6.8% F500S (TTC→TCC)  opuA → glycine betaine ABC transporter ATP‑binding protein/permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,159,5230TC 6.8% 213.6 / 3.0 73F500S (TTC→TCC) opuAglycine betaine ABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base C (3/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GTGCGCACCGTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCACGCCGCTGTTGGATGTCATGCA  >  NC_002947/3159375‑3159670
                                                                                                                                                    |                                                                                                                                                   
gTGCGCACCGTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTc                                                                                                                                                    <  2:477602/150‑1 (MQ=255)
        cGTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCTGCGCGCTGCTGGTCCGCTCGACTGCGTAAGCCCTGGGGTATATCGGCCTGTTCGGTTTcgg                                                                                                                                            <  2:2078357/150‑3 (MQ=255)
        cGTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTg                                                                                                                                            <  1:465691/150‑1 (MQ=255)
        cGTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTg                                                                                                                                            <  1:1411469/150‑1 (MQ=255)
         gTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTAGCTGGCCTGGCGCTCTGCGTGCGTGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTgg                                                                                                                                           >  2:73260/1‑150 (MQ=255)
         gTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTgg                                                                                                                                           <  1:2254121/150‑1 (MQ=255)
         gTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTgg                                                                                                                                           <  1:86014/150‑1 (MQ=255)
                         aaGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCCGGGAGCGTACGATTTcca                                                                                                                           >  2:1246617/1‑150 (MQ=255)
                          aGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTccac                                                                                                                          <  2:1491312/150‑1 (MQ=255)
                             tGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGGTCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCAt                                                                                                                       >  1:1914585/1‑150 (MQ=255)
                                       tggATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTc                                                                                                             >  2:1114650/1‑150 (MQ=255)
                                       tggATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCCGGGAGCGTACGATTTCCACCATCGCTCTGGTc                                                                                                             <  1:1246617/150‑1 (MQ=255)
                                           ttCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGtt                                                                                                                                                     <  2:556169/106‑1 (MQ=255)
                                           ttCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGtt                                                                                                                                                     >  1:556169/1‑106 (MQ=255)
                                           ttCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGCGCGTTCGCTTTGCTCTCTCGGTCTGCTCGgcc                                                                                                         >  2:1662877/1‑149 (MQ=255)
                                               gTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTTCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTc                                                                                                     <  1:856183/150‑1 (MQ=255)
                                               gTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTc                                                                                                     >  1:179511/1‑150 (MQ=255)
                                               gTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTc                                                                                                     >  1:68856/1‑150 (MQ=255)
                                                 ggCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGc                                                                                                                                       >  2:689902/1‑114 (MQ=255)
                                                 ggCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGGTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGc                                                                                                                                       <  1:689902/114‑1 (MQ=255)
                                                    cGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGtgc                                                                                                >  1:2073367/1‑150 (MQ=255)
                                                     ggTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGtgct                                                                                               <  1:873698/150‑1 (MQ=255)
                                                      gTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGtgctg                                                                                              <  2:1357349/150‑1 (MQ=255)
                                                      gTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGtgctg                                                                                              <  2:1671337/150‑1 (MQ=255)
                                                             cgcTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCgc                                                                                       <  1:1114650/150‑1 (MQ=255)
                                                                ttGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCgctgc                                                                                    >  2:203245/1‑150 (MQ=255)
                                                                ttGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCgctgc                                                                                    <  2:1890971/150‑1 (MQ=255)
                                                                        gTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGc                                                                            >  1:2016462/1‑150 (MQ=255)
                                                                         tCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCGGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCa                                                                           <  2:507691/150‑1 (MQ=255)
                                                                           tACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCAtt                                                                         >  2:1826783/1‑150 (MQ=255)
                                                                            aCCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTc                                                                        <  1:1826783/150‑1 (MQ=255)
                                                                             cctggcctggCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTccacc                                                                                                                         >  1:1461426/1‑100 (MQ=255)
                                                                             cctggcctggCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTccacc                                                                                                                         <  2:1461426/100‑1 (MQ=255)
                                                                             cctggcctggCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTcc                                                                       >  2:1495870/1‑150 (MQ=255)
                                                                              ctggcctggcGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTc                                                                                         <  1:51643/131‑1 (MQ=255)
                                                                              ctggcctggcGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTc                                                                                         >  2:51643/1‑131 (MQ=255)
                                                                              ctggcctggcGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTccc                                                                      <  1:1434455/150‑1 (MQ=255)
                                                                               tggcctggcGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCa                                                                     >  2:1769756/1‑150 (MQ=255)
                                                                                  cctggcGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCg                                                                  >  1:554770/1‑150 (MQ=255)
                                                                                    tggcGCTCCGCCTGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCAc                                                                    <  2:241826/146‑1 (MQ=255)
                                                                                    tggcGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCAc                                                                    >  1:241826/1‑146 (MQ=255)
                                                                                        gcTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTCCGGTTTCTgg                                                                                                                                           >  1:2108218/1‑71 (MQ=255)
                                                                                        gcTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTCCGGTTTCTgg                                                                                                                                           <  2:2108218/71‑1 (MQ=255)
                                                                                        gcTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTCCGGTTTCTgg                                                                                                                                           <  2:2108139/71‑1 (MQ=255)
                                                                                        gcTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTCCGGTTTCTgg                                                                                                                                           >  1:2108139/1‑71 (MQ=255)
                                                                                         cTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGTCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCAttct                                                           >  2:972077/1‑150 (MQ=255)
                                                                                            ctctggCGAGATCCGTGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATtcttct                                                        <  2:132628/146‑1 (MQ=255)
                                                                                            cgccgGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGtgctg                                                                                              >  2:753780/1‑112 (MQ=255)
                                                                                            cgccgGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGtgctg                                                                                              <  1:753780/112‑1 (MQ=255)
                                                                                                 gcgcGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTgg                                                      <  2:680861/147‑1 (MQ=255)
                                                                                                 gcgcGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTgg                                                      >  1:680861/1‑147 (MQ=255)
                                                                                                 gcgcGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTgagat                                                   >  1:132628/1‑146 (MQ=255)
                                                                                                  cgcgATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCaa                                                  >  1:1539315/1‑150 (MQ=255)
                                                                                                      aTCCTGGCTCTGACGACTGCGTTAGGCCTGGGGTATATCGGCCGGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCgg                                              <  1:2015205/150‑1 (MQ=255)
                                                                                                         cTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCggaaa                                           >  2:1799866/1‑147 (MQ=255)
                                                                                                         cTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGgcg                                           >  2:378671/1‑150 (MQ=255)
                                                                                                           ggCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCt                                         <  1:203245/150‑1 (MQ=255)
                                                                                                             cccTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGACCGGTTTCTGGGAGCGTACGATTTATACCATCGTTCTTGTGTTCGCGTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTgg                                       >  2:852653/1‑150 (MQ=255)
                                                                                                                 gacgacTGCGTTAGCCCTGGGGTATATCGTCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCg                                   <  1:972077/150‑1 (MQ=255)
                                                                                                                 gacgacTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGTGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCg                                   <  1:1662877/150‑1 (MQ=255)
                                                                                                                 gacgacTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCg                                   <  1:1769756/150‑1 (MQ=255)
                                                                                                                         cGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTgg                                                                                                               <  2:838394/66‑1 (MQ=255)
                                                                                                                         cGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTgg                                                                                                               >  1:838394/1‑66 (MQ=255)
                                                                                                                           ttAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTc                         >  1:1058827/1‑150 (MQ=255)
                                                                                                                             aGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCAc                       >  1:1511298/1‑150 (MQ=255)
                                                                                                                              gCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCAcg                      <  2:1539315/150‑1 (MQ=255)
                                                                                                                                ccTGGGGTATATCGGGCTGTAAAGTTTGTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACTGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCAcgcc                    >  2:735400/1‑150 (MQ=255)
                                                                                                                                       tataTCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCTCTTGGTCACGCCGCTGTTg             <  1:73260/150‑1 (MQ=255)
                                                                                                                                          atCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCACGCCGCTGTTGGAt          >  1:2329002/1‑150 (MQ=255)
                                                                                                                                            cGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCACGCCGCTGTTGGATGt        <  2:1914585/150‑1 (MQ=255)
                                                                                                                                              gCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTAAAGTGCATCTGTCGCTGCGTAGCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCACGCCGCTGTTGGATGTCa      >  1:1453552/1‑150 (MQ=255)
                                                                                                                                               ccTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATtcttct                                                        <  2:32200/99‑1 (MQ=255)
                                                                                                                                               ccTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATtcttct                                                        >  1:32200/1‑99 (MQ=255)
                                                                                                                                                 tGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTCCTGGCCGAGCGGCCGCTGGCTCGGCGCTTGGTCACGCCGCTGTTGGATGTCATGc   <  1:1799866/150‑1 (MQ=255)
                                                                                                                                                  gTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCACGCCGCTGTTGGATGTCATGCa  <  2:554770/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGCGCACCGTCATCAATGGCGTCGAAGGTGCGCTGGGCTGGATTCCGTGGCCGGTACCGGCGCTTGCATTGGTCTACCTGGCCTGGCGCTCCGCCGGCGCGATCCTGGCCCTGACGACTGCGTTAGCCCTGGGGTATATCGGCCTGTTCGGTTTCTGGGAGCGTACGATTTCCACCATCGCTCTGGTCGGCGCTTCGGTGCTGCTGTCGCTGCTGATCGGCATTCCCACCGGCATTCTTCTGGCCAAGCGGGCGCTGGCTCGGCGCTTGGTCACGCCGCTGTTGGATGTCATGCA  >  NC_002947/3159375‑3159670

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: