Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,041,026 C→A 11.0% 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,0260CA11.0% 1125.6 / 67.7 307D459Y (GAC→TAC) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (154/119);  new base A (27/7);  total (181/126)
Fisher's exact test for biased strand distribution p-value = 9.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

CGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGATTT  >  CP009273/3040877‑3041175
                                                                                                                                                     |                                                                                                                                                     
cGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTc                                                                                                                                                       <  1:3858282/150‑1 (MQ=255)
   gATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTg                                                                                                                                                    >  1:4125107/1‑150 (MQ=255)
         cGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCAc                                                                                                                                              <  1:2833857/150‑1 (MQ=255)
          gtcTGATTCAGAGCCAGATCTTTACGCTCGAGCGAGTTCATATAGCGCATCTTTTCGGTTTCGCTGTGGTAGAGATTAAACATCGGATGTGTGAGGATTTCGTCTTCGCGAAGCATCGCTGTCTGGATAGAGCGGCTGTCGTGAGCCACg                                                                                                                                             <  1:3501850/148‑1 (MQ=255)
             tGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCt                                                                                                                                          >  1:1582731/1‑150 (MQ=255)
               aTTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCttt                                                                                                                                        <  1:177397/150‑1 (MQ=255)
               aTTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCttt                                                                                                                                        <  1:2015556/150‑1 (MQ=255)
                 tCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGt                                                                                                                                      >  2:2559467/1‑150 (MQ=255)
                   aGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTcc                                                                                                                                    >  1:4160911/1‑150 (MQ=255)
                    gcgcCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCa                                                                                                                                   >  2:1239849/1‑150 (MQ=255)
                    gcgcCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCa                                                                                                                                   >  2:3543659/1‑150 (MQ=255)
                      gcCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGc                                                                                                                                 >  1:1960721/1‑150 (MQ=255)
                         aGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCgtg                                                                                                                              <  1:1120836/150‑1 (MQ=255)
                         aGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCgtg                                                                                                                              >  2:3409/1‑150 (MQ=255)
                         aGATCATTACGCACCAGCGAGTGCATATAGCGCATCATTTCGGATTCGCTGTGGTAGCGATAAAACACCGGATGTGAGAGTATTTCGTCGTCTCGCAGCATCGCTGGCTGGATAGAGCGGCAGTAGTGAGCCACGTCTTTGTCAAgagag                                                                                                                              >  1:1994023/1‑146 (MQ=255)
                           aTCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTc                                                                                                                            <  1:54785/150‑1 (MQ=255)
                             cTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGa                                                                                                                          >  1:938527/1‑150 (MQ=255)
                              tttACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTTAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGAt                                                                                                                         <  2:2898991/150‑1 (MQ=255)
                               ttACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATg                                                                                                                        <  1:16969/150‑1 (MQ=255)
                                 aCGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTc                                                                                                                      >  1:2224669/1‑150 (MQ=255)
                                 aCGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTc                                                                                                                      >  1:646790/1‑150 (MQ=255)
                                 aCGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTc                                                                                                                      <  2:722142/150‑1 (MQ=255)
                                    cTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAg                                                                                                                   >  2:3133660/1‑150 (MQ=255)
                                     tCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAgg                                                                                                                  >  1:646771/1‑150 (MQ=255)
                                            aGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGt                                                                                                           >  2:4037155/1‑150 (MQ=255)
                                            aGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGt                                                                                                           >  2:3528361/1‑150 (MQ=255)
                                                cATATAGCGCATCATTTCGGTTTCTCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCTGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATc                                                                                                       <  1:309026/150‑1 (MQ=255)
                                                cATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATc                                                                                                       >  1:1135978/1‑150 (MQ=255)
                                                 atatAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCg                                                                                                      <  2:3348481/150‑1 (MQ=255)
                                                  tataGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGc                                                                                                     >  2:4603749/1‑150 (MQ=255)
                                                    taGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGccc                                                                                                   >  2:4461528/1‑150 (MQ=255)
                                                     aGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCca                                                                                                  >  1:3188132/1‑150 (MQ=255)
                                                      gcgcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcag                                                                                                 <  1:3009822/150‑1 (MQ=255)
                                                       cgcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagc                                                                                                >  1:2425209/1‑150 (MQ=255)
                                                          atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                             >  2:2468560/1‑150 (MQ=255)
                                                          atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                             <  2:2048941/150‑1 (MQ=255)
                                                          atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                             >  2:3085008/1‑150 (MQ=255)
                                                          atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                             >  2:3889073/1‑150 (MQ=255)
                                                          atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                             >  2:4315710/1‑150 (MQ=255)
                                                            catTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCAc                                                                                           >  2:637602/1‑150 (MQ=255)
                                                            catTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCAc                                                                                           >  1:869729/1‑150 (MQ=255)
                                                             atTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACg                                                                                          <  2:275222/150‑1 (MQ=255)
                                                               ttCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGtt                                                                                        >  1:841491/1‑150 (MQ=255)
                                                                 cGGTTTCGCTGTGGTAGCGATTAAACACCGGAGGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGa                                                                                      >  2:279020/1‑150 (MQ=255)
                                                                       cGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCt                                                                                >  1:3709653/1‑150 (MQ=255)
                                                                          tgtgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCa                                                                             >  2:1091526/1‑150 (MQ=255)
                                                                            tgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCAtt                                                                           >  2:1372119/1‑150 (MQ=255)
                                                                             ggTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTa                                                                          <  1:670838/150‑1 (MQ=255)
                                                                              gTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTAc                                                                         <  2:2007303/150‑1 (MQ=255)
                                                                               tAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACg                                                                        >  1:2822130/1‑150 (MQ=255)
                                                                                 gCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGtt                                                                      <  1:138472/150‑1 (MQ=255)
                                                                                   gATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGtttt                                                                    >  2:4460160/1‑150 (MQ=255)
                                                                                     ttAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCa                                                                  <  1:2559467/150‑1 (MQ=255)
                                                                                     ttAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCa                                                                  >  2:1681247/1‑150 (MQ=255)
                                                                                      tAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAc                                                                 >  2:2929276/1‑150 (MQ=255)
                                                                                       aaaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcg                                                                <  1:4635944/150‑1 (MQ=255)
                                                                                       aaaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcg                                                                >  2:2104872/1‑150 (MQ=255)
                                                                                        aaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcgc                                                               <  1:3684184/150‑1 (MQ=255)
                                                                                            ccGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTgg                                                           <  2:1023489/150‑1 (MQ=255)
                                                                                             cGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgt                                                          >  1:3423349/1‑150 (MQ=255)
                                                                                              ggATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgtt                                                         >  2:3657066/1‑150 (MQ=255)
                                                                                              ggATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgtt                                                         <  1:4158223/150‑1 (MQ=255)
                                                                                              ggATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgtt                                                         >  1:2916224/1‑150 (MQ=255)
                                                                                              ggATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgtt                                                         <  2:591946/150‑1 (MQ=255)
                                                                                               gATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttg                                                        <  2:1204190/150‑1 (MQ=255)
                                                                                               gATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttg                                                        <  1:711850/150‑1 (MQ=255)
                                                                                                aTGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgt                                                       >  1:2290953/1‑150 (MQ=255)
                                                                                                 tGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgtt                                                      >  2:4106959/1‑150 (MQ=255)
                                                                                                    gTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtca                                                   <  2:3332855/150‑1 (MQ=255)
                                                                                                     tGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtcat                                                  <  1:289105/150‑1 (MQ=255)
                                                                                                     tGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtcat                                                  >  2:2162899/1‑150 (MQ=255)
                                                                                                         gATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAg                                              >  1:833400/1‑150 (MQ=255)
                                                                                                          aTTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAgg                                             >  1:3372556/1‑150 (MQ=255)
                                                                                                           tttCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAggg                                            <  2:4125107/150‑1 (MQ=255)
                                                                                                            ttCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGt                                           <  2:563476/150‑1 (MQ=255)
                                                                                                              cgtcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGa                                         >  1:4166496/1‑150 (MQ=255)
                                                                                                                 cgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATcc                                      >  1:2611011/1‑150 (MQ=255)
                                                                                                                  gtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGGTGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATccc                                     >  1:1455990/1‑150 (MQ=255)
                                                                                                                    cgcgCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCCGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCaa                                   >  2:1303941/1‑150 (MQ=255)
                                                                                                                     gcgcAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAc                                  <  2:4578139/150‑1 (MQ=255)
                                                                                                                           caTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGtt                            >  2:1965955/1‑150 (MQ=255)
                                                                                                                            aTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGtcc                           <  1:2760836/150‑3 (MQ=255)
                                                                                                                            aTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTc                           >  1:2103935/1‑150 (MQ=255)
                                                                                                                            aTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTc                           >  1:3078608/1‑150 (MQ=255)
                                                                                                                              cGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAg                         <  1:4449480/150‑1 (MQ=255)
                                                                                                                               gCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGa                        >  1:3975478/1‑150 (MQ=255)
                                                                                                                                 aGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAat                      >  2:1212981/1‑150 (MQ=255)
                                                                                                                                   gCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAatat                    >  2:4218135/1‑150 (MQ=255)
                                                                                                                                   gCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAatat                    >  2:826122/1‑150 (MQ=255)
                                                                                                                                   gCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAatat                    >  2:1589409/1‑150 (MQ=255)
                                                                                                                                    cTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATc                   <  1:2373332/150‑1 (MQ=255)
                                                                                                                                         tAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTAc              >  1:1553787/1‑150 (MQ=255)
                                                                                                                                         tAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTAc              <  2:4525125/150‑1 (MQ=255)
                                                                                                                                          agagCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACg             >  2:546212/1‑150 (MQ=255)
                                                                                                                                          agagCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACg             >  2:3291316/1‑150 (MQ=255)
                                                                                                                                               ggCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGt        >  2:223153/1‑150 (MQ=255)
                                                                                                                                                gCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGtt       >  2:837814/1‑150 (MQ=255)
                                                                                                                                                gCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGtt       <  1:4579444/150‑1 (MQ=255)
                                                                                                                                                  tGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGa     >  1:1233524/1‑150 (MQ=255)
                                                                                                                                                    tCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGAtt   <  1:2667429/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGATTT  >  CP009273/3040877‑3041175

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: