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,8090GT95.6% 577.8 / 5.7 221E352* (GAA→TAA) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base G (12/2);  new base T (122/85);  total (134/87)
Fisher's exact test for biased strand distribution p-value = 5.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 221 total aligned reads displayed.

TCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATA  >  CP009273/4214664‑4214953
                                                                                                                                                  |                                                                                                                                                
 ctcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAg                                                                                                                                               >  1:2480153/1‑148 (MQ=255)
 ctcAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAg                                                                                                                                               >  1:2927766/1‑148 (MQ=255)
   cAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGAt                                                                                                                                             <  2:1594675/148‑1 (MQ=255)
    aaTATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGaaa                                                                                                                                            >  1:2430776/1‑146 (MQ=255)
    aaTATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATa                                                                                                                                            >  2:3157706/1‑148 (MQ=255)
     atatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAg                                                                                                                                           <  1:1658261/148‑1 (MQ=255)
          gtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGt                                                                                                                                      <  1:3191059/148‑1 (MQ=255)
          gtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGt                                                                                                                                      <  1:3413957/148‑1 (MQ=255)
                  ctgctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAAc                                                                                                                              >  2:3048718/1‑148 (MQ=255)
                  ctgctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAAc                                                                                                                              >  2:1848596/1‑148 (MQ=255)
                  ctgctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAAc                                                                                                                              >  1:356936/1‑148 (MQ=255)
                  ctgctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGAACATTGGGGAAGATAGCCTGTTTGTGAAc                                                                                                                              >  1:2655940/1‑148 (MQ=255)
                    gctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGt                                                                                                                            <  1:1039365/148‑1 (MQ=255)
                    gctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGGTAATATAGCCTGTTTTTGaaact                                                                                                                            >  1:783950/1‑145 (MQ=255)
                     ctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTg                                                                                                                           >  2:3327551/1‑148 (MQ=255)
                       gtgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTggg                                                                                                                         <  2:2813850/148‑1 (MQ=255)
                        tgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGt                                                                                                                        >  2:2234757/1‑148 (MQ=255)
                           caccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGaa                                                                                                                     <  2:1006710/148‑1 (MQ=255)
                           caccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGaa                                                                                                                     <  1:950035/148‑1 (MQ=255)
                             ccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACg                                                                                                                   >  2:2778857/1‑148 (MQ=255)
                               acGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCa                                                                                                                 <  2:2572001/148‑1 (MQ=255)
                                 gCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCAcc                                                                                                               >  1:3437118/1‑148 (MQ=255)
                                    acaacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAAc                                                                                                            >  1:3094563/1‑148 (MQ=255)
                                    acaacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGGACAATGGCGAATATAGCCAGTTTTTGAAAGTGGGGGAACGCACCCAc                                                                                                            >  1:1200988/1‑148 (MQ=255)
                                     caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           >  1:521368/1‑148 (MQ=255)
                                     caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           >  1:1455975/1‑148 (MQ=255)
                                     caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACg                                                                                                           >  1:2933015/1‑148 (MQ=255)
                                     caacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGGGGGTGAACGCACCAACg                                                                                                           >  1:521122/1‑148 (MQ=255)
                                        caTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCa                                                                                                        >  1:424491/1‑148 (MQ=255)
                                         atatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGTACGCACCAACGTCAc                                                                                                       >  1:1025424/1‑148 (MQ=255)
                                          tatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCAcc                                                                                                      >  2:2901034/1‑148 (MQ=255)
                                          tatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCTAAGATAGCCTGTTTGTGAACGTGGGTGTAAGCACCAACGTCAcc                                                                                                      >  1:2985517/1‑148 (MQ=255)
                                           atTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                                                                     <  1:2638082/148‑1 (MQ=255)
                                           atTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                                                                     >  1:2794343/1‑148 (MQ=255)
                                             tGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGt                                                                                                   >  1:2043635/1‑148 (MQ=255)
                                             tGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGt                                                                                                   <  2:1076000/148‑1 (MQ=255)
                                             tGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGt                                                                                                   >  2:2009252/1‑148 (MQ=255)
                                              gCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGtt                                                                                                  <  1:2315563/148‑1 (MQ=255)
                                               cAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTc                                                                                                 >  2:1196068/1‑148 (MQ=255)
                                               cAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTc                                                                                                 >  2:913888/1‑148 (MQ=255)
                                               cAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTc                                                                                                 >  1:1669080/1‑148 (MQ=255)
                                                  cGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATATCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGc                                                                                              >  1:3002159/1‑148 (MQ=255)
                                                    ctGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTa                                                                                            <  2:2522151/147‑1 (MQ=255)
                                                    aTGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCACACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTa                                                                                            <  2:3383551/148‑1 (MQ=255)
                                                    aTGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTa                                                                                            >  2:441983/1‑148 (MQ=255)
                                                    aTGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCCGCCTGGAGCCGCGGAAAAATGGCTAATATATCCTGTTTTTGAACGTGGGTGTAAGCACCAACCTCACCGGTTccccga                                                                                            >  1:1249551/1‑144 (MQ=255)
                                                      gAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAAACATTGCGAAGATAACCATTTTTTGAACCTGGGGGAACCCACCAACGTCACCCGTTACGCTaat                                                                                          >  1:1556535/1‑147 (MQ=255)
                                                      gAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCTAAGATAGCCTGTTTTTGAACGTGGGTGAACGCACCAACCTCACCGGTTCCGCGAAg                                                                                          >  1:1382978/1‑148 (MQ=255)
                                                               aGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgcgc                                                                                 <  1:1934938/148‑1 (MQ=255)
                                                                 tAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCt                                                                               >  1:2780546/1‑148 (MQ=255)
                                                                  aGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTg                                                                              <  2:3289618/148‑1 (MQ=255)
                                                                   gAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGa                                                                             <  1:2347595/148‑1 (MQ=255)
                                                                    aaGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGAt                                                                            <  1:2161125/148‑1 (MQ=255)
                                                                      ggATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                                                          >  2:1252731/1‑148 (MQ=255)
                                                                       gATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGGTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCaa                                                                         >  1:1585865/1‑148 (MQ=255)
                                                                       gATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCaa                                                                         >  1:2605305/1‑148 (MQ=255)
                                                                         ttAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCGAAGTTCAAGCGCCTGATCAaag                                                                       >  2:2665766/1‑148 (MQ=255)
                                                                          tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                      <  2:1339325/148‑1 (MQ=255)
                                                                          tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaaga                                                                      <  2:1338733/148‑1 (MQ=255)
                                                                            gcgcCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaagaag                                                                    <  2:809311/148‑1 (MQ=255)
                                                                            gcgcCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCGGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCGAAGTTCAAGCGCCGGATCAaagaag                                                                    >  1:1589276/1‑148 (MQ=255)
                                                                               ccgcGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAgaga                                                                 <  2:763983/148‑1 (MQ=255)
                                                                               ccgcGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAgaga                                                                 <  2:1699543/148‑1 (MQ=255)
                                                                                cgcgCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGaa                                                                <  1:2028318/148‑1 (MQ=255)
                                                                                 gcgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGaaa                                                               <  1:436305/148‑1 (MQ=255)
                                                                                   gcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATa                                                             <  2:545681/148‑1 (MQ=255)
                                                                                   gcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATa                                                             <  2:356936/148‑1 (MQ=255)
                                                                                     aaaCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACa                                                           <  1:1013029/148‑1 (MQ=255)
                                                                                     aaaCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACa                                                           <  2:2082067/148‑1 (MQ=255)
                                                                                      aaCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTAAGATATCCTGTTTGTGAACGGGGGTGAACGCACCAACGTCACCGGTTCCGCGAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAg                                                          >  1:1962535/1‑148 (MQ=255)
                                                                                      aaCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAg                                                          >  2:2904344/1‑148 (MQ=255)
                                                                                       aCTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTTGTTGTGAACGTGGGTGTACGCACCAACGTCACCGGTTCCGCTAAGTTCAATCGCGTGATCAAAGTAGAGAAATACAGc                                                         >  2:590470/1‑148 (MQ=255)
                                                                                        cTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTTAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCg                                                        <  1:3102685/148‑1 (MQ=255)
                                                                                        cTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCg                                                        <  1:1850753/148‑1 (MQ=255)
                                                                                              gAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcgc                                                  >  2:1248514/1‑148 (MQ=255)
                                                                                                 aTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTgg                                               <  2:859826/148‑1 (MQ=255)
                                                                                                 aTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCGAAGTTCAAGCGCCTGATCAAAGAAGAGAAATAAAGCGAGGCGCTgg                                               >  2:1509177/1‑148 (MQ=255)
                                                                                                  ttCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTTGCTAAGATAGCCAGTTTTTGAAAGTGGGTGAAAGCACCCACGTCACCCGTTGCGCGAATGTAAAGCGCCAGGTCAAACAAGAGTAAAACACCGAGGCGCTGGa                                              >  1:2413949/1‑148 (MQ=255)
                                                                                                  ttCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGa                                              >  1:3194460/1‑148 (MQ=255)
                                                                                                  ttCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGa                                              >  1:1448928/1‑148 (MQ=255)
                                                                                                  ttCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAAAATTTGCTAAGATAGCCTGTTTTTGTACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAATCGCCGGATCAAATAAGAGAAAAACAGCGAGGCGCGGGa                                              >  1:1363832/1‑148 (MQ=255)
                                                                                                   tCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAATGGCTAAGATAGCCGGTTTGTTAAACTGGGTGAACGCCCCCACCTCCCCCGTTCCGCGAAGATAAAGCGCCCGATAAAATAAGAGAAAAAAACAGAGGCGCTGGAt                                             >  1:1947530/1‑148 (MQ=255)
                                                                                                    cccGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATg                                            <  1:256615/148‑1 (MQ=255)
                                                                                                      cGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTc                                          >  2:2875562/1‑148 (MQ=255)
                                                                                                      cGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTc                                          >  2:3121675/1‑148 (MQ=255)
                                                                                                      cGTAGCCTGCCGTTTGTCCGGCCATGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTTTGAACGTGGGTGAACGCACCAACGTCACCGGTTACGCTAAGTTAAAGCCCCTGATAAAAGAAGAGAAATAAAGCGAGGCGCTGGATGTc                                          >  2:58441/1‑148 (MQ=255)
                                                                                                           ccTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgcgcg                                     >  2:2175022/1‑148 (MQ=255)
                                                                                                                cGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAAc                                >  2:1100434/1‑148 (MQ=255)
                                                                                                                 gTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACa                               <  2:2480153/148‑1 (MQ=255)
                                                                                                                 gTTTGTCCGGCCTGGAGCCGCTGAACATTGGCTAAGAGAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACa                               <  2:476024/148‑1 (MQ=255)
                                                                                                                       ccGGCCTGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTAAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGa                         >  1:724262/1‑148 (MQ=255)
                                                                                                                             tGGAGCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACgg                   >  2:2339936/1‑148 (MQ=255)
                                                                                                                                 gCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGTTGTCGCGCGTCAACAGGTGGAAAACGgcgcg               >  2:587498/1‑148 (MQ=255)
                                                                                                                                 gCCGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgcg               >  2:3076884/1‑148 (MQ=255)
                                                                                                                                  ccGCTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGgcgcgc              >  1:85061/1‑148 (MQ=255)
                                                                                                                                     cTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGa           >  2:3433515/1‑148 (MQ=255)
                                                                                                                                     cTGAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGa           <  1:1708017/148‑1 (MQ=255)
                                                                                                                                       gAACATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGAtt         >  1:1150409/1‑148 (MQ=255)
                                                                                                                                         aaCATTGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTa        <  2:1937214/148‑1 (MQ=255)
                                                                                                                                              tGGCTAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGat   >  1:1873378/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATA  >  CP009273/4214664‑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: