Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,476,198 G→T 6.4% R85L (CGG→CTG)  PP_1290 → polysaccharide deacetylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,476,1980GT 6.4% 294.5 / 5.9 93R85L (CGG→CTG) PP_1290polysaccharide deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base G (48/39);  new base T (2/4);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTGCGTTCTATA  >  NC_002947/1476050‑1476347
                                                                                                                                                    |                                                                                                                                                     
gACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCgg                                                                                                                                                      >  2:1166861/1‑150 (MQ=255)
      ccGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGgcgcgc                                                                                                                                                >  1:987153/1‑150 (MQ=255)
      ccGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGgcgcgc                                                                                                                                                >  1:926554/1‑150 (MQ=255)
      ccGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGgcgcgc                                                                                                                                                >  2:68762/1‑150 (MQ=255)
       cGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCCGAATCTGgcgcacc                                                                                                                                               <  2:175015/150‑4 (MQ=255)
          ctcTGTTCGATGTGGCTTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCAtgt                                                                                                                                            <  2:876434/148‑1 (MQ=255)
          cgcTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCCGAATCTGgcgcaccttc                                                                                                                                            <  2:652742/150‑7 (MQ=255)
          cgcTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCAtgt                                                                                                                                            <  1:473183/150‑1 (MQ=255)
            cTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTgg                                                                                                                                          <  1:1065415/150‑1 (MQ=255)
              gTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGGGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCa                                                                                                                                        <  1:68762/150‑1 (MQ=255)
              gTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCa                                                                                                                                        >  2:2418443/1‑150 (MQ=255)
               ttCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAATCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCCGAATCTGgcgca                                                                                                                                                 >  1:768760/1‑139 (MQ=255)
               ttCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCCGAATCTGgcgca                                                                                                                                                 <  2:768760/140‑2 (MQ=255)
               ttCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCCGAATCTGgcgca                                                                                                                                                 >  1:1996499/1‑139 (MQ=255)
               ttCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCCGAATCTGgcgca                                                                                                                                                 <  2:1996499/140‑2 (MQ=255)
               ttCGATGTGGCCTCTCGGGCAGAAATCCGTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGAGGCAGATCAATCTGCAGATGGTGCGCCAGATTCGGGCGCGCACGTGGCAg                                                                                                                                       >  2:1249554/1‑150 (MQ=255)
                        gCCTCTTGGGCAGAAATCCTTTCGTTCGCTCAAGCGATGCATTAAAGCGAGTTGCTGGATGACGCAAGCCTCCTGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTgg                                                                                                                              <  1:1739400/150‑1 (MQ=255)
                        gCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTgg                                                                                                                              <  2:317176/150‑1 (MQ=255)
                         ccTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTggg                                                                                                                             >  1:893099/1‑150 (MQ=255)
                          ctctCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGa                                                                                                                            >  1:1201768/1‑150 (MQ=255)
                          ctctCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGa                                                                                                                            <  1:1166861/150‑1 (MQ=255)
                           tctcGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAg                                                                                                                           <  2:2070442/150‑1 (MQ=255)
                               gggCAGAAATCCTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTa                                                                                                                      <  1:2146901/150‑1 (MQ=255)
                                   aGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGa                                                                                                                            <  1:985692/141‑1 (MQ=255)
                                   aGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGa                                                                                                                            >  2:985692/1‑141 (MQ=255)
                                   aGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAgg                                                                                                                          >  2:1274281/1‑143 (MQ=255)
                                   aGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAgg                                                                                                                          <  1:1274281/143‑1 (MQ=255)
                                    gAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCaa                                                                                                                  <  1:441828/150‑1 (MQ=255)
                                     aaaTCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAAc                                                                                                                 >  2:1169113/1‑150 (MQ=255)
                                       aTCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAg                                                                                                               <  1:2324263/150‑1 (MQ=255)
                                           tttCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgc                                                                                                           >  1:509360/1‑150 (MQ=255)
                                           tttCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgc                                                                                                           >  2:1042586/1‑150 (MQ=255)
                                           tttCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGAATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgc                                                                                                           >  1:508859/1‑150 (MQ=255)
                                               gttcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCCGATCCCGCAGCAGAAGGTGCGCCAGAGTCGGTCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgcagcg                                                                                                       >  2:268984/1‑150 (MQ=255)
                                               gttcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGTTGCGCCAGATTCg                                                                                                                                                       <  2:1542455/102‑1 (MQ=255)
                                               gttcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTgg                                                                                                                                          <  2:516708/115‑1 (MQ=255)
                                               gttcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTgg                                                                                                                                          >  1:516708/1‑115 (MQ=255)
                                               gttcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTGTCAGCAGgagcagcg                                                                                                       >  1:1075203/1‑150 (MQ=255)
                                               gttcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgcagcg                                                                                                       >  2:1617165/1‑150 (MQ=255)
                                               ggtcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCg                                                                                                                                                       >  1:1542455/3‑102 (MQ=255)
                                                ttcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgcagcgc                                                                                                      >  2:2153765/1‑150 (MQ=255)
                                                 tcgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgcagcgca                                                                                                     <  2:742298/150‑1 (MQ=255)
                                                  cgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGcgcagcgcag                                                                                                    <  1:1567526/150‑1 (MQ=255)
                                                  cgCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCAACGGGGGCGGCAGGACCATTGGTTGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGcgcagcgcag                                                                                                    >  1:2123707/1‑150 (MQ=255)
                                                   gCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGt                                                                                                   >  1:339907/1‑150 (MQ=255)
                                                    cTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGtt                                                                                                  <  2:963777/150‑1 (MQ=255)
                                                     tCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTtg                                                                                                 >  1:2259109/1‑150 (MQ=255)
                                                     tCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTtg                                                                                                 >  2:2439029/1‑150 (MQ=255)
                                                        aGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGa                                                                                              <  2:759857/150‑1 (MQ=255)
                                                           gcTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCa                                                                                           >  1:1421118/1‑150 (MQ=255)
                                                           gcTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCa                                                                                           <  2:2082032/150‑1 (MQ=255)
                                                            cTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGGAAGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAg                                                                                          >  2:2298070/1‑150 (MQ=255)
                                                            cTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAg                                                                                          <  1:1684972/150‑1 (MQ=255)
                                                               cACGAAAGCGAGGTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGAAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGAt                                                                                       >  1:2267482/1‑150 (MQ=255)
                                                                aCGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATg                                                                                      <  2:1421118/150‑1 (MQ=255)
                                                                    aaGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCa                                                                                                                                        >  2:732800/1‑96 (MQ=255)
                                                                    aaGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCa                                                                                                                                        <  1:732800/96‑1 (MQ=255)
                                                                       cGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCAtt                                                                               <  2:987153/150‑1 (MQ=255)
                                                                          gTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATg                                                                                      <  2:311791/140‑1 (MQ=255)
                                                                          gTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATg                                                                                      >  1:311791/1‑140 (MQ=255)
                                                                          gTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGAAGATCCATCTGGAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTg                                                                            <  2:666210/150‑1 (MQ=255)
                                                                           ttGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGa                                                                                              <  1:101123/131‑1 (MQ=255)
                                                                           ttGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGa                                                                                              >  2:101123/1‑131 (MQ=255)
                                                                                    gACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGtt                                                                  >  2:1749701/1‑150 (MQ=255)
                                                                                       gaaaGCATGGCAGCCAAGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACg                                                               >  2:1731209/3‑150 (MQ=255)
                                                                                        ccaGCCGGGCCGCCAGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCAAGATTCGGGCGCGAAAGTGACAGCGGTTGTGGGAGGGGTAGCAACAGGCGCAGCGCAGTTGTGAGCAGGATGAAAAATACTGCTACCCGGTTACgg                                                              >  1:1306930/3‑149 (MQ=255)
                                                                                               ggcggcCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgca                                                                                                          >  2:810344/1‑99 (MQ=255)
                                                                                               ggcggcCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgca                                                                                                          <  1:810344/99‑1 (MQ=255)
                                                                                                     ccGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGtt                                                 <  1:2153765/150‑1 (MQ=255)
                                                                                                          tGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGTGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGca                                            <  2:1201768/150‑1 (MQ=255)
                                                                                                          tGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGca                                            <  1:1617165/150‑1 (MQ=255)
                                                                                                             ggTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGcacac                                         >  1:1760092/1‑150 (MQ=255)
                                                                                                               ttGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACaa                                       >  2:1337692/1‑150 (MQ=255)
                                                                                                                 gCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAgc                                     <  1:1169113/150‑1 (MQ=255)
                                                                                                                   ggCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTtgt                                                                                                >  1:2205830/1‑89 (MQ=255)
                                                                                                                   ggCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTtgt                                                                                                <  2:2205830/89‑1 (MQ=255)
                                                                                                                   ggCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAgcgg                                   >  2:383371/1‑150 (MQ=255)
                                                                                                                      aGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGGGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGa                                <  2:202417/150‑1 (MQ=255)
                                                                                                                       gATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTTGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGAc                               <  2:1306930/150‑1 (MQ=255)
                                                                                                                          cAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGTGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGtt                            >  2:2199256/1‑150 (MQ=255)
                                                                                                                           aaTCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCATGATGCATCATGCTTCTACCGGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTc                           <  2:926554/150‑1 (MQ=255)
                                                                                                                            aTCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCGGCGGTTGGGGTAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGGGGCGACGTTCg                          >  1:1871798/1‑150 (MQ=255)
                                                                                                                            aTCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCg                          >  2:286464/1‑150 (MQ=255)
                                                                                                                             tCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGc                         >  2:29388/1‑150 (MQ=255)
                                                                                                                             tCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGc                         >  1:2228557/1‑150 (MQ=255)
                                                                                                                               tGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTg                       >  2:956131/1‑150 (MQ=255)
                                                                                                                                 cAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGcc                     >  2:132300/1‑150 (MQ=255)
                                                                                                                                 cAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTACTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGAGGCGACGTTCGCTGcc                     >  2:1216307/1‑150 (MQ=255)
                                                                                                                                   gAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGGGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGa                   >  2:214570/1‑150 (MQ=255)
                                                                                                                                          ttccAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGt            <  1:132300/148‑1 (MQ=255)
                                                                                                                                           gcCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTg           >  2:1016700/1‑150 (MQ=255)
                                                                                                                                            ccAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTGc          <  1:383371/150‑1 (MQ=255)
                                                                                                                                             cAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGGTGCGCACACAAGCGGCGACGTTCCCTGCAGATCTCgtgccg         >  2:1715122/1‑145 (MQ=255)
                                                                                                                                               gATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTGCGtt       <  2:1760092/150‑1 (MQ=255)
                                                                                                                                               gATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTATGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTGCGtt       <  1:1749701/150‑1 (MQ=255)
                                                                                                                                                  tCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGAGGCGATGCGCAGGCAGGCGGCGAGGTTCGCGGCCGATGTCGGTGCGTTCta    >  1:1982715/1‑150 (MQ=255)
                                                                                                                                                    gggCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTGCGTTCtata  <  2:2259109/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTGCGTTCTATA  >  NC_002947/1476050‑1476347

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