Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,214,809 G→T 100% E352* (GAA→TAA)  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,8090GT100.0% 582.6 / ‑6.7 201E352* (GAA→TAA) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base T (103/92);  total (109/92)
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 201 total aligned reads displayed.

TCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATA  >  CP009273/4214666‑4214953
                                                                                                                                               |                                                                                                                                                
 cAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGAt                                                                                                                                             >  1:1835317/1‑148 (MQ=255)
    tatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTCAGATAGc                                                                                                                                          >  1:3286819/1‑148 (MQ=255)
    tatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGc                                                                                                                                          >  2:767807/1‑148 (MQ=255)
    tatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGc                                                                                                                                          >  2:412107/1‑148 (MQ=255)
    tatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTACCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGc                                                                                                                                          >  2:1265157/1‑148 (MQ=255)
     atCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGcc                                                                                                                                         <  2:899019/148‑1 (MQ=255)
       cgtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTg                                                                                                                                       >  2:1777231/1‑148 (MQ=255)
        gtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGt                                                                                                                                      <  1:2213443/148‑1 (MQ=255)
         tcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTGCGGCCTGGAGCCGCGGAACAATTGCTAAGATAGCCTGtt                                                                                                                                     >  1:1347708/1‑148 (MQ=255)
          cggcggCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGttt                                                                                                                                    <  1:424580/148‑1 (MQ=255)
              ggcTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGa                                                                                                                                <  1:885117/148‑1 (MQ=255)
               gctgctGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGaa                                                                                                                               <  2:3294361/148‑1 (MQ=255)
                ctgctgTGACACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAAc                                                                                                                              >  2:2836715/1‑148 (MQ=255)
                   ctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACCTTGGCTAAGATAGCCTGTTTGTGAACGTg                                                                                                                           >  2:2255900/1‑148 (MQ=255)
                   ctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTg                                                                                                                           <  1:1114007/148‑1 (MQ=255)
                     gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATCGCCTGTTTGTGAACGTggg                                                                                                                         >  2:2634101/1‑148 (MQ=255)
                     gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTggg                                                                                                                         >  1:1507264/1‑148 (MQ=255)
                     gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTg                                                                                                                           >  2:570907/1‑146 (MQ=255)
                     gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTg                                                                                                                           <  1:570907/146‑1 (MQ=255)
                        gCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATATCCTGTTTGTGAACGTGGGTGa                                                                                                                      <  2:1599570/148‑1 (MQ=255)
                           ccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACg                                                                                                                   >  2:340028/1‑148 (MQ=255)
                            cacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGc                                                                                                                  >  1:819343/1‑148 (MQ=255)
                            cacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGc                                                                                                                  >  1:1148151/1‑148 (MQ=255)
                            cacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGc                                                                                                                  >  2:501538/1‑148 (MQ=255)
                             acGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCa                                                                                                                 <  2:1359678/148‑1 (MQ=255)
                                ccACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCa                                                                                                              <  2:3250851/148‑1 (MQ=255)
                                ccACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCa                                                                                                              <  2:1993287/148‑1 (MQ=255)
                                   caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           >  2:2293326/1‑148 (MQ=255)
                                   caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTGTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           >  1:1040793/1‑148 (MQ=255)
                                       atatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACCTccc                                                                                                       >  1:3663234/1‑146 (MQ=255)
                                         atTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTTTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                                                                     >  2:1920410/1‑148 (MQ=255)
                                         atTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                                                                     <  2:1842515/148‑1 (MQ=255)
                                         atTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                                                                     <  2:1845804/148‑1 (MQ=255)
                                            gCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGATATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGtt                                                                                                  <  1:2328035/148‑1 (MQ=255)
                                             cAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTc                                                                                                 <  1:1940823/148‑1 (MQ=255)
                                               gCGATGCGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCg                                                                                               <  2:1523468/148‑1 (MQ=255)
                                                cGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGATTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGc                                                                                              >  2:191674/1‑148 (MQ=255)
                                                  aTGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTa                                                                                            >  2:2937920/1‑148 (MQ=255)
                                                     aGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGt                                                                                         >  2:2807239/1‑148 (MQ=255)
                                                     aGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCTAAGATAGCCAGTTTTTGAACGTGGGTGAACGCACCAACCTCACCGGTTCCGCTAAGt                                                                                         >  1:503250/1‑148 (MQ=255)
                                                      gTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGtt                                                                                        <  2:128992/148‑1 (MQ=255)
                                                       tCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGACCGCACCAACGTCACCGGTTCCGCTAAGTTc                                                                                       >  2:1570275/1‑148 (MQ=255)
                                                       tCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTc                                                                                       >  1:3012407/1‑148 (MQ=255)
                                                          tGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAg                                                                                    <  1:1492984/148‑1 (MQ=255)
                                                             aGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgcgc                                                                                 >  1:3091725/1‑148 (MQ=255)
                                                              gTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGcc                                                                                <  1:3513430/148‑1 (MQ=255)
                                                               tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               >  2:2406807/1‑148 (MQ=255)
                                                               tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               >  2:1055481/1‑148 (MQ=255)
                                                               tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               <  1:3257302/148‑1 (MQ=255)
                                                                  aaGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATATCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGAt                                                                            >  1:899671/1‑148 (MQ=255)
                                                                    ggATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                                                          >  2:2082306/1‑148 (MQ=255)
                                                                    ggATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                                                          >  2:1469864/1‑148 (MQ=255)
                                                                    ggATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                                                          >  2:3674413/1‑148 (MQ=255)
                                                                     gATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTGTGTCCGGCCTGGAGCCGCTGAACCTTGGCTAAGATAGCCTGTTTGTGAACGGGGGTGAACGCACCAAAGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCaa                                                                         >  1:835526/1‑148 (MQ=255)
                                                                        tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                      >  2:2829153/1‑148 (MQ=255)
                                                                           cgccgcGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTTTCCGGCCTGGAGCCGCTGTACAATTGCAAAGAAAGCCTGTTTTTGAACTTGGGTGAACGCACCAACATCAACGGTTACCCGAAGTTCAAGCGCCCGATCAAagaaga                                                                   >  1:1218191/1‑148 (MQ=255)
                                                                           cgccgcGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTTTGAACGTGGGTGAACGCACCAACCTCACCGGTTACGCTAACTTCAATCGCCTGATCAAataaga                                                                   >  1:1215589/1‑148 (MQ=255)
                                                                             ccgcGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAgaga                                                                 <  2:3325409/148‑1 (MQ=255)
                                                                                 gcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATa                                                             <  2:1958965/148‑1 (MQ=255)
                                                                                 gcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATa                                                             <  1:636785/148‑1 (MQ=255)
                                                                                     aCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGc                                                         >  2:451260/1‑148 (MQ=255)
                                                                                       tGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGa                                                       <  2:3012407/148‑1 (MQ=255)
                                                                                       tGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGa                                                       <  1:3094899/148‑1 (MQ=255)
                                                                                        gCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGGTAAGATAGCCGGTTTTTGTAAGTGGGTGTAAGCCCCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCGGATCAAATAATATAAATAAATCGAg                                                      >  1:1750648/1‑148 (MQ=255)
                                                                                            gAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGTACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcgc                                                  >  1:3139580/1‑148 (MQ=255)
                                                                                            gAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcgc                                                  <  2:3695402/148‑1 (MQ=255)
                                                                                            gAAATTCCCGTAGCCTGCCGTTTGTCCGGCCGGGAGCCGCTGAAAATTGGGGAAGATAGCCTGTTTTTGTACGTGGGTGAACGCACCCACGTCACCGGTTGCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcgc                                                  >  1:1988572/1‑148 (MQ=255)
                                                                                              aaTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTg                                                >  1:2525094/1‑148 (MQ=255)
                                                                                               aTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTgg                                               <  1:767807/148‑1 (MQ=255)
                                                                                                ttCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGa                                              <  2:2185592/148‑1 (MQ=255)
                                                                                                  cccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATg                                            <  2:2661939/148‑1 (MQ=255)
                                                                                                  cccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATg                                            <  2:1148151/148‑1 (MQ=255)
                                                                                                  cccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATg                                            <  2:819343/148‑1 (MQ=255)
                                                                                                      tAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGCTGTcgc                                        >  2:2852663/1‑148 (MQ=255)
                                                                                                      tAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgc                                        >  2:2852722/1‑148 (MQ=255)
                                                                                                      tAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgc                                        >  2:3427216/1‑148 (MQ=255)
                                                                                                       aGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgcg                                       >  2:3355563/1‑148 (MQ=255)
                                                                                                            gCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCACGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCa                                  <  1:1612427/148‑1 (MQ=255)
                                                                                                             ccGTTTGTCCGGCCTGGTGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCaa                                 <  1:3124431/148‑1 (MQ=255)
                                                                                                               gTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACa                               >  1:959154/1‑148 (MQ=255)
                                                                                                               gTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACa                               >  2:2293849/1‑148 (MQ=255)
                                                                                                                 ttGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAgg                             >  2:2552139/1‑148 (MQ=255)
                                                                                                                 ttGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAgg                             <  1:1253454/148‑1 (MQ=255)
                                                                                                                  tGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGt                            <  2:2682580/148‑1 (MQ=255)
                                                                                                                        gCCTGGAGCCGCTGAACATTGGCTATGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGaaaa                      <  2:2202234/148‑1 (MQ=255)
                                                                                                                         ccTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAATAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAAc                     >  2:1772004/1‑148 (MQ=255)
                                                                                                                           tGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACgg                   <  1:2082303/148‑1 (MQ=255)
                                                                                                                             gAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcg                 >  1:2666181/1‑148 (MQ=255)
                                                                                                                                  gCTGAACATGGGCTAAGATAGCCTGTTTGTGAAAGTGGGTGAACGCACCAACGTCACCGGTTCCTCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATCCAGCGTGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAg            <  2:683099/148‑1 (MQ=255)
                                                                                                                                   cTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGa           <  2:3170898/148‑1 (MQ=255)
                                                                                                                                     gAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCTCAGAtt         <  2:136019/148‑1 (MQ=255)
                                                                                                                                     gAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAtt         >  2:3112723/1‑148 (MQ=255)
                                                                                                                                     gAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAtt         <  1:2432818/148‑1 (MQ=255)
                                                                                                                                     gAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAtt         <  2:3642217/148‑1 (MQ=255)
                                                                                                                                      aaCATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTa        <  2:1835317/148‑1 (MQ=255)
                                                                                                                                      aaCATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTa        >  1:588339/1‑148 (MQ=255)
                                                                                                                                       aCATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTAt       <  1:2028177/148‑1 (MQ=255)
                                                                                                                                         aTTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAACTACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCg     >  2:2988716/1‑148 (MQ=255)
                                                                                                                                           tGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCTTCAACAGGTGGAAAACGGCGCGCAGATTATCGat   <  2:288703/148‑1 (MQ=255)
                                                                                                                                            ggCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGATAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGata  <  1:382764/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATA  >  CP009273/4214666‑4214953

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