Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,592,990 A→C 8.6% intergenic (+28/+82) ypfN → / ← ypfH putative membrane protein, UPF0370 family/palmitoyl‑CoA esterase activity, uncertain physiological substrate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,592,9900AC 8.6% 280.6 / 10.7 105intergenic (+28/+82)ypfN/ypfHputative membrane protein, UPF0370 family/palmitoyl‑CoA esterase activity, uncertain physiological substrate
Reads supporting (aligned to +/- strand):  ref base A (36/59);  new base C (7/3);  total (43/62)
Fisher's exact test for biased strand distribution p-value = 8.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

TCAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTT  >  NC_000913/2592841‑2593139
                                                                                                                                                     |                                                                                                                                                     
tcAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATaa                                                                                                                                                       >  2:1376993/1‑150 (MQ=255)
 cAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAg                                                                                                                                                      >  2:567211/1‑150 (MQ=255)
  aaGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAgg                                                                                                                                                     <  2:1687369/150‑1 (MQ=255)
    gcgcGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCg                                                                                                                                                   <  2:1324469/150‑1 (MQ=255)
     cgcgTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGt                                                                                                                                                  <  2:430832/150‑1 (MQ=255)
          cgacactAACAAATTTCTCGCCAACAATCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGGAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAAcg                                                                                                                                             <  1:1219214/144‑1 (MQ=255)
                tAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCAt                                                                                                                                       >  1:1561728/1‑150 (MQ=255)
                 aaGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATc                                                                                                                                      >  2:1771067/1‑150 (MQ=255)
                  aGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                                                     <  1:1636298/150‑1 (MQ=255)
                   gAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCg                                                                                                                                    <  2:829009/150‑1 (MQ=255)
                    aaaTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                                                   <  1:1772731/150‑1 (MQ=255)
                    aaaTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                                                   <  1:833390/150‑1 (MQ=255)
                      aTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCa                                                                                                                                 <  2:407367/150‑1 (MQ=255)
                        ttCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCACCCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAAt                                                                                                                               <  1:1151547/150‑1 (MQ=255)
                         tctcGCCAACAAGCCTGAGCTTCCCGCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATg                                                                                                                              <  2:263543/150‑1 (MQ=255)
                            cGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTg                                                                                                                           >  2:1059393/1‑150 (MQ=255)
                            cGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTg                                                                                                                           <  2:542620/150‑1 (MQ=255)
                             gCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGa                                                                                                                          <  2:1711361/150‑1 (MQ=255)
                             gCCAACAAGCCGGAGCTTCCACCGCATCGTGATTTCAACGATAAGTGTGACGATGACGACGACTGGCCGAAAAATGATCAACCGAATAAGTAAGCTGAGTGCTCAACAGTAGTCCGGATAAGTCGTTAACGCATCATCAGGCAATGGTGa                                                                                                                          >  2:355210/1‑150 (MQ=255)
                              ccAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGaa                                                                                                                         >  2:755746/1‑150 (MQ=255)
                                aacaaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACg                                                                                                                       >  2:547982/1‑150 (MQ=255)
                                   aaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATg                                                                                                                    <  1:1397611/150‑1 (MQ=255)
                                      ccGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATACGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCt                                                                                                                 >  1:1148528/1‑150 (MQ=255)
                                       cGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTg                                                                                                                <  2:496244/150‑1 (MQ=255)
                                        ggAGCTTCCCCCGAATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACAGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGa                                                                                                               >  1:577743/1‑150 (MQ=255)
                                          aGCTTCCCCCGCATGGTGATTTCCTCGATAAGGGGGACGATGACGACGACTGGCCGAAAATGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATg                                                                                                             <  1:891673/150‑1 (MQ=38)
                                             ttCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGa                                                                                                          >  1:563642/1‑150 (MQ=255)
                                                 cccGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCt                                                                                                      >  1:564289/1‑150 (MQ=255)
                                                  ccGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTc                                                                                                     >  1:1382046/1‑150 (MQ=255)
                                                   cGCCTCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCa                                                                                                    <  1:317380/150‑1 (MQ=255)
                                                   cGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGCCTGGCCGAAAAAGGATCAACCGAAGAAGAAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCa                                                                                                    <  2:1209481/150‑1 (MQ=255)
                                                      aTCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCAcgc                                                                                                 <  2:750121/150‑1 (MQ=255)
                                                       tCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCAcgcg                                                                                                <  1:1380329/150‑1 (MQ=255)
                                                         gaggttacAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGCTCAACAGAAGAAGTAAGCTGAGTGCTCATCCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTc                                                                                              <  1:339738/143‑1 (MQ=38)
                                                         gTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTc                                                                                              <  1:1280589/150‑1 (MQ=255)
                                                         gTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATTGTGAACGATGCCTGATGCGACGCTCACGCGTc                                                                                              <  2:21449/150‑1 (MQ=255)
                                                         gTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTc                                                                                              <  1:1081498/150‑1 (MQ=255)
                                                         gTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTc                                                                                              >  1:93691/1‑150 (MQ=255)
                                                          tGATTTCAACGATAAGTGGGACGATGCCGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCt                                                                                             >  1:977599/1‑150 (MQ=255)
                                                              ttCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATc                                                                                         <  1:1552567/150‑1 (MQ=255)
                                                                 aaCGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAgg                                                                                      <  2:1391258/150‑1 (MQ=255)
                                                                  aCGATAAGTGGGACGATGACGACGACTGGCCGAACAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGc                                                                                     <  1:1446793/150‑1 (MQ=255)
                                                                    gATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCa                                                                                   <  1:1468959/150‑1 (MQ=255)
                                                                    gATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCa                                                                                   <  2:742932/150‑1 (MQ=255)
                                                                    gATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATAAGACGCGTGAGCGTCGCATCAGGcat                                                                                   <  1:623419/150‑3 (MQ=37)
                                                                     attaGTGGGACGTTGACTACGACGGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGATCACGCGTCTTATCAGGCCAg                                                                                  <  2:926725/147‑1 (MQ=37)
                                                                      tAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAgg                                                                                 <  2:1697216/150‑1 (MQ=255)
                                                                      tAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATAAGACGCGTGAGCGTCGCATCAGGcatcg                                                                                 <  2:1008708/150‑5 (MQ=37)
                                                                      tAAGTGGGACGATGACGACGACTGGCCGAAAAAGGAGCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAgg                                                                                 <  2:498336/150‑1 (MQ=255)
                                                                      aaaGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGACGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAgg                                                                                 <  1:259350/149‑1 (MQ=255)
                                                                        aGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTTAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACTCTCACGCGTCTTATCAGGCCAGGTa                                                                               <  1:1392668/150‑1 (MQ=255)
                                                                         gTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTaa                                                                              >  2:1510954/1‑150 (MQ=255)
                                                                            ggACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCt                                                                           >  2:188615/1‑150 (MQ=255)
                                                                               cGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTc                                                                        >  2:1671729/1‑150 (MQ=255)
                                                                                gATTACGACGACTGGCCGAAAAAGGATCAACCGAAGAACTAAGCTGAGTGCTCTACCGGAGGCCGGATACGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCt                                                                       <  2:315729/150‑1 (MQ=38)
                                                                                   tacgactacTGGCAGAAAACGTATCAACCGAAGAAGTATGCTGAGTGCTCAACAGTCGGCCGGATAAGGCGTGAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGAAGCTCTCGCGTCTTAGCAGGCCAGGTAATCTGTCtttt                                                                    <  2:1093439/149‑1 (MQ=37)
                                                                                   gacgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCtttt                                                                    <  2:657246/150‑1 (MQ=255)
                                                                                     cgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTTGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAc                                                                  <  1:1255134/150‑1 (MQ=255)
                                                                                     cgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAc                                                                  <  1:1822010/150‑1 (MQ=255)
                                                                                     cgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAc                                                                  >  1:1218159/1‑150 (MQ=255)
                                                                                              gCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCa                                                         <  1:1314000/150‑1 (MQ=255)
                                                                                               ccGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCaa                                                        >  2:799623/1‑150 (MQ=255)
                                                                                                cGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAAt                                                       <  2:156056/150‑1 (MQ=255)
                                                                                                 gAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATc                                                      <  2:135101/150‑1 (MQ=255)
                                                                                                  aaaaaGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATca                                                     <  1:1059393/150‑1 (MQ=255)
                                                                                                   aaaaGGATCACCCGAAGAAGTCAGCTGAGTGCTCAACCGTCGGCCGGATACGTCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATcac                                                    <  2:1221389/150‑1 (MQ=38)
                                                                                                   aaaaGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATcac                                                    >  2:708904/1‑150 (MQ=255)
                                                                                                      aGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatc                                                 <  2:1486485/150‑1 (MQ=255)
                                                                                                       ggATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCCTCCGGCAATGGTGAACGCTGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatca                                                >  1:335563/1‑150 (MQ=255)
                                                                                                       ggATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAGTCAcatca                                                >  1:1329392/1‑150 (MQ=255)
                                                                                                       ggATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatca                                                >  1:1617284/1‑150 (MQ=255)
                                                                                                        gATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGGCTTTTACATCATTTCAATCAcatcat                                               >  2:464107/1‑150 (MQ=255)
                                                                                                         aTCCACCGAAGAAGTAAGCGGAGGGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatcatc                                              <  2:143072/150‑1 (MQ=255)
                                                                                                         aTCAACCGAAGAAGTAAGCTGAGTGCTCACCCGTCGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatcatc                                              <  2:918933/150‑1 (MQ=255)
                                                                                                         aTCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTCAACGATGCCTGATGCGACGCTCCCGCGTCTTAGCAGGCCAGTTAAGCTGTCTTTTACATCATTTCAAGCCcaacatc                                              >  1:350685/1‑150 (MQ=25)
                                                                                                         aTCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCCGCAATGTTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCATTTACATCATTTCAATCAcatcatc                                              >  1:115794/1‑150 (MQ=255)
                                                                                                            aaCCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCCTCCGGCAATGGTGAACGATGCCTGATGCGACGCTCCCGCGTCTTATCAGGCCAGGTAATCTGTCTGTTACATCATTTCAATCACatcatcatc                                           >  1:141760/1‑150 (MQ=255)
                                                                                                                gaagaaGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCa                                       <  1:1336522/150‑1 (MQ=255)
                                                                                                                gaagaaGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCa                                       >  2:413036/1‑150 (MQ=255)
                                                                                                                 aagaagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTTAACGATGCCTGATGCGACGCTCACGCTTATTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCTCCAg                                      <  1:1137142/150‑1 (MQ=255)
                                                                                                                         aGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCg                              <  1:1404289/150‑1 (MQ=255)
                                                                                                                          gCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGAAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCCGGCCAGGTAATCTGTCGTTTACATCATTTCAATCCCATCATCATCGCCAGGCTTACCGc                             >  2:84235/1‑150 (MQ=255)
                                                                                                                              aGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACt                         <  1:956210/150‑1 (MQ=255)
                                                                                                                                 gCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAg                      >  1:1789382/1‑150 (MQ=255)
                                                                                                                                  cTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAgc                     >  1:261292/1‑150 (MQ=255)
                                                                                                                                    cAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAgcgc                   >  2:30366/1‑150 (MQ=255)
                                                                                                                                      aCCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCtt                 <  1:1771067/150‑1 (MQ=255)
                                                                                                                                       ccctAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACAGCCACTTAGCGCTtc                >  2:1574016/4‑150 (MQ=255)
                                                                                                                                         cattgCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAAGCCCATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtcgt              >  1:296896/5‑150 (MQ=255)
                                                                                                                                          tAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtcgtc             <  2:92209/150‑1 (MQ=255)
                                                                                                                                          attgCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATTCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtcgtc             >  2:767653/4‑150 (MQ=255)
                                                                                                                                          attgCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtcgtc             >  1:1028380/4‑150 (MQ=255)
                                                                                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCACTCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaa          >  1:1195104/1‑150 (MQ=255)
                                                                                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCCTCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaa          >  1:1626131/1‑150 (MQ=255)
                                                                                                                                             gCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACAACATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaa          <  2:1382046/150‑1 (MQ=255)
                                                                                                                                             gCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATACCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaa          <  1:1562615/150‑1 (MQ=255)
                                                                                                                                              ccGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaaa         >  1:1769483/1‑150 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaaa         >  2:1506027/1‑150 (MQ=255)
                                                                                                                                               cGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAaat        >  1:600814/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGtt  <  2:764686/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTT  >  NC_000913/2592841‑2593139

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