Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,476,192 A→C 5.6% Q83P (CAG→CCG)  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,1920AC 5.6% 285.1 / 6.4 99Q83P (CAG→CCG) PP_1290polysaccharide deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base A (50/43);  new base C (2/4);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GCTGGACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTGCG  >  NC_002947/1476046‑1476340
                                                                                                                                                  |                                                                                                                                                    
gCTGGACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTAGCTGCTCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGAt                                                                                                                                                   <  1:1636908/150‑1 (MQ=255)
  tGGACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTc                                                                                                                                                 <  1:395007/150‑1 (MQ=255)
  tGGACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTc                                                                                                                                                 <  1:902645/150‑1 (MQ=255)
  tGGACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTc                                                                                                                                                 <  1:2102743/150‑1 (MQ=255)
    gACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCgg                                                                                                                                               >  2:1166861/1‑150 (MQ=255)
          ccGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGgcgcgc                                                                                                                                         >  1:926554/1‑150 (MQ=255)
          ccGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGgcgcgc                                                                                                                                         >  2:68762/1‑150 (MQ=255)
          ccGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGgcgcgc                                                                                                                                         >  1:987153/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)
               gcTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTc                                                                                                                                                 <  1:2416001/137‑1 (MQ=255)
               gcTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTc                                                                                                                                                 >  2:2416001/1‑137 (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                                                                                                                                          >  1:1996499/1‑139 (MQ=255)
                   ttCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCCGAATCTGgcgca                                                                                                                                          <  2:1996499/140‑2 (MQ=255)
                   ttCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCCGAATCTGgcgca                                                                                                                                          <  2:768760/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:1166861/150‑1 (MQ=255)
                              ctctCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGa                                                                                                                     >  1:1201768/1‑150 (MQ=255)
                               tctcGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAg                                                                                                                    <  2:2070442/150‑1 (MQ=255)
                                   gggCAGAAATCC‑TTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTa                                                                                                               <  1:2146901/150‑1 (MQ=255)
                                       aGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGa                                                                                                                     <  1:985692/141‑1 (MQ=255)
                                       aGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGa                                                                                                                     >  2:985692/1‑141 (MQ=255)
                                       aGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAgg                                                                                                                   <  1:1274281/143‑1 (MQ=255)
                                       aGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAgg                                                                                                                   >  2:1274281/1‑143 (MQ=255)
                                        gAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCaa                                                                                                           <  1:441828/150‑1 (MQ=255)
                                         aaaTCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAAc                                                                                                          >  2:1169113/1‑150 (MQ=255)
                                           aTCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAg                                                                                                        <  1:2324263/150‑1 (MQ=255)
                                               tttCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgc                                                                                                    >  2:1042586/1‑150 (MQ=255)
                                               tttCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGgcgc                                                                                                    >  1:509360/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                                                                                    <  2:2082032/150‑1 (MQ=255)
                                                               gcTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCa                                                                                    >  1:1421118/1‑150 (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                                                                                                                                 <  1:732800/96‑1 (MQ=255)
                                                                        aaGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCa                                                                                                                                 >  2:732800/1‑96 (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)
                                                                                                                                                  |                                                                                                                                                    
GCTGGACAGCCCGGCGCTGTTCGATGTGGCCTCTCGGGCAGAAATCCTTTCGTTCGCTCAAGCGCTGCACGAAAGCGAGTTGCTGGATGACGCAAGCCTGGCGGCCCGGCTGGGGTTGCGGCAGATCAATCTGCAGAAGGTGCGCCAGATTCGGGCGCGCATGTGGCAGCGGTTGTGGGAGGGGTATCAACAGGCGCAGCGCAGTTGTGAGCAGGATGCATCATTCTGCTACCCGGTTACGTCCATGGCGGAGTTGCGCACACAAGCGGCGACGTTCGCTGCCGATGTCGGTGCG  >  NC_002947/1476046‑1476340

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