Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,816,611 C→A intergenic (‑62/‑65) rph ← / → yicC ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH/UPF0701 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,816,6110CA97.7% 317.4 / ‑3.4 89intergenic (‑62/‑65)rph/yicCribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH/UPF0701 family protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (41/45);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAA  >  NC_000913/3816465‑3816758
                                                                                                                                                  |                                                                                                                                                   
gCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTAt                                                                                                                                                 <  3:1081442/151‑1 (MQ=255)
    cGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTg                                                                                                                                             <  3:1432721/151‑1 (MQ=255)
    cGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTg                                                                                                                                             <  3:929598/151‑1 (MQ=255)
     gAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc                                                                                                                                               >  3:1206858/1‑148 (MQ=255)
     gAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTg                                                                                                                                             >  2:1377359/1‑150 (MQ=255)
       gCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGAc                                                                                                                                           >  2:558428/1‑150 (MQ=255)
             tGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCAt                                                                                                                                    >  3:989834/1‑151 (MQ=255)
              gCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATc                                                                                                                                   <  3:499370/151‑1 (MQ=255)
              gCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATc                                                                                                                                   >  2:90172/1‑151 (MQ=255)
              gCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATc                                                                                                                                   <  3:1250808/151‑1 (MQ=255)
               cATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCg                                                                                                                                  <  3:1605791/151‑1 (MQ=255)
               cATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCg                                                                                                                                  >  3:1188200/1‑151 (MQ=255)
               cATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCg                                                                                                                                  <  3:1099796/151‑1 (MQ=255)
                aTGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGa                                                                                                                                 <  3:383607/151‑1 (MQ=255)
                  gTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATg                                                                                                                               <  2:1528619/151‑1 (MQ=255)
                  gTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATg                                                                                                                               <  3:1592209/151‑1 (MQ=255)
                     ttGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGc                                                                                                                             <  2:1181353/150‑1 (MQ=255)
                     ttGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGc                                                                                                                             <  3:550213/150‑1 (MQ=255)
                       gTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCtat                                                                                                                          <  2:1009856/151‑1 (MQ=255)
                        tataGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCtat                                                                                                                          >  2:1231232/1‑150 (MQ=255)
                         ataGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCtata                                                                                                                         <  3:382006/150‑1 (MQ=255)
                         ataGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATaa                                                                                                                        <  2:929148/151‑1 (MQ=255)
                           aGTTGCGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAAt                                                                                                                       <  2:1327666/150‑1 (MQ=255)
                           aGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAAt                                                                                                                       <  3:7540/150‑1 (MQ=255)
                           aGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATc                                                                                                                      >  2:518536/1‑151 (MQ=255)
                            gTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATcc                                                                                                                     >  3:101309/1‑151 (MQ=255)
                              tACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCtt                                                                                                                   >  2:335667/1‑151 (MQ=255)
                               aCGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTc                                                                                                                  >  2:1079496/1‑151 (MQ=255)
                               aCGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTc                                                                                                                  <  2:888176/151‑1 (MQ=255)
                                  aGTCAGGGAAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGATTGAAACAAGTGTGGCTGCGCATTATACGGACTTCCGGCGGGTTTTAATGTCCTGAGAAGCCATCGATGGGTATAATTCtt                                                                                                                   >  3:1274928/1‑147 (MQ=255)
                                   gTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCAcc                                                                                                              >  3:795044/1‑151 (MQ=255)
                                     cAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACct                                                                                                             <  3:1411479/150‑1 (MQ=255)
                                     cAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACct                                                                                                             <  3:419831/150‑1 (MQ=255)
                                      aGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACctc                                                                                                            >  4:10134/1‑150 (MQ=255)
                                          tAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCtt                                                                                                       >  2:921111/1‑151 (MQ=255)
                                           aaCGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCttt                                                                                                      >  3:1183628/1‑151 (MQ=255)
                                           aaCGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCtt                                                                                                       <  2:1066749/150‑1 (MQ=255)
                                                gACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATaa                                                                                                 <  3:991173/151‑1 (MQ=255)
                                                     aCCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAgg                                                                                            <  3:586140/151‑1 (MQ=255)
                                                     aCCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAgg                                                                                            >  2:526998/1‑151 (MQ=255)
                                                      ccTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAgg                                                                                            <  3:1354938/150‑1 (MQ=255)
                                                         gATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACg                                                                                        >  2:1037637/1‑151 (MQ=255)
                                                            tattGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTc                                                                                      <  3:13930/150‑1 (MQ=255)
                                                             attGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCt                                                                                     <  3:15209/150‑1 (MQ=255)
                                                             attGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTa                                                                                    <  3:397121/151‑1 (MQ=255)
                                                                gCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATg                                                                                  >  2:1602143/1‑150 (MQ=255)
                                                                 cTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGa                                                                                 >  2:1269855/1‑150 (MQ=255)
                                                                  tACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATc                                                                               >  2:936363/1‑151 (MQ=255)
                                                                  tACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATc                                                                               >  3:1389789/1‑151 (MQ=255)
                                                                   aCGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATc                                                                               <  3:673627/150‑1 (MQ=255)
                                                                    cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCg                                                                             <  3:546911/151‑1 (MQ=255)
                                                                     ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGc                                                                            <  2:1439545/151‑1 (MQ=255)
                                                                      gCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGc                                                                            >  3:1316470/1‑150 (MQ=255)
                                                                         tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGt                                                                         >  3:228835/1‑150 (MQ=255)
                                                                                 gCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGcc                                                                <  2:195355/151‑1 (MQ=255)
                                                                                  cATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCt                                                               >  2:1165572/1‑151 (MQ=255)
                                                                                   atatTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTTTTACTATCCTGACAAGGCATCGATGGCTATAATCCTTTCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCt                                                               >  3:1294061/1‑150 (MQ=255)
                                                                                     atTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACg                                                            <  3:1232748/151‑1 (MQ=255)
                                                                                     atTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGAAACGTCTATGATCCGCAGTATGACCGCCTAc                                                             <  3:181477/150‑1 (MQ=255)
                                                                                       tGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGcc                                                          >  2:1040110/1‑151 (MQ=255)
                                                                                       tGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGcc                                                          >  2:936806/1‑151 (MQ=255)
                                                                                               ccGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTg                                                   <  3:663466/150‑1 (MQ=255)
                                                                                                cGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGaa                                                 <  2:691052/151‑1 (MQ=255)
                                                                                                cGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGa                                                  >  2:1109085/1‑150 (MQ=255)
                                                                                                 ggCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGaaa                                                <  3:805007/151‑1 (MQ=255)
                                                                                                      gAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCa                                             >  2:355173/1‑149 (MQ=255)
                                                                                                      gAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAg                                           <  2:1550537/151‑1 (MQ=255)
                                                                                                         aCAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATGCTGTCGAGGCATGGATGGCTGTAATTCTTCCATCTCTCCTTTTATAAACAGGAAGGGCTATGAACCGCAGGATGACCGCCTAGGCCCGGCGAGAAATCAAg                                           >  3:96175/1‑148 (MQ=255)
                                                                                                         aCAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAggg                                         >  3:1240589/1‑150 (MQ=255)
                                                                                                          cAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTg                                       <  2:194250/151‑1 (MQ=255)
                                                                                                          cAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTg                                       >  3:281506/1‑151 (MQ=255)
                                                                                                            aaTGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGa                                      >  3:1548966/1‑150 (MQ=255)
                                                                                                            aaTGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGa                                      >  2:221257/1‑150 (MQ=255)
                                                                                                             aTGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAAt                                    <  3:375827/151‑1 (MQ=255)
                                                                                                               gtgGCTGCGCATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATg                                   <  2:760847/150‑1 (MQ=255)
                                                                                                                       gcATTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAgcgc                          >  2:1491996/1‑151 (MQ=255)
                                                                                                                         aTTATACGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAgcgc                          >  3:1215404/1‑149 (MQ=255)
                                                                                                                            ataCGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAAcc                      <  3:1326980/150‑1 (MQ=255)
                                                                                                                            ataCGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCt                     >  2:577538/1‑151 (MQ=255)
                                                                                                                            ataCGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCt                     <  2:415172/151‑1 (MQ=255)
                                                                                                                               cGGACTTCCGGCGGTTATTTCTATCCTGACAAGGCATCGATGGCTACAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCCCCTACGCCCGGCGTGAAATAAAGGGTGGATGCGGCAGCGCAACCt                     >  3:632019/1‑148 (MQ=255)
                                                                                                                               cGGACTTCCGGCGGTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTGTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTggg                  <  3:741624/151‑1 (MQ=255)
                                                                                                                                ggACTTCCGGCGGTTATTCCTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTggg                  >  3:694339/1‑150 (MQ=255)
                                                                                                                                       cggcggTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGGGAATGGGGGAGCGCAACCTGGGAAATgcgc          >  2:955289/1‑151 (MQ=255)
                                                                                                                                       cggcggTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGGGAATGGGGGAGCGCAACCTGGGAAATgcg           >  2:89519/1‑150 (MQ=255)
                                                                                                                                       cggcggTTATTACTATCCTGACAAAGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGGGAATGGGGGAGCGCAACCTGGGAAATgcg           >  3:1030615/1‑150 (MQ=255)
                                                                                                                                        ggcggTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCt         <  3:1584062/151‑1 (MQ=255)
                                                                                                                                         gcggTTATTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTc        >  3:1325592/1‑151 (MQ=255)
                                                                                                                                               aTTACTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTaaa  <  3:837187/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGCATCGATGGCTATAATCCTTCCACCTCTCCTTTTATAAACAGGAACGTCTATGATCCGCAGTATGACCGCCTACGCCCGGCGTGAAATCAAGGGTGAATGGGGGAGCGCAACCTGGGAAATGCGCTCGGTAAA  >  NC_000913/3816465‑3816758

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: ‑