Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,214,563 A→T 100% T270S (ACC→TCC)  metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,214,5630AT99.2% 895.6 / ‑2.4 259T270S (ACC→TCC) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (112/145);  total (112/147)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 259 total aligned reads displayed.

GCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATAT  >  CP009273/4214417‑4214708
                                                                                                                                                  |                                                                                                                                                 
gcACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTcc                                                                                                                                                 >  2:3436907/1‑149 (MQ=255)
gcACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTcc                                                                                                                                                 <  2:3646259/149‑1 (MQ=255)
gcACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTcc                                                                                                                                                 <  2:2062759/149‑1 (MQ=255)
  aCGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCcgc                                                                                                                                               >  1:2924581/1‑149 (MQ=255)
   cGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCcgcg                                                                                                                                              >  2:325212/1‑149 (MQ=255)
   cGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCcgcg                                                                                                                                              >  2:4663729/1‑149 (MQ=255)
    gatctcCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCgcgc                                                                                                                                             >  2:4064401/3‑149 (MQ=255)
        tcCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCAccc                                                                                                                                         <  2:4289899/149‑1 (MQ=255)
             gCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACg                                                                                                                                    >  2:2500981/1‑149 (MQ=255)
              cAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGc                                                                                                                                   >  1:4047116/1‑149 (MQ=255)
               aGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGcc                                                                                                                                  <  2:1099361/149‑1 (MQ=255)
                gACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCg                                                                                                                                 >  1:1856958/1‑149 (MQ=255)
                gACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCg                                                                                                                                 >  1:1011389/1‑149 (MQ=255)
                   caccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGc                                                                                                                              <  1:4220754/149‑1 (MQ=255)
                   caccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGc                                                                                                                              >  1:1562052/1‑149 (MQ=255)
                   caccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGc                                                                                                                              <  1:1238602/149‑1 (MQ=255)
                    accGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCt                                                                                                                             <  1:1135282/149‑1 (MQ=255)
                       gAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTAcc                                                                                                                          >  2:1268415/1‑149 (MQ=255)
                       gAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTATGACCATTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGCAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCCGGCTAcc                                                                                                                          >  2:4554814/1‑149 (MQ=255)
                        aaGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTAccc                                                                                                                         >  1:459266/1‑149 (MQ=255)
                         aGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCa                                                                                                                        >  2:986388/1‑149 (MQ=255)
                          gCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCaa                                                                                                                       >  1:478047/1‑149 (MQ=255)
                                 aCAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCttt                                                                                                                <  1:3021101/149‑1 (MQ=255)
                                 aCAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCttt                                                                                                                <  1:1767755/149‑1 (MQ=255)
                                    aCTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGCACGCCGGGCTACCCAACGCCTTTGGt                                                                                                             <  1:3685971/149‑1 (MQ=255)
                                      tCATTGCGCCATGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCGCGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGa                                                                                                           <  2:2427309/149‑1 (MQ=255)
                                        aTTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGt                                                                                                         >  2:3541797/1‑149 (MQ=255)
                                            cgcCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGa                                                                                                     <  1:202548/149‑1 (MQ=255)
                                             gcCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAt                                                                                                    <  2:3316739/149‑1 (MQ=255)
                                              ccACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtc                                                                                                   <  2:1517130/149‑1 (MQ=255)
                                              ccACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtc                                                                                                   >  1:382669/1‑149 (MQ=255)
                                               cACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtct                                                                                                  <  2:1907101/149‑1 (MQ=255)
                                               cACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtct                                                                                                  <  1:4563006/149‑1 (MQ=255)
                                                aCGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtctc                                                                                                 <  1:3582981/149‑1 (MQ=255)
                                                aCGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtctc                                                                                                 <  1:667739/149‑1 (MQ=255)
                                                  gCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGa                                                                                               <  2:1600923/149‑1 (MQ=255)
                                                  gCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGa                                                                                               <  1:2338373/149‑1 (MQ=255)
                                                     gAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGc                                                                                            <  2:809247/149‑1 (MQ=255)
                                                      aaGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGcc                                                                                           <  1:1201619/149‑1 (MQ=255)
                                                      aaGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGcc                                                                                           <  2:1201831/149‑1 (MQ=255)
                                                         ctctGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGac                                                                                        >  1:3643385/1‑149 (MQ=255)
                                                          tctGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGaca                                                                                       >  1:3524045/1‑149 (MQ=255)
                                                           ctGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGacac                                                                                      <  2:1628789/149‑1 (MQ=255)
                                                             gACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGa                                                                                    >  1:844868/1‑149 (MQ=255)
                                                               ccTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATg                                                                                  <  1:3697897/149‑1 (MQ=255)
                                                                   tGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCaa                                                                              <  1:3541797/149‑1 (MQ=255)
                                                                      ccTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAAc                                                                           >  2:3200375/1‑149 (MQ=255)
                                                                       cTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACa                                                                          <  2:4118046/149‑1 (MQ=255)
                                                                          aaCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGAt                                                                       <  1:679815/149‑1 (MQ=255)
                                                                           aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATa                                                                      >  2:2340256/1‑149 (MQ=255)
                                                                               tgCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTg                                                                  <  1:18240/149‑1 (MQ=255)
                                                                                gcgcTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGGGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGa                                                                 >  1:4272982/1‑149 (MQ=255)
                                                                                 cgcTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGaa                                                                <  2:1041354/149‑1 (MQ=255)
                                                                                  gcTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAAt                                                               >  1:4286891/1‑149 (MQ=255)
                                                                                   cTGGGGCCAGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATg                                                              <  1:2415593/149‑1 (MQ=255)
                                                                                    tGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATgg                                                             <  2:3578680/149‑1 (MQ=255)
                                                                                     ggggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATggg                                                            <  2:2867647/149‑1 (MQ=255)
                                                                                     ggggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCGTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATggg                                                            >  1:972965/1‑149 (MQ=255)
                                                                                       ggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGgcg                                                          >  2:3783705/1‑149 (MQ=255)
                                                                                          ccGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCaa                                                       >  2:4334596/1‑149 (MQ=255)
                                                                                            gATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGc                                                     <  2:3986529/149‑1 (MQ=255)
                                                                                              tGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGACTGGGCGCAAGCgg                                                   >  2:2974121/1‑149 (MQ=255)
                                                                                              tGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCgg                                                   >  1:1695302/1‑149 (MQ=255)
                                                                                                aaCTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGt                                                 <  2:384644/149‑1 (MQ=255)
                                                                                                   tGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGtttt                                              <  1:1057789/149‑1 (MQ=255)
                                                                                                    gcgcCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTtc                                             <  1:223468/149‑1 (MQ=255)
                                                                                                       ccAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCa                                          >  2:309933/1‑149 (MQ=255)
                                                                                                       ccAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCa                                          >  2:2702182/1‑149 (MQ=255)
                                                                                                        cAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCaa                                         >  2:3658786/1‑149 (MQ=255)
                                                                                                         aGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                                        <  2:3724775/149‑1 (MQ=255)
                                                                                                           tACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAatat                                      <  2:1376282/149‑1 (MQ=255)
                                                                                                            aCGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtc                                     <  2:4685976/149‑1 (MQ=255)
                                                                                                              gTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgt                                   >  1:2123263/1‑149 (MQ=255)
                                                                                                              gTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgt                                   <  1:501862/149‑1 (MQ=255)
                                                                                                              gTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgt                                   <  1:282554/149‑1 (MQ=255)
                                                                                                               tGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgtc                                  <  2:873346/149‑1 (MQ=255)
                                                                                                               tGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgtc                                  <  2:4645554/149‑1 (MQ=255)
                                                                                                                gCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgtcg                                 <  1:4453310/149‑1 (MQ=255)
                                                                                                                  aGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTcggc                               <  1:3308507/149‑1 (MQ=255)
                                                                                                                   ggAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTcggcg                              >  1:2649241/1‑149 (MQ=255)
                                                                                                                    gAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTcggcgg                             >  1:3273144/1‑149 (MQ=255)
                                                                                                                     aGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCggcggc                            <  2:2924581/149‑1 (MQ=255)
                                                                                                                     aGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCggcggc                            <  1:2652287/149‑1 (MQ=255)
                                                                                                                          tCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGctgctg                       >  2:1453812/1‑149 (MQ=255)
                                                                                                                          tCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGctgctg                       <  2:2836163/149‑1 (MQ=255)
                                                                                                                            aCGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCtgtg                     <  2:2325521/149‑1 (MQ=255)
                                                                                                                            aCGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCtgtg                     <  2:300280/149‑1 (MQ=255)
                                                                                                                              ggATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGc                   <  1:3892843/149‑1 (MQ=255)
                                                                                                                               gATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGca                  <  1:1739906/149‑1 (MQ=255)
                                                                                                                               gATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGACAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGca                  >  1:3199309/1‑149 (MQ=255)
                                                                                                                               gATTGAGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATCTCGTCGGCGGCTGCTGTGGca                  <  1:4129863/149‑1 (MQ=255)
                                                                                                                                aTTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGcac                 >  1:91999/1‑149 (MQ=255)
                                                                                                                                  tGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGcacca               <  2:1886486/149‑1 (MQ=255)
                                                                                                                                   gCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGcaccac              >  1:3325313/1‑149 (MQ=255)
                                                                                                                                   gCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGcaccac              >  2:517639/1‑149 (MQ=255)
                                                                                                                                      gAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGcc           >  2:2475747/1‑149 (MQ=255)
                                                                                                                                          gCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCacaa       >  1:4639571/1‑149 (MQ=255)
                                                                                                                                           cTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCACTATCGTCGGCGGCTGCTGTGGCACCACGCCacaac      >  1:132146/1‑149 (MQ=255)
                                                                                                                                           cTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCacaac      <  2:2647541/149‑1 (MQ=255)
                                                                                                                                               gTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACatat  >  1:435412/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATAT  >  CP009273/4214417‑4214708

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