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% 1793.7 / 34.5 595D681Y (GAT→TAT) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (188/320);  new base T (82/3);  total (270/325)
Fisher's exact test for biased strand distribution p-value = 9.10e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 601 total aligned reads displayed.

TGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACG  >  NC_000913/4014944‑4015242
                                                                                                                                                     |                                                                                                                                                       
 ggCGAAAGATGACACACAACTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                                        <  2:6903689/150‑1 (MQ=255)
 ggCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                                        <  2:2696262/150‑1 (MQ=255)
 ggCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                                        <  2:1966032/150‑1 (MQ=255)
  gCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTAt                                                                                                                                                       >  1:7586093/1‑150 (MQ=255)
          tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCaaa                                                                                                                                              <  1:938802/150‑1 (MQ=255)
              acacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATACCAAAtaa                                                                                                                                           >  1:909897/1‑148 (MQ=255)
               cacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGaaa                                                                                                                                         <  2:1111431/150‑1 (MQ=255)
                acaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAAt                                                                                                                                        <  2:7306690/150‑1 (MQ=255)
                acaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAAt                                                                                                                                        >  1:2124564/1‑150 (MQ=255)
                   aaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                                        <  1:1028947/132‑1 (MQ=255)
                       cacaCTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCt                                                                                                                                 <  1:1799123/150‑1 (MQ=255)
                        acacTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTAAGAGTTTTATTATACAAATTAAATAGGTCCTg                                                                                                                                >  1:5545270/1‑150 (MQ=255)
                        acacTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTTATTATACAAAATAAAACGGTCCTg                                                                                                                                >  2:850918/1‑150 (MQ=255)
                            tCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGt                                                                                                                            <  1:2343743/150‑1 (MQ=255)
                              acaTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTATACAAATGAAATAGGTCCTGGCGTCt                                                                                                                          >  1:7342943/1‑150 (MQ=255)
                               caTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCTGACGTCATCACCATCGAAACCTCGCGTTCCGTCTTGGAGTTGCTGGAGTCGTTTGTTGTGTTTGTTTTTTGAATTGTAATAGGTGGTGGCGtgtt                                                                                                                         >  1:5142243/1‑147 (MQ=255)
                               caTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTa                                                                                                                         <  1:8348530/150‑1 (MQ=255)
                                 tgtgTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATg                                                                                                                       <  1:78749/150‑1 (MQ=255)
                                  gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                                                                                                     <  2:7967107/150‑1 (MQ=255)
                                    gtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACa                                                                                                                   <  2:7227278/150‑1 (MQ=255)
                                    gtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACa                                                                                                                   <  1:1433605/150‑1 (MQ=255)
                                     ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGGCTATGACAt                                                                                                                  >  2:8390109/1‑150 (MQ=255)
                                     ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAAACGGTCCAGGCGTCTATGACAt                                                                                                                  >  1:8575535/1‑150 (MQ=255)
                                       aTTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTc                                                                                                                <  2:2172792/150‑1 (MQ=255)
                                        ttGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAATCGGTCCTGGCGTCTATGAAAttaa                                                                                                               >  1:8892592/1‑148 (MQ=255)
                                          gCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTGATTATACAAATGAAATCGGTCCTGGCGGCTATGACATTCACt                                                                                                             >  1:2863915/1‑150 (MQ=255)
                                           cGAGTTCAACGACATCACGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTc                                                                                                            <  1:6035037/150‑1 (MQ=255)
                                             aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                          <  2:6075473/150‑1 (MQ=255)
                                             aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                          <  1:4632613/150‑1 (MQ=255)
                                               ttCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTATACAAAAGAAATCCGTCCTGGCGTCTATGAAAATCACACGCCa                                                                                                        >  2:348763/1‑150 (MQ=255)
                                                  aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTTAGTATCCAAATTAAATAGGTCCTGGCGTCTATGAAAATCACACGCCaaaa                                                                                                     >  2:4128952/1‑149 (MQ=255)
                                                  aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAAc                                                                                                     >  2:5816570/1‑150 (MQ=255)
                                                    cgatctcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGt                                                                                                   <  1:7481730/145‑1 (MQ=255)
                                                     gACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTGATTATCCAAATGAAATAGGTCCTTGCGTCTATGACATTCACACGCCAAACGTa                                                                                                  >  2:2248562/1‑150 (MQ=255)
                                                         tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATTAAATAGGTCCTGGCGTCTATGAAAATCACACGCCAAACCTaacca                                                                                              >  1:5913899/1‑146 (MQ=255)
                                                           atGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAATGAAATCGGTCCAGGCGTCTATGACATTCACTCGCCAAACCTAccccgc                                                                                            >  2:8188820/1‑146 (MQ=255)
                                                              gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAATAGTTTTTTTATGCAAAATAAAAAGGTACATGCGTAAATGAAATTAACACGCCAAAAGTACCGAGCCTg                                                                                         >  1:8188943/1‑150 (MQ=255)
                                                              gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTATCCAAAAGAAAACGGTCCAGGCGTCTATGACATTAACTCGCCAAACGTACCCAGCggg                                                                                         >  1:4801156/1‑148 (MQ=255)
                                                               aTTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTgg                                                                                        <  2:2343347/150‑1 (MQ=255)
                                                                ttCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATACAAATGAAAAAGGTCCTGGCGTCTATGACATTCACACGCCAAACGGAACGAGCGTGGa                                                                                       >  1:6467082/1‑150 (MQ=255)
                                                                ttCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAATGAAAACGGTCCAGGCGTCTATGACAATAAAAACCCAAACGTACCCAGCGTGGa                                                                                       >  1:8472106/1‑150 (MQ=255)
                                                                 tCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGaa                                                                                      <  2:5692544/150‑1 (MQ=255)
                                                                   gATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATg                                                                                    <  2:447789/150‑1 (MQ=255)
                                                                    agtGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATgg                                                                                   <  1:122909/148‑1 (MQ=255)
                                                                    aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAAGAAATAGGTCCTGGCGTCGATGACATTCACTCGCCCAACGTACCGAGCGTGGAATgg                                                                                   >  2:7446616/1‑150 (MQ=255)
                                                                      tGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCAGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGAt                                                                                 <  2:3702217/150‑1 (MQ=255)
                                                                       gcggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGAtt                                                                                <  2:1728457/150‑1 (MQ=255)
                                                                        cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTTAGAGTTTTTTTATCAAAAATAAAAAGGTCCTGGCGTCGATTACATTAAATCGCCAAACCTAACGAGCGTGGTATGGGttt                                                                               >  1:1558517/1‑149 (MQ=255)
                                                                        cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTg                                                                               <  2:3632013/150‑1 (MQ=255)
                                                                        cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTg                                                                               <  2:6346525/150‑1 (MQ=255)
                                                                         ggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGa                                                                              >  2:4732020/1‑150 (MQ=255)
                                                                          gcgcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGaa                                                                             <  1:299699/150‑1 (MQ=255)
                                                                              tGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTTTTTATAAAAAAGAAAAAGGTCCGGGCGTCGAATAAAATCACACCCCAAAAATAACGAGCCTGGGAAGGATTGAAGcca                                                                         >  1:3957267/1‑149 (MQ=255)
                                                                               ggACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCtt                                                                        <  1:4324945/150‑1 (MQ=255)
                                                                               ggACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCtt                                                                        >  2:6173608/1‑150 (MQ=255)
                                                                                 aCGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTATCCAAAATAAATCGGGCCTGGCGTCTATGACATTCACTCGCCCAACGTAACCAGCGTGGCAAGGATTGAAGCCTTGc                                                                      >  2:3346262/1‑150 (MQ=255)
                                                                                 aCGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGc                                                                      <  1:6768645/150‑1 (MQ=255)
                                                                                  cGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTTTTTATACAAATTAAAACGGTCCTGGCGGCTATTACAATAACTAGCCAAAAGTACCCAGCGTGGTAAGGATTTAAGCCGTGCt                                                                     >  1:6938583/1‑150 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAATAGGTCCTGGCGTCTATTACATTAAATAGCCAAACGGACCGAGCGTGGAATGGATTGAATCCTTTCTgta                                                                  >  2:4935761/1‑148 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  <  1:3827069/150‑1 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGGCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  >  2:5292248/1‑150 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTTTTTTTACCAAAGTAAACGGTCCTGGCGGCTATGGCATTCACTCGCCAAACGTACCGCGCGGGGGATGGGTTGTAGCCTTGCTgga                                                                  >  1:1325843/1‑148 (MQ=255)
                                                                                            atcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTATACAAAAGAAAACGGTCCAGGCGGCTATGGCAATAAAACGCCAAACGGAACCAGCGGGGGATGGGTTGTAGCCGTGCTGAAGCAAGCg                                                           >  1:7787476/1‑150 (MQ=255)
                                                                                            atcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCg                                                           <  2:4624594/150‑1 (MQ=255)
                                                                                            atcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCg                                                           <  2:4294736/150‑1 (MQ=255)
                                                                                              caCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTTAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGc                                                         >  2:9304987/1‑150 (MQ=255)
                                                                                                ccATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTTTTATACAAAAGAAAAAGGTCCAGGCGTGTATGTACATAACACGCCAAACCTACCGAGCCTGGGATGGATTTTAGCCATGCTGAAGAAATCGGCaa                                                       >  1:6258363/1‑150 (MQ=255)
                                                                                                  aTCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAATGAAATCGGTCCAGGCGTCTATGACATTCAATCGCCAAACGTACCGAGCGTGGAAAGGATTGAAGCCATGCTGAAGAAAGCGGCaaaa                                                     >  2:2935452/1‑150 (MQ=255)
                                                                                                   tCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAAc                                                    <  2:8990243/150‑1 (MQ=255)
                                                                                                     gAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGc                                                  <  1:6098756/150‑1 (MQ=255)
                                                                                                      aaaCCTCGCGTTCCGACATGGAGGTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGACATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCa                                                 <  1:4156196/150‑1 (MQ=255)
                                                                                                           tCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTATCCAAATGAAATAGGTCCTGGCGTCTATGACATTAACACGCCAAACCTACCGAGCGTGGAATGGATTGAAGCCGTGCTGAATAAAGCGGCCAAACGCAttacc                                            >  2:3560357/1‑147 (MQ=255)
                                                                                                           tCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCg                                            <  1:185715/150‑1 (MQ=255)
                                                                                                           tCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCg                                            <  1:5307791/150‑1 (MQ=255)
                                                                                                            cgcgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCgg                                           <  1:3291432/150‑1 (MQ=255)
                                                                                                              cgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCa                                         <  2:4414353/150‑1 (MQ=255)
                                                                                                                ttCCGACATGGAGTTGCTGGAGTCGTTTTTAGAGTTTGATTATACAAATGAAAAAGGTCCTGGCGTCTATGACATTAACACGCCAAAAGTACCGAGCGTGGAATTGATTGAAGCCTTGCTGAAGAAATCGGCAAAACGCAATACGGCaga                                       >  2:5533637/1‑150 (MQ=255)
                                                                                                                ttCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCCCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAACAAAGCGGCAAAACGCATTCCGGCaga                                       <  1:2589252/150‑1 (MQ=255)
                                                                                                                 tCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCagag                                      <  2:2504572/150‑1 (MQ=255)
                                                                                                                  ccGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATTGATTGAAGCCTTGCTGAAGAAAGCGGCCAAACGCATTCCGGCAGAgc                                     >  2:8288589/1‑150 (MQ=255)
                                                                                                                  ccGACATGGAGTTGCTGGAGTCGTGTGCAGAGTTTGTTTATCCAAAAGGAATCGGTCCTGGCGTCTATGACATTCACTAGCCCAACGAACCGAGCGTGGAATGGATTGTAGCCGTGCGGAAGAAAGCGGCAAAACGCCTTCCGGCAGAgc                                     >  2:8834022/1‑150 (MQ=255)
                                                                                                                   cGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCAGGCGTCTATGACATTCACTCGCCAAACCTAACGAGCGTGGAATGGATTTAAGCCTTGCTGAACAAAGCGGCAAAACGCATTCCCGCAGAgcg                                    >  2:8092825/1‑150 (MQ=255)
                                                                                                                    gACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAgcgc                                   <  2:6675237/150‑1 (MQ=255)
                                                                                                                          gAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTgg                             <  2:2075276/150‑1 (MQ=255)
                                                                                                                          gAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTgg                             <  2:6480033/150‑1 (MQ=255)
                                                                                                                                     gTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa                  <  2:6256682/150‑1 (MQ=255)
                                                                                                                                     gTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa                  <  1:3529625/150‑1 (MQ=255)
                                                                                                                                         tttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGCTTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgt              >  2:7063452/1‑150 (MQ=255)
                                                                                                                                         tttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgt              >  2:155685/1‑150 (MQ=255)
                                                                                                                                           tGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTgg            <  1:3643530/150‑1 (MQ=255)
                                                                                                                                            gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc           <  2:1308569/150‑1 (MQ=255)
                                                                                                                                             aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          <  2:1807789/150‑1 (MQ=255)
                                                                                                                                             aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          >  1:7157526/1‑150 (MQ=255)
                                                                                                                                              agagTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCt         <  2:3908056/150‑1 (MQ=255)
                                                                                                                                              agagTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCt         <  2:7652310/150‑1 (MQ=255)
                                                                                                                                                  tttGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  1:5678222/1‑150 (MQ=255)
                                                                                                                                                  tttGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  2:860083/1‑150 (MQ=255)
                                                                                                                                                  tttGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     <  1:4501697/150‑1 (MQ=255)
                                                                                                                                                  tttGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  2:6573002/1‑150 (MQ=255)
                                                                                                                                                   ttGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaaa    <  2:3807079/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: