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,5630AT96.8% 600.5 / 7.4 186T270S (ACC→TCC) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base T (87/93);  total (90/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 186 total aligned reads displayed.

CACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACAT  >  CP009273/4214418‑4214706
                                                                                                                                                 |                                                                                                                                               
cACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTcc                                                                                                                                               >  1:3447163/1‑148 (MQ=255)
       tcCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCAcc                                                                                                                                        <  2:2378152/148‑1 (MQ=255)
        ccGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCAccc                                                                                                                                       >  2:813307/1‑148 (MQ=255)
        ccGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCAccc                                                                                                                                       >  2:2289833/1‑148 (MQ=255)
           ggCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGaa                                                                                                                                    >  1:967250/1‑148 (MQ=255)
              aGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGc                                                                                                                                 <  2:1383926/148‑1 (MQ=255)
                accaccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCg                                                                                                                               >  2:329012/1‑148 (MQ=255)
                accaccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCg                                                                                                                               >  2:580798/1‑148 (MQ=255)
                  caccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCggg                                                                                                                             <  2:2063175/148‑1 (MQ=255)
                  caccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCggg                                                                                                                             <  2:343392/148‑1 (MQ=255)
                   accGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGc                                                                                                                            >  2:2752282/1‑148 (MQ=255)
                    ccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCt                                                                                                                           >  1:2734165/1‑148 (MQ=255)
                    ccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCt                                                                                                                           >  1:3149919/1‑148 (MQ=255)
                    ccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCt                                                                                                                           >  1:2468935/1‑148 (MQ=255)
                     cGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTa                                                                                                                          <  1:1874183/148‑1 (MQ=255)
                       aaGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTAcc                                                                                                                        <  1:1983695/148‑1 (MQ=255)
                         gCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAATGCCGGGCTACCCa                                                                                                                      >  2:1514962/1‑148 (MQ=255)
                          cATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCaa                                                                                                                     <  1:2310810/148‑1 (MQ=255)
                          cATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCaa                                                                                                                     >  2:998796/1‑148 (MQ=255)
                            ttttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACg                                                                                                                   >  2:2217547/1‑148 (MQ=255)
                            ttttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACg                                                                                                                   >  2:2217814/1‑148 (MQ=255)
                            ttttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACg                                                                                                                   >  2:495005/1‑148 (MQ=255)
                             tttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTCCCCAACGc                                                                                                                  >  1:1003173/1‑148 (MQ=255)
                             tttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGc                                                                                                                  <  2:1601711/148‑1 (MQ=255)
                             tttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGc                                                                                                                  >  1:2765708/1‑148 (MQ=255)
                             tttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGc                                                                                                                  >  1:2779210/1‑148 (MQ=255)
                              ttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGcc                                                                                                                 >  2:1102530/1‑148 (MQ=255)
                              ttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGcc                                                                                                                 <  1:2518696/148‑1 (MQ=255)
                                 cAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCttt                                                                                                              >  1:1261767/1‑148 (MQ=255)
                                  aaCTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCtgtg                                                                                                             >  2:2532110/1‑148 (MQ=255)
                                  aaCTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTg                                                                                                             >  1:2825074/1‑148 (MQ=255)
                                   aCTCATTGCGCCACGCCGGAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTgg                                                                                                            <  2:2952352/148‑1 (MQ=255)
                                    cTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGt                                                                                                           >  2:2639849/1‑148 (MQ=255)
                                        ttGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGt                                                                                                       <  2:2479321/148‑1 (MQ=255)
                                        ttGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGt                                                                                                       <  2:1951387/148‑1 (MQ=255)
                                              cACGCCGAAGCTCTGACCTTTGGCCTTAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtc                                                                                                 <  1:2284993/148‑1 (MQ=255)
                                                  ccGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGa                                                                                             <  1:2708150/148‑1 (MQ=255)
                                                    gAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACg                                                                                           >  2:1591914/1‑148 (MQ=255)
                                                    gAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACg                                                                                           <  1:239099/148‑1 (MQ=255)
                                                      aGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGcc                                                                                         >  1:1186776/1‑148 (MQ=255)
                                                       gCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCg                                                                                        >  2:2352358/1‑148 (MQ=255)
                                                        ctctGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGa                                                                                       >  2:1169399/1‑148 (MQ=255)
                                                          ctGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGaca                                                                                     >  1:3295772/1‑148 (MQ=255)
                                                          ctGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGaca                                                                                     >  2:1350579/1‑148 (MQ=255)
                                                             aCCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGa                                                                                  >  2:664746/1‑148 (MQ=255)
                                                             aCCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGa                                                                                  >  2:665616/1‑148 (MQ=255)
                                                                tttGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATgg                                                                               >  1:2210991/1‑148 (MQ=255)
                                                                    gCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCaaa                                                                           <  1:2221259/148‑1 (MQ=255)
                                                                    gCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCaaa                                                                           <  1:3076198/148‑1 (MQ=255)
                                                                      cTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAAc                                                                         <  1:2469345/148‑1 (MQ=255)
                                                                       tGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACa                                                                        <  1:2960267/148‑1 (MQ=255)
                                                                        gAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAg                                                                       >  1:66070/1‑148 (MQ=255)
                                                                          aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCGCGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGAt                                                                     <  2:604468/148‑1 (MQ=255)
                                                                          aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCTACGCCGACACGATGGCAAAACAGAt                                                                     >  1:1013092/1‑148 (MQ=255)
                                                                          aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGAt                                                                     >  1:3250862/1‑148 (MQ=255)
                                                                          aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGAt                                                                     <  2:1705570/148‑1 (MQ=255)
                                                                          aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGAt                                                                     >  1:3251035/1‑148 (MQ=255)
                                                                           cTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATa                                                                    <  1:704557/148‑1 (MQ=255)
                                                                             gtgCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACg                                                                  <  1:1921533/148‑1 (MQ=255)
                                                                                 gcTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGaa                                                              <  2:312385/148‑1 (MQ=255)
                                                                                  cTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCGGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAAt                                                             <  1:2661004/148‑1 (MQ=255)
                                                                                   tGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATg                                                            >  2:3402525/1‑148 (MQ=255)
                                                                                     gggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATggg                                                          <  2:1016838/148‑1 (MQ=255)
                                                                                     gggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATggg                                                          <  1:2877981/148‑1 (MQ=255)
                                                                                         ccGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCa                                                      <  1:2835850/148‑1 (MQ=255)
                                                                                           gATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAg                                                    >  1:429164/1‑148 (MQ=255)
                                                                                           gATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAg                                                    >  1:428764/1‑148 (MQ=255)
                                                                                            aTGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGc                                                   <  1:1481899/148‑1 (MQ=255)
                                                                                                  tGCGCCAGTACGTGCAGGAGCTGTCTCGGATTGCGGAATGCTACGTCTCCGCGCACCCGATCGCCGGGCTACCCAACGCCTTTGGTGCGTACGATCTCGACGCCGACTCGATGGCAAAACATATACGTGAATGGGCGCAAGCGGGttt                                             <  2:1796152/148‑1 (MQ=255)
                                                                                                  tGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGttt                                             <  1:1658604/148‑1 (MQ=255)
                                                                                                     gcCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTtct                                          <  1:329012/148‑1 (MQ=255)
                                                                                                      ccAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTtctc                                         >  2:3301236/1‑148 (MQ=255)
                                                                                                       cAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCa                                        <  2:2409398/148‑1 (MQ=255)
                                                                                                       cAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCa                                        >  2:2529917/1‑148 (MQ=255)
                                                                                                        aGTACGTGCAGGAGCTGTCACGGAGTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCaa                                       <  2:3450898/148‑1 (MQ=255)
                                                                                                         gTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                                      <  2:1881808/148‑1 (MQ=255)
                                                                                                         gTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                                      <  1:1344729/148‑1 (MQ=255)
                                                                                                          tACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAata                                     >  2:1648407/1‑148 (MQ=255)
                                                                                                          tACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAata                                     <  1:3375202/148‑1 (MQ=255)
                                                                                                          tACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAata                                     <  1:620141/148‑1 (MQ=255)
                                                                                                             gTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcg                                  <  2:906993/148‑1 (MQ=255)
                                                                                                             gTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcg                                  <  2:967250/148‑1 (MQ=255)
                                                                                                              tGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAATGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgt                                 <  1:1514962/148‑1 (MQ=255)
                                                                                                              tGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGACTGGGCGCACGCGGGTTTTCTCAATAtcgt                                 >  2:1011691/1‑148 (MQ=255)
                                                                                                               gCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgtc                                <  2:2760542/148‑1 (MQ=255)
                                                                                                               gCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATAtcgtc                                <  1:1185600/148‑1 (MQ=255)
                                                                                                                     gCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCggcggc                          <  1:2987379/148‑1 (MQ=255)
                                                                                                                      cTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgct                         <  1:255701/148‑1 (MQ=255)
                                                                                                                       tGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgctg                        <  1:2126438/148‑1 (MQ=255)
                                                                                                                         tCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgctgct                      >  2:1589611/1‑148 (MQ=255)
                                                                                                                          cACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGctgctg                     <  1:3455330/148‑1 (MQ=255)
                                                                                                                            cGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCtgtg                   >  2:347309/1‑148 (MQ=255)
                                                                                                                              gATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGc                 >  2:1165397/1‑148 (MQ=255)
                                                                                                                                  gCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGANACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGcacca             <  1:1806137/148‑1 (MQ=255)
                                                                                                                                    ggAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACg           <  1:46961/148‑1 (MQ=255)
                                                                                                                                      aaTGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAatat                                    <  2:356282/121‑1 (MQ=255)
                                                                                                                                          cTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCacaa     <  2:602205/148‑1 (MQ=255)
                                                                                                                                           tACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCacaac    <  1:467568/148‑1 (MQ=255)
                                                                                                                                            aCGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCacaaca   <  2:3149919/148‑1 (MQ=255)
                                                                                                                                             cGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACat  <  2:2468959/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CACGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACAT  >  CP009273/4214418‑4214706

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: