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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,0850AT 7.8% 1101.0 / 12.0 365E678V (GAA→GTA) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (119/204);  new base T (25/2);  total (159/206)
Fisher's exact test for biased strand distribution p-value = 8.82e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 379 total aligned reads displayed.

GCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCC  >  NC_000913/4014937‑4015234
                                                                                                                                                    |                                                                                                                                                     
 ccgccgTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAag                                                                                                                                                     <  2:1148771/150‑1 (MQ=255)
  cgccgTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAaga                                                                                                                                                    <  2:5188265/150‑1 (MQ=255)
                aTGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGATTGAAGAGTTTGATTATCCaa                                                                                                                                      <  2:5151803/150‑1 (MQ=255)
                 tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCaaa                                                                                                                                     >  1:4654542/1‑150 (MQ=255)
                 tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCaaa                                                                                                                                     <  1:1918632/150‑1 (MQ=255)
                 tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCaaa                                                                                                                                     <  2:4087779/150‑1 (MQ=255)
                 tGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGAGTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTGTGATTATCCaac                                                                                                                                     >  1:5805756/1‑149 (MQ=255)
                  gACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAAt                                                                                                                                    >  2:4620794/1‑150 (MQ=255)
                     acacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGaa                                                                                                                                 <  2:3166110/150‑1 (MQ=255)
                     acacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAAGaa                                                                                                                                 >  2:3835853/1‑150 (MQ=255)
                     acacaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGATCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGaa                                                                                                                                 >  1:3445967/1‑150 (MQ=255)
                       acaAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGGGGTTGATTATCCCAATGaaag                                                                                                                               >  1:3410946/1‑149 (MQ=255)
                        caAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATc                                                                                                                              <  1:3672213/150‑1 (MQ=255)
                         aaaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCg                                                                                                                             <  2:3303662/150‑1 (MQ=255)
                         aaaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCg                                                                                                                             <  1:3009054/150‑1 (MQ=255)
                         aaaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCg                                                                                                                             >  2:2417592/1‑150 (MQ=255)
                           aTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGt                                                                                                                           <  1:5677650/150‑1 (MQ=255)
                            tCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTc                                                                                                                          <  2:1402015/150‑1 (MQ=255)
                             ccACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTTGAGTTTTTTTTTCCAAATTAAATAGGtag                                                                                                                         >  2:2060706/1‑148 (MQ=255)
                             ccACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCCAATGAAATCGGTcc                                                                                                                         >  2:3169152/1‑150 (MQ=255)
                              cacaCTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGTTTATAAAAAAGAAATccgacca                                                                                                                        >  2:3003795/1‑144 (MQ=255)
                                    cacaTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTc                                                                                                                  <  1:647933/150‑1 (MQ=255)
                                     acaTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAATACGTCCTGGCGTcg                                                                                                                 >  2:741592/1‑149 (MQ=255)
                                     acaTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCt                                                                                                                 <  2:3587164/150‑1 (MQ=255)
                                     acaTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCt                                                                                                                 <  2:3724683/150‑1 (MQ=255)
                                      caTGTGTTAGTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTa                                                                                                                <  1:135804/150‑1 (MQ=255)
                                        tgtgTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAAAAGGTCCTGGCGTCTATg                                                                                                              >  2:1887673/1‑150 (MQ=255)
                                         gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                                                                                             <  2:1026732/150‑1 (MQ=255)
                                         gtgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                                                                                             <  1:4369689/150‑1 (MQ=255)
                                            ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAAAGAAATCCGTCCGGGCGTCGATGACaa                                                                                                          >  2:147349/1‑149 (MQ=255)
                                            ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTGGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACAt                                                                                                          >  1:2259152/1‑150 (MQ=255)
                                            ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTGATTATCCCAAAGTAATCGGTCCTGGCGTCTATGACAt                                                                                                          >  1:1395497/1‑150 (MQ=255)
                                               ttGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCa                                                                                                       <  2:483012/150‑1 (MQ=255)
                                                tGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTATAGGTTGATTATGCAAATTAAATACGTGCTGGCGTCGATGAAAATCaa                                                                                                      >  1:2903708/1‑149 (MQ=255)
                                                  cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTc                                                                                                    <  1:1701696/150‑1 (MQ=255)
                                                  cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTc                                                                                                    <  1:1160410/150‑1 (MQ=255)
                                                  cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTAGATTATCCAAATGAAATAGGTCCTGGCGTCTATGACATTCACTc                                                                                                    >  1:3536699/1‑150 (MQ=255)
                                                  cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGTAGCGTTTGATTATACAAATGAAATAGGTCCTGGCGTCTATTACATTCAcgc                                                                                                    >  1:1377925/1‑148 (MQ=255)
                                                   gAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCg                                                                                                   <  2:2495787/150‑1 (MQ=255)
                                                    aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                  <  2:5423468/150‑1 (MQ=255)
                                                    aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                                                                                                  <  2:5223846/150‑1 (MQ=255)
                                                      ttCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAATCGGTCCTGGCGTCTATGACAATCACACGCCa                                                                                                >  2:558432/1‑150 (MQ=255)
                                                      ttCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCCAAAGAAATCGGTCCGGGCGGCGATGACAATCACGCGCCa                                                                                                >  1:5184471/1‑150 (MQ=255)
                                                         aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTATAGTTTGATTATCCCAATGAAATAGGTCCGGGCGTCGATGACATTCACTCGCCCAAc                                                                                             >  1:2314278/1‑150 (MQ=255)
                                                         aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTTATTATCCAAATGAAATAGGTCCTGGCGGCTATGAAATTAACGCGCCAAAc                                                                                             >  1:274692/1‑150 (MQ=255)
                                                         aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAAc                                                                                             >  2:5300382/1‑150 (MQ=255)
                                                         aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGATTTTGTTTATCCAAAAGAAATAGGTCCTGGCGTCTATGAAAATAACTCGCCAAAc                                                                                             >  2:596257/1‑150 (MQ=255)
                                                          aCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTATCCAAATGAAATCGGTCCAGGCGTCTATGACATTAACTCGCCAAACg                                                                                            >  2:4029079/1‑150 (MQ=255)
                                                            gACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGGAGAGTTTGATTATACAAATGAAATCGGTCCTGGCGTCTATGACATTAACTAGCCAAACGTa                                                                                          >  1:266286/1‑150 (MQ=255)
                                                                 catGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAg                                                                                     <  1:284373/150‑1 (MQ=255)
                                                                    ggATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGGCTATGACATTCACTCGCCAAACGTACCGAGCGt                                                                                  >  2:4389074/1‑150 (MQ=255)
                                                                     gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTg                                                                                 <  1:941371/150‑1 (MQ=255)
                                                                        tCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGaa                                                                              <  2:2802408/150‑1 (MQ=255)
                                                                        tCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGaa                                                                              <  1:3303935/150‑1 (MQ=255)
                                                                        tCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGTTTATCCCAATGCAATCGGTCCTGGCGTCTATTACAGTCACTCGCCCAACGGACCGCGCGTGGac                                                                              >  1:3834361/1‑149 (MQ=255)
                                                                           aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTTTTATACAAAAGAAATCGGTCCTGGCGTCTATGACATTAACTCGCCAAACCTAACCAGCGTGGAATgg                                                                           >  1:1870075/1‑150 (MQ=255)
                                                                              gcggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAAACGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGAtt                                                                        >  2:661656/1‑150 (MQ=255)
                                                                                 gcgcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTTTTTATACAAAATAAAAAGGTCCTGGCGTCGATGAAAATCACACGCCAAAACAACCCAGCGTGGAATGGATTGaa                                                                     >  2:2532336/1‑150 (MQ=255)
                                                                                     tGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTTAAGAGTTTTATTATACAAATTAAATCGGTCCTGGCGTCTATGACATTAAATCGCCAAACGTACCGAGCGTGGAATGGATTGAAGcca                                                                 >  2:4531752/1‑149 (MQ=255)
                                                                                       gACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTTTTTATACAAAATAAATCGGTCCGGGCGGCGATTACATTCAAACGCCAAACCGAACCAGCGTGGGATGGATTGAAGCCtag                                                               >  1:4029494/1‑148 (MQ=255)
                                                                                       gACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACAATCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTg                                                               >  2:3378224/1‑150 (MQ=255)
                                                                                       gACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAATAGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTg                                                               >  2:2332352/1‑150 (MQ=255)
                                                                                         cGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCt                                                             <  1:2417592/150‑1 (MQ=255)
                                                                                         cGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCt                                                             <  1:2593918/150‑1 (MQ=255)
                                                                                         cGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAATAGTTTGATTATCCAAATTAAATAGGTCCTGGCGTCTATGACAATAACACGCCAAACGTAACGAGCGTGGAATTGATTTAATCCTTgat                                                             >  1:4211028/1‑148 (MQ=255)
                                                                                           cAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTTGTGTTTGATTATGCAAATGAAATAGGTGCTGGCGTGTATTTCATTAACTCGCCAAACGTACCGAGCGTTGTATTGATTGTATCCTTGGTga                                                           >  2:321559/1‑150 (MQ=255)
                                                                                                    tcaCCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGTAGAGTTTGATTATCCAAATTAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCgg                                                  <  1:5566735/150‑1 (MQ=255)
                                                                                                          tCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAAc                                            <  1:4339239/150‑1 (MQ=255)
                                                                                                           cGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGGAGAGTTTGTTTTTCCCAATGAAATAGGTCCAGGCGTCTATGACAATCACACGCCAAACCTACCGAGCGTGGAATGGGTTGAAGCCGTGCAGAAGAAAGCGGGAAAACg                                           >  1:3269100/1‑150 (MQ=255)
                                                                                                            gAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGGTTGGTTATCCCCATGAAATCGGTCCGGGCGGCGATGCCCTTCACGCGCCAAACCGACCGAGCGTGGAATGGATTGAAGCCGTGCTGAAGAAAGCGGCACAACGc                                          >  2:863841/1‑150 (MQ=255)
                                                                                                             aaaCCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAATCGGTCCAGGCGTCTATGACATTCACTCGCCAAACGTACCCAGCGTGGAATGGATTGAAGCCTTGCTGAATAAAGCGGCAAAACGCa                                         >  2:5087983/1‑150 (MQ=255)
                                                                                                              aaCCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCAt                                        <  2:4341740/150‑1 (MQ=255)
                                                                                                                ccTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTc                                      <  1:4053681/150‑1 (MQ=255)
                                                                                                                 cTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTcc                                     <  1:4932319/150‑1 (MQ=255)
                                                                                                                    gcgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGc                                  >  2:177039/1‑150 (MQ=255)
                                                                                                                    gcgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATACAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGc                                  >  2:345966/1‑150 (MQ=255)
                                                                                                                    gcgTTCCGACATGGAGTTGCTGGAGTCGTGTGAAGAGTTTGATTATCCAAAGGAAAACGGTCCTGGCGTCGAAGACAGTCACTCCCCAAACCTACCGAGCGTGGAATGGATTGAAGCCATGCAGAACAAAGCGGCAAAACGCATTCCCGc                                  >  2:250307/1‑150 (MQ=255)
                                                                                                                     cgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCa                                 <  2:2892950/150‑1 (MQ=255)
                                                                                                                     cgTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCa                                 <  2:2723032/150‑1 (MQ=255)
                                                                                                                          cGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATGCAAATGAAATAGGTGCTGGCGTCTATGTCATTCACTCGCCAATCGTACCGAGGGTGGAATGGATTGTAGCCTTGGTGAAGAAATAGGCAAAACGCATTCCGGCATAgcg                            >  2:4722207/1‑150 (MQ=255)
                                                                                                                          cGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGGCTATGACATTCAATCGCCAAACGTACCGAGCGTGGAATGGGTTGAAGCCTTGCAGAAGAAAGCGGCAAAAAGCAATACGGCAGAgcg                            >  2:1280753/1‑150 (MQ=255)
                                                                                                                             cATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCt                         <  1:677413/150‑1 (MQ=255)
                                                                                                                              aTGGAGTTGCTGGAGTCGTTTGAAGAGTTTTATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCtg                        >  2:908812/1‑150 (MQ=255)
                                                                                                                                 gAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTgg                     <  1:1658184/150‑1 (MQ=255)
                                                                                                                                   gTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGt                   <  1:1594428/150‑1 (MQ=255)
                                                                                                                                            gTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa          <  1:2165361/150‑1 (MQ=255)
                                                                                                                                            gTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGa          <  1:5382441/150‑1 (MQ=255)
                                                                                                                                               gTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtg       >  2:1094244/1‑150 (MQ=255)
                                                                                                                                                tttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgt      <  2:80542/150‑1 (MQ=255)
                                                                                                                                                tttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgt      >  1:3052756/1‑150 (MQ=255)
                                                                                                                                                 ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg     >  2:2135698/1‑150 (MQ=255)
                                                                                                                                                 ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg     >  1:1074188/1‑150 (MQ=255)
                                                                                                                                                 ttGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACtgtg     >  2:5731528/1‑150 (MQ=255)
                                                                                                                                                   gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc   <  1:2523843/150‑1 (MQ=255)
                                                                                                                                                   gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc   <  1:819127/150‑1 (MQ=255)
                                                                                                                                                   gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc   <  2:1531490/150‑1 (MQ=255)
                                                                                                                                                   gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc   <  2:2259152/150‑1 (MQ=255)
                                                                                                                                                   gAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGc   <  1:101138/150‑1 (MQ=255)
                                                                                                                                                    aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc  >  1:5438872/1‑150 (MQ=255)
                                                                                                                                                    aaGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGcc  <  2:1647516/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGCCGTGGCGAAAGATGACACACAAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGAATGGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCC  >  NC_000913/4014937‑4015234

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