Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,040,479 T→G 13.1% Q641P (CAG→CCG)  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,4790TG13.1% 753.4 / 120.8 328Q641P (CAG→CCG) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base T (148/134);  new base G (18/25);  total (167/161)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 330 total aligned reads displayed.

CCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCG  >  CP009273/3040332‑3040628
                                                                                                                                                   |                                                                                                                                                     
ccAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGAACGGTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATCCAGGAGAGGTTATCACCCGCCTGt                                                                                                                                                     >  2:2175893/1‑150 (MQ=255)
       gggCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCtt                                                                                                                                              <  1:3254899/150‑1 (MQ=255)
         gCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTc                                                                                                                                            <  1:2447257/150‑1 (MQ=255)
          cGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTcg                                                                                                                                           >  1:564764/1‑150 (MQ=255)
          cGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTAcg                                                                                                                                           >  2:2544169/1‑150 (MQ=255)
           gTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTcgc                                                                                                                                          >  2:593180/1‑150 (MQ=255)
            ttCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTcgcg                                                                                                                                         >  1:2787194/1‑150 (MQ=255)
              cATGTTCTCGCAATCAAGTTAAACCTGACCGCAGAACTGATGCACGACTTCACACACTTCACGGATCGTTTATTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTcgcgcg                                                                                                                                       <  1:4185652/150‑1 (MQ=255)
               aTGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCgcgcgc                                                                                                                                      >  2:3698357/1‑150 (MQ=255)
                   tCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGAt                                                                                                                                  <  2:1582950/150‑1 (MQ=255)
                    cgcgCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAtc                                                                                                                                 <  1:597288/150‑1 (MQ=255)
                     gcgcCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAtca                                                                                                                                <  1:3129315/150‑1 (MQ=255)
                       gcCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAtcagt                                                                                                                              <  2:1771598/150‑1 (MQ=255)
                        ccATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAtcagtc                                                                                                                             >  2:1976048/1‑150 (MQ=255)
                         cATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAtcagtca                                                                                                                            >  1:2511458/1‑150 (MQ=255)
                           tCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGa                                                                                                                          <  2:4996481/150‑1 (MQ=255)
                            cacggTAAACCGGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATGGTTTCTTCATAGACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATAAGTCAGAt                                                                                                                         <  1:4983684/147‑1 (MQ=255)
                            cAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGAt                                                                                                                         <  2:4038309/150‑1 (MQ=255)
                             aaGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATc                                                                                                                        <  1:3818649/150‑1 (MQ=255)
                               ggTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCTCGCGCAGATCAGTCAGATCGa                                                                                                                      >  1:4497275/1‑150 (MQ=255)
                                  aaaCCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGt                                                                                                                   >  1:3370783/1‑150 (MQ=255)
                                  aaaCCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGt                                                                                                                   >  2:3224723/1‑150 (MQ=255)
                                   aaCCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGtt                                                                                                                  >  2:1202811/1‑150 (MQ=255)
                                     ccTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGc                                                                                                                <  2:3831587/150‑1 (MQ=255)
                                     ccTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGc                                                                                                                >  1:355573/1‑150 (MQ=255)
                                      cTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGcc                                                                                                               <  1:1098673/150‑1 (MQ=255)
                                          ccgccgAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGttt                                                                                                           <  2:3117575/150‑1 (MQ=255)
                                           cgccgAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGtttt                                                                                                          <  1:1842696/150‑1 (MQ=255)
                                           cgccgAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGtttt                                                                                                          <  2:4114424/150‑1 (MQ=255)
                                               gAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCCGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATc                                                                                                      <  1:4699906/150‑1 (MQ=255)
                                                atcTGATGCTCGAATTCACACACTTCACTGTTCTTTTCTTCATACTCGTCGTGGGTAGTAGGATAAGTTACGATGATTCTGGAGATGTTATAGCCGGCCTGTTGCGCTTTCGCGCGATGATCAGTGTGATCGATGTTGCCGTTTTTATca                                                                                                     <  2:5155561/148‑1 (MQ=255)
                                                 aCTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATcac                                                                                                    <  2:1239020/150‑1 (MQ=255)
                                                  cTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATcaca                                                                                                   <  2:74682/150‑1 (MQ=255)
                                                   tGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGAGCAGTCAGATCGATGTTGCCGTTTTTATCacac                                                                                                  >  1:292011/1‑150 (MQ=255)
                                                     aTGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGc                                                                                                >  2:4704191/1‑150 (MQ=255)
                                                       gCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCca                                                                                              >  1:4779394/1‑150 (MQ=255)
                                                        cACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCcac                                                                                             <  1:4345030/150‑1 (MQ=255)
                                                           gACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAac                                                                                          >  1:1209966/1‑150 (MQ=255)
                                                            aCTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAacc                                                                                         <  2:3021254/150‑1 (MQ=255)
                                                                 acacacTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAaccaccac                                                                                    >  1:3746734/1‑150 (MQ=255)
                                                                  cacacTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAaccaccacc                                                                                   <  2:5239925/150‑1 (MQ=255)
                                                                   acacTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATAGCCCGCCACTACCGCTATCGCGCGCAGATCAGACAGATCGATTTTGCCGTTATTATCACACGCCACAaccaccacca                                                                                  >  1:4990355/1‑149 (MQ=255)
                                                                      cTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCa                                                                               >  2:2211779/1‑150 (MQ=255)
                                                                        tCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCAtt                                                                             >  2:641158/1‑150 (MQ=255)
                                                                           cGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCt                                                                          >  2:2826836/1‑150 (MQ=255)
                                                                           cGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCt                                                                          <  2:1340753/150‑1 (MQ=255)
                                                                            ggATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTg                                                                         >  1:1586275/1‑150 (MQ=255)
                                                                              aTCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGcc                                                                       >  1:5075824/1‑150 (MQ=255)
                                                                                cGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCat                                                                     >  1:2680711/1‑150 (MQ=255)
                                                                                     cttcATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATAtgtg                                                                <  2:367331/150‑1 (MQ=255)
                                                                                         aTACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCAGAATTCCTGCCATATGTGCAGa                                                            >  1:2471043/1‑150 (MQ=255)
                                                                                             acGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCg                                                        <  1:5100982/150‑1 (MQ=255)
                                                                                               gCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCggg                                                      >  1:2992967/1‑150 (MQ=255)
                                                                                               gCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCggg                                                      <  1:2872476/150‑1 (MQ=255)
                                                                                                  gTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGtt                                                   >  1:221948/1‑150 (MQ=255)
                                                                                                   tGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTa                                                  >  2:3736259/1‑150 (MQ=255)
                                                                                                   tGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTa                                                  <  1:3604562/150‑1 (MQ=255)
                                                                                                    gggTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAg                                                 >  1:3961012/1‑150 (MQ=255)
                                                                                                    gggTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAg                                                 <  2:1473617/150‑1 (MQ=255)
                                                                                                       tAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTc                                              >  1:1182986/1‑150 (MQ=255)
                                                                                                        aGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTcc                                             >  2:4638678/1‑150 (MQ=255)
                                                                                                        aGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTcc                                             >  2:4899631/1‑150 (MQ=255)
                                                                                                          aaGGATAAGTCACCCTGATACAGGAGAGGGTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATGAGGCCGATCGATGTTGCCGTTTTTATCACCCGCCACAACCACCCCCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGt                                           <  2:1534266/150‑1 (MQ=255)
                                                                                                          aaGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCAGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGt                                           >  1:4913815/1‑150 (MQ=255)
                                                                                                           aGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTg                                          <  2:4446296/150‑1 (MQ=255)
                                                                                                            ggATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTgc                                         >  2:2047481/1‑150 (MQ=255)
                                                                                                            ggATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTgc                                         <  1:86881/150‑1 (MQ=255)
                                                                                                              aTAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTgcgc                                       >  1:3329797/1‑150 (MQ=255)
                                                                                                                aaGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAg                                     >  1:1868231/1‑150 (MQ=255)
                                                                                                                aaGTCACCATGATACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAg                                     >  2:4170232/1‑150 (MQ=255)
                                                                                                                 aGTCACCATGATACAGGAGAGGTTATCGCCCGCCGGTTCTGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGa                                    <  1:127060/150‑1 (MQ=255)
                                                                                                                     aCCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGcc                                <  2:3541897/150‑1 (MQ=255)
                                                                                                                     aCCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGcc                                >  1:4881405/1‑150 (MQ=255)
                                                                                                                      ccATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCg                               >  1:4570151/1‑150 (MQ=255)
                                                                                                                      ccATGATACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCg                               >  1:150991/1‑150 (MQ=255)
                                                                                                                       cATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCgg                              <  1:4073720/150‑1 (MQ=255)
                                                                                                                           aTACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATc                          <  1:2485174/150‑1 (MQ=255)
                                                                                                                           aTACAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATc                          <  1:636672/150‑1 (MQ=255)
                                                                                                                              cAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGATTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAgg                       <  2:1675960/150‑1 (MQ=255)
                                                                                                                              cAGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAgg                       <  2:3949352/150‑1 (MQ=255)
                                                                                                                               aGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGc                      >  2:3490781/1‑150 (MQ=255)
                                                                                                                               aGGAGAGGTTATCGCCCGCCGGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCAATCCTGCCATATGTGCATAAGCGGGGTTAGATCCGTGAGCAGAAGCCGGGATCAGGc                      >  2:231839/1‑150 (MQ=255)
                                                                                                                                ggAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCa                     <  2:4667803/150‑1 (MQ=255)
                                                                                                                                ggAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCa                     >  1:2422163/1‑150 (MQ=255)
                                                                                                                                ggAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCa                     >  2:2469353/1‑150 (MQ=255)
                                                                                                                                 gagaGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAg                    <  2:4779394/150‑1 (MQ=255)
                                                                                                                                    agGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGata                 >  1:5199059/1‑150 (MQ=255)
                                                                                                                                     gggtATCGCCCGCCGGGTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACCCCCACCACCTGCATGCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGatat                <  2:409011/147‑1 (MQ=255)
                                                                                                                                     ggTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGatat                <  2:3327391/150‑1 (MQ=255)
                                                                                                                                     ggTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGatat                <  1:5160069/150‑1 (MQ=255)
                                                                                                                                      gTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATc               >  1:1022391/1‑150 (MQ=255)
                                                                                                                                        tATCGCCCGCCTGTTCCGCATTCGCGCGCAGATAAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCCTGCCTGCCATATGTGCAGAAGCGGGGGTAGTTCAGTGCGCAGAAGCCGGGATCAGGCAGTTATcgc             >  2:3905804/1‑150 (MQ=255)
                                                                                                                                         aTCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACCCGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATcgcg            <  2:4035877/150‑1 (MQ=255)
                                                                                                                                         aTCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATcgcg            <  1:3101764/150‑1 (MQ=255)
                                                                                                                                           cgcccgccTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGAt          <  1:3068113/150‑1 (MQ=255)
                                                                                                                                            gcccgccTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATg         <  2:4316849/150‑1 (MQ=255)
                                                                                                                                                  cTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTc   >  2:4056799/1‑150 (MQ=255)
                                                                                                                                                  cTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTc   <  2:2957205/150‑1 (MQ=255)
                                                                                                                                                  cTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGACATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTc   <  1:335875/150‑1 (MQ=255)
                                                                                                                                                   tGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCg  <  2:4536093/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAACCTGACCGCCGAACTGATGCACGACTTCACACACTTCACGGATCGTTTCTTCATACACGCCGTGGGTAGAAGGATAAGTCACCATGATACAGGAGAGGTTATCGCCCGCCTGTTCCGCTTTCGCGCGCAGATCAGTCAGATCGATGTTGCCGTTTTTATCACACGCCACAACCACCACCTGCATTCCTGCCATATGTGCAGAAGCGGGGTTAGTTCCGTGCGCAGAAGCCGGGATCAGGCAGATATCGCGATGCCCTTCG  >  CP009273/3040332‑3040628

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