Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,585,827 A→G G220G (GGT→GGC)  ydeP ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,585,8270AG99.0% 382.6 / ‑4.4 100G220G (GGT→GGC) ydePputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (48/51);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 101 total aligned reads displayed.

AGTCTCAGAGTTCGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTACCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGAAAGTTATTGCTC  >  NC_000913/1585680‑1585974
                                                                                                                                                   |                                                                                                                                                   
aGTCTCAGAGTTCGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGcc                                                                                                                                                   >  1:228477/1‑150 (MQ=255)
      agagTTCGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGtt                                                                                                                                             >  3:314804/1‑150 (MQ=255)
      agagTTCGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTa                                                                                                                                            <  4:1177983/151‑1 (MQ=255)
         gTTCGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGc                                                                                                                                         <  4:1223813/151‑1 (MQ=255)
            cGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGccc                                                                                                                                       >  2:58815/1‑150 (MQ=255)
            cGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGccc                                                                                                                                       >  2:1574370/1‑150 (MQ=255)
            cGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGccc                                                                                                                                       >  4:12442/1‑150 (MQ=255)
            cGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCa                                                                                                                                      <  2:903179/151‑1 (MQ=255)
              gcaccATATCAAACGGGGTTTGGGGAGAGGTAAAGGGCTACAGGCAACGTTCCTGGAGAGGATGGATGGCTATCATAATCACTCCCCGATTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTACCAGGGTTATGCCCa                                                                                                                                      >  2:875840/5‑149 (MQ=255)
               cAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATg                                                                                                                                   <  2:776161/151‑1 (MQ=255)
                 gCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGc                                                                                                                                  >  4:1332685/1‑150 (MQ=255)
                 gCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCa                                                                                                                                 <  2:1327200/151‑1 (MQ=255)
                   aTTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCaaa                                                                                                                               <  2:1591795/151‑1 (MQ=255)
                   aTTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCaa                                                                                                                                >  2:1115367/1‑150 (MQ=255)
                    tttCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAAt                                                                                                                              >  2:412674/1‑151 (MQ=255)
                    tttCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAAt                                                                                                                              >  4:1294454/1‑151 (MQ=255)
                     ttCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTTATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGTTTATGCCCAATGCAAAt                                                                                                                              <  2:1307751/150‑1 (MQ=255)
                     ttCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAAt                                                                                                                              <  4:1460688/150‑1 (MQ=255)
                      tCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATg                                                                                                                             <  2:65491/150‑1 (MQ=255)
                       cAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATgc                                                                                                                            <  2:1334711/150‑2 (MQ=255)
                       cAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGa                                                                                                                            <  2:1663510/150‑1 (MQ=255)
                       cAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGAc                                                                                                                           >  2:172824/1‑151 (MQ=255)
                        aaaCGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACt                                                                                                                          <  3:82982/151‑1 (MQ=255)
                           cGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTaaa                                                                                                                       <  2:41627/151‑1 (MQ=255)
                            gggTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAAt                                                                                                                      >  2:518596/1‑151 (MQ=255)
                             ggTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATc                                                                                                                     <  4:628159/151‑1 (MQ=255)
                               ttttGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGc                                                                                                                   >  4:1086460/1‑151 (MQ=255)
                                 ttGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCa                                                                                                                  <  3:232257/150‑1 (MQ=255)
                                   gcggtgcggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACtt                                                                                                               <  4:77058/151‑1 (MQ=255)
                                    cggtgcggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTtc                                                                                                              <  4:748779/151‑1 (MQ=255)
                                     ggtgcggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTtct                                                                                                             <  2:1076215/151‑1 (MQ=255)
                                       tgcggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTtctc                                                                                                            <  2:1144199/150‑1 (MQ=255)
                                        gcggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCaa                                                                                                          <  4:1543115/151‑1 (MQ=255)
                                        gcggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCa                                                                                                           <  4:1116779/150‑1 (MQ=255)
                                        gcggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCa                                                                                                           <  4:768165/150‑1 (MQ=255)
                                         cggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCaaa                                                                                                         <  2:1475129/151‑1 (MQ=255)
                                         cggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCaaa                                                                                                         <  2:819564/151‑1 (MQ=255)
                                         cggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCaa                                                                                                          >  4:1533104/1‑150 (MQ=255)
                                          ggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCaa                                                                                                          >  4:520338/1‑149 (MQ=255)
                                          ggtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAg                                                                                                        <  1:441500/151‑1 (MQ=255)
                                           gtAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGt                                                                                                       >  1:465426/1‑151 (MQ=255)
                                              aaTCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCtt                                                                                                    <  4:1072184/151‑1 (MQ=255)
                                               aTCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTc                                                                                                   <  4:1058485/151‑1 (MQ=255)
                                               aTCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTc                                                                                                   <  2:829578/151‑1 (MQ=255)
                                                 cGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTcc                                                                                                  >  4:1119544/1‑150 (MQ=255)
                                                 cGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCa                                                                                                 <  2:61885/151‑1 (MQ=255)
                                                  gCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTcc                                                                                                  >  3:22149/1‑149 (MQ=255)
                                                  gCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTcc                                                                                                  >  2:1147121/1‑149 (MQ=255)
                                                      cAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAac                                                                                            <  4:1359229/151‑1 (MQ=255)
                                                       aGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAaca                                                                                           <  4:121817/151‑1 (MQ=255)
                                                        ggCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAacac                                                                                          <  2:895318/151‑1 (MQ=255)
                                                         gCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACg                                                                                         <  2:551336/151‑1 (MQ=255)
                                                         gCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACg                                                                                         >  2:186034/1‑151 (MQ=255)
                                                            aCGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTc                                                                                      >  4:1679967/1‑151 (MQ=255)
                                                                tCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCt                                                                                   <  2:1449831/150‑1 (MQ=255)
                                                                 ccTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCttt                                                                                 >  4:860798/1‑151 (MQ=255)
                                                                    gTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTAcc                                                                              <  4:195193/151‑1 (MQ=255)
                                                                        agGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTaca                                                                          <  2:506677/151‑1 (MQ=255)
                                                                             tGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGGGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGa                                                                      >  4:1352865/1‑150 (MQ=255)
                                                                              gATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGAt                                                                     >  2:632932/1‑150 (MQ=255)
                                                                               aTGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATAc                                                                   >  4:1252418/1‑151 (MQ=255)
                                                                                 ggCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACt                                                                  <  2:1523281/150‑1 (MQ=255)
                                                                                  gCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTc                                                                 >  2:1434561/1‑150 (MQ=255)
                                                                                    gATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCt                                                              <  4:785518/151‑1 (MQ=255)
                                                                                     aTCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTg                                                             >  4:199240/1‑151 (MQ=255)
                                                                                        aTTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGcc                                                           >  1:200138/1‑150 (MQ=255)
                                                                                            tCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAAcc                                                      >  4:246169/1‑151 (MQ=255)
                                                                                            tCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAAc                                                       <  4:1230707/150‑1 (MQ=255)
                                                                                             cGCTCCCCGTTTCACTAAAGCGCGCAACGGAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAAccc                                                     >  2:808682/1‑151 (MQ=255)
                                                                                              gCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAAccc                                                     >  2:598357/1‑150 (MQ=255)
                                                                                              gCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAAccc                                                     >  4:1618728/1‑150 (MQ=255)
                                                                                               cTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCAc                                                   >  2:1213074/1‑151 (MQ=255)
                                                                                                tCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACg                                                  >  2:795087/1‑151 (MQ=255)
                                                                                                       ttCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTc                                            <  3:381323/150‑1 (MQ=255)
                                                                                                        tCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCgg                                          >  4:506487/1‑151 (MQ=255)
                                                                                                          aCTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGtt                                        <  4:1463350/151‑1 (MQ=255)
                                                                                                              aaGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCAt                                     >  4:514177/1‑150 (MQ=255)
                                                                                                              aaGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATg                                    <  2:505785/151‑1 (MQ=255)
                                                                                                                 cgcgcAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGc                                  >  2:1371338/1‑150 (MQ=255)
                                                                                                                  gcgcAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGcac                                <  1:353998/151‑1 (MQ=255)
                                                                                                                   cgcAACGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGcac                                >  2:558239/1‑150 (MQ=255)
                                                                                                                      aaCGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGt                            >  1:396115/1‑151 (MQ=255)
                                                                                                                      aaCGAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGt                            >  1:408038/1‑151 (MQ=255)
                                                                                                                         gAAGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTgg                         <  4:156864/151‑1 (MQ=255)
                                                                                                                          aaGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGa                        <  4:983147/151‑1 (MQ=255)
                                                                                                                          aaGTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGa                        <  2:1103544/151‑1 (MQ=255)
                                                                                                                            gTCAGCATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGc                      >  3:462239/1‑151 (MQ=255)
                                                                                                                                 cATGCGAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCg                 >  1:240045/1‑151 (MQ=255)
                                                                                                                                 cATGCGAGGGTGGGTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCACTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTc                  >  2:29675/1‑150 (MQ=255)
                                                                                                                                      gAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGaaa            <  3:50588/151‑1 (MQ=255)
                                                                                                                                      gAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGaa             >  2:1050792/1‑150 (MQ=255)
                                                                                                                                      gAGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGaa             >  2:1630877/1‑150 (MQ=255)
                                                                                                                                       aGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGaa             >  2:415104/1‑149 (MQ=255)
                                                                                                                                       aGGGTGGTTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGAAAg           <  4:1250301/151‑1 (MQ=255)
                                                                                                                                         ggtggtTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGAAAGtt         >  4:463453/1‑151 (MQ=255)
                                                                                                                                          gtggtTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGTTTCATGGCGCATGTTGGAGCAGTCGGGAAAGTTa        <  4:1148579/151‑1 (MQ=255)
                                                                                                                                            ggtTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGAAAGTTAtt      <  4:1662104/151‑1 (MQ=255)
                                                                                                                                             gtTTGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGAAAGTTATTg     >  1:69407/1‑151 (MQ=255)
                                                                                                                                               ttGTGCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGAAAGTTATTGc    <  4:413266/150‑1 (MQ=255)
                                                                                                                                                tgtgCCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGAAAGTTATTGCTc  >  3:178671/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTCTCAGAGTTCGTCAGCATTTCAAACGGGTTTTGCGGTGCGGTAAATCGCTCCAGGCCACGTTCCTGTAGAGGATTGATGGCGATCATTTTCGCTCCCCGTTTCACTAAAGCGCGCAACGAAGTCAGCATGCGAGGGTGGTTTGTACCAGGGTTATGCCCAATGCAAATGACTAAATCGCACTTCTCAAAGTCTTCCAGCAACACGGTCCCTTTACCTACACCGATACTCGCTGCCAAACCCACGCTTGTCGGTTCATGGCACATGTTGGAGCAGTCGGGAAAGTTATTGCTC  >  NC_000913/1585680‑1585974

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑