Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,093 G→T 10.3% 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.3% 1805.5 / 30.1 576D681Y (GAT→TAT) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (156/339);  new base T (79/2);  total (235/341)
Fisher's exact test for biased strand distribution p-value = 2.54e-32
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 583 total aligned reads displayed.

TGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACG  >  NC_000913/4014944‑4015242
                                                                                                                                                     |                                                                                                                                                     
 ggCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                                      <  1:7718488/150‑1 (MQ=255)
 ggCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                                      <  1:984413/150‑1 (MQ=255)
   cGAAAGATGACACACACCTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGAtt                                                                                                                                                    <  2:7403122/150‑1 (MQ=255)
   cGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGAtt                                                                                                                                                    <  2:3416893/150‑1 (MQ=255)
    gAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGTGTTTGAttt                                                                                                                                                   >  1:3385765/1‑149 (MQ=255)
    gAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGGGTTTGTTTa                                                                                                                                                   >  1:5402950/1‑150 (MQ=255)
    gAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTa                                                                                                                                                   >  2:4542325/1‑150 (MQ=255)
    gAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTa                                                                                                                                                   >  1:1363095/1‑150 (MQ=255)
          tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCaaa                                                                                                                                             <  2:3984134/150‑1 (MQ=255)
              acacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGaa                                                                                                                                         <  1:5538609/150‑1 (MQ=255)
                acaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAATAGTTTGATTATACAAATGaaaa                                                                                                                                       >  2:3298663/1‑149 (MQ=255)
                acaAATCCACACTCACATGTGGTATTGCGAGTTCAACGACATCATGGATTCGAGTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGGGTGAAGAGGTTGATTATCCAAAGGAAAt                                                                                                                                       >  2:7514411/1‑150 (MQ=255)
                 caAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATc                                                                                                                                      <  1:7588903/150‑1 (MQ=255)
                  aaaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCg                                                                                                                                     <  1:192436/150‑1 (MQ=255)
                   aaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCgg                                                                                                                                    <  2:6878733/150‑1 (MQ=255)
                      ccACACTCACATTTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACTTGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAATCGGTcc                                                                                                                                 <  1:3842369/150‑1 (MQ=255)
                       cacaCTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGGTTGAATATCCAAAAGAAAACCGTccg                                                                                                                                >  2:1687553/1‑149 (MQ=255)
                        acacTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTg                                                                                                                               <  2:3136856/150‑1 (MQ=255)
                            tCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGt                                                                                                                           <  1:1753021/150‑1 (MQ=255)
                            tCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCACATGAAATCGGTGCAGGCGt                                                                                                                           >  1:4451085/1‑150 (MQ=255)
                            tCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATACAAATGAAATAGGTCCTGGCGt                                                                                                                           >  1:1254490/1‑150 (MQ=255)
                               caTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTa                                                                                                                        <  2:6435283/150‑1 (MQ=255)
                                  gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                                                                                                     <  1:3083083/150‑1 (MQ=255)
                                  gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                                                                                                     <  2:3757825/150‑1 (MQ=255)
                                  gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                                                                                                     <  2:6216636/150‑1 (MQ=255)
                                    gtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACa                                                                                                                   <  2:5163642/150‑1 (MQ=255)
                                     ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACCTGGAGTTGCTGGAGTCGTTTGTCGGGTTTTTTTATCCAAATGTAATCGGTCCGGGCGTCGATGGCAt                                                                                                                  >  1:712438/1‑150 (MQ=255)
                                     ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTTTTATACAAAAGAAAACGGTCCTGGCGTCTATGACAt                                                                                                                  >  1:4174933/1‑150 (MQ=255)
                                         tGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCAc                                                                                                              <  1:580173/150‑1 (MQ=255)
                                         tGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGATTTTGATTATCCAAAAGAAATCTGTACTGGCGTCTATTGCATTCAc                                                                                                              >  2:1553453/1‑150 (MQ=255)
                                          gCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAATAGGTCCTGGCGTCTATGACATTAAca                                                                                                             >  1:8084471/1‑149 (MQ=255)
                                           cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTTTTTATCCAAATTAAATCGGTCCTGGCGTCTATGACATTCACTc                                                                                                            >  2:6336792/1‑150 (MQ=255)
                                             aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                          <  1:5735188/150‑1 (MQ=255)
                                             aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                          <  1:3298690/150‑1 (MQ=255)
                                             aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                          <  2:1131643/150‑1 (MQ=255)
                                             aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                          <  1:7900455/150‑1 (MQ=255)
                                                tCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaa                                                                                                       <  2:1866101/150‑1 (MQ=255)
                                                   aCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGGGTTTGATTATCCAAATTAAAACGGTCCTGGCGTCTATGACATTCACTCGCCAAACg                                                                                                    >  1:4783562/1‑150 (MQ=255)
                                                         tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGa                                                                                              <  2:4183278/150‑1 (MQ=255)
                                                         tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTGTTTATACAAAAGAAATAGGGCCTGGCGTCTGTGACATTAACACGCCAAACGTAccca                                                                                              >  2:3502648/1‑148 (MQ=255)
                                                           atGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTTATTATCCAAATGAAATCGGTCCAGGCGTCGAATACATTCACTCGCCCAACGTACCGAGc                                                                                            >  1:7823238/1‑150 (MQ=255)
                                                              gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGTGTTTTTTTATGCAAAAGAAAACGGTCCGGGCGTCGATGACATTAACTCGCCAAACGTACCGAGCGTg                                                                                         >  2:300123/1‑150 (MQ=255)
                                                              gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAATAGTTTTATTATACAAAATAAAAAGGTCCAGGCGTATATGACAATAACACGCCAAACCTAACCAGCGTg                                                                                         >  1:5049650/1‑150 (MQ=255)
                                                              gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGGTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTg                                                                                         <  2:3067138/150‑1 (MQ=255)
                                                               aTTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTAAGAGTTTGATTATACAAATGTAATAGGTGCTGGCGTGTATGACATTCAATAGCCAAACGTAACGAGCGTga                                                                                        >  2:1337980/1‑149 (MQ=255)
                                                               aTTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATACAAATGAAAACGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTgg                                                                                        >  2:4724664/1‑150 (MQ=255)
                                                                 tCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGaa                                                                                      <  1:5905958/150‑1 (MQ=255)
                                                                 tCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCAAAATGAAATCGGTCCTGGCGTCTATGACAATAACAAACCAAACGTACCGAGCGTGGaa                                                                                      >  1:6827557/1‑150 (MQ=255)
                                                                    aTTGCGGCGCTGGACGCAGTCGTCATCACCATCGAAACCTCGCGTTCCGACTTGGAGTTGCTGGAGTCGTTTGATGAGTTTGATTATCCTAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATgg                                                                                   <  2:8146275/150‑1 (MQ=255)
                                                                    aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATACAAATGAAATAGGTCCTGGCGTCTATTTAATTAACTCGCCAAAAGTACCCAGCGTGGTATgg                                                                                   >  2:7717610/1‑150 (MQ=255)
                                                                    aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATgg                                                                                   <  2:4818518/150‑1 (MQ=255)
                                                                    aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATgg                                                                                   <  1:2704710/150‑1 (MQ=255)
                                                                    aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATgg                                                                                   <  2:2304231/150‑1 (MQ=255)
                                                                    aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATgg                                                                                   <  1:4392256/150‑1 (MQ=255)
                                                                       gcggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGAtt                                                                                <  2:2756278/150‑1 (MQ=255)
                                                                         ggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAATAGGTCCTGGCGTCTATGACATTCACTCGCCACACGTACCGAGCGTGGGATGGATTGa                                                                              >  1:5182484/1‑150 (MQ=255)
                                                                         ggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATAGGTCCTGGCGTCTATGTCATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGa                                                                              >  2:6730514/1‑150 (MQ=255)
                                                                         ggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGATGGTTTTTTTATGCAAATGAAAAAGGGCCAGGCGTCTATGACAATCACACGCCCAACGTACCCAGCGTGGGATGGGTTGa                                                                              >  1:6252170/1‑150 (MQ=255)
                                                                         ggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGCCGTTGCGTATCCCAATGCAATCCGGCCGGGCGGCGCTGACCGTCACACGCCCAACGGACCGCGCGTGGAAGGGaggga                                                                              >  2:7053567/1‑146 (MQ=255)
                                                                          gcgcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGaa                                                                             <  1:4247553/150‑1 (MQ=255)
                                                                          gcgcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGaa                                                                             <  2:5864023/150‑1 (MQ=255)
                                                                              tGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATACAAAAGAAATAGGGGCTTGCGTCTAATACAATAACTCGCCAAAAGGACCGCGCGTGGAAAGGATTGAAGCCt                                                                         >  1:3982589/1‑150 (MQ=255)
                                                                              tGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTTTTTTCCAAAATAAAAAGGTCCAGGCGTCTATGACAATCACACGCCAAACCTACCCAGCGTGGAAAGGATTGTAGcca                                                                         >  1:1589738/1‑149 (MQ=255)
                                                                                  cGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTTCCGAGCGTGGAATGGATTGAAGCCTTGCt                                                                     >  2:5416306/1‑150 (MQ=255)
                                                                                   gCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTGGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTg                                                                    <  2:6683801/150‑1 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  <  1:2444677/150‑1 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  <  1:4471355/150‑1 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  <  2:6006384/150‑1 (MQ=255)
                                                                                                 cATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCaaa                                                      >  2:550959/1‑150 (MQ=255)
                                                                                                  aTCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTGTTTATCCAAATGAAATCGGTCCGGGCGGCGATGACAATCACTCGCCAAAAGTAACGAGCGGGGTATGGTGTGAAGCCTTGCGGCAGGTAGCggggcaa                                                     >  2:7619942/1‑145 (MQ=255)
                                                                                                         ccTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTc                                              <  2:7220040/150‑1 (MQ=255)
                                                                                                          cTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTcc                                             <  1:7195886/150‑1 (MQ=255)
                                                                                                          cTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTcc                                             <  2:5029076/150‑1 (MQ=255)
                                                                                                           tCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCGTGCTGAAGAAAGCGGCAAAACGCATTCCg                                            <  2:1059446/150‑1 (MQ=255)
                                                                                                           tCGCGTTCCGACATGGAGTTGCAGGACTCGATTGAAGAGTTTGAATATCCAAAACAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCg                                            <  2:5572577/150‑1 (MQ=255)
                                                                                                               gTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATACAAAAGAAAAAGGTCCAGGCGGCGATGACAATCACACGCCAAACGTACCGAGCGTGGAAAGGATTGAAGCCGTGCTGTAGAAAGCGGCAAAACGCATTCCGGCag                                        >  2:3596091/1‑150 (MQ=255)
                                                                                                                  ccGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAATAGGTCCTGGCGGCGATGACATTCACTCGCCAAACGTACCGGGCGTGGAATGGCTTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAgc                                     >  2:7962472/1‑150 (MQ=255)
                                                                                                                      cATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCt                                 <  1:1017550/150‑1 (MQ=255)
                                                                                                                       aTGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATAGGTCCTGGCGTCTATGACATTCAAACGCCAAACGTAACGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCtg                                >  2:5714091/1‑150 (MQ=255)
                                                                                                                         ggAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCtgtg                              <  1:2155231/150‑1 (MQ=255)
                                                                                                                          gAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCTCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCTCCTGTgg                             <  2:2570736/150‑1 (MQ=255)
                                                                                                                                 tGGAGTCGTTTGAATAGTTTTTTTATACAAAAGAAAACGGTCCTGGCGTCTATGACATTAACTCGCCAAAAGTACCGAGCGTGGAATGGATTTAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGATCGCCTGTGGGTCAAcc                      >  2:2033737/1‑150 (MQ=255)
                                                                                                                                     gTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa                  <  2:4459905/150‑1 (MQ=255)
                                                                                                                                     gTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa                  <  1:2937024/150‑1 (MQ=255)
                                                                                                                                     gTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa                  <  1:6358954/150‑1 (MQ=255)
                                                                                                                                        gTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtg               >  1:8252080/1‑150 (MQ=255)
                                                                                                                                        gTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtg               >  2:5297880/1‑150 (MQ=255)
                                                                                                                                         tttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTgg                             >  2:6559668/1‑135 (MQ=255)
                                                                                                                                         tttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTgg                             <  1:6559668/135‑1 (MQ=255)
                                                                                                                                          ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg             <  1:6027418/150‑1 (MQ=255)
                                                                                                                                          ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg             <  2:3837354/150‑1 (MQ=255)
                                                                                                                                            gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc           <  2:2873833/150‑1 (MQ=255)
                                                                                                                                             aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          <  1:4682130/150‑1 (MQ=255)
                                                                                                                                             aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          >  2:86537/1‑150 (MQ=255)
                                                                                                                                             aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          >  2:2633843/1‑150 (MQ=255)
                                                                                                                                                 gTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACAGTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaa      <  1:2500543/150‑1 (MQ=255)
                                                                                                                                                  tttGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGTGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  1:5841083/1‑150 (MQ=255)
                                                                                                                                                  tttGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  1:1095533/1‑150 (MQ=255)
                                                                                                                                                  tttGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  1:3647010/1‑150 (MQ=255)
                                                                                                                                                    tGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAAc   <  2:1427483/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACG  >  NC_000913/4014944‑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: