Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,085 A→T 9.5% E678V (GAA→GTA)  metE → 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,0850AT 9.5% 1978.9 / 28.7 671E678V (GAA→GTA) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (192/393);  new base T (57/2);  total (276/395)
Fisher's exact test for biased strand distribution p-value = 2.89e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 684 total aligned reads displayed.

CGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCC  >  NC_000913/4014936‑4015234
                                                                                                                                                     |                                                                                                                                                     
cgccgcCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGa                                                                                                                                                       <  2:912471/150‑1 (MQ=255)
 gccgccGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGaa                                                                                                                                                      <  1:9249486/150‑1 (MQ=255)
  ccgccgTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTag                                                                                                                                                     >  1:4935995/1‑150 (MQ=255)
    gccgTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAagag                                                                                                                                                   <  1:8989133/150‑1 (MQ=255)
     ccgTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGt                                                                                                                                                  <  1:176080/150‑1 (MQ=255)
       gTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGttt                                                                                                                                                <  1:9816701/150‑1 (MQ=255)
         ggCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                              <  2:9639900/150‑1 (MQ=255)
         ggCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                              <  2:1410475/150‑1 (MQ=255)
         ggCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                              <  1:1020623/150‑1 (MQ=255)
          gCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGCAGTGTTTGAt                                                                                                                                             >  2:1691883/1‑150 (MQ=255)
            gAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTGAGAGTTttatta                                                                                                                                           >  2:9277067/1‑150 (MQ=255)
            gAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTa                                                                                                                                           >  2:7045563/1‑150 (MQ=255)
               aGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATcc                                                                                                                                        <  2:10583307/150‑1 (MQ=255)
                  tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCaaa                                                                                                                                     >  1:7567483/1‑150 (MQ=255)
                   gACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAATAGTTTTATTATACAAAt                                                                                                                                    >  2:10469857/1‑150 (MQ=255)
                   gACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAt                                                                                                                                    <  1:7811782/150‑1 (MQ=255)
                   gACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAt                                                                                                                                    <  2:4443731/150‑1 (MQ=255)
                   gACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAt                                                                                                                                    <  2:10351839/150‑1 (MQ=255)
                   gACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGGGGCGGGGTTGATTATCCCAAt                                                                                                                                    >  2:2951650/1‑150 (MQ=255)
                    acacacAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATg                                                                                                                                   >  2:2494724/1‑150 (MQ=255)
                        acaACTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAAt                                                                                                                               <  1:1954714/150‑1 (MQ=255)
                        acaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATACAAATGAAAt                                                                                                                               >  1:2667168/1‑150 (MQ=255)
                            aTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATAGGt                                                                                                                           >  1:6086125/1‑150 (MQ=255)
                             tCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATACAAATGAAATCGGTc                                                                                                                          >  1:2770045/1‑150 (MQ=255)
                             tCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTc                                                                                                                          <  2:1888206/150‑1 (MQ=255)
                              ccACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTcc                                                                                                                         <  1:7739666/150‑1 (MQ=255)
                               cacaCTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGaa                                                                                                                                                      >  1:2957019/1‑120 (MQ=255)
                                acacTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTATGCAAATTTAATAGGTCCTg                                                                                                                       >  1:10100736/1‑150 (MQ=255)
                                        aTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAATGAAATAGGTCCTGGCGTCTAt                                                                                                               >  2:8669172/1‑150 (MQ=255)
                                          gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                                                                                             <  1:42013/150‑1 (MQ=255)
                                          gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAAGAAAACCGTCCAGGCCTCAATGa                                                                                                             >  2:3626202/1‑150 (MQ=255)
                                            gtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACa                                                                                                           <  1:3966032/150‑1 (MQ=255)
                                                ttGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAATAGGTCCTGGCGTCTATGACATTCa                                                                                                       >  2:6062059/1‑150 (MQ=255)
                                                 tGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCAc                                                                                                      <  2:7735516/150‑1 (MQ=255)
                                                   cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAATCGGTCCTGGCGTCTATGAAATTCAcac                                                                                                    >  2:1502048/1‑148 (MQ=255)
                                                    gAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATTAAATCGGTCCTGGCGGCGATGACATTAACTCg                                                                                                   >  1:3442853/1‑150 (MQ=255)
                                                     aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                  <  2:9201282/150‑1 (MQ=255)
                                                      gTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTGATTATACAAAAGAAATCGGTACTGGCGTCTATGACAATCACACGcc                                                                                                 >  1:250058/1‑150 (MQ=255)
                                                       ttCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATCCAAAATAAATAGGTCCTGGCGGCGATGACATTCACTCGCCa                                                                                                >  1:10426017/1‑150 (MQ=255)
                                                         cAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaaa                                                                                              <  2:3973651/150‑1 (MQ=255)
                                                         cAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaaa                                                                                              >  2:5769833/1‑150 (MQ=255)
                                                          aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTGATTATACAAATGAAAACGGTCCTGGCGTCTATTACATTCACTCGCCAAAc                                                                                             >  1:1281860/1‑150 (MQ=255)
                                                               catcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAAGAAAACGGTCCAGGCGTCTATGACAATCACACGCCAAACGTAcc                                                                                        >  2:7948499/1‑150 (MQ=255)
                                                                atcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCg                                                                                       <  1:4495737/150‑1 (MQ=255)
                                                                 tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGa                                                                                      <  2:4216242/150‑1 (MQ=255)
                                                                 tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTATACAAATTAAATCGGTCCTGGCGTCTATGAAAATAACTCGCCAAACGTACCGa                                                                                      >  1:10562623/1‑150 (MQ=255)
                                                                   atGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGGGGAGAGGTTGATTATCCAAAAGAAATCCGTCCAGGCGGCGATGACCATCACACGCCCAACCGACCgcgc                                                                                    >  2:3770133/1‑150 (MQ=255)
                                                                   atGGACTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGc                                                                                    >  2:2364034/1‑150 (MQ=255)
                                                                     ggATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGt                                                                                  <  1:8734065/150‑1 (MQ=255)
                                                                      gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTg                                                                                 <  2:10032547/150‑1 (MQ=255)
                                                                           gATTGCGGCGCTGGACGCAGACGTCATCACCATCGACACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACAGTCACTCGCCAAACGTACCGAGCGTGGAATg                                                                            <  2:10699907/150‑1 (MQ=255)
                                                                           gATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATg                                                                            <  1:9379897/150‑1 (MQ=255)
                                                                             ttGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGa                                                                          <  2:4966472/150‑1 (MQ=255)
                                                                             ttGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGa                                                                          <  2:7517670/150‑1 (MQ=255)
                                                                             ttGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGa                                                                          <  2:586354/150‑1 (MQ=255)
                                                                             ttGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTTAAGAGTTTGATTATCCAAAAGAAAACGGTCCTGGCGTCTATGACATTCACACGCCAAACCGACCCAGCGGGGAAAggc                                                                          >  2:3384151/1‑149 (MQ=255)
                                                                               gcggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGGTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGAtt                                                                        <  1:7461544/150‑1 (MQ=255)
                                                                                cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTg                                                                       <  2:8433331/150‑1 (MQ=255)
                                                                                  gcgcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGaa                                                                     >  2:7356331/1‑150 (MQ=255)
                                                                                       ggACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGTGTTTTATTATCCAAATGTAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAAAGGATTTAAGCCtt                                                                >  2:7312481/1‑150 (MQ=255)
                                                                                       ggACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCtt                                                                <  1:10763405/150‑1 (MQ=255)
                                                                                          cGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCt                                                             <  2:6483668/150‑1 (MQ=255)
                                                                                            cAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTga                                                           <  1:5967924/150‑1 (MQ=255)
                                                                                             aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                          <  1:2175210/150‑1 (MQ=255)
                                                                                             aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                          <  2:515730/150‑1 (MQ=255)
                                                                                              gACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaag                                                         <  1:3244163/150‑1 (MQ=255)
                                                                                              gACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCCAGCGTGGAATGGATTGAAGCCTTGCTgaag                                                         >  2:1202119/1‑150 (MQ=255)
                                                                                                 gTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAATAGTTTGTTTATCCAAAAGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGaaa                                                      >  2:5444099/1‑150 (MQ=255)
                                                                                                  tcatcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAg                                                     <  1:10000487/150‑1 (MQ=255)
                                                                                                  tcatcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAg                                                     <  1:5684659/150‑1 (MQ=255)
                                                                                                    atcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCg                                                   <  1:5769833/150‑1 (MQ=255)
                                                                                                     tcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCgg                                                  <  2:9424980/150‑1 (MQ=255)
                                                                                                     tcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTGTTTATACAAAAGAAAAAGGGCCTGGCGGCGATGACAATCACACGCCAAAACGAACGAGCGTGGAAATGATTGGAGCCATGCTGGAGAAAGCgg                                                  >  1:5198555/1‑150 (MQ=255)
                                                                                                       aCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                              >  2:6802622/1‑56 (MQ=255)
                                                                                                       aCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                              <  1:6802622/56‑1 (MQ=255)
                                                                                                         cATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCGATGACATTCACTCGCCAAACGTAACGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCaaa                                              >  2:10805383/1‑150 (MQ=255)
                                                                                                           tCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTTATTATACAAATGAAATAGGTCCTGGCGTCGATGACAATCACACGCCAAACGTACCGAGCGTGGAAATGATTGAAGCCTTGCTGAAGAAAGCGGCAAAAc                                            >  2:2619260/1‑150 (MQ=255)
                                                                                                           tCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTAACCAAATTAAAAAGGTCCTGGCGGCTATGACATTCACACGCCAAACCTACCGAGCGTGGAAAGGGTTGAAGCCGTGCTGAAGAAAGCGGCaaaaa                                            >  2:4126271/1‑149 (MQ=255)
                                                                                                           tCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGGGAAGAGTTTTATTATACAAATGAAAAAGGTCCGGGCGGGTATGACAATCACACGCCCAAACTACCGAGCGTGGGATGGTTTGTAGCCATGCTGAAGAAATCGGCCAAAc                                            >  1:8056272/1‑150 (MQ=255)
                                                                                                            cGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACAATCACACGCCAAACCTACCGAGCGTGGAATGGCTTGAAGCCGTGCTGAAGAAATCGGCaaaaag                                           >  2:4705796/1‑148 (MQ=255)
                                                                                                                  cTCGCGTTCCGACCTGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTcc                                     <  2:6436299/150‑1 (MQ=255)
                                                                                                                   tCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCg                                    <  1:7318570/150‑1 (MQ=255)
                                                                                                                    cgcgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCgg                                   <  1:702641/150‑1 (MQ=255)
                                                                                                                        ttCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCaga                               <  2:3139913/150‑1 (MQ=255)
                                                                                                                                  gAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTgg                     <  1:6465021/150‑1 (MQ=255)
                                                                                                                                       gCTGGAGTCGTTTTAAGAGTTTTATTATACAAATTAAAAAGGTCCAGGCGTCTATGAAAATAAAACGCCAAACGTACCCAGCGTGGAAAGGATTTAAACCATTCAGAAGAAAGCGGCAAAAACCAATACGGCAGAGCGCCTGTGGGTCaa                >  2:10255463/1‑150 (MQ=255)
                                                                                                                                         tGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAAcc              <  2:5681629/150‑1 (MQ=255)
                                                                                                                                          ggAGTCGTTTGAAGTGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAAccc             <  1:8452714/150‑1 (MQ=255)
                                                                                                                                           gAGTCGTTTGAAGAGTTTTATTATACAAATGAAATAGGTCCTGGCGTCTATGACAATCACACGCCAAACGTACCCAGCGTGGAATGGATTTAAGCCTTTCTTAAGAAAGCGGCAAAAAGCAATACGGCACAGCGCCAGTGGGTCAACCCg            >  2:6173742/1‑150 (MQ=255)
                                                                                                                                              tCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGAc         <  2:9711903/150‑1 (MQ=255)
                                                                                                                                               cGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACt        <  2:8147909/150‑1 (MQ=255)
                                                                                                                                                gTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtg       >  2:8860182/1‑150 (MQ=255)
                                                                                                                                                 tttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgt      <  1:8295723/150‑1 (MQ=255)
                                                                                                                                                  ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg     <  2:2085168/150‑1 (MQ=255)
                                                                                                                                                  ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg     >  1:308606/1‑150 (MQ=255)
                                                                                                                                                  ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg     <  1:7948499/150‑1 (MQ=255)
                                                                                                                                                  ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGAGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGTATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg     <  1:5090732/150‑1 (MQ=255)
                                                                                                                                                   tGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTgg    >  2:769132/1‑150 (MQ=255)
                                                                                                                                                    gAAGAGTTTGATTATCCACATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc   <  1:3770133/150‑1 (MQ=255)
                                                                                                                                                    gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc   <  1:1448766/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCC  >  NC_000913/4014936‑4015234

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