Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,041,026 C→A 16.7% D459Y (GAC→TAC)  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,041,0260CA16.7% 959.8 / 121.7 306D459Y (GAC→TAC) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (149/106);  new base A (27/24);  total (176/130)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 310 total aligned reads displayed.

CGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGAT  >  CP009273/3040877‑3041173
                                                                                                                                                     |                                                                                                                                                   
cGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTc                                                                                                                                                     <  1:3968937/150‑1 (MQ=255)
cGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTc                                                                                                                                                     <  2:28844/150‑1 (MQ=255)
      catcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGc                                                                                                                                               >  2:1748417/1‑150 (MQ=255)
      catcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGc                                                                                                                                               >  2:563125/1‑150 (MQ=255)
       atcGGCTGATTATGCGCCAGTTCTTTTCGCTCCATCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGcc                                                                                                                                              <  2:4105893/150‑1 (MQ=255)
        tcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCa                                                                                                                                             >  2:2073375/1‑150 (MQ=255)
         cGCCTGATACAGCGCCAGATCTTTACGCACCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGACGCGCAGCATCGCAGGCTGGATAGAGCTGCTGTCGTGAGCCAc                                                                                                                                            >  1:1538116/1‑150 (MQ=255)
           ccTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGt                                                                                                                                          >  2:3915004/1‑150 (MQ=255)
           ccTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGt                                                                                                                                          <  1:3166043/150‑1 (MQ=255)
           ccTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGt                                                                                                                                          >  2:956223/1‑150 (MQ=255)
           ccTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCGGGATAGAGCGGCTGTAGTGAGCCACGt                                                                                                                                          >  2:4050497/1‑150 (MQ=255)
               aTTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCttt                                                                                                                                      <  2:3133572/150‑1 (MQ=255)
                 tCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCGGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGt                                                                                                                                    >  2:1816367/1‑150 (MQ=255)
                 tCAGCGCCAGATCTTTACGCTCCAGAGAGTTCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGTATAGAGCGGCTGTAGTGAGCCACGTCTTTGt                                                                                                                                    <  2:4462449/150‑1 (MQ=255)
                   aGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTcc                                                                                                                                  >  2:2968888/1‑150 (MQ=255)
                         aGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCgtg                                                                                                                            >  1:2347774/1‑150 (MQ=255)
                           aTCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTc                                                                                                                          <  1:2620402/150‑1 (MQ=255)
                               ttACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATg                                                                                                                      >  2:1205364/1‑150 (MQ=255)
                               ttACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATg                                                                                                                      <  2:832012/150‑1 (MQ=255)
                                  cGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTcc                                                                                                                   >  2:337764/1‑150 (MQ=255)
                                    cTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAg                                                                                                                 >  1:4396218/1‑150 (MQ=255)
                                       cAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGcc                                                                                                              >  1:2689539/1‑150 (MQ=255)
                                        aGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCg                                                                                                             >  1:3817716/1‑150 (MQ=255)
                                        aGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCg                                                                                                             <  2:3526852/150‑1 (MQ=255)
                                         gCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGt                                                                                                            >  2:1577564/1‑150 (MQ=255)
                                           gAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTgg                                                                                                          <  2:2924232/150‑1 (MQ=255)
                                             gTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGtt                                                                                                        <  2:2631064/150‑1 (MQ=255)
                                                    taGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGccc                                                                                                 >  2:3793333/1‑150 (MQ=255)
                                                     aGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCca                                                                                                >  2:3315242/1‑150 (MQ=255)
                                                      gcgcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcag                                                                                               <  2:2404544/150‑1 (MQ=255)
                                                       cgcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagc                                                                                              >  1:5075/1‑150 (MQ=255)
                                                         catcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagcag                                                                                            <  2:4593584/150‑1 (MQ=255)
                                                          atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                           <  2:4500381/150‑1 (MQ=255)
                                                            catTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCAc                                                                                         >  1:3608823/1‑150 (MQ=255)
                                                            catTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCAc                                                                                         <  1:756998/150‑1 (MQ=255)
                                                                  ggTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGaa                                                                                   >  1:34338/1‑150 (MQ=255)
                                                                  ggTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGaa                                                                                   >  1:2793634/1‑150 (MQ=255)
                                                                   gTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGaaa                                                                                  >  2:2523940/1‑150 (MQ=255)
                                                                      tCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGc                                                                               >  1:736887/1‑150 (MQ=255)
                                                                          tgtgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCa                                                                           >  1:2076796/1‑150 (MQ=255)
                                                                          tgtgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCa                                                                           >  1:2076962/1‑150 (MQ=255)
                                                                          tgtgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCa                                                                           >  1:2301912/1‑150 (MQ=255)
                                                                           gtgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCAt                                                                          >  1:2868212/1‑150 (MQ=255)
                                                                            tgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCAtt                                                                         <  2:4262936/150‑1 (MQ=255)
                                                                              gTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTAc                                                                       <  1:373599/150‑1 (MQ=255)
                                                                               tAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTGCg                                                                      >  1:4248140/1‑150 (MQ=255)
                                                                                aGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGt                                                                     >  2:2150747/1‑150 (MQ=255)
                                                                                aGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGt                                                                     >  1:1862063/1‑150 (MQ=255)
                                                                                aGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGt                                                                     >  2:1458536/1‑150 (MQ=255)
                                                                                aGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGt                                                                     >  1:1992023/1‑150 (MQ=255)
                                                                                   gATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGtttt                                                                  >  2:1880513/1‑150 (MQ=255)
                                                                                    aTTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTc                                                                 >  2:1905179/1‑150 (MQ=255)
                                                                                     ttAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCa                                                                >  1:2576443/1‑150 (MQ=255)
                                                                                     ttAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCa                                                                <  2:1970715/150‑1 (MQ=255)
                                                                                      tAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAc                                                               <  1:2939341/150‑1 (MQ=255)
                                                                                      tAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAc                                                               >  1:1446163/1‑150 (MQ=255)
                                                                                      tAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAc                                                               <  1:1012714/150‑1 (MQ=255)
                                                                                       aaaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCTATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcg                                                              >  2:4317752/1‑150 (MQ=255)
                                                                                       aaaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcg                                                              <  1:3664716/150‑1 (MQ=255)
                                                                                       aaaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcg                                                              >  2:233359/1‑150 (MQ=255)
                                                                                       aaaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGACTTTGTCCAGCGTGTCGATGTCCAGGCCGAGGTTATCTCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcg                                                              >  2:2673984/1‑150 (MQ=255)
                                                                                          cacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGt                                                           <  2:2763939/150‑1 (MQ=255)
                                                                                          cacCGGATGGGTGAGGATTTCGTCGTCGCGCAACATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATAACGTTTTCACGCGt                                                           >  1:4063958/1‑150 (MQ=255)
                                                                                            ccGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTgg                                                         <  1:3320138/150‑1 (MQ=255)
                                                                                              ggATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgtt                                                       <  1:112455/150‑1 (MQ=255)
                                                                                               gATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAGAAGCTGCATTACGTTTTCACGCGTGgttg                                                      <  2:1759859/150‑1 (MQ=255)
                                                                                               gATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttg                                                      <  2:72098/150‑1 (MQ=255)
                                                                                                 tGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgtt                                                    >  1:2468790/1‑150 (MQ=255)
                                                                                                 tGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgtt                                                    >  1:3068693/1‑150 (MQ=255)
                                                                                                 tGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgtt                                                    <  2:1001326/150‑1 (MQ=255)
                                                                                                   ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                  <  2:1059056/150‑1 (MQ=255)
                                                                                                   ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                  >  2:3289286/1‑150 (MQ=255)
                                                                                                             tcgtcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTg                                        >  1:1922124/1‑150 (MQ=255)
                                                                                                               gtcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGAt                                      >  1:22910/1‑150 (MQ=255)
                                                                                                               gtcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGAt                                      <  2:3604896/150‑1 (MQ=255)
                                                                                                               gtcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGAt                                      <  2:45590/150‑1 (MQ=255)
                                                                                                                 cgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATcc                                    >  1:40983/1‑150 (MQ=255)
                                                                                                                 cgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATcc                                    >  2:3516524/1‑150 (MQ=255)
                                                                                                                     gcgcAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAc                                <  1:4503988/150‑1 (MQ=255)
                                                                                                                          gcaTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGt                           >  1:4550081/1‑150 (MQ=255)
                                                                                                                            aTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTc                         <  1:1761046/150‑1 (MQ=255)
                                                                                                                             tCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCa                        >  2:1427599/1‑150 (MQ=255)
                                                                                                                              cGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAg                       <  1:4363175/150‑1 (MQ=255)
                                                                                                                                 aGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAat                    >  1:3489235/1‑150 (MQ=255)
                                                                                                                                    cTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATc                 <  2:4477931/150‑1 (MQ=255)
                                                                                                                                     tGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCg                >  1:1583918/1‑150 (MQ=255)
                                                                                                                                     tGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCg                >  2:279220/1‑150 (MQ=255)
                                                                                                                                     tGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCg                >  1:1605206/1‑150 (MQ=255)
                                                                                                                                       gATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCACCAAGGGTGATCCCAACCGCGTTCAGAATATCGCt              <  2:3422988/150‑1 (MQ=255)
                                                                                                                                          agagCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACg           >  2:492025/1‑150 (MQ=255)
                                                                                                                                          agagCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACg           >  2:147878/1‑150 (MQ=255)
                                                                                                                                            agCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCa         >  2:3954944/1‑150 (MQ=255)
                                                                                                                                            agCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCa         <  2:131742/150‑1 (MQ=255)
                                                                                                                                             gCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAg        <  1:2885380/150‑1 (MQ=255)
                                                                                                                                               ggCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCGTCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGt      <  2:2312180/150‑1 (MQ=255)
                                                                                                                                               ggCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGt      >  2:2212342/1‑150 (MQ=255)
                                                                                                                                               ggCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGt      <  2:1347564/150‑1 (MQ=255)
                                                                                                                                                  tGTCGTGAGTCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGa   >  1:1333721/1‑150 (MQ=255)
                                                                                                                                                  tGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGa   <  1:609165/150‑1 (MQ=255)
                                                                                                                                                  tGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGa   <  1:1207813/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGAT  >  CP009273/3040877‑3041173

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: ‑