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.

TGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACG  >  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)
          tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCaaa                                                                                                                                              <  1:938802/150‑1 (MQ=255)
              acacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATACCAAAtaa                                                                                                                                           >  1:909897/1‑148 (MQ=255)
               cacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGaaa                                                                                                                                         <  2:1111431/150‑1 (MQ=255)
                acaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAAt                                                                                                                                        <  2:7306690/150‑1 (MQ=255)
                acaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAAt                                                                                                                                        >  1:2124564/1‑150 (MQ=255)
                   aaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                                                                                                        <  1:1028947/132‑1 (MQ=255)
                       cacaCTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCt                                                                                                                                 <  1:1799123/150‑1 (MQ=255)
                        acacTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTAAGAGTTTTATTAT‑ACAAATTAAATAGGTCCTg                                                                                                                                >  1:5545270/1‑150 (MQ=255)
                        acacTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTTATTAT‑ACAAAATAAAACGGTCCTg                                                                                                                                >  2:850918/1‑150 (MQ=255)
                            tCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGt                                                                                                                            <  1:2343743/150‑1 (MQ=255)
                              acaTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTAT‑ACAAATGAAATAGGTCCTGGCGTCt                                                                                                                          >  1:7342943/1‑150 (MQ=255)
                               caTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCTGACGTCATCACCATCGAAACCTCGCGTTCCGTCTTGGAGTTGCTGGAGTCGTTTGTTGTGTTTGTTTTT‑TGAATTGTAATAGGTGGTGGCGtgtt                                                                                                                         >  1:5142243/1‑147 (MQ=255)
                               caTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTa                                                                                                                         <  1:8348530/150‑1 (MQ=255)
                                 tgtgTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATg                                                                                                                       <  1:78749/150‑1 (MQ=255)
                                  gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑a                                                                                                                     <  2:7967107/150‑1 (MQ=255)
                                    gtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACa                                                                                                                   <  2:7227278/150‑1 (MQ=255)
                                    gtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACa                                                                                                                   <  1:1433605/150‑1 (MQ=255)
                                     ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGGCTATG‑ACAt                                                                                                                  >  2:8390109/1‑150 (MQ=255)
                                     ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑ACAAATGAAAACGGTCCAGGCGTCTATG‑ACAt                                                                                                                  >  1:8575535/1‑150 (MQ=255)
                                       aTTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTc                                                                                                                <  2:2172792/150‑1 (MQ=255)
                                        ttGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑ACAAATGAAATCGGTCCTGGCGTCTATG‑AAAttaa                                                                                                               >  1:8892592/1‑148 (MQ=255)
                                          gCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTGATTAT‑ACAAATGAAATCGGTCCTGGCGGCTATG‑ACATTCACt                                                                                                             >  1:2863915/1‑150 (MQ=255)
                                           cGAGTTCAACGACATCACGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTc                                                                                                            <  1:6035037/150‑1 (MQ=255)
                                             aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGc                                                                                                          <  2:6075473/150‑1 (MQ=255)
                                             aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGc                                                                                                          <  1:4632613/150‑1 (MQ=255)
                                               ttCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTAT‑ACAAAAGAAATCCGTCCTGGCGTCTATG‑AAAATCACACGCCa                                                                                                        >  2:348763/1‑150 (MQ=255)
                                                  aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTTAGTAT‑CCAAATTAAATAGGTCCTGGCGTCTATG‑AAAATCACACGCCaaaa                                                                                                     >  2:4128952/1‑149 (MQ=255)
                                                  aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAAc                                                                                                     >  2:5816570/1‑150 (MQ=255)
                                                    cgatctcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGt                                                                                                   <  1:7481730/145‑1 (MQ=255)
                                                     gACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTGATTAT‑CCAAATGAAATAGGTCCTTGCGTCTATG‑ACATTCACACGCCAAACGTa                                                                                                  >  2:2248562/1‑150 (MQ=255)
                                                         tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑ACAAATTAAATAGGTCCTGGCGTCTATG‑AAAATCACACGCCAAACCTaacca                                                                                              >  1:5913899/1‑146 (MQ=255)
                                                           atGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCAGGCGTCTATG‑ACATTCACTCGCCAAACCTAccccgc                                                                                            >  2:8188820/1‑146 (MQ=255)
                                                              gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAATAGTTTTTTTAT‑GCAAAATAAAAAGGTACATGCGTAAATG‑AAATTAACACGCCAAAAGTACCGAGCCTg                                                                                         >  1:8188943/1‑150 (MQ=255)
                                                              gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTATTAT‑CCAAAAGAAAACGGTCCAGGCGTCTATG‑ACATTAACTCGCCAAACGTACCCAGCggg                                                                                         >  1:4801156/1‑148 (MQ=255)
                                                               aTTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTgg                                                                                        <  2:2343347/150‑1 (MQ=255)
                                                                ttCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTAT‑ACAAATGAAAAAGGTCCTGGCGTCTATG‑ACATTCACACGCCAAACGGAACGAGCGTGGa                                                                                       >  1:6467082/1‑150 (MQ=255)
                                                                ttCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTAT‑CCAAATGAAAACGGTCCAGGCGTCTATG‑ACAATAAAAACCCAAACGTACCCAGCGTGGa                                                                                       >  1:8472106/1‑150 (MQ=255)
                                                                 tCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGaa                                                                                      <  2:5692544/150‑1 (MQ=255)
                                                                   gATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATg                                                                                    <  2:447789/150‑1 (MQ=255)
                                                                    agtGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATgg                                                                                   <  1:122909/148‑1 (MQ=255)
                                                                    aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAAAGAAATAGGTCCTGGCGTCGATG‑ACATTCACTCGCCCAACGTACCGAGCGTGGAATgg                                                                                   >  2:7446616/1‑150 (MQ=255)
                                                                      tGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCAGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGAt                                                                                 <  2:3702217/150‑1 (MQ=255)
                                                                       gcggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGAtt                                                                                <  2:1728457/150‑1 (MQ=255)
                                                                        cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTTAGAGTTTTTTTAT‑CAAAAATAAAAAGGTCCTGGCGTCGATT‑ACATTAAATCGCCAAACCTAACGAGCGTGGTATGGGttt                                                                               >  1:1558517/1‑149 (MQ=255)
                                                                        cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTg                                                                               <  2:3632013/150‑1 (MQ=255)
                                                                        cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTg                                                                               <  2:6346525/150‑1 (MQ=255)
                                                                         ggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGa                                                                              >  2:4732020/1‑150 (MQ=255)
                                                                          gcgcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGaa                                                                             <  1:299699/150‑1 (MQ=255)
                                                                              tGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTTTTTAT‑AAAAAAGAAAAAGGTCCGGGCGTCGAAT‑AAAATCACACCCCAAAAATAACGAGCCTGGGAAGGATTGAAGcca                                                                         >  1:3957267/1‑149 (MQ=255)
                                                                               ggACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCtt                                                                        <  1:4324945/150‑1 (MQ=255)
                                                                               ggACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCtt                                                                        >  2:6173608/1‑150 (MQ=255)
                                                                                 aCGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTAT‑CCAAAATAAATCGGGCCTGGCGTCTATG‑ACATTCACTCGCCCAACGTAACCAGCGTGGCAAGGATTGAAGCCTTGc                                                                      >  2:3346262/1‑150 (MQ=255)
                                                                                 aCGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGc                                                                      <  1:6768645/150‑1 (MQ=255)
                                                                                  cGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTTTTTAT‑ACAAATTAAAACGGTCCTGGCGGCTATT‑ACAATAACTAGCCAAAAGTACCCAGCGTGGTAAGGATTTAAGCCGTGCt                                                                     >  1:6938583/1‑150 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTAT‑CCAAATGAAATAGGTCCTGGCGTCTATT‑ACATTAAATAGCCAAACGGACCGAGCGTGGAATGGATTGAATCCTTTCTgta                                                                  >  2:4935761/1‑148 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  <  1:3827069/150‑1 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGGCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTgaa                                                                  >  2:5292248/1‑150 (MQ=255)
                                                                                     aGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGAGTTTTTTTTT‑ACCAAAGTAAACGGTCCTGGCGGCTATG‑GCATTCACTCGCCAAACGTACCGCGCGGGGGATGGGTTGTAGCCTTGCTgga                                                                  >  1:1325843/1‑148 (MQ=255)
                                                                                            atcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTAT‑ACAAAAGAAAACGGTCCAGGCGGCTATG‑GCAATAAAACGCCAAACGGAACCAGCGGGGGATGGGTTGTAGCCGTGCTGAAGCAAGCg                                                           >  1:7787476/1‑150 (MQ=255)
                                                                                            atcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCg                                                           <  2:4624594/150‑1 (MQ=255)
                                                                                            atcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCg                                                           <  2:4294736/150‑1 (MQ=255)
                                                                                              caCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTTTAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGc                                                         >  2:9304987/1‑150 (MQ=255)
                                                                                                ccATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTTTTAT‑ACAAAAGAAAAAGGTCCAGGCGTGTATGTACA‑TAACACGCCAAACCTACCGAGCCTGGGATGGATTTTAGCCATGCTGAAGAAATCGGCaa                                                       >  1:6258363/1‑150 (MQ=255)
                                                                                                  aTCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCAGGCGTCTATG‑ACATTCAATCGCCAAACGTACCGAGCGTGGAAAGGATTGAAGCCATGCTGAAGAAAGCGGCaaaa                                                     >  2:2935452/1‑150 (MQ=255)
                                                                                                   tCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAAc                                                    <  2:8990243/150‑1 (MQ=255)
                                                                                                     gAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGc                                                  <  1:6098756/150‑1 (MQ=255)
                                                                                                      aaaCCTCGCGTTCCGACATGGAGGTGCTGGAGTCGTTTGAAGAGTTTTATTAT‑CCAAATGACATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCa                                                 <  1:4156196/150‑1 (MQ=255)
                                                                                                           tCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTAT‑CCAAATGAAATAGGTCCTGGCGTCTATG‑ACATTAACACGCCAAACCTACCGAGCGTGGAATGGATTGAAGCCGTGCTGAATAAAGCGGCCAAACGCAttacc                                            >  2:3560357/1‑147 (MQ=255)
                                                                                                           tCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCg                                            <  1:185715/150‑1 (MQ=255)
                                                                                                           tCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCg                                            <  1:5307791/150‑1 (MQ=255)
                                                                                                            cgcgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCgg                                           <  1:3291432/150‑1 (MQ=255)
                                                                                                              cgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCa                                         <  2:4414353/150‑1 (MQ=255)
                                                                                                                ttCCGACATGGAGTTGCTGGAGTCGTTTTTAGAGTTTGATTAT‑ACAAATGAAAAAGGTCCTGGCGTCTATG‑ACATTAACACGCCAAAAGTACCGAGCGTGGAATTGATTGAAGCCTTGCTGAAGAAATCGGCAAAACGCAATACGGCaga                                       >  2:5533637/1‑150 (MQ=255)
                                                                                                                ttCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCCCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAACAAAGCGGCAAAACGCATTCCGGCaga                                       <  1:2589252/150‑1 (MQ=255)
                                                                                                                 tCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCagag                                      <  2:2504572/150‑1 (MQ=255)
                                                                                                                  ccGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATTGATTGAAGCCTTGCTGAAGAAAGCGGCCAAACGCATTCCGGCAGAgc                                     >  2:8288589/1‑150 (MQ=255)
                                                                                                                  ccGACATGGAGTTGCTGGAGTCGTGTGCAGAGTTTGTTTAT‑CCAAAAGGAATCGGTCCTGGCGTCTATG‑ACATTCACTAGCCCAACGAACCGAGCGTGGAATGGATTGTAGCCGTGCGGAAGAAAGCGGCAAAACGCCTTCCGGCAGAgc                                     >  2:8834022/1‑150 (MQ=255)
                                                                                                                   cGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCAGGCGTCTATG‑ACATTCACTCGCCAAACCTAACGAGCGTGGAATGGATTTAAGCCTTGCTGAACAAAGCGGCAAAACGCATTCCCGCAGAgcg                                    >  2:8092825/1‑150 (MQ=255)
                                                                                                                    gACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAgcgc                                   <  2:6675237/150‑1 (MQ=255)
                                                                                                                          gAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTgg                             <  2:2075276/150‑1 (MQ=255)
                                                                                                                          gAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTgg                             <  2:6480033/150‑1 (MQ=255)
                                                                                                                                     gTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa                  <  2:6256682/150‑1 (MQ=255)
                                                                                                                                     gTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa                  <  1:3529625/150‑1 (MQ=255)
                                                                                                                                         tttGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGCTTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgt              >  2:7063452/1‑150 (MQ=255)
                                                                                                                                         tttGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgt              >  2:155685/1‑150 (MQ=255)
                                                                                                                                           tGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTgg            <  1:3643530/150‑1 (MQ=255)
                                                                                                                                            gAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc           <  2:1308569/150‑1 (MQ=255)
                                                                                                                                             aaGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          <  2:1807789/150‑1 (MQ=255)
                                                                                                                                             aaGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc          >  1:7157526/1‑150 (MQ=255)
                                                                                                                                              agagTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCt         <  2:3908056/150‑1 (MQ=255)
                                                                                                                                              agagTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCt         <  2:7652310/150‑1 (MQ=255)
                                                                                                                                                  tttGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  1:5678222/1‑150 (MQ=255)
                                                                                                                                                  tttGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  2:860083/1‑150 (MQ=255)
                                                                                                                                                  tttGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     <  1:4501697/150‑1 (MQ=255)
                                                                                                                                                  tttGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaa     >  2:6573002/1‑150 (MQ=255)
                                                                                                                                                   ttGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGaaaa    <  2:3807079/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTAT‑CCAAATGAAATCGGTCCTGGCGTCTATG‑ACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACG  >  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: ‑