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.5% 614.4 / ‑4.6 185T270S (ACC→TCC) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (93/91);  total (94/91)
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.99e-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 185 total aligned reads displayed.

CGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATA  >  CP009273/4214420‑4214707
                                                                                                                                               |                                                                                                                                                
cGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCcgc                                                                                                                                              >  1:3165652/1‑148 (MQ=255)
 gCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCcgcg                                                                                                                                             >  1:1698337/1‑148 (MQ=255)
        gggCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGa                                                                                                                                      <  2:3122791/148‑1 (MQ=255)
        gggCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGa                                                                                                                                      <  2:3540720/148‑1 (MQ=255)
         ggCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGaa                                                                                                                                     <  2:2736108/148‑1 (MQ=255)
         ggCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGaa                                                                                                                                     <  1:2403256/148‑1 (MQ=255)
              accaccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCg                                                                                                                                >  1:1007992/1‑148 (MQ=255)
                caccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCggg                                                                                                                              >  1:309308/1‑148 (MQ=255)
                  ccGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCt                                                                                                                            >  1:2581242/1‑148 (MQ=255)
                     aaGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTAcc                                                                                                                         <  2:2567456/148‑1 (MQ=255)
                     aaGCATTTTACAACTCATGGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTAcc                                                                                                                         <  2:2398029/148‑1 (MQ=255)
                      aGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTAccc                                                                                                                        >  1:1179623/1‑148 (MQ=255)
                      aGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTAccc                                                                                                                        >  2:1406592/1‑148 (MQ=255)
                        cATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCaa                                                                                                                      >  1:1350467/1‑148 (MQ=255)
                        cATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCaa                                                                                                                      >  1:2625128/1‑148 (MQ=255)
                          ttttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACg                                                                                                                    <  2:1866453/148‑1 (MQ=255)
                           tttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGc                                                                                                                   >  1:2302365/1‑148 (MQ=255)
                           tttACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGc                                                                                                                   >  1:2774659/1‑148 (MQ=255)
                             tACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCt                                                                                                                 <  2:42113/148‑1 (MQ=255)
                              aCAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCtt                                                                                                                >  1:3461374/1‑148 (MQ=255)
                                 aCTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTgg                                                                                                             >  1:400393/1‑148 (MQ=255)
                                     aTTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAg                                                                                                         <  1:2484769/148‑1 (MQ=255)
                                       tGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTa                                                                                                       <  1:2379579/148‑1 (MQ=255)
                                       tGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTa                                                                                                       >  2:2952030/1‑148 (MQ=255)
                                          gcCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGa                                                                                                    <  1:1057206/148‑1 (MQ=255)
                                          gcCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGa                                                                                                    <  1:1057154/148‑1 (MQ=255)
                                           ccACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAt                                                                                                   <  1:921872/148‑1 (MQ=255)
                                           ccACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAt                                                                                                   <  2:799957/148‑1 (MQ=255)
                                             aCGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtct                                                                                                 <  1:343014/148‑1 (MQ=255)
                                             aCGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtct                                                                                                 >  1:510271/1‑148 (MQ=255)
                                             aCGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGAtct                                                                                                 >  1:2303781/1‑148 (MQ=255)
                                               gCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCg                                                                                               <  1:323164/148‑1 (MQ=255)
                                                ccGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGa                                                                                              <  1:330436/148‑1 (MQ=255)
                                                  gAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACg                                                                                            >  2:1931702/1‑148 (MQ=255)
                                                   aaGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGc                                                                                           >  2:3322321/1‑148 (MQ=255)
                                                     gCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCg                                                                                         >  2:1897065/1‑148 (MQ=255)
                                                           aCCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCATCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGa                                                                                   <  2:3176200/148‑1 (MQ=255)
                                                           aCCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGa                                                                                   >  2:3294508/1‑148 (MQ=255)
                                                            ccTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGAt                                                                                  >  2:2900127/1‑148 (MQ=255)
                                                             cTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATg                                                                                 >  2:842032/1‑148 (MQ=255)
                                                             cTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATg                                                                                 >  2:2664635/1‑148 (MQ=255)
                                                               ttGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGc                                                                               >  1:728623/1‑148 (MQ=255)
                                                                tGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCa                                                                              <  1:2017563/148‑1 (MQ=255)
                                                                  gCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCaaa                                                                            >  1:1857890/1‑148 (MQ=255)
                                                                  gCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCaaa                                                                            <  2:3451710/148‑1 (MQ=255)
                                                                  gCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCaaa                                                                            >  2:1285766/1‑148 (MQ=255)
                                                                    cTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAAc                                                                          <  2:2735726/148‑1 (MQ=255)
                                                                      gAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAg                                                                        <  1:337739/148‑1 (MQ=255)
                                                                      gAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAg                                                                        >  2:2523532/1‑148 (MQ=255)
                                                                        aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGAt                                                                      <  1:2739356/148‑1 (MQ=255)
                                                                        aCTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGAt                                                                      <  1:1825642/148‑1 (MQ=255)
                                                                         cTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATa                                                                     >  1:1832526/1‑148 (MQ=255)
                                                                           gtgCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACg                                                                   >  1:763352/1‑148 (MQ=255)
                                                                           gtgCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACg                                                                   <  2:102786/148‑1 (MQ=255)
                                                                            tgCGCTGGGGCCCGGTGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGt                                                                  <  1:525200/148‑1 (MQ=255)
                                                                                cTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAAt                                                              <  2:1403399/148‑1 (MQ=255)
                                                                                cTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAAt                                                              <  1:1081135/148‑1 (MQ=255)
                                                                                   gggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATggg                                                           >  2:1165155/1‑148 (MQ=255)
                                                                                   gggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATggg                                                           >  1:3483341/1‑148 (MQ=255)
                                                                                   gggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATggg                                                           >  2:2986865/1‑148 (MQ=255)
                                                                                   gggCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATggg                                                           >  2:1772283/1‑148 (MQ=255)
                                                                                          aTGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGc                                                    >  1:3037379/1‑148 (MQ=255)
                                                                                             aaCTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCggg                                                 >  1:2952656/1‑148 (MQ=255)
                                                                                             aaCTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCggg                                                 <  2:1698337/148‑1 (MQ=255)
                                                                                             aaCTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCggg                                                 <  2:1698332/148‑1 (MQ=255)
                                                                                                  cgcCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTtc                                            >  2:1964011/1‑148 (MQ=255)
                                                                                                    ccAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTtctc                                          >  2:1612363/1‑148 (MQ=255)
                                                                                                     cAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCa                                         >  1:3152155/1‑148 (MQ=255)
                                                                                                     cAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCa                                         >  1:2247671/1‑148 (MQ=255)
                                                                                                      aGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCaa                                        <  2:2150835/148‑1 (MQ=255)
                                                                                                       gTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                                       <  2:311356/148‑1 (MQ=255)
                                                                                                       gTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                                       >  1:1657199/1‑148 (MQ=255)
                                                                                                       gTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAat                                       <  1:2117136/148‑1 (MQ=255)
                                                                                                                  aGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTcggcgg                            <  2:1533049/148‑1 (MQ=255)
                                                                                                                  aGCTGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTcggcgg                            >  1:3552489/1‑148 (MQ=255)
                                                                                                                     tGTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgctg                         <  1:2878957/148‑1 (MQ=255)
                                                                                                                      gTCACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGgctgc                        <  2:499414/148‑1 (MQ=255)
                                                                                                                        cACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGctgctg                      >  2:3365122/1‑148 (MQ=255)
                                                                                                                        cACGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGctgctg                      <  2:3008958/148‑1 (MQ=255)
                                                                                                                         aCGGATTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCtgt                     <  1:758971/148‑1 (MQ=255)
                                                                                                                             aTTGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGca                 >  1:280101/1‑148 (MQ=255)
                                                                                                                               tGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGcacc               <  2:1422714/148‑1 (MQ=255)
                                                                                                                               tGCGGAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGcacc               <  2:46126/148‑1 (MQ=255)
                                                                                                                                   gAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGc           <  1:2655811/148‑1 (MQ=255)
                                                                                                                                   gAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGc           >  1:2736672/1‑148 (MQ=255)
                                                                                                                                   gAATGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGc           >  1:2836708/1‑148 (MQ=255)
                                                                                                                                     aTGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCca         >  1:1842108/1‑148 (MQ=255)
                                                                                                                                     aTGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCca         <  1:1405471/148‑1 (MQ=255)
                                                                                                                                     aTGCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCca         <  1:1406592/148‑1 (MQ=255)
                                                                                                                                       gCTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCcaca       >  2:3256692/1‑148 (MQ=255)
                                                                                                                                        cTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCacaa      <  2:1501786/148‑1 (MQ=255)
                                                                                                                                        cTACGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCacaa      <  1:187349/148‑1 (MQ=255)
                                                                                                                                          aCGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCacaaca    >  2:2266883/1‑148 (MQ=255)
                                                                                                                                           cGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACat   <  1:2965338/148‑1 (MQ=255)
                                                                                                                                           cGTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACat   <  1:1490518/148‑1 (MQ=255)
                                                                                                                                            gTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACata  <  1:2952030/148‑1 (MQ=255)
                                                                                                                                            gTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACata  <  2:1549374/148‑1 (MQ=255)
                                                                                                                                            gTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACata  <  1:1772283/148‑1 (MQ=255)
                                                                                                                                            gTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACata  <  1:253580/148‑1 (MQ=255)
                                                                                                                                            gTCTCCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACata  <  2:264571/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CGCTCTCCGGGCAGACCACCGAAGCATTTTACAACTCATTGCGCCACGCCGAAGCTCTGACCTTTGGCCTGAACTGTGCGCTGGGGCCCGATGAACTGCGCCAGTACGTGCAGGAGCTGTCACGGATTGCGGAATGCTACGTCACCGCGCACCCGAACGCCGGGCTACCCAACGCCTTTGGTGAGTACGATCTCGACGCCGACACGATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATA  >  CP009273/4214420‑4214707

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: ‑