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,8090GT97.2% 509.6 / ‑0.5 190E352* (GAA→TAA) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (10/0);  new base T (93/87);  total (103/87)
Fisher's exact test for biased strand distribution p-value = 2.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 190 total aligned reads displayed.

CTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATA  >  CP009273/4214665‑4214953
                                                                                                                                                |                                                                                                                                                
  cAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGAt                                                                                                                                             <  1:1612363/148‑1 (MQ=255)
  cAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGAt                                                                                                                                             >  2:1798083/1‑148 (MQ=255)
     tatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGc                                                                                                                                          <  2:1832526/148‑1 (MQ=255)
       tcgtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCt                                                                                                                                        <  2:280101/148‑1 (MQ=255)
              cggcTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTtgtg                                                                                                                                 <  1:448728/148‑1 (MQ=255)
                  tgctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACg                                                                                                                             >  1:2529965/1‑148 (MQ=255)
                    ctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTg                                                                                                                           <  2:3045643/148‑1 (MQ=255)
                       tgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGt                                                                                                                        <  2:2625128/148‑1 (MQ=255)
                        ggCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTg                                                                                                                       <  1:1771357/148‑1 (MQ=255)
                          caccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTAAGATAGCCAGTTTTTGTACGTGGGTGaa                                                                                                                     >  1:1589412/1‑148 (MQ=255)
                           accacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAAc                                                                                                                    >  1:1349298/1‑148 (MQ=255)
                            ccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGGTTGTGAACGTGGGTGAACg                                                                                                                   >  1:52874/1‑148 (MQ=255)
                            ccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGAACATTGGGGAAGATATCCTGTTTGTGAAAGTGGGTGTACg                                                                                                                   >  1:1530213/1‑148 (MQ=255)
                             cacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGc                                                                                                                  >  1:2863464/1‑148 (MQ=255)
                             cacGCCACAACATAGTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGc                                                                                                                  <  1:1964011/148‑1 (MQ=255)
                               cGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCAc                                                                                                                >  2:2160511/1‑148 (MQ=255)
                                 ccACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCa                                                                                                              <  2:2952656/148‑1 (MQ=255)
                                  cacaACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCaa                                                                                                             <  1:1631116/148‑1 (MQ=255)
                                  cacaACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCaa                                                                                                             >  2:856125/1‑148 (MQ=255)
                                    caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTAAGATAGCCTGTTTGTGAAAGTGGGTGAACGCACCAACg                                                                                                           >  1:2812243/1‑148 (MQ=255)
                                    caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           <  1:3256692/148‑1 (MQ=255)
                                    caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           >  2:778862/1‑148 (MQ=255)
                                    caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           <  1:2266883/148‑1 (MQ=255)
                                        atatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCAc                                                                                                       >  1:2309999/1‑148 (MQ=255)
                                         tatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCAcc                                                                                                      >  2:1258609/1‑148 (MQ=255)
                                         tatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCAcc                                                                                                      <  1:2226545/148‑1 (MQ=255)
                                            tGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGt                                                                                                   <  2:1553532/148‑1 (MQ=255)
                                             gCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGtt                                                                                                  <  2:2045374/148‑1 (MQ=255)
                                                 cGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATATCCTGTTTGTGAACGTGGGTGAACGCACCCACGTCACCCGGTCCGc                                                                                              >  1:3005084/1‑148 (MQ=255)
                                                 cGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGc                                                                                              >  1:2068771/1‑148 (MQ=255)
                                                 cGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGc                                                                                              >  1:1190796/1‑148 (MQ=255)
                                                 cGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATATCCTGTTTGTGTACGTGGGTGAACGCACCCACCGCACCCGTTCCGc                                                                                              >  1:677064/1‑148 (MQ=255)
                                                  gATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCt                                                                                             <  2:1179623/148‑1 (MQ=255)
                                                   ctGAGGCGTGCAGTAGAAGGATTAGCGCCGCGCACACTGCCGGAAATTCGCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTa                                                                                            <  2:1855764/147‑1 (MQ=255)
                                                   aTGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTa                                                                                            <  1:2900127/148‑1 (MQ=255)
                                                   aTGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTa                                                                                            >  1:1915785/1‑148 (MQ=255)
                                                     gAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCGAAGATATCCTGTTTTTGAACGTGGGTGAACGCACCCACCTCACCCGTTCCGCTAAg                                                                                          >  1:2169178/1‑148 (MQ=255)
                                                     gAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATATCCTGTTTGTGAACGTGGGTGAACGCCCCAACGTCACCGGTTCCGCTAAg                                                                                          >  1:2582093/1‑148 (MQ=255)
                                                     gAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAg                                                                                          >  1:2582023/1‑148 (MQ=255)
                                                         cGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCa                                                                                      <  2:2303781/148‑1 (MQ=255)
                                                          gTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCaa                                                                                     <  2:460768/148‑1 (MQ=255)
                                                            gCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGGTCCGCTAAGTTCAAgc                                                                                   >  1:1996466/1‑148 (MQ=255)
                                                              aGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgcgc                                                                                 <  2:1657199/148‑1 (MQ=255)
                                                               gTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGcc                                                                                <  2:1751767/148‑1 (MQ=255)
                                                               gTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGcc                                                                                <  2:1343721/148‑1 (MQ=255)
                                                                tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               >  2:253201/1‑148 (MQ=255)
                                                                tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCCGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               >  2:2414133/1‑148 (MQ=255)
                                                                 aGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTg                                                                              >  1:2360616/1‑148 (MQ=255)
                                                                   aaGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGAt                                                                            >  1:1044803/1‑148 (MQ=255)
                                                                     ggATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                                                          >  1:2090052/1‑148 (MQ=255)
                                                                         tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                      <  2:3381383/148‑1 (MQ=255)
                                                                          aGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaagaa                                                                     <  1:1349989/148‑1 (MQ=255)
                                                                           gcgcCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaagaag                                                                    >  1:1482010/1‑148 (MQ=255)
                                                                           gcgcCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaagaag                                                                    <  1:1399472/148‑1 (MQ=255)
                                                                              ccgcGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAgaga                                                                 <  1:1547594/148‑1 (MQ=255)
                                                                                 cgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAAt                                                              >  2:280647/1‑148 (MQ=255)
                                                                                 cgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAAt                                                              >  2:279644/1‑148 (MQ=255)
                                                                                 cgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAAt                                                              <  1:679514/148‑1 (MQ=255)
                                                                                 cgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAAt                                                              >  2:3410069/1‑148 (MQ=255)
                                                                                   cAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATAc                                                            <  2:3037379/148‑1 (MQ=255)
                                                                                   cAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCGAAGATAGCCTGTTTGTGTACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATAc                                                            >  1:3487866/1‑148 (MQ=255)
                                                                                      aCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGc                                                         >  2:1594733/1‑148 (MQ=255)
                                                                                       cTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCg                                                        >  2:223437/1‑148 (MQ=255)
                                                                                          ccGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAgg                                                     <  2:2836450/148‑1 (MQ=255)
                                                                                            ggAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcg                                                   <  1:3157978/148‑1 (MQ=255)
                                                                                             gAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcgc                                                  >  2:1398967/1‑148 (MQ=255)
                                                                                             gAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcgc                                                  >  2:690810/1‑148 (MQ=255)
                                                                                              aaaTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCt                                                 >  2:380785/1‑148 (MQ=255)
                                                                                                   cccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATg                                            <  1:1311220/148‑1 (MQ=255)
                                                                                                   cccGGAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATg                                            <  2:2736672/148‑1 (MQ=255)
                                                                                                    ccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGt                                           >  1:3545691/1‑148 (MQ=255)
                                                                                                     cGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAAAAATTGGGAACATATCCTTTTTGTGAAAGTTGGTGAACGCACCAAACGCACCCGTTCCGCTAAGTTAAAGCGCCTGATCAAAGAAAAGAAATACAGCGAGGCGCTgggtttg                                          >  1:80726/1‑143 (MQ=255)
                                                                                                          ccTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgcgcg                                     >  2:1924445/1‑148 (MQ=255)
                                                                                                          ccTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgcgcg                                     <  1:2664635/148‑1 (MQ=255)
                                                                                                           cTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGt                                    <  1:2883066/148‑1 (MQ=255)
                                                                                                           cTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGt                                    <  2:953098/148‑1 (MQ=255)
                                                                                                           cTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGt                                    <  1:1201747/148‑1 (MQ=255)
                                                                                                            tGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCGAAGATATCCTGTTTGTTAACGTGGGGGAACGCACCAACCTCACCCGTTCCGCGAAGTTCAAACCCCGGATAAAATAAGAGAAAAAAACCGAGGCGCTGGATGTCGCGCGTc                                   >  1:140285/1‑148 (MQ=255)
                                                                                                              ccGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCCCCCGTTCCGCGAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACATCGAGGCGCTGGATGTCGCGCGTCaa                                 >  1:3364792/1‑148 (MQ=255)
                                                                                                              ccGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCaa                                 <  2:2169178/148‑1 (MQ=255)
                                                                                                               cGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAAc                                <  1:1123361/148‑1 (MQ=255)
                                                                                                                  ttGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAgg                             >  2:357696/1‑148 (MQ=255)
                                                                                                                            tGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACgg                   <  1:279644/148‑1 (MQ=255)
                                                                                                                            tGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACgg                   <  1:280647/148‑1 (MQ=255)
                                                                                                                            tGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACgg                   <  1:3410069/148‑1 (MQ=255)
                                                                                                                             ggAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgc                  >  2:3319801/1‑148 (MQ=255)
                                                                                                                              gAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcg                 >  1:663623/1‑148 (MQ=255)
                                                                                                                               aGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgc                >  2:2747121/1‑148 (MQ=255)
                                                                                                                               aGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgc                >  2:2747124/1‑148 (MQ=255)
                                                                                                                                gCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgcg               >  2:2455466/1‑148 (MQ=255)
                                                                                                                                gCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgcg               >  2:748293/1‑148 (MQ=255)
                                                                                                                                   gCTGAACTTTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAg            <  1:2414152/148‑1 (MQ=255)
                                                                                                                                   gCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACAAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAg            <  1:3189819/148‑1 (MQ=255)
                                                                                                                                     tGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAt          >  1:1604404/1‑148 (MQ=255)
                                                                                                                                         cATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATc      >  1:1619692/1‑148 (MQ=255)
                                                                                                                                           ttGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGa    >  2:2378513/1‑148 (MQ=255)
                                                                                                                                            tGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGat   >  1:916765/1‑148 (MQ=255)
                                                                                                                                            tGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGat   >  1:3195047/1‑148 (MQ=255)
                                                                                                                                             ggCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATCCAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGata  <  2:146945/148‑1 (MQ=255)
                                                                                                                                             ggCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGata  <  2:1731495/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATA  >  CP009273/4214665‑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: