Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,040,655 C→A 6.8% A582A (GCG→GCT gcvP ← glycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,040,6550CA 6.8% 806.0 / 22.0 218A582A (GCG→GCTgcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (104/99);  new base A (8/7);  total (112/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 220 total aligned reads displayed.

ATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTG  >  CP009273/3040508‑3040801
                                                                                                                                                   |                                                                                                                                                  
aTCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATc                                                                                                                                                    >  2:3130869/1‑148 (MQ=255)
aTCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATc                                                                                                                                                    >  2:2328965/1‑148 (MQ=255)
 tCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCg                                                                                                                                                   >  2:2509358/1‑148 (MQ=255)
 tCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCg                                                                                                                                                   >  2:726965/1‑148 (MQ=255)
  cGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGc                                                                                                                                                  <  2:3105199/148‑1 (MQ=255)
  cGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTACTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGc                                                                                                                                                  >  1:1978306/1‑148 (MQ=255)
  cGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCAGTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGc                                                                                                                                                  <  2:1615599/148‑1 (MQ=255)
   gATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGcc                                                                                                                                                 <  1:1395522/148‑1 (MQ=255)
    aTGTTGCCGTTTTTATCACACGCCACAACCACCACCAGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAATCCGGGATCAGGCAGATATCGCGATGCCCGTCGTTGCGGCTTTCATGATAATGACGAATCGCca                                                                                                                                                >  1:2200398/1‑148 (MQ=255)
        tGCCGTTTTTATCACACCCCACAACCACCACCAGCAATCCTGCCAGATGTTCAGAAGCGGGGTTAGTTACGTGCGCAGAAACCGGGATCAGGCAGATATCGCGAATCCCTTCGTTGCGGCTTTCATGATAAAGACCAAACGCcaccag                                                                                                                                            >  1:746393/1‑148 (MQ=255)
         gCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAgg                                                                                                                                           <  1:939378/148‑1 (MQ=255)
         gCCGTTTTTATCACACGCCACAACCACCACCTGCAGTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATAGCCAGCAgg                                                                                                                                           <  2:2211271/148‑1 (MQ=255)
             tttttATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCcg                                                                                                                                       >  2:1454494/1‑148 (MQ=255)
             tttttATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCcg                                                                                                                                       >  2:1454354/1‑148 (MQ=255)
                ttATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGt                                                                                                                                    >  1:1067644/1‑148 (MQ=255)
                ttATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGt                                                                                                                                    >  1:2842757/1‑148 (MQ=255)
                   tCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGGTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTAtt                                                                                                                                 <  2:329024/148‑1 (MQ=255)
                       acGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGACTCGCCAGCAGGCCCGCGTATTCGcc                                                                                                                             >  2:2799523/1‑148 (MQ=255)
                        cGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGccc                                                                                                                            >  2:3237121/1‑148 (MQ=255)
                          ccACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCtg                                                                                                                          >  2:1776811/1‑148 (MQ=255)
                          ccACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCtg                                                                                                                          <  2:1928789/148‑1 (MQ=255)
                              aaCCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTgcg                                                                                                                      >  2:1488429/1‑148 (MQ=255)
                                    caccTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCagag                                                                                                                <  1:1818331/148‑1 (MQ=255)
                                      ccTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGtt                                                                                                              <  1:1924095/148‑1 (MQ=255)
                                      ccTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGtt                                                                                                              <  1:630231/148‑1 (MQ=255)
                                      ccTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGtt                                                                                                              >  1:2492458/1‑148 (MQ=255)
                                        tGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCg                                                                                                            >  2:352373/1‑148 (MQ=255)
                                         gCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCgg                                                                                                           >  2:1002302/1‑148 (MQ=255)
                                         gCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCgg                                                                                                           <  2:713088/148‑1 (MQ=255)
                                             tCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGc                                                                                                       >  1:3130005/1‑148 (MQ=255)
                                             tCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGc                                                                                                       >  1:1419237/1‑148 (MQ=255)
                                             tCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGc                                                                                                       >  1:1121137/1‑148 (MQ=255)
                                               cTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCAt                                                                                                     <  1:3236126/148‑1 (MQ=255)
                                                  ccATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACa                                                                                                  >  1:2727923/1‑148 (MQ=255)
                                                  ccATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACa                                                                                                  >  1:2278952/1‑148 (MQ=255)
                                                   cATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACaa                                                                                                 <  1:3000871/148‑1 (MQ=255)
                                                         tgCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGcg                                                                                           >  2:238663/1‑148 (MQ=255)
                                                              aaGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTa                                                                                      >  2:1384929/1‑148 (MQ=255)
                                                               aGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTaa                                                                                     <  1:1434983/148‑1 (MQ=255)
                                                               aGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTaa                                                                                     <  1:482390/148‑1 (MQ=255)
                                                               aGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTaa                                                                                     <  1:2035476/148‑1 (MQ=255)
                                                               aGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCGGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTaa                                                                                     <  1:239596/148‑1 (MQ=255)
                                                                gCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAAc                                                                                    <  2:1415216/148‑1 (MQ=255)
                                                                 cGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAAcc                                                                                   >  2:3214129/1‑148 (MQ=255)
                                                                 cGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAAcc                                                                                   >  2:1550808/1‑148 (MQ=255)
                                                                  ggggTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATAGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCg                                                                                  >  2:2809662/1‑148 (MQ=255)
                                                                   gggTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCgg                                                                                 >  1:2941287/1‑148 (MQ=255)
                                                                    ggTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGt                                                                                >  2:1242883/1‑148 (MQ=255)
                                                                     gTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTc                                                                               >  1:1503891/1‑148 (MQ=255)
                                                                     gTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTc                                                                               <  1:563899/148‑1 (MQ=255)
                                                                     gTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTc                                                                               <  2:1400202/148‑1 (MQ=255)
                                                                      ttAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCa                                                                              <  2:1985990/148‑1 (MQ=255)
                                                                      ttAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCa                                                                              <  1:2758778/148‑1 (MQ=255)
                                                                       tAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAg                                                                             <  2:1748414/148‑1 (MQ=255)
                                                                          ttCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGttt                                                                          <  1:3475841/148‑1 (MQ=255)
                                                                            ccGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCa                                                                        <  2:412437/148‑1 (MQ=255)
                                                                                     gAAGCCGGGATCAGGCAGATCTCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGt                                                               >  1:2863662/1‑148 (MQ=255)
                                                                                     gAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGt                                                               <  2:2448816/148‑1 (MQ=255)
                                                                                       aGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTca                                                             <  2:3447919/148‑1 (MQ=255)
                                                                                       aGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTca                                                             >  1:3216820/1‑148 (MQ=255)
                                                                                       aGCCGGGATCAGGCAGAGATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTACCCGGTCAGTGTCACCAGCCAGTca                                                             >  2:1854312/1‑148 (MQ=255)
                                                                                            ggATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTcagccag                                                        <  1:3287424/148‑1 (MQ=255)
                                                                                               tCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTg                                                     >  2:592841/1‑148 (MQ=255)
                                                                                                cAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTACCCGGTCAGTTTCACCAGCCAGTCAGCCAGCTgc                                                    >  2:3331102/1‑148 (MQ=255)
                                                                                                   gCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCa                                                 <  2:2737730/148‑1 (MQ=255)
                                                                                                       atatCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAatca                                             >  2:1089444/1‑148 (MQ=255)
                                                                                                       atatCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATAGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAatca                                             >  2:3100144/1‑148 (MQ=255)
                                                                                                           cgcgATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATctg                                         <  1:2828378/148‑1 (MQ=255)
                                                                                                           cgcgATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATctg                                         <  1:858034/148‑1 (MQ=255)
                                                                                                            gcgATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATctgc                                        <  1:1758727/148‑1 (MQ=255)
                                                                                                            gcgATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATctgc                                        <  1:1550808/148‑1 (MQ=255)
                                                                                                             cgATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATctgct                                       >  2:3078464/1‑148 (MQ=255)
                                                                                                               aTGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGa                                     <  1:1786767/148‑1 (MQ=255)
                                                                                                               aTGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGa                                     >  1:109399/1‑148 (MQ=255)
                                                                                                               aTGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGa                                     <  2:215046/148‑1 (MQ=255)
                                                                                                                  cccTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATaa                                  <  2:2164072/148‑1 (MQ=255)
                                                                                                                  cccTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATaa                                  <  1:2058732/148‑1 (MQ=255)
                                                                                                                  cccTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATaa                                  <  2:2164140/148‑1 (MQ=255)
                                                                                                                    cTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAAcc                                <  1:2499625/148‑1 (MQ=255)
                                                                                                                     ttCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCt                               >  2:1140709/1‑148 (MQ=255)
                                                                                                                     ttCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCt                               >  2:1141529/1‑148 (MQ=255)
                                                                                                                      tCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCtt                              >  1:374631/1‑148 (MQ=255)
                                                                                                                           gCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGcc                         <  2:2593888/148‑1 (MQ=255)
                                                                                                                              gCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGc                      >  1:1090667/1‑148 (MQ=255)
                                                                                                                              gCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGc                      >  1:3322336/1‑148 (MQ=255)
                                                                                                                              gCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGc                      >  1:1773752/1‑148 (MQ=255)
                                                                                                                               cTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCt                     <  1:814933/148‑1 (MQ=255)
                                                                                                                               cTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCt                     <  2:2544837/148‑1 (MQ=255)
                                                                                                                                 ttCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTcc                   >  2:1236429/1‑148 (MQ=255)
                                                                                                                                 ttCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTcc                   >  1:2443282/1‑148 (MQ=255)
                                                                                                                                   cATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCcgg                 >  1:429099/1‑148 (MQ=255)
                                                                                                                                     tGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCcggcg               <  1:1701752/148‑1 (MQ=255)
                                                                                                                                     tGATAATGACGAATAGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCcggcg               <  1:1002942/148‑1 (MQ=255)
                                                                                                                                            gACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGaa        >  1:3325930/1‑148 (MQ=255)
                                                                                                                                             aCGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAAc       >  2:1134120/1‑148 (MQ=255)
                                                                                                                                               gAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACgg     <  2:903718/148‑1 (MQ=255)
                                                                                                                                                aaTCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACggg    <  2:2842828/148‑1 (MQ=255)
                                                                                                                                                aaTCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACggg    <  2:1067644/148‑1 (MQ=255)
                                                                                                                                                aaTCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACggg    <  1:318552/148‑1 (MQ=255)
                                                                                                                                                  tCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTg  >  2:2931007/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCGTTGCGGCTTTCATGATAATGACGAATCGCCAGCAGGCCCGCGTATTCGCCCTGTGCGCCAGAGTTCGGCTGCATACAAACGGCGTCGTAACCGGTCAGTTTCACCAGCCAGTCAGCCAGCTGCGCAATCATCTGCTGATAACCTTCGGCCTGCTCCGGCGGGCAGAACGGGTG  >  CP009273/3040508‑3040801

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: