Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 683,664 C→T 5.5% G64S (GGC→AGC)  PP_0585 ← MerR family metal‑sensitive transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947683,6640CT 5.5% 262.8 / 3.7 90G64S (GGC→AGC) PP_0585MerR family metal‑sensitive transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base T (3/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CTGACAGTGCGAAACCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGG  >  NC_002947/683515‑683813
                                                                                                                                                     |                                                                                                                                                     
cTGACAGTGCGAAACCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGTTCCTCGTTGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGcc                                                                                                                                                       <  2:233433/150‑1 (MQ=255)
   aCAGTGCGAAACCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGAc                                                                                                                                                    <  2:1884494/150‑1 (MQ=255)
     aGTGCGAAACCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGAcct                                                                                                                                                  >  2:1062883/1‑150 (MQ=255)
         cGAAACCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACctcctc                                                                                                                                              <  2:2329509/150‑1 (MQ=255)
           aaaCCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGa                                                                                                                                            <  1:188623/150‑1 (MQ=255)
             aCCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACTTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAgg                                                                                                                                          <  2:481937/150‑1 (MQ=255)
             aCCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAgg                                                                                                                                          <  1:1060451/150‑1 (MQ=255)
                   tCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAg                                                                                                                                    >  2:1967577/1‑150 (MQ=255)
                       cGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCa                                                                                                                                <  1:1062883/150‑1 (MQ=255)
                               tCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCgg                                                                                                                        <  2:1474861/150‑1 (MQ=255)
                                 gcgcAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGa                                                                                                                      >  1:772355/1‑150 (MQ=255)
                                     aGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAgcgc                                                                                                                  <  1:1025645/150‑1 (MQ=255)
                                         tCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGaaaa                                                                                                                                     <  1:880974/127‑1 (MQ=255)
                                         tCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTtga                                                                                                              <  2:772355/150‑1 (MQ=255)
                                         tCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTtga                                                                                                              <  2:1764761/150‑1 (MQ=255)
                                         tCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTtga                                                                                                              >  2:1093566/1‑150 (MQ=255)
                                         tCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCAGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGAAAACGCCCAGGTCGCGCGAGCGGTtga                                                                                                              >  1:506050/1‑150 (MQ=255)
                                         tCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTAACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGaaaa                                                                                                                                     >  2:880974/1‑127 (MQ=255)
                                          cACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTtgat                                                                                                             <  2:141988/150‑1 (MQ=255)
                                           aCCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTtgctg                                                                                                            >  1:1943782/1‑150 (MQ=255)
                                              aGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCTACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAg                                                                                                         <  1:2187401/150‑1 (MQ=255)
                                                  cctcGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCa                                                                                                     >  2:640901/1‑150 (MQ=255)
                                                   ctcGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTTCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAg                                                                                                         <  2:1438100/145‑1 (MQ=255)
                                                   ctcGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTTCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAg                                                                                                         >  1:1438100/1‑145 (MQ=255)
                                                   ctcGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAg                                                                                                    >  1:605391/1‑150 (MQ=255)
                                                           gcCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAgg                                                                                                                                          >  2:1487532/1‑104 (MQ=255)
                                                           gcCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAgg                                                                                                                                          <  1:1487532/104‑1 (MQ=255)
                                                                   agttcaTCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTgc                                                                                    <  1:1967577/150‑1 (MQ=255)
                                                                    gttcaTCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTGGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTgcg                                                                                   >  1:1323943/1‑150 (MQ=255)
                                                                     ttcaTCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTgcgc                                                                                  >  1:1689623/1‑150 (MQ=255)
                                                                        aTCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTg                                                                               >  2:611688/1‑150 (MQ=255)
                                                                        aTCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTg                                                                               >  1:117402/1‑150 (MQ=255)
                                                                        aTCAATGTGCTGCATGGCCAGCACCTTCACATCGGCGCTGGCACGTTGTCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGATGGCAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTg                                                                               <  1:2080414/150‑1 (MQ=255)
                                                                           gATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTa                                                                            >  2:2249328/1‑150 (MQ=255)
                                                                             tgtgGTTCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAg                                                                                                                                           <  1:2090050/85‑1 (MQ=255)
                                                                             tgtgGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAg                                                                                                                                           >  2:2090050/1‑85 (MQ=255)
                                                                               tggtgCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAgg                                                                        >  2:626208/1‑150 (MQ=255)
                                                                                ggtgCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGc                                                                       <  1:640901/150‑1 (MQ=255)
                                                                                 gtgCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCg                                                                      >  1:1385868/1‑150 (MQ=255)
                                                                                   gCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGt                                                                    >  1:51595/1‑150 (MQ=255)
                                                                                    cATGGCCAGCGCCTTCACATCGGCGCTTGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTa                                                                   <  1:611688/150‑1 (MQ=255)
                                                                                      tGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAACGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCa                                                                                                     >  1:411846/1‑114 (MQ=255)
                                                                                      tGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCa                                                                                                     <  2:411846/114‑1 (MQ=255)
                                                                                             cgcCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTc                                                          <  2:51595/150‑1 (MQ=255)
                                                                                                 ttCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTg                                                      <  2:605391/150‑1 (MQ=255)
                                                                                                  tCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGGCCTCCTCGAGGGAAAANCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTg                                                                                     <  1:1178723/118‑1 (MQ=255)
                                                                                                  tCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGGCCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTg                                                                                     >  2:1178723/1‑118 (MQ=255)
                                                                                                    acaTCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCgg                                                   >  2:1749875/1‑149 (MQ=255)
                                                                                                    acaTCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGa                                                   <  2:1385868/150‑1 (MQ=255)
                                                                                                      aTCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTg                                                                               <  1:678367/120‑1 (MQ=255)
                                                                                                      aTCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGt                                                 >  1:1898081/1‑150 (MQ=255)
                                                                                                      aTCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATAAAGGCCAGGTTGTGCAAGTCTTGCGCCTg                                                                               >  2:678367/1‑120 (MQ=255)
                                                                                                      aTCGGCGCTCGCACGTTGGCGGTCCTGCCACAAGGTCAGTAACTTGCCGACCTCCTCAAGGGAAAAGCCCAGATCGCGGGAGCGTTTGATGAAGGCGAGGTTGTGCAAGTCCTCAGCCTGGCACCGCCGGTAGCAGCTGTAGCTGCgtgt                                                 >  1:2301804/1‑150 (MQ=255)
                                                                                                         ggagcTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGc                                              >  1:2400816/4‑150 (MQ=255)
                                                                                                         ggCGCTGGCAGGTTGGCGCTCCTGCCACAAGGTCAGCAACTTGCCGACCACCTCGTGGGCCAGCACTCGGTCCCCTTAGCGCTTGATGAACTCCAAGTTGTGATCGTCTTTTTCCGCTTACAGTCGGTAGCCGCTGTCGCTGCGAGTGGc                                              <  1:1749875/150‑1 (MQ=255)
                                                                                                          gcgcTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAg                                                                                                         >  2:666649/1‑90 (MQ=255)
                                                                                                          gcgcTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAg                                                                                                         <  1:666649/90‑1 (MQ=255)
                                                                                                               ggCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGttt                                        <  2:2400816/150‑1 (MQ=255)
                                                                                                                 gtcGTTGGCGGTCCTGCCACAAGGTCAGCGACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGa                                      <  2:1214691/148‑1 (MQ=255)
                                                                                                                 cACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGCGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGa                                      <  1:376566/150‑1 (MQ=255)
                                                                                                                    gTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCgg                                            <  2:1616397/141‑1 (MQ=255)
                                                                                                                    gTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCgg                                            >  1:1616397/1‑141 (MQ=255)
                                                                                                                     ttGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGGGGGACAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAg                                  >  2:2011329/1‑150 (MQ=255)
                                                                                                                     ttGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAg                                  >  1:1291784/1‑150 (MQ=255)
                                                                                                                         cGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCg                              <  1:1359046/150‑1 (MQ=255)
                                                                                                                            tCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATg                           >  1:625615/1‑150 (MQ=255)
                                                                                                                             ccTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATgg                          <  1:2249328/150‑1 (MQ=255)
                                                                                                                              cTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCt                                                             >  1:1491509/1‑114 (MQ=255)
                                                                                                                              cTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCt                                                             <  2:1491509/114‑1 (MQ=255)
                                                                                                                               ctcgaggaGGTCGGCAAGTTGCTGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGTGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGAt                        >  2:315116/8‑150 (MQ=255)
                                                                                                                                  cacaAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGcc                                                                                                      <  2:2021603/69‑1 (MQ=255)
                                                                                                                                  cacaAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGcc                                                                                                      >  1:2021603/1‑69 (MQ=255)
                                                                                                                                   acaAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGAt                            >  1:2345274/1‑142 (MQ=255)
                                                                                                                                   acaAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGAt                            <  2:2345274/142‑1 (MQ=255)
                                                                                                                                    caAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCgta                   >  2:834458/1‑150 (MQ=255)
                                                                                                                                      aGNTCAGCAACTTGCCGACCTCCTCGANGGAAAAGCCCAGNTCGCGGNANCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGCTTCgtagt                 >  1:1604032/1‑150 (MQ=255)
                                                                                                                                        gtcggCAAGTTGCTGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCgg                                            <  1:1349715/117‑1 (MQ=255)
                                                                                                                                        gTCGGCAAGTTGCTGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCgg                                            >  2:1349715/1‑121 (MQ=255)
                                                                                                                                        gTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCg                                                                      <  2:445094/95‑1 (MQ=255)
                                                                                                                                        gTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCg                                                                      >  1:445094/1‑95 (MQ=255)
                                                                                                                                           ggCAAGTTGCTGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCgta                   >  1:1117546/2‑143 (MQ=255)
                                                                                                                                           ggCAAGTTGCTGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCgta                   <  2:1117546/142‑1 (MQ=255)
                                                                                                                                            gCAACTTGCCGACCTCCTCGAGGGAAAAGCCCATGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGa           >  2:2250264/1‑150 (MQ=255)
                                                                                                                                            gCAACTTGCCGACCTCCTCGAGGGAAAAGCCCATGTCGCGGGAGCGCTCGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGa           >  2:2250484/1‑150 (MQ=255)
                                                                                                                                            gCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTgc                                                                                    >  1:22840/1‑77 (MQ=255)
                                                                                                                                            gCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTgc                                                                                    <  2:22840/77‑1 (MQ=255)
                                                                                                                                             cAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGat          <  2:1291784/150‑1 (MQ=255)
                                                                                                                                               aCTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGatca        >  2:920361/1‑150 (MQ=255)
                                                                                                                                                 ttGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTc                      >  1:917777/1‑134 (MQ=255)
                                                                                                                                                 ttGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTc                      <  2:917777/134‑1 (MQ=255)
                                                                                                                                                 ttGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGGTGGAGCAGTCCGATTGATTCGTAGTAACGGataatc      >  2:2310273/1‑150 (MQ=255)
                                                                                                                                                    ccGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTg   <  1:1093566/150‑1 (MQ=255)
                                                                                                                                                     cGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTgg  >  1:922414/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGACAGTGCGAAACCAGCTCGCCGAGGGTGTCGCGCAGGCTCACCAGCTCCTCGATGCGCCGGTTCAGTTCATCGATGTGGTGCATGGCCAGCGCCTTCACATCGGCGCTGGCACGTTGGCGGTCCTGCCACAAGGTCAGCAACTTGCCGACCTCCTCGAGGGAAAAGCCCAGGTCGCGGGAGCGCTTGATGAAGGCCAGGTTGTGCAAGTCTTGCGCCTGGTACAGGCGGTAGCCGCTGTCGCTGCGAGTGGCGGGTTTGAGCAGGCCGATGGATTCGTAGTAACGGATCATCTTGG  >  NC_002947/683515‑683813

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