Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,232,804 A→T 6.1% intergenic (+76/‑74) cdd → / → sanA cytidine/deoxycytidine deaminase/DUF218 domain‑containing protein SanA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,232,8040AT 6.1% 238.9 / 3.1 75intergenic (+76/‑74)cdd/sanAcytidine/deoxycytidine deaminase/DUF218 domain‑containing protein SanA
Reads supporting (aligned to +/- strand):  ref base A (31/39);  new base T (2/3);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGTGGGATGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTCGGCCTCGAT  >  NC_000913/2232656‑2232949
                                                                                                                                                    |                                                                                                                                                 
tCAGTGGGATGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAg                                                                                                                                                  >  2:1398268/1‑150 (MQ=255)
tCAGTGGGATGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTACGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAg                                                                                                                                                  <  1:1313083/150‑1 (MQ=255)
  aGTGGGATGGCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGtt                                                                                                                                                <  1:368880/150‑1 (MQ=255)
  aGTGGGATGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGtt                                                                                                                                                >  1:1406458/1‑150 (MQ=255)
   gTGGGATGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGttt                                                                                                                                               >  2:212625/1‑150 (MQ=255)
         tGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGtt                                                                                                                                         >  2:46831/1‑150 (MQ=255)
         tGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGtt                                                                                                                                         >  1:98317/1‑150 (MQ=255)
         tGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGCATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGGCGGATGCGGCGTGAACGCCTTGTCCCGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACTGTTTGCCggg                                                                                                                                         <  2:493921/150‑3 (MQ=255)
          gCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACTGTTTGCCggga                                                                                                                                        <  2:368308/150‑4 (MQ=255)
              ccTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTg                                                                                                                                    <  1:151720/150‑1 (MQ=255)
               cTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTgc                                                                                                                                   >  2:1772544/1‑150 (MQ=255)
                 ccGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTgcgc                                                                                                                                 >  1:1306042/1‑150 (MQ=255)
                 ccGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTgcgc                                                                                                                                 <  1:1148354/150‑1 (MQ=255)
                 ccGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTgcgc                                                                                                                                 <  2:1133271/150‑1 (MQ=255)
                  cGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTTGTGTCTAAATTACCGGCAAACAGTTTGTCGTTTCTTgatca                                                                                                                                <  2:1136973/150‑5 (MQ=255)
                  cGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTc                                                                                                                                       <  2:833730/143‑1 (MQ=255)
                  cGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTc                                                                                                                                       >  1:833730/1‑143 (MQ=255)
                  cGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCa                                                                                                                                >  2:273436/1‑150 (MQ=255)
                  cGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCa                                                                                                                                >  2:1605390/1‑150 (MQ=255)
                       cGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAAcc                                                                                                                           >  1:1006213/1‑150 (MQ=255)
                           gAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCg                                                                                                                       >  2:1148735/1‑150 (MQ=255)
                             aaGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCggg                                                                                                                     >  2:656365/1‑150 (MQ=255)
                              aGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACTGTTTGCCggg                                                                                                                                         >  2:1410612/1‑127 (MQ=255)
                              aGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACTGTTTGCCggg                                                                                                                                         <  1:1410612/129‑3 (MQ=255)
                              aGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGt                                                                                                                    >  1:1604705/1‑150 (MQ=255)
                                ctctCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTaa                                                                                                                  <  2:1876544/150‑1 (MQ=255)
                                 tctcGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTaaa                                                                                                                 <  1:1474454/150‑1 (MQ=255)
                                 tctcGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCATCTCCTGGTGTCGAAATTCACGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTaaa                                                                                                                 <  1:1956244/150‑1 (MQ=255)
                                    cGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTa                                                                                                              <  1:1640867/150‑1 (MQ=255)
                                     ggCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAg                                                                                                             <  2:1406458/150‑1 (MQ=255)
                                       cTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGcc                                                                                                           <  1:1604915/150‑1 (MQ=255)
                                          tCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGa                                                                                                        <  2:1315766/150‑1 (MQ=255)
                                          tCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGa                                                                                                        <  2:305430/150‑1 (MQ=255)
                                            acaGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATg                                                                                                      <  1:46831/150‑1 (MQ=255)
                                                   cGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAAttt                                                                                               >  2:168110/1‑150 (MQ=255)
                                                        gAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCtt                                                                                          <  2:1306042/150‑1 (MQ=255)
                                                            ccTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGAt                                                                                      >  1:442145/2‑150 (MQ=255)
                                                             cTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATc                                                                                     <  2:98317/150‑1 (MQ=255)
                                                              ttCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCg                                                                                    >  1:1140861/1‑150 (MQ=255)
                                                              ttCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCg                                                                                    <  1:1252671/150‑1 (MQ=255)
                                                              ttCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCg                                                                                    >  1:1898616/1‑150 (MQ=255)
                                                              ttCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCg                                                                                    <  2:1483718/150‑1 (MQ=255)
                                                                ctcGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAg                                                                                  <  2:416975/150‑1 (MQ=255)
                                                                ctcGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAg                                                                                  >  2:1163310/1‑150 (MQ=255)
                                                                         gCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCa                                                                         <  1:1986506/150‑1 (MQ=255)
                                                                         gCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCa                                                                         <  1:1518356/150‑1 (MQ=255)
                                                                           cTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATg                                                                       <  1:1148735/150‑1 (MQ=255)
                                                                              gTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTa                                                                    >  2:99167/1‑150 (MQ=255)
                                                                                gCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTaaa                                                                  >  1:3105/1‑150 (MQ=255)
                                                                                 ccGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAg                                                                 <  1:1605390/150‑1 (MQ=255)
                                                                                     aTGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGcgcg                                                             <  2:1604705/150‑1 (MQ=255)
                                                                                        cggcgTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATg                                                                       >  2:328258/1‑137 (MQ=255)
                                                                                        cggcgTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATg                                                                       <  1:328258/137‑1 (MQ=255)
                                                                                        cggcgTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCgtgt                                                          >  1:483189/1‑150 (MQ=255)
                                                                                           cgTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTcc                                                       <  1:273436/150‑1 (MQ=255)
                                                                                              gAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCa                                                    <  1:168110/150‑1 (MQ=255)
                                                                                               aaCGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAGATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAg                                                   >  2:774397/1‑150 (MQ=255)
                                                                                               aaCGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAg                                                   <  2:1898616/150‑1 (MQ=255)
                                                                                                   ccTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTg                                               <  1:668871/150‑1 (MQ=255)
                                                                                                         ccGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGcc                                                 <  1:1966444/142‑1 (MQ=255)
                                                                                                         ccGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGcc                                                 >  2:1966444/1‑142 (MQ=255)
                                                                                                           ggCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCt                                       <  1:1007630/150‑1 (MQ=255)
                                                                                                                gCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCg                                  >  1:768490/1‑150 (MQ=255)
                                                                                                                      ccTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTtgc                            <  1:1469310/150‑1 (MQ=255)
                                                                                                                              ggTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTg                    <  1:656365/150‑1 (MQ=255)
                                                                                                                                  tCGAGATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACtg                <  1:774397/150‑1 (MQ=255)
                                                                                                                                   cGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACtgt               >  1:1516438/1‑150 (MQ=255)
                                                                                                                                    gAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACAAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGAGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGATGCTGTTGACtttg              >  2:665112/1‑150 (MQ=255)
                                                                                                                                       aTTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTc           >  2:1180557/1‑150 (MQ=255)
                                                                                                                                       aTTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTc           <  1:212625/150‑1 (MQ=255)
                                                                                                                                       aTTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTc           >  1:1407030/1‑150 (MQ=255)
                                                                                                                                          aacGGCAAACTGTTTGCCGTTTCTTCCGCAAAACCAGCTGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATCTTAAAAAGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTCGGc        >  2:706788/3‑150 (MQ=255)
                                                                                                                                            cGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTCGGCCt      <  2:733608/150‑1 (MQ=255)
                                                                                                                                              gCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTCGGCCTCg    <  2:406296/150‑1 (MQ=255)
                                                                                                                                                aaaCAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTAATCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTCGGCCTCGAt  >  1:1135056/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCAGTGGGATGCCACCTCCGCAACGCTGAAAGCTCTCGGCTGTCACAGTATCGACCGAGTGCTTCTCGCTTAAGCCTGGTGCCGGATGCGGCGTGAACGCCTTGTCCGGCTTGCCAGCCCTCTCCTGGTGTCGAAATTCCCGGCAAACAGTTTGCCGTTTCTTGCGCAAAACCAGCGGGTAAAGTAGCCTGATGGAAATTTTCCTTAGATCGAGTCTCCTGCATGTTAAAGCGCGTGTTCCTCAGCCTGTTAGTCCTGATCGGCTTGCTGCTGTTGACTGTGCTCGGCCTCGAT  >  NC_000913/2232656‑2232949

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: