Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 636,069 A→G 5.2% R167R (CGT→CGC ybdN ← putative PAPS reductase/DUF3440 domain‑containing protein YbdN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913636,0690AG 5.2% 239.9 / 2.6 77R167R (CGT→CGCybdNputative PAPS reductase/DUF3440 domain‑containing protein YbdN
Reads supporting (aligned to +/- strand):  ref base A (32/41);  new base G (2/2);  total (34/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 = 6.95e-01

TCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTC  >  NC_000913/635929‑636218
                                                                                                                                             |                                                                                                                                                     
acAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCg                                                                                                                                              <  2:652481/149‑1 (MQ=255)
 cAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGTGCGGATgagtc                                                                                                                                             <  2:963380/150‑5 (MQ=255)
 cAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGTGCGGATgagtc                                                                                                                                             <  2:910618/150‑5 (MQ=255)
 cAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACg                                                                                                                                                     >  2:952801/1‑142 (MQ=255)
 cAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACg                                                                                                                                                     <  1:952801/142‑1 (MQ=255)
   gCCACTTTCCAGTCGTAAATTTTGTAAATTTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCATATCGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGatc                                                                                                                                           <  2:450006/150‑1 (MQ=255)
   gCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGatc                                                                                                                                           >  1:188936/1‑150 (MQ=255)
      aCTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACAGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGatcatc                                                                                                                                        >  2:64731/1‑150 (MQ=255)
       cTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGatcatca                                                                                                                                       <  2:1375144/150‑1 (MQ=255)
              gTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATcgccgc                                                                                                                                >  1:1433570/1‑150 (MQ=255)
                          gTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTtg                                                                                                                    <  2:128142/150‑1 (MQ=255)
                            aaaTGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTtgtg                                                                                                                  <  2:187249/150‑1 (MQ=255)
                             aaTGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGa                                                                                                                 >  2:809115/1‑150 (MQ=255)
                                gTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGaaaa                                                                                                              >  2:697127/1‑150 (MQ=255)
                                gTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGaaaa                                                                                                              <  2:1424490/150‑1 (MQ=255)
                                    cAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCATGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCAt                                                                                                          >  1:127536/1‑150 (MQ=255)
                                      gCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTc                                                                                                        <  2:292950/150‑1 (MQ=255)
                                      gCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAACCATTCg                                                                                                       <  2:369501/150‑1 (MQ=255)
                                        tATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAAATGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCgg                                                                                                      <  2:1223049/150‑1 (MQ=255)
                                         aTGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGc                                                                                                     >  2:1108391/1‑150 (MQ=255)
                                         aTGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGc                                                                                                     <  2:381003/150‑1 (MQ=255)
                                            aCCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGATGTTTATTTAAACTGGCGATTGCGACAAAACGGTTGTAGGAATCATCCGCACGGATGCCGATCATCATCCCCCCCGGACGTTTTTGTGAAAACCAATCGGCaaa                                                                                                  >  2:1191840/1‑150 (MQ=255)
                                             ccGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAAc                                                                                                 >  2:545166/1‑150 (MQ=255)
                                             ccGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAAc                                                                                                 <  2:1387996/150‑1 (MQ=255)
                                              cGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTGGGCAAACt                                                                                                <  2:411446/150‑1 (MQ=255)
                                                  tGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACg                                                                                            <  1:135198/150‑1 (MQ=255)
                                                   ggTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGt                                                                                           >  2:17069/1‑150 (MQ=255)
                                                      gCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACa                                                                                        <  1:529524/150‑1 (MQ=255)
                                                       cGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACaa                                                                                       <  2:1433570/150‑1 (MQ=255)
                                                         gCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATAGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGAACAAAt                                                                                     >  2:323773/1‑150 (MQ=255)
                                                          ccGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAAtt                                                                                    >  2:1418404/1‑150 (MQ=255)
                                                             tAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGtt                                                                                 <  1:1563494/150‑1 (MQ=255)
                                                                tCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGa                                                                              >  2:729938/1‑150 (MQ=255)
                                                                  cAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAg                                                                            <  2:251264/150‑1 (MQ=255)
                                                                      gTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCa                                                                        <  1:806782/150‑1 (MQ=255)
                                                                      gTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCa                                                                        >  2:1263408/1‑150 (MQ=255)
                                                                       tttATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGAAGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCAt                                                                       >  1:563985/1‑150 (MQ=255)
                                                                        ttATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATg                                                                      <  1:697127/150‑1 (MQ=255)
                                                                          aTCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGTATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTGACGTACAAATTGTTCGAAGGTCATGcc                                                                    <  1:739352/150‑1 (MQ=255)
                                                                          aTCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGcc                                                                    <  2:1065338/150‑1 (MQ=255)
                                                                          aTCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGcc                                                                    <  1:1071808/150‑1 (MQ=255)
                                                                            cgtcgGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGTTTGTAGGACTCATCCGCAGGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAAGCATTCGGCAAACTCACGTAGAAATTGTTCGAAGGTCATGCctg                                                                  <  2:585508/150‑1 (MQ=255)
                                                                                ggCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCctggctg                                                              <  2:992832/150‑1 (MQ=255)
                                                                                 gtaaatagtttttTTTTTAAACGGGCGCTGGCCCACAAAACGGTTTTTGTACTAATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCctggctg                                                              <  1:1554192/143‑1 (MQ=255)
                                                                                 gCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGcc                                                                    <  2:591422/143‑1 (MQ=255)
                                                                                 gCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGcc                                                                    >  1:591422/1‑143 (MQ=255)
                                                                                   aaaaCGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTa                                                           <  1:321359/150‑1 (MQ=255)
                                                                                      aCGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGGTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTaaaa                                                        <  2:932939/150‑1 (MQ=255)
                                                                                       cGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAg                                                       <  1:266598/150‑1 (MQ=255)
                                                                                        gttgttTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGcc                                                                    >  1:1627519/1‑136 (MQ=255)
                                                                                        gttgttTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGcc                                                                    <  2:1627519/136‑1 (MQ=255)
                                                                                           gttTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGa                                                   >  1:892360/1‑150 (MQ=255)
                                                                                            tttatttaAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGaa                                                  <  1:1263408/150‑1 (MQ=255)
                                                                                                tttaAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAg                                              >  1:1179586/1‑150 (MQ=255)
                                                                                                 ttaAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGt                                             >  1:1114869/1‑150 (MQ=255)
                                                                                                  taAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTc                                            >  2:1340463/1‑150 (MQ=255)
                                                                                                    aaCTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGTCGATCATCATCGCCGCCGGACGTTTTTCTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAg                                          <  2:683251/150‑1 (MQ=255)
                                                                                                    aaCTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAg                                          <  1:545166/150‑1 (MQ=255)
                                                                                                    aaCTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAg                                          <  1:261198/150‑1 (MQ=255)
                                                                                                      cTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGANAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAggg                                        >  2:279306/1‑150 (MQ=255)
                                                                                                      cTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAggg                                        >  2:273682/1‑150 (MQ=255)
                                                                                                      cTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAggg                                        <  2:188936/150‑1 (MQ=255)
                                                                                                       tggcgatggcgaCAAAACTGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGGTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGt                                       <  2:1020672/150‑1 (MQ=255)
                                                                                                               gcgaCAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATc                               <  1:1190222/150‑1 (MQ=255)
                                                                                                                 cgaCAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCg                              >  1:1372852/1‑150 (MQ=255)
                                                                                                                  gaCAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGc                             >  1:1357575/1‑150 (MQ=255)
                                                                                                                     aaaaCGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATc                          <  1:131831/150‑1 (MQ=255)
                                                                                                                         gggTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTgg                      >  1:1480139/2‑150 (MQ=255)
                                                                                                                         cGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTgg                      >  2:1197851/1‑150 (MQ=255)
                                                                                                                               tgatcggcATCCGTGCGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGt                >  1:1366243/8‑150 (MQ=255)
                                                                                                                                  gACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGa             <  1:64731/150‑1 (MQ=255)
                                                                                                                                   aCTCATCCGCACGGATGCCGATCTTCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGAc            >  2:254673/1‑150 (MQ=255)
                                                                                                                                     gcATCCGTGCGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGt          >  2:1460052/2‑150 (MQ=255)
                                                                                                                                        tCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTAcc       <  2:1169104/150‑1 (MQ=255)
                                                                                                                                         ccGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTAccc      >  2:1597161/1‑150 (MQ=255)
                                                                                                                                         ccGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTCGACGTAccc      >  1:141565/1‑150 (MQ=255)
                                                                                                                                            cACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCAtt   >  1:660761/1‑150 (MQ=255)
                                                                                                                                             aCGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTc  >  2:1151121/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
TCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGGTTGTAGGACTCATCCGCACGGATGCCGATCATCATCGCCGCCGGACGTTTTTGTGAAAACCATTCGGCAAACTCACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTC  >  NC_000913/635929‑636218

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