Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,272,822 T→C 6.0% D24D (GAT→GAC PP_3743 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,272,8220TC 6.0% 351.7 / 8.0 116D24D (GAT→GACPP_3743hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (56/53);  new base C (2/5);  total (58/58)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

Only 100 of 120 total aligned reads displayed.

AACTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTATGTGGAGTCGGCGC  >  NC_002947/4272674‑4272970
                                                                                                                                                    |                                                                                                                                                    
aaCTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTGCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGATCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATg                                                                                                                                                     <  1:1965487/150‑1 (MQ=255)
  cTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCg                                                                                                                                                   >  2:518954/1‑150 (MQ=255)
     ccccAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGc                                                                                                                                                >  1:1384219/1‑150 (MQ=255)
      cccAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCg                                                                                                                                               >  1:1102391/1‑150 (MQ=255)
      cccAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCg                                                                                                                                               >  1:898309/1‑150 (MQ=255)
        cAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCATCGACGCGGGccag                                                                                                                                             <  2:2390261/150‑4 (MQ=255)
            gCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTCGCCCGCATCGACGCGGGccagggca                                                                                                                                         <  2:1950916/150‑8 (MQ=255)
            gCCCCCGCCGTCTTAGATACAGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTCGCCCGCATCGACGCGGGccagggca                                                                                                                                         <  2:771675/150‑8 (MQ=255)
              ccccGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCACGTCGATGCGGGCGACGTCGcc                                                                                                                                       >  2:1199252/1‑150 (MQ=255)
                        ttAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACCTACAACACCTACTGCTCAATTGCGCACTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCg                                                                                                                             >  1:1485909/1‑150 (MQ=255)
                        ttAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCg                                                                                                                             >  2:357699/1‑150 (MQ=255)
                        ttAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCg                                                                                                                             >  2:358021/1‑150 (MQ=255)
                          aGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGAt                                                                                                                                                      >  1:1260287/1‑123 (MQ=255)
                              aCTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGt                                                                                                                                           >  1:276813/1‑130 (MQ=255)
                              aCTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGt                                                                                                                                           <  2:276813/130‑1 (MQ=255)
                                  aGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGc                                                                                                                   <  2:1708316/150‑1 (MQ=255)
                                  aGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGc                                                                                                                   >  1:1008688/1‑150 (MQ=255)
                                    tGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGcc                                                                                                                              >  2:1338538/1‑137 (MQ=255)
                                      gTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCg                                                                                                               >  2:1057344/1‑150 (MQ=255)
                                      gTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCCGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCg                                                                                                               >  2:771423/1‑150 (MQ=255)
                                           tttCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTa                                                                                                          <  1:554928/150‑1 (MQ=255)
                                            ttCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTAt                                                                                                         >  1:1656723/1‑150 (MQ=255)
                                                 tgtgCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                                    <  1:2057068/150‑1 (MQ=255)
                                                 tgtgCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                                    <  2:1448829/150‑1 (MQ=255)
                                                 tgtgCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGAAAATTCTGTCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                                    <  2:1797473/150‑1 (MQ=255)
                                                 tgtgCTGCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                                    >  1:772956/1‑150 (MQ=255)
                                                   tgCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCt                                                                                                  >  2:1837662/1‑150 (MQ=255)
                                                     cTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTc                                                                                                >  1:897336/1‑150 (MQ=255)
                                                      ttCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCgg                                                                                                                                                  >  2:41368/1‑99 (MQ=255)
                                                      ttCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCgg                                                                                                                                                  <  1:41368/99‑1 (MQ=255)
                                                      ttCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCa                                                                                               >  1:1472568/1‑150 (MQ=255)
                                                           gCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgac                                                                                          <  2:1118259/150‑1 (MQ=255)
                                                           gCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgac                                                                                          >  1:1762235/1‑150 (MQ=255)
                                                           gCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCCGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgac                                                                                          <  1:771423/150‑1 (MQ=255)
                                                              aCTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgacgac                                                                                       <  1:1837662/150‑1 (MQ=255)
                                                                tGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGc                                                                                                                                                    <  2:2292379/87‑1 (MQ=255)
                                                                tGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGc                                                                                                                                                    >  1:2292379/1‑87 (MQ=255)
                                                                tGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACaa                                                                                     >  2:105859/1‑150 (MQ=255)
                                                                  gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGTTGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAAc                                                                                   <  1:743365/150‑1 (MQ=255)
                                                                  gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTc                                                                                                                >  2:1950352/1‑121 (MQ=255)
                                                                  gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                                    >  1:1387568/1‑133 (MQ=255)
                                                                  gagaACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAAc                                                                                   <  2:1762235/150‑1 (MQ=255)
                                                                     aaCACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGa                                                                                <  2:800418/150‑1 (MQ=255)
                                                                      acacGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGat                                                                               <  2:1310780/150‑1 (MQ=255)
                                                                      acacGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTACTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGat                                                                               <  2:1596799/150‑1 (MQ=255)
                                                                       cacGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGatc                                                                              <  1:583117/150‑1 (MQ=255)
                                                                       cacGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGatc                                                                              >  1:513651/1‑150 (MQ=255)
                                                                          gACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGatcatc                                                                           <  1:2373984/150‑1 (MQ=255)
                                                                          gACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTCGCCCGCATCGACGCGGGccag                                                                                                                                             <  1:1556131/84‑4 (MQ=38)
                                                                            cATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCaa                                                                                              <  1:1247339/129‑1 (MQ=255)
                                                                            cATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATACGACCTTCaa                                                                                              >  2:1247339/1‑129 (MQ=255)
                                                                                 gcCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGCTGATCCCTGccc                                                                    >  1:1566710/1‑150 (MQ=255)
                                                                                  ccACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCg                                                                   >  2:1290536/1‑150 (MQ=255)
                                                                                   cACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGt                                                                  >  2:1437881/1‑150 (MQ=255)
                                                                                     ccTACTGCTCGCCCCTTCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGACACCACGGCAGCCATTCCCCCTTCGGTGTTTTCCGTCGTCAAGTACGACAAACGGATCATCCCTGCCCGTGa                                                                <  1:1749215/150‑1 (MQ=255)
                                                                                     ccTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGgac                                                                                          <  2:2132374/124‑1 (MQ=255)
                                                                                     ccTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGa                                                                <  2:1472568/150‑1 (MQ=255)
                                                                                      cTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGAc                                                               <  1:1057344/150‑1 (MQ=255)
                                                                                         cTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                                    <  2:1340987/110‑1 (MQ=255)
                                                                                         cTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGAc                                                                                                    >  1:1340987/1‑110 (MQ=255)
                                                                                            cTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGa                                                                                >  1:817472/1‑127 (MQ=255)
                                                                                            cTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGa                                                                                <  2:817472/127‑1 (MQ=255)
                                                                                            cTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGa                                                                                >  2:283441/1‑127 (MQ=255)
                                                                                               gCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCa                                                      <  2:1384219/150‑1 (MQ=255)
                                                                                               gCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCa                                                      >  2:2094968/1‑150 (MQ=255)
                                                                                                 cccTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGga                                                                                           >  1:448123/1‑111 (MQ=255)
                                                                                                 cccTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGc                                                    <  2:1743562/150‑1 (MQ=255)
                                                                                                     gcgcTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTtc                                                <  2:405006/150‑1 (MQ=255)
                                                                                                        ctgctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAg                                                                                             <  1:2030963/102‑1 (MQ=255)
                                                                                                        ctgctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAg                                                                                             >  2:2030963/1‑102 (MQ=255)
                                                                                                        ctgctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTtcgtc                                             <  1:1544953/150‑1 (MQ=255)
                                                                                                          gctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGc                                           <  2:897336/150‑1 (MQ=255)
                                                                                                          gctgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGc                                           <  1:518954/150‑1 (MQ=255)
                                                                                                            tgctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCa                                         >  1:801144/1‑150 (MQ=255)
                                                                                                             gctCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAg                                        >  2:2088464/1‑150 (MQ=255)
                                                                                                                 cTGATGCCTGCCGGTGCTGCCCTGGCCCGAGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAgcggc                                    >  2:646927/1‑150 (MQ=255)
                                                                                                                  tGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAgcggcg                                   >  1:424777/1‑150 (MQ=255)
                                                                                                                   gATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAg                                                     >  2:785421/1‑131 (MQ=255)
                                                                                                                   gATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAg                                                     <  1:785421/131‑1 (MQ=255)
                                                                                                                    aTGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGt                                 >  1:882289/1‑150 (MQ=255)
                                                                                                                      gcctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGc                                                                                                              <  2:2376838/71‑1 (MQ=255)
                                                                                                                      gcctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGc                                                                                                              >  1:2376838/1‑71 (MQ=255)
                                                                                                                      acctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTtc                                                <  1:1609914/132‑1 (MQ=255)
                                                                                                                      acctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTtc                                                >  2:1609914/2‑133 (MQ=255)
                                                                                                                        ctgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCa                             >  2:851096/1‑150 (MQ=255)
                                                                                                                         tgccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTANCAt                            <  2:1509772/150‑1 (MQ=255)
                                                                                                                          gccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATc                           <  2:424777/150‑1 (MQ=255)
                                                                                                                          gccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATc                           <  1:358021/150‑1 (MQ=255)
                                                                                                                          gccGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATc                           <  1:357699/150‑1 (MQ=255)
                                                                                                                                  tgcCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGACGGTGCCATCACGGCCCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGgcgc                   <  2:62976/150‑1 (MQ=255)
                                                                                                                                    cccTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATcc                                                                          <  1:775062/93‑1 (MQ=255)
                                                                                                                                    cccTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCt                 >  1:2069515/1‑150 (MQ=255)
                                                                                                                                       tcgCCCGCATCGACGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTAtg              >  2:154809/3‑150 (MQ=255)
                                                                                                                                        ggCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTAtgt             <  1:2094968/150‑1 (MQ=255)
                                                                                                                                         gCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTAtgtg            >  1:2202312/1‑150 (MQ=255)
                                                                                                                                                tCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTATGTGGAGTCgg     <  1:59089/150‑1 (MQ=255)
                                                                                                                                                   aTNCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCNCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTATGTGGAGTCGgcgc  >  1:1792594/1‑150 (MQ=255)
                                                                                                                                                   aTGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTATGTGGAGTCGgcgc  <  2:237464/150‑1 (MQ=255)
                                                                                                                                                    tGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGc                                           >  1:1805618/1‑108 (MQ=255)
                                                                                                                                                    tGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGc                                           <  2:1805618/108‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACTACCCCAACGCCCCCGCCGTCTTAGATACTGAGTGGTTTATTTCCCTGTGCTTCACGCTACTGGAGAACACGACATGCGCCACCTACTGCTCGCCCCTGCGCTGCTGCTCCTGATGCCTGCCGGTGCTGCCCTGGCCCGCGTCGATGCGGGCGACGTCGCCACCTCGGCCGGTATTTCCGCTTCGCTGTATTCGACCTTCAAGGACGACAAACGGATCATCCCTGCCCGTGACGAACTTTCAGCGTTCGTCGCCAGCGGCGGTACCATCCGCGGCGCCTATGTGGAGTCGGCGC  >  NC_002947/4272674‑4272970

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