Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,926,556 C→A 7.3% D118Y (GAC→TAC) ‡ PP_1726 ← ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,926,5560CA 7.3% 310.1 / 9.0 96D118Y (GAC→TAC) ‡PP_1726ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (53/36);  new base A (4/3);  total (57/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CTTGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGGTGT  >  NC_002947/1926410‑1926702
                                                                                                                                                  |                                                                                                                                                  
cTTGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCa                                                                                                                                                 >  1:2382477/1‑150 (MQ=255)
  tGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTCGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAgg                                                                                                                                               >  1:1465445/1‑150 (MQ=255)
      cTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgcc                                                                                                                                           <  2:1646674/150‑2 (MQ=255)
        tGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttg                                                                                                                                         <  2:1072383/150‑1 (MQ=255)
        tGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttg                                                                                                                                         <  2:2138777/150‑1 (MQ=255)
           cTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgcc                                                                                                                                           >  1:1091430/1‑144 (MQ=255)
           cTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgcc                                                                                                                                           <  2:1091430/145‑2 (MQ=255)
                tcCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTaa                                                                                                                                 >  2:2197024/1‑150 (MQ=255)
                tcCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgccggaat                                                                                                                                      <  1:824269/145‑7 (MQ=255)
                tcCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGGACCAGGTgccggaat                                                                                                                                      >  2:824269/1‑139 (MQ=255)
                 cctggTCATGGCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAg                                                                                                                                <  1:646815/147‑1 (MQ=255)
                  cAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggctt                                                                                                                                     >  1:730374/1‑144 (MQ=255)
                  cAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggctt                                                                                                                                     <  2:730374/144‑1 (MQ=255)
                   aGGTCATGCTAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACGTGGTCCCAGGTGCTTGGCTTGTAAggc                                                                                                                              <  1:1781054/150‑1 (MQ=255)
                    ggTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggct                                                                                                                             <  2:1701906/150‑1 (MQ=255)
                      tCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctgg                                                                                                                           <  1:921402/150‑1 (MQ=255)
                          gCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctggc                                                                                                                          <  2:369447/147‑1 (MQ=255)
                          gCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctggc                                                                                                                          >  1:369447/1‑147 (MQ=255)
                          gCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctggc                                                                                                                          <  1:2185888/147‑1 (MQ=255)
                          gCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctggc                                                                                                                          >  2:2185888/1‑147 (MQ=255)
                          gCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGa                                                                                                                       >  2:1273068/1‑150 (MQ=255)
                            cAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCATTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACg                                                                                                                     >  1:1496667/1‑150 (MQ=255)
                                  ttCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                               <  2:2382477/150‑1 (MQ=255)
                                    cGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggc                                                                                                                                       <  1:1824590/124‑1 (MQ=255)
                                    cGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggc                                                                                                                                       >  2:1824590/1‑124 (MQ=255)
                                        cGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCAcc                                                                                                         >  1:837460/1‑150 (MQ=255)
                                        cGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCAcc                                                                                                         >  2:968798/1‑150 (MQ=255)
                                          ttCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACcgc                                                                                                       <  1:1273068/150‑1 (MQ=255)
                                             ctctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgctt                                                                                                                                          >  1:1802300/1‑112 (MQ=255)
                                             ctctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgctt                                                                                                                                          <  2:1802300/112‑1 (MQ=255)
                                             ctctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                               >  1:1883769/1‑139 (MQ=255)
                                             ctctTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                               <  2:1883769/139‑1 (MQ=255)
                                                  ccTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGc                                                                                               >  2:201185/1‑150 (MQ=255)
                                                  ccTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGc                                                                                               >  1:1568596/1‑150 (MQ=255)
                                                  ccTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGc                                                                                               <  2:1246093/150‑1 (MQ=255)
                                                   cTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggctt                                                                                                                                     >  2:1912359/1‑111 (MQ=255)
                                                   cTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTgcttggctt                                                                                                                                     <  1:1912359/111‑1 (MQ=255)
                                                    ttGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggct                                                                                                                             <  2:1408758/118‑1 (MQ=255)
                                                    ttGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggct                                                                                                                             >  1:1408758/1‑118 (MQ=255)
                                                    ttGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGAc                                                                                                                      <  1:766705/125‑1 (MQ=255)
                                                    ttGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGAc                                                                                                                      >  2:766705/1‑125 (MQ=255)
                                                     tGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGa                                                                                            >  2:130883/1‑150 (MQ=255)
                                                       aCCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTccc                                                                                                                                                  <  1:2000947/94‑1 (MQ=255)
                                                       aCCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTccc                                                                                                                                                  >  2:2000947/1‑94 (MQ=255)
                                                         cAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAggc                                                                                        >  2:1753529/1‑150 (MQ=255)
                                                         cAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAggc                                                                                        <  2:1568596/150‑1 (MQ=255)
                                                          aGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGCAggcg                                                                                       <  2:1140394/150‑1 (MQ=255)
                                                          aGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGCATTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAggcg                                                                                       <  1:1757782/150‑1 (MQ=255)
                                                            gTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggct                                                                                                                             <  2:1941536/110‑1 (MQ=255)
                                                            gTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggct                                                                                                                             >  1:1941536/1‑110 (MQ=255)
                                                                ttgttgATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCAccc                                                                                 >  1:303509/1‑150 (MQ=255)
                                                                      atgatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAg                                                                                                                                                >  2:993788/1‑81 (MQ=255)
                                                                      atgatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAg                                                                                                                                                <  1:993788/81‑1 (MQ=255)
                                                                        gatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctggc                                                                                                                          >  1:1098344/1‑101 (MQ=255)
                                                                        gatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctggc                                                                                                                          <  2:1098344/101‑1 (MQ=255)
                                                                        gatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                               >  2:1248920/1‑112 (MQ=255)
                                                                        gatgaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTg                                                                                                               <  1:1248920/112‑1 (MQ=255)
                                                                           gaAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGCCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGAt                                                                      >  1:2289757/1‑150 (MQ=255)
                                                                            aaaGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCGCGGATg                                                                     >  2:1111642/1‑150 (MQ=255)
                                                                            aaaGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTACTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTACGGCTTGCTGACGCCCGTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATg                                                                     >  1:91214/1‑150 (MQ=255)
                                                                              aGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTCTCACTGGCCGATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTTGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTc                                                                   <  2:319206/150‑1 (MQ=255)
                                                                               gCGATGGGGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCg                                                                  >  2:705822/1‑150 (MQ=255)
                                                                                   tGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGgcgc                                                              >  2:1637914/1‑150 (MQ=255)
                                                                                          ggTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctgg                                                                                                                           >  1:1780597/1‑82 (MQ=255)
                                                                                          ggTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAggctgg                                                                                                                           <  2:1780597/82‑1 (MQ=255)
                                                                                                  ccAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTtct                                               >  2:1744458/1‑150 (MQ=255)
                                                                                                   cAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCa                                                                                                           >  1:1382300/1‑89 (MQ=255)
                                                                                                   cAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCa                                                                                                           <  2:1382300/89‑1 (MQ=255)
                                                                                                   cAGCGCCCAGTGGCCGTCCTTGGCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTtctc                                              <  1:1012792/150‑1 (MQ=255)
                                                                                                           aGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGa                                      <  1:201185/150‑1 (MQ=255)
                                                                                                            gTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGaa                                     <  1:968798/150‑1 (MQ=255)
                                                                                                                  gtccttgtccttGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGc                               <  1:705822/150‑1 (MQ=255)
                                                                                                                  gtccttgtccttGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGc                               >  2:125908/1‑150 (MQ=255)
                                                                                                                  gtccttgtccttGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGc                               >  1:1584371/1‑150 (MQ=255)
                                                                                                                     cttgtccttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAg                            >  1:415315/1‑150 (MQ=255)
                                                                                                                      ttgtccttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGt                           >  1:466493/1‑150 (MQ=255)
                                                                                                                           cttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTg                      >  1:869222/1‑150 (MQ=255)
                                                                                                                           cttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTg                      >  2:204167/1‑150 (MQ=255)
                                                                                                                           cttgGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCAACCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTg                      <  1:1111642/150‑1 (MQ=255)
                                                                                                                               gCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGcc                  >  1:1820804/1‑150 (MQ=255)
                                                                                                                                cccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAggcggc                                                                                     >  1:205145/1‑82 (MQ=255)
                                                                                                                                cccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAggcggc                                                                                     <  2:205145/82‑1 (MQ=255)
                                                                                                                                cccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCg                 >  1:1850805/1‑150 (MQ=255)
                                                                                                                                cccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCg                 >  1:1422202/1‑150 (MQ=255)
                                                                                                                                cccATTCCGGCACCTGGTCCCAGGTGCTTGGCGTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAAGCGGCACCCACGTCACCGATGTCTGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCg                 >  2:1076118/1‑150 (MQ=255)
                                                                                                                                 ccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGa                >  2:1939913/1‑150 (MQ=255)
                                                                                                                                 ccATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGa                >  2:2261001/1‑150 (MQ=255)
                                                                                                                                  cATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAg               >  2:4135/1‑150 (MQ=255)
                                                                                                                                   aTTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGc              >  2:1550903/1‑150 (MQ=255)
                                                                                                                                   aTTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCGGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGc              >  2:1392786/1‑150 (MQ=255)
                                                                                                                                      caagCACCTGGGACCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCa           >  2:42716/4‑150 (MQ=255)
                                                                                                                                       cGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTg                                 >  1:1263813/1‑127 (MQ=255)
                                                                                                                                       cGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTg                                 <  2:1263813/127‑1 (MQ=255)
                                                                                                                                       aagCACCTGGGACCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCAt          >  2:2350884/3‑150 (MQ=255)
                                                                                                                                             cTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGgt    >  2:1346935/1‑150 (MQ=255)
                                                                                                                                             cTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGgt    >  1:1941003/1‑150 (MQ=255)
                                                                                                                                              tGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGgtg   >  2:1412610/1‑150 (MQ=255)
                                                                                                                                               ggTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGgtgt  <  2:1496667/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTTGTACTTGCCTTTCTCCAGGTCATGCCAGGTTTTCGGCCGTTCCTCTTCCTTGACCAGGTCCTTGTTGATGATGAAAGCGATGGTGCCGGTATAGGCCAGCGCCCAGTGGCCGTCCTTGTCCTTGGCCCATTCCGGCACCTGGTCCCAGGTGCTTGGCTTGTAAGGCTGGCTGACGCCCTTGGTCACCGCGATCGGGCCGAAGGCGGCACCCACGTCACCGATGTCGGCGCTGGCGTTGTCCTTCTCGGCTTCGAACTTGGCCAGTTCCTGGGCCGAGCTCATGTCGGTGT  >  NC_002947/1926410‑1926702

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: