Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,041,026 C→A 6.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,0260CA 6.0% 891.8 / 21.0 237D459Y (GAC→TAC) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (113/109);  new base A (4/11);  total (117/120)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

GGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGATTT  >  CP009273/3040878‑3041175
                                                                                                                                                    |                                                                                                                                                     
gggATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCg                                                                                                                                                      >  2:417610/1‑150 (MQ=255)
 ggATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGt                                                                                                                                                     >  2:940244/1‑150 (MQ=255)
   atcatcGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGa                                                                                                                                                   <  2:1006494/150‑1 (MQ=255)
        cGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCAc                                                                                                                                              <  1:3384592/150‑1 (MQ=255)
        cGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCAc                                                                                                                                              >  1:3963354/1‑150 (MQ=255)
        cGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCAc                                                                                                                                              <  2:2719592/150‑1 (MQ=255)
        cGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCAc                                                                                                                                              <  2:3420355/150‑1 (MQ=255)
        cGCCTGAGTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCAc                                                                                                                                              <  2:2719731/150‑1 (MQ=255)
           cTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTc                                                                                                                                           >  1:4080915/1‑150 (MQ=255)
             gATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCtt                                                                                                                                         >  2:2083155/1‑150 (MQ=255)
              aTTCAGCGCCATAACTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTAGCTGTGGTAGCGATTAAAAACCGGAGGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCttt                                                                                                                                        >  2:3852074/1‑150 (MQ=255)
              aTTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCttt                                                                                                                                        >  2:3852064/1‑150 (MQ=255)
                  aGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTcc                                                                                                                                    >  1:2513000/1‑150 (MQ=255)
                        aGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCgtg                                                                                                                              >  2:1231192/1‑150 (MQ=255)
                        aGATCTTTACGCTCCAGCGACTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCgtg                                                                                                                              <  2:3768432/150‑1 (MQ=255)
                         gATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCgtgt                                                                                                                             >  1:976633/1‑150 (MQ=255)
                           tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                           >  1:1237029/1‑150 (MQ=255)
                           tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                           >  1:2721184/1‑150 (MQ=255)
                           tCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCg                                                                                                                           >  1:2721226/1‑150 (MQ=255)
                              ttACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATg                                                                                                                        <  2:4372269/150‑1 (MQ=255)
                                   cTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAg                                                                                                                   <  1:588909/150‑1 (MQ=255)
                                      cAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGcc                                                                                                                <  2:1483235/150‑1 (MQ=255)
                                         cGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTg                                                                                                             <  1:3420011/150‑1 (MQ=255)
                                           aGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGt                                                                                                           >  2:125555/1‑150 (MQ=255)
                                           aGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGt                                                                                                           <  1:1389002/150‑1 (MQ=255)
                                            gTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGtt                                                                                                          <  2:2766552/150‑1 (MQ=255)
                                            gTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGtt                                                                                                          <  2:2765306/150‑1 (MQ=255)
                                            gTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGtt                                                                                                          <  2:2519797/150‑1 (MQ=255)
                                            gTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGtt                                                                                                          <  2:3872721/150‑1 (MQ=255)
                                            gTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGtt                                                                                                          >  2:4356067/1‑150 (MQ=255)
                                                atatAGCGCATCATTTCGGTTTCGCTGTGGTAGCTATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCg                                                                                                      >  2:3112900/1‑150 (MQ=255)
                                                atatAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCg                                                                                                      >  2:891201/1‑150 (MQ=255)
                                                    aGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCca                                                                                                  <  1:1765718/150‑1 (MQ=255)
                                                    aGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCca                                                                                                  >  2:2836101/1‑150 (MQ=255)
                                                     gcgcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcag                                                                                                 <  1:1448682/150‑1 (MQ=255)
                                                     gcgcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcag                                                                                                 <  2:3498368/150‑1 (MQ=255)
                                                      cgcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGGCGTGAGCCACGTCTTGGTCCAGCGTGTCGAGGTCCAGGCCGGGGTTATCGCCcagc                                                                                                >  2:3336185/1‑150 (MQ=255)
                                                       gcATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCcagca                                                                                               <  2:2663277/150‑1 (MQ=255)
                                                         atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                             <  2:1887864/150‑1 (MQ=255)
                                                         atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                             >  1:770915/1‑150 (MQ=255)
                                                         atcatTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCagcagc                                                                                             >  1:398778/1‑150 (MQ=255)
                                                            atTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACg                                                                                          <  1:2672266/150‑1 (MQ=255)
                                                                cGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGa                                                                                      <  1:469034/150‑1 (MQ=255)
                                                                  gTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGaaa                                                                                    >  1:1889962/1‑150 (MQ=255)
                                                                    ttCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAg                                                                                  <  2:1389228/150‑1 (MQ=255)
                                                                     tCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGc                                                                                 >  1:150779/1‑150 (MQ=255)
                                                                     tCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGc                                                                                 >  1:1022555/1‑150 (MQ=255)
                                                                     tCGCTGTGGTAGAGATTAAACACCGGATTGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGTCTGGATTGAGCGGCTGTAGTGATCCACTTCTTTGTCCAGCGTGTAGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAATAGc                                                                                 <  2:1201308/150‑1 (MQ=255)
                                                                       gCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTg                                                                               >  2:2597862/1‑150 (MQ=255)
                                                                        cTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGc                                                                              <  2:1737745/150‑1 (MQ=255)
                                                                          gtgGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCAt                                                                            <  1:1991315/150‑1 (MQ=255)
                                                                             gTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTAc                                                                         <  2:1092811/150‑1 (MQ=255)
                                                                                  gATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGtttt                                                                    >  1:2226855/1‑150 (MQ=255)
                                                                                  gATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGtttt                                                                    >  1:3643824/1‑150 (MQ=255)
                                                                                    ttAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCa                                                                  <  1:3088598/150‑1 (MQ=255)
                                                                                       aaCACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcgc                                                               >  2:4171022/1‑150 (MQ=255)
                                                                                        acacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcgcg                                                              >  2:827772/1‑150 (MQ=255)
                                                                                        acacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCAcgcg                                                              <  1:2920395/150‑1 (MQ=255)
                                                                                         cacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGt                                                             >  1:113762/1‑150 (MQ=255)
                                                                                         cacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGt                                                             <  1:3737480/150‑1 (MQ=255)
                                                                                         cacCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGt                                                             >  2:3813011/1‑150 (MQ=255)
                                                                                         cacCGGAGGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGt                                                             <  1:230763/150‑1 (MQ=255)
                                                                                           ccGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTgg                                                           <  2:1857779/150‑1 (MQ=255)
                                                                                              gATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttg                                                        <  1:2997281/150‑1 (MQ=255)
                                                                                                tGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGgttgtt                                                      >  2:217454/1‑150 (MQ=255)
                                                                                                 gggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGttt                                                     >  2:293988/1‑150 (MQ=255)
                                                                                                 gggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGttt                                                     <  2:2774424/150‑1 (MQ=255)
                                                                                                  ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                    <  1:1799091/150‑1 (MQ=255)
                                                                                                  ggTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtc                                                    <  2:3149245/150‑1 (MQ=255)
                                                                                                   gTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtca                                                   >  2:414106/1‑150 (MQ=255)
                                                                                                      aGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtcatca                                                <  2:2513000/150‑1 (MQ=255)
                                                                                                      aGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTtcatca                                                <  2:6201/150‑1 (MQ=255)
                                                                                                          tttCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAggg                                            <  1:2891102/150‑1 (MQ=255)
                                                                                                               tcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATc                                       <  2:491671/150‑1 (MQ=255)
                                                                                                               tcgtcgCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATc                                       <  2:1925708/150‑1 (MQ=255)
                                                                                                                    gcgcAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGACTCTTCTTTGTACAGCGTGTAGATGTACTGTCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCTAc                                  <  1:416487/150‑1 (MQ=255)
                                                                                                                     cgcAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAAcc                                 <  2:1187478/150‑1 (MQ=255)
                                                                                                                      gcagcaTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACcg                                >  2:2677850/1‑150 (MQ=255)
                                                                                                                        agcaTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACcgcg                              >  1:859211/1‑150 (MQ=255)
                                                                                                                          caTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGtt                            <  2:3545651/150‑1 (MQ=255)
                                                                                                                          caTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGtt                            >  1:559067/1‑150 (MQ=255)
                                                                                                                           aTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTc                           >  1:2981870/1‑150 (MQ=255)
                                                                                                                           aTCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTc                           >  1:2982793/1‑150 (MQ=255)
                                                                                                                            tCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCa                          <  1:747686/150‑1 (MQ=255)
                                                                                                                             cGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAg                         <  2:3982569/150‑1 (MQ=255)
                                                                                                                                  gCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAatat                    <  2:472414/150‑1 (MQ=255)
                                                                                                                                  gCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAatat                    <  1:1190756/150‑1 (MQ=255)
                                                                                                                                    tGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCg                  <  1:3236336/150‑1 (MQ=255)
                                                                                                                                    tGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCg                  >  2:3474513/1‑150 (MQ=255)
                                                                                                                                     ggATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGc                 >  2:2995875/1‑150 (MQ=255)
                                                                                                                                     ggATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGc                 >  2:1103479/1‑150 (MQ=255)
                                                                                                                                       aTAGAGCGGCTGTAGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTa               <  2:1642488/150‑1 (MQ=255)
                                                                                                                                         agagCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACg             <  2:3781356/150‑1 (MQ=255)
                                                                                                                                            gCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAg          <  2:1253520/150‑1 (MQ=255)
                                                                                                                                            gCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAg          <  2:2395741/150‑1 (MQ=255)
                                                                                                                                                cTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTg      <  2:4080915/150‑1 (MQ=255)
                                                                                                                                                cTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTg      <  2:753875/150‑1 (MQ=255)
                                                                                                                                                  gTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGAt    <  1:2068925/150‑1 (MQ=255)
                                                                                                                                                   tCGTGAGCCACGTCTTTGTCGAGCGTGTCGATGTCCAGGCCGTGGTTTTCGCCCAGCAGCACGTTGAATAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGGTTCGCAGGTTGAtt   <  2:1998289/150‑1 (MQ=255)
                                                                                                                                                    cGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGAttt  >  1:3452622/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGATCATCGCCTGATTCAGCGCCAGATCTTTACGCTCCAGCGAGTGCATATAGCGCATCATTTCGGTTTCGCTGTGGTAGCGATTAAACACCGGATGGGTGAGGATTTCGTCGTCGCGCAGCATCGCAGGCTGGATAGAGCGGCTGTCGTGAGCCACGTCTTTGTCCAGCGTGTCGATGTCCAGGCCGTGGTTATCGCCCAGCAGCACGTTGAAAAGCTGCATTACGTTTTCACGCGTGGTTGTTTCATCAAGGGTGATCCCAACCGCGTTCAGAATATCGCTACGCAGGTTGATTT  >  CP009273/3040878‑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: