Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,093 G→T 10.7% D681Y (GAT→TAT)  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,0930GT10.7% 1151.8 / 19.4 377D681Y (GAT→TAT) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (126/197);  new base T (51/3);  total (177/200)
Fisher's exact test for biased strand distribution p-value = 1.25e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 378 total aligned reads displayed.

GGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACG  >  NC_000913/4014945‑4015242
                                                                                                                                                    |                                                                                                                                                     
 gCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGAt                                                                                                                                                     <  2:6040913/150‑1 (MQ=255)
  cGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGAtt                                                                                                                                                    <  1:305870/150‑1 (MQ=255)
   gAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGAttt                                                                                                                                                   >  1:264550/1‑149 (MQ=255)
   gAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTa                                                                                                                                                   <  2:5727613/150‑1 (MQ=255)
     aaGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATc                                                                                                                                                 <  2:1851720/150‑1 (MQ=255)
      aGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATcc                                                                                                                                                <  1:4398643/150‑1 (MQ=255)
       gATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCa                                                                                                                                               <  2:832370/150‑1 (MQ=255)
         tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCaaa                                                                                                                                             <  2:2566237/150‑1 (MQ=255)
          gACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGCTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTGATTATCCCAAt                                                                                                                                            >  2:5458836/1‑150 (MQ=255)
            cacacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTTATTATCAAAATGa                                                                                                                                          >  1:5427586/1‑150 (MQ=255)
               acaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAAt                                                                                                                                       >  2:4839238/1‑150 (MQ=255)
                 aaaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCg                                                                                                                                     <  2:6765514/150‑1 (MQ=255)
                 aaaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAATCg                                                                                                                                     >  2:1095988/1‑150 (MQ=255)
                  aaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCgg                                                                                                                                    <  1:2848956/150‑1 (MQ=255)
                   aTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGt                                                                                                                                   <  1:2271883/150‑1 (MQ=255)
                     ccACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAATAAATCGGTcc                                                                                                                                 >  2:6664744/1‑150 (MQ=255)
                      cacaCTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCt                                                                                                                                <  2:1970500/150‑1 (MQ=255)
                         acTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATCCAAATGAAATCGGTCCAGGc                                                                                                                             >  2:6473344/1‑150 (MQ=255)
                         acTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTATACAAATGAAATCGGTCCTGGc                                                                                                                             >  1:1003804/1‑150 (MQ=255)
                           tCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTTATTATACAAATGAAAACGGTCCTGGCgg                                                                                                                           >  1:3422731/1‑149 (MQ=255)
                           tCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTTATTATCCAAATGAAATCGGGCCTGGCgg                                                                                                                           >  1:5693778/1‑149 (MQ=255)
                             acaTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCt                                                                                                                         >  1:6568938/1‑150 (MQ=255)
                              caTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTAAGATTTTTATTATACAAAAGAAATAGGTCCAGGCGTCTa                                                                                                                        >  2:2685723/1‑150 (MQ=255)
                               aTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGGCGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTAt                                                                                                                       <  2:1614088/150‑1 (MQ=255)
                                 gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTTATTATACAAATGAAAACGGTACTGGCGTCTAtta                                                                                                                     >  2:2931981/1‑148 (MQ=255)
                                 gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                                                                                                     <  2:2387569/150‑1 (MQ=255)
                                   gtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAAGAAAACGGTCCTGGCGGCTATGACa                                                                                                                   >  2:775663/1‑150 (MQ=255)
                                    ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACAt                                                                                                                  <  1:4152236/150‑1 (MQ=255)
                                      aTTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTTTGACATTc                                                                                                                <  1:3503409/150‑1 (MQ=255)
                                         gCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATACAAATGAAATCGGTCCAGGCGGCTATGACAATCAca                                                                                                             >  2:4183901/1‑149 (MQ=255)
                                          cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCGGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTc                                                                                                            <  2:4516377/150‑1 (MQ=255)
                                             gTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGCAGAGTTTGATAAACCAAATGAAATCGGTCCAGGCGTCTATGACAATACCTCGcc                                                                                                         >  1:4722650/1‑150 (MQ=255)
                                                cAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaaa                                                                                                      <  1:3610694/150‑1 (MQ=255)
                                                 aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTTATTATCCAAATGAAATCGGTCCTGGCGGCGATGACATTAACGCGCCCAAc                                                                                                     >  1:2969399/1‑150 (MQ=255)
                                                 aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAAc                                                                                                     >  2:2819387/1‑150 (MQ=255)
                                                 aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAAc                                                                                                     <  2:5077889/150‑1 (MQ=255)
                                                 aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGAATGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAAc                                                                                                     <  2:2102926/150‑1 (MQ=255)
                                                   cGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGt                                                                                                   <  1:6304512/150‑1 (MQ=255)
                                                    gACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTa                                                                                                  <  1:6046732/150‑1 (MQ=255)
                                                    gACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGGTGGAAGAGTTTGATTATCCAAATTAAAAAGGGCCTGGCGGCTATGGCATTCACACGCCAAACGtt                                                                                                  >  2:6369122/1‑149 (MQ=255)
                                                      catcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTAcc                                                                                                <  2:1039176/150‑1 (MQ=255)
                                                        tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAATAGTTTGATTATCCAAAAGAAATACGTCCTGGCGTCTATGACATTCACACGCCAAACCTaacca                                                                                              >  2:3184092/1‑146 (MQ=255)
                                                           tGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGGCTATGACAATCACTCGCCAAACGTACCGAGCg                                                                                           >  2:1383881/1‑150 (MQ=255)
                                                            ggATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGt                                                                                          <  2:3061341/150‑1 (MQ=255)
                                                                tCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGaa                                                                                      <  2:2763299/150‑1 (MQ=255)
                                                                  gATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATg                                                                                    <  2:3758083/150‑1 (MQ=255)
                                                                  gATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAAAGAAATCGGTCCAGGCGGCGATGTCATTAACGCGCCAAACGTACCGCGCGTGGacag                                                                                    >  2:5142632/1‑147 (MQ=255)
                                                                    ttGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGa                                                                                  <  1:962837/150‑1 (MQ=255)
                                                                       cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTg                                                                               <  2:4698809/150‑1 (MQ=255)
                                                                       cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTg                                                                               <  1:3973563/150‑1 (MQ=255)
                                                                       cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTg                                                                               <  1:4830286/150‑1 (MQ=255)
                                                                        ggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGGTTATACAAAAGAAATAGGTCCAGGCGTCTATGACAGTAAAACGCCAAAAGTACCCAGCGTGGCATGGAttta                                                                              >  1:3354489/1‑148 (MQ=255)
                                                                           gcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGc                                                                           <  2:5877964/150‑1 (MQ=255)
                                                                            cTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTTTTTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGTATGGATTGAAGcc                                                                          >  2:3261614/1‑150 (MQ=255)
                                                                             tGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTTGAGTGGTTTTTTATCCAAAAGAAATAGGGCCAGGGGGCGATGGCCATCAATCGCCAAACGTTCCGAGCGTGGGATGGATTGTAGCCt                                                                         >  1:2586949/1‑150 (MQ=255)
                                                                             tGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTTGAGTTTTATTATACAAATTAAATAGGTCCAGGCGGCTATGAAAATCACACGCCAAACGTACCCAGCGTGGAATGGATTGAATccg                                                                         >  1:146175/1‑149 (MQ=255)
                                                                                 cGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTTTTTATACAAATTAAATCGGTCCAGGCGTCGATGACAATCACACGCCAAACGTTACGGGCGGGGAATGGTTTGAAGCCTTGCt                                                                     >  1:6131825/1‑150 (MQ=255)
                                                                                    aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  <  2:5095091/150‑1 (MQ=255)
                                                                                    aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  <  1:910240/150‑1 (MQ=255)
                                                                                    aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  <  1:2685723/150‑1 (MQ=255)
                                                                                      aCGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaaga                                                                <  2:147076/150‑1 (MQ=255)
                                                                                        gTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGGGGGTTTGGTTATGCCAATGAAAAAGGGCCGGGCGGCGTTGGACGTCACGAGCCCAACGTCCCGCGCGAGGGATGGCTTGAATCCGTGCTGAAGaaa                                                              >  1:598417/1‑150 (MQ=255)
                                                                                         tcatcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAg                                                             <  2:4916347/150‑1 (MQ=255)
                                                                                          catcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGc                                                            <  2:6189065/150‑1 (MQ=255)
                                                                                           atcaCCATCGACACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCg                                                           <  1:5365465/150‑1 (MQ=255)
                                                                                           atcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTATCCAAATGAAATAGGTCCTGGCGTCTATGGCATCCACTCGCCAAACGTAACGAGCGTGGAAAGGATTGAAGCCTTGCTGAAGAAATCg                                                           >  2:3001223/1‑150 (MQ=255)
                                                                                            tcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCgg                                                          <  1:6640096/150‑1 (MQ=255)
                                                                                            tcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCgg                                                          <  2:3680499/150‑1 (MQ=255)
                                                                                              aCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTGTCCAAATGAAATAGGTCCTGGCGTCTATTGCATACACTCGCCAAACGTAACGAGCGTGGCATGGGTTGAAGCCATGCTGAAGCAAGCGGCa                                                        >  2:1802362/1‑150 (MQ=255)
                                                                                               ccATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCaa                                                       <  1:1658700/150‑1 (MQ=255)
                                                                                               ccATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGGCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCaa                                                       >  2:4357522/1‑150 (MQ=255)
                                                                                                 aTCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATACAAATGAAATCGGTCCTGGCCTCTATGACATTCACACGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCaaaa                                                     >  2:4795698/1‑150 (MQ=255)
                                                                                                  tCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATACAAAAGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAAc                                                    >  2:3331293/1‑150 (MQ=255)
                                                                                                  tCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAAc                                                    <  2:1232654/150‑1 (MQ=255)
                                                                                                   cGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACg                                                   <  1:3787512/150‑1 (MQ=255)
                                                                                                           cgcgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCgg                                           <  2:6005897/150‑1 (MQ=255)
                                                                                                             cgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCa                                         <  1:6411882/150‑1 (MQ=255)
                                                                                                             cgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCa                                         <  1:3702883/150‑1 (MQ=255)
                                                                                                               ttCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCaga                                       <  1:2180879/150‑1 (MQ=255)
                                                                                                               ttCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCaga                                       <  1:1910627/150‑1 (MQ=255)
                                                                                                                  cGACATGGAGTTGCTGGAGTCGGGTGAAGAGTTTGATTATCCAAAAGAAATAGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAAAGGATTGAAGCCTTTCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAgcg                                    >  2:5022893/1‑150 (MQ=255)
                                                                                                                   gACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATTACAATAACTCGCCAAACGTACCGAGCGTGGAATTGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCAATCCGGCAGAgcgc                                   >  2:348664/1‑150 (MQ=255)
                                                                                                                    aCATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTACTGGCGTCTATGACAATCAAACGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGcccc                                  >  2:232809/1‑150 (MQ=255)
                                                                                                                        ggAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCtgtg                              <  1:455095/150‑1 (MQ=255)
                                                                                                                         gAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGTAAAAGGTCCTGGCGTCTATTACATTCACAAGCCAAACGTAACCAGCGTGGAAAGGATTTAAGCCATGCTGAAGAAAGCGGCAAAAAGCAATACGGCAGAGCGCCGGTgg                             >  2:1728744/1‑150 (MQ=255)
                                                                                                                          aGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTggg                            <  1:1363467/150‑1 (MQ=255)
                                                                                                                           gTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATTAAATAGGTCCTGGCGTCTATGACATTCACTAGCCAAAAGTACCCAGCGTGGAATGGATTGAAGCCTTTCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGt                           >  2:5336792/1‑150 (MQ=255)
                                                                                                                           gTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGt                           <  1:510653/150‑1 (MQ=255)
                                                                                                                           gTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGt                           <  1:3273645/150‑1 (MQ=255)
                                                                                                                              gCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCaa                        <  2:4968739/150‑1 (MQ=255)
                                                                                                                               cTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAAACGGTCCTGGCGTCTATGACATTCACACGCCAAAAGTAACGAGCGTGGAATGGATTGAAGCCATTCTGAAGAAATCGGCAAAACGCAATACGGCAGAGCGCCAGTGGGTCAAc                       >  2:965071/1‑150 (MQ=255)
                                                                                                                               cTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAAc                       <  1:2324890/150‑1 (MQ=255)
                                                                                                                                tGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAAcc                      <  1:27840/150‑1 (MQ=255)
                                                                                                                                    gTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa                  <  1:2930992/150‑1 (MQ=255)
                                                                                                                                       gTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtg               >  1:77592/1‑150 (MQ=255)
                                                                                                                                           gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc           <  1:5073598/150‑1 (MQ=255)
                                                                                                                                           gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc           <  1:1750442/150‑1 (MQ=255)
                                                                                                                                            aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          >  2:5283405/1‑150 (MQ=255)
                                                                                                                                            aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          >  2:1949740/1‑150 (MQ=255)
                                                                                                                                                gTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaa      <  2:3860039/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACG  >  NC_000913/4014945‑4015242

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