Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,214,915 C→T A387V (GCG→GTG)  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,9150CT99.7% 1420.2 / NA 344A387V (GCG→GTG) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (172/171);  minor base G (1/0);  total (173/171)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 349 total aligned reads displayed.

GTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTC  >  CP009273/4214767‑4215062
                                                                                                                                                    |                                                                                                                                                   
 tAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTcgtgcg                                                                                                                                                  >  2:943967/1‑151 (MQ=255)
  aGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGt                                                                                                                                                 <  2:477941/151‑1 (MQ=255)
  aGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGt                                                                                                                                                 <  2:3126597/151‑1 (MQ=255)
   gCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTc                                                                                                                                                >  1:3962924/1‑151 (MQ=255)
       gCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATACCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGacaaca                                                                                                                                            <  2:3174237/151‑1 (MQ=255)
        ccGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAg                                                                                                                                           >  1:4036874/1‑151 (MQ=255)
         cGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAgg                                                                                                                                          <  1:1562517/151‑1 (MQ=255)
          gTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGt                                                                                                                                         >  1:2570293/1‑151 (MQ=255)
           tttGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTg                                                                                                                                        <  2:456166/151‑1 (MQ=255)
             tGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGa                                                                                                                                      <  1:2060538/151‑1 (MQ=255)
               tCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGaaa                                                                                                                                    <  1:4928203/151‑1 (MQ=255)
                 cGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAAc                                                                                                                                  <  1:567202/151‑1 (MQ=255)
                 cGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAAc                                                                                                                                  <  1:677542/151‑1 (MQ=255)
                 cGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAAc                                                                                                                                  <  2:2465311/151‑1 (MQ=255)
                  ggCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACg                                                                                                                                 <  1:739385/151‑1 (MQ=255)
                     cTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCCCCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGgcg                                                                                                                              <  2:3677107/151‑1 (MQ=255)
                       ggAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGAGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGCAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGgcgcg                                                                                                                            <  2:2024503/151‑1 (MQ=255)
                            cGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGAt                                                                                                                       <  1:1135026/151‑1 (MQ=255)
                             gCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAATGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGagt                                                                                                                      >  1:4502401/1‑149 (MQ=255)
                              cTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTa                                                                                                                     >  1:3141055/1‑151 (MQ=255)
                              cTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTa                                                                                                                     <  1:1890354/151‑1 (MQ=255)
                                 aaCATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCg                                                                                                                  <  1:1240603/151‑1 (MQ=255)
                                   cATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGat                                                                                                                <  1:286195/151‑1 (MQ=255)
                                    aTTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCTCCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGTGGCGCTGGATGTCGTGCGTCAACTGGTGGAAAACGGCGCGCAGATTATCGata                                                                                                               >  2:5183748/1‑151 (MQ=255)
                                     ttGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGatat                                                                                                              <  2:196415/151‑1 (MQ=255)
                                     ttGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGatat                                                                                                              >  1:1307120/1‑151 (MQ=255)
                                     ttGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGatat                                                                                                              >  2:3281053/1‑151 (MQ=255)
                                       ggCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCa                                                                                                            <  1:2795846/151‑1 (MQ=255)
                                       ggCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCa                                                                                                            >  2:4032364/1‑151 (MQ=255)
                                        gCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCaa                                                                                                           >  1:2224245/1‑151 (MQ=255)
                                        gCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCaa                                                                                                           >  2:3506281/1‑151 (MQ=255)
                                           aaGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACAt                                                                                                        >  1:4610712/1‑151 (MQ=255)
                                           aaGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACAt                                                                                                        <  1:999703/151‑1 (MQ=255)
                                            aGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATg                                                                                                       <  1:2974344/151‑1 (MQ=255)
                                            aGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGAACAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATg                                                                                                       <  2:4110374/151‑1 (MQ=255)
                                               tAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGAt                                                                                                    <  2:3829830/151‑1 (MQ=255)
                                                  ccTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACCTGGATGaa                                                                                                 >  2:2773751/1‑151 (MQ=255)
                                                      tttGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGa                                                                                             <  1:169466/151‑1 (MQ=255)
                                                      tttGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGa                                                                                             <  2:173831/151‑1 (MQ=255)
                                                       ttGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGAt                                                                                            <  1:1830049/151‑1 (MQ=255)
                                                           gAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACCTGGATGAAGGGATGCTc                                                                                        >  2:1638297/1‑151 (MQ=255)
                                                               gTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATg                                                                                    <  1:374525/151‑1 (MQ=255)
                                                                tGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGc                                                                                   >  2:2161799/1‑151 (MQ=255)
                                                                       aCGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAgcg                                                                            <  2:4493105/151‑1 (MQ=255)
                                                                           aCCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGa                                                                        >  2:1889518/1‑151 (MQ=255)
                                                                           aCCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGa                                                                        <  1:3828523/151‑1 (MQ=255)
                                                                            ccAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGAt                                                                       <  1:1032300/151‑1 (MQ=255)
                                                                             cAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATg                                                                      >  2:3053901/1‑151 (MQ=255)
                                                                             cAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGAGTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATg                                                                      >  1:3624262/1‑151 (MQ=255)
                                                                              aaCGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATgg                                                                     >  2:988259/1‑151 (MQ=255)
                                                                               aCGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGt                                                                    <  1:2595274/151‑1 (MQ=255)
                                                                                cGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTg                                                                   <  2:4276284/151‑1 (MQ=255)
                                                                                 gTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGc                                                                  >  2:2003633/1‑151 (MQ=255)
                                                                                  tCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCg                                                                 <  2:2066740/151‑1 (MQ=255)
                                                                                    aCCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGtt                                                               >  1:5232026/1‑151 (MQ=255)
                                                                                      cGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGtttt                                                             >  1:939375/1‑151 (MQ=255)
                                                                                          tCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCa                                                         >  2:1159550/1‑151 (MQ=255)
                                                                                          tCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCa                                                         <  2:4991207/151‑1 (MQ=255)
                                                                                          tCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCa                                                         >  1:5201280/1‑151 (MQ=255)
                                                                                           ccGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCaa                                                        <  2:1607160/151‑1 (MQ=255)
                                                                                           ccGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCaa                                                        >  2:2292275/1‑151 (MQ=255)
                                                                                               tAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTg                                                    >  2:1224790/1‑151 (MQ=255)
                                                                                                  gTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGAtt                                                 >  1:4317296/1‑151 (MQ=255)
                                                                                                   ttCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTg                                                >  1:2283267/1‑151 (MQ=255)
                                                                                                   ttCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTg                                                >  2:4663175/1‑151 (MQ=255)
                                                                                                    tCAAGCGCCTGATCAAAGTAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGc                                               >  1:5278177/1‑151 (MQ=255)
                                                                                                    tCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGc                                               >  1:86035/1‑151 (MQ=255)
                                                                                                     cAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGcc                                              <  1:3350944/151‑1 (MQ=255)
                                                                                                        gcgcCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGt                                           <  2:152189/151‑1 (MQ=255)
                                                                                                        gcgcCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGt                                           <  2:4206151/151‑1 (MQ=255)
                                                                                                            cTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAAc                                       >  2:3411601/1‑151 (MQ=255)
                                                                                                            cTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTATTCTCAATCTGATTGCCGGTGAAc                                       >  2:52659/1‑151 (MQ=255)
                                                                                                               aTCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATCTCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCgg                                    >  2:5100681/1‑151 (MQ=255)
                                                                                                               aTCAAAGAAGAGAAAGACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCgg                                    <  2:1375032/151‑1 (MQ=255)
                                                                                                                 cAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGat                                  >  1:2637525/1‑151 (MQ=255)
                                                                                                                   aagaagAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGatat                                <  1:502428/151‑1 (MQ=255)
                                                                                                                   aagaagAGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGatat                                <  1:1508858/151‑1 (MQ=255)
                                                                                                                    agaagaGCAAGACAGCGCGGCGCTGGATGTCGTGCGTCAACCGGTGGAAAACGGCGCGCAGATTATCGATATCTACATGGATGAAGGGATGCTAGATGCCGCAGCGGCTATTGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGCATAtc                               <  2:845535/151‑1 (MQ=255)
                                                                                                                    agaagaGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtc                               <  2:3964026/151‑1 (MQ=255)
                                                                                                                     gaagaGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcg                              <  2:617009/151‑1 (MQ=255)
                                                                                                                     gaagaGAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcg                              <  1:1788252/151‑1 (MQ=255)
                                                                                                                       agagAAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcgct                            <  1:4852683/151‑1 (MQ=255)
                                                                                                                        gagaAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcgctc                           >  2:4876648/1‑151 (MQ=255)
                                                                                                                        gagaAATACAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcgctc                           <  1:1591429/151‑1 (MQ=255)
                                                                                                                         agaAATACAGCGAGGCGCTGGATGTCGTGCGTCCACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACCTGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATAtcgctcg                          <  2:899901/151‑1 (MQ=255)
                                                                                                                               aCAGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATCTCGCTCGCGTGcc                    >  1:2356882/1‑151 (MQ=255)
                                                                                                                                 aGCGAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGTCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCCATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGa                  >  1:3549986/1‑151 (MQ=255)
                                                                                                                                 aGCGAGGCGCTGGATGTCGTGCGTCAACAGGAGGAACACGGCGCGCAGATTATCGATATCAACATGGATGAAGTGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCAGATTGCCGGTGAACCTGATATAGCTCGCGTGCCGa                  >  1:3240247/1‑151 (MQ=255)
                                                                                                                                    gAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACCTGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTa               <  2:1777550/151‑1 (MQ=255)
                                                                                                                                    gAGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTa               >  1:3085265/1‑151 (MQ=255)
                                                                                                                                     aGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTAt              <  1:4999761/151‑1 (MQ=255)
                                                                                                                                     aGGCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTAt              >  2:2717419/1‑151 (MQ=255)
                                                                                                                                      ggCGCTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGAGTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATg             <  2:3447606/151‑1 (MQ=255)
                                                                                                                                       gcgcTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGa            >  1:1121195/1‑151 (MQ=255)
                                                                                                                                         gcTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATc          >  2:1319861/1‑151 (MQ=255)
                                                                                                                                          cTGGATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCg         <  1:5131040/151‑1 (MQ=255)
                                                                                                                                             gATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGAct      >  2:768517/1‑151 (MQ=255)
                                                                                                                                             gATGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACAGGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGAct      >  2:5076060/1‑151 (MQ=255)
                                                                                                                                             gATGCCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGAct      <  2:259318/151‑1 (MQ=255)
                                                                                                                                               tGTCGTGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCACCATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGActcc    <  2:3255764/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGTCGCGCGTCAACAGGTGGAAAACGGCGCGCAGATTATCGATATCAACATGGATGAAGGGATGCTCGATGCCGAAGCGGCGATGGTGCGTTTTCTCAATCTGATTGCCGGTGAACCGGATATCGCTCGCGTGCCGATTATGATCGACTCCTC  >  CP009273/4214767‑4215062

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