Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,822,975 C→A 53.8% A339A (GCG→GCT astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,822,9750CA53.8% 33.9 / inf 223A339A (GCG→GCTastDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (59/44);  new base A (70/50);  total (129/94)
Fisher's exact test for biased strand distribution p-value = 8.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 227 total aligned reads displayed.

AAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTG  >  CP009273/1822826‑1823124
                                                                                                                                                     |                                                                                                                                                     
  aGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGc                                                                                                                                                     >  1:1351630/1‑150 (MQ=255)
  aGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGc                                                                                                                                                     <  2:546229/150‑1 (MQ=255)
  aGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGc                                                                                                                                                     <  2:1845401/150‑1 (MQ=255)
   gTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCt                                                                                                                                                    <  1:4538998/150‑1 (MQ=255)
   gTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCt                                                                                                                                                    >  1:3441635/1‑150 (MQ=255)
     aTCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTc                                                                                                                                                  <  2:1932177/150‑1 (MQ=255)
         tAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGt                                                                                                                                              <  2:852893/150‑1 (MQ=255)
          aaCGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGtt                                                                                                                                             >  1:1871223/1‑150 (MQ=255)
                aGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGcc                                                                                                                                       >  1:2370991/1‑150 (MQ=255)
                   cgcgCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATg                                                                                                                                    >  1:3669158/1‑150 (MQ=255)
                    gcgcAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGc                                                                                                                                   >  2:4562613/1‑150 (MQ=255)
                      gcAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAg                                                                                                                                 >  2:1244250/1‑150 (MQ=255)
                      gcAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAg                                                                                                                                 >  2:1244067/1‑150 (MQ=255)
                        aaTAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTa                                                                                                                               >  2:2433301/1‑150 (MQ=255)
                        aaTAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTa                                                                                                                               >  1:2070041/1‑150 (MQ=255)
                         aTAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTaa                                                                                                                              >  1:2168125/1‑150 (MQ=255)
                          tAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAac                                                                                                                             <  1:3576912/150‑1 (MQ=255)
                                   gAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACctgctg                                                                                                                    <  2:4585096/150‑1 (MQ=255)
                                     acacCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCtgtg                                                                                                                  >  1:232418/1‑150 (MQ=255)
                                      cacCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCACAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGc                                                                                                                 >  2:3433039/1‑150 (MQ=255)
                                        ccTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCgg                                                                                                               >  1:2219768/1‑150 (MQ=255)
                                        ccTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCgg                                                                                                               >  1:2905169/1‑150 (MQ=255)
                                        ccTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCgg                                                                                                               >  2:161319/1‑150 (MQ=255)
                                         ctctTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGc                                                                                                              <  2:589808/150‑1 (MQ=255)
                                          tctTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGcc                                                                                                             <  2:1420497/150‑1 (MQ=255)
                                            ttCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTg                                                                                                           <  2:2911339/150‑1 (MQ=255)
                                            ttCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTg                                                                                                           >  2:3144833/1‑150 (MQ=255)
                                             tCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGt                                                                                                          >  2:1034553/1‑150 (MQ=255)
                                                  tGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGa                                                                                                     <  1:85253/150‑1 (MQ=255)
                                                  tGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGa                                                                                                     <  1:1054217/150‑1 (MQ=255)
                                                    gTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGaaa                                                                                                   <  2:3112976/150‑1 (MQ=255)
                                                      aCGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATc                                                                                                 <  2:4270557/150‑1 (MQ=255)
                                                        gCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAg                                                                                               <  2:2294597/150‑1 (MQ=255)
                                                         ccAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGc                                                                                              >  1:355830/1‑150 (MQ=255)
                                                          cAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGcc                                                                                             >  1:576506/1‑150 (MQ=255)
                                                           aGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGccc                                                                                            >  1:1118817/1‑150 (MQ=255)
                                                           aGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGccc                                                                                            >  1:2291256/1‑150 (MQ=255)
                                                            gCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCg                                                                                           <  2:3833687/150‑1 (MQ=255)
                                                             cAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGc                                                                                          >  2:1254977/1‑150 (MQ=255)
                                                              aaCGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGcc                                                                                         >  2:3721135/1‑150 (MQ=255)
                                                              aaCGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGcc                                                                                         >  1:2905783/1‑150 (MQ=255)
                                                               aCGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCa                                                                                        <  1:441198/150‑1 (MQ=255)
                                                               aCGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCa                                                                                        >  2:3333868/1‑150 (MQ=255)
                                                                 gCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAAt                                                                                      >  1:142341/1‑150 (MQ=255)
                                                                 gCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAAt                                                                                      >  1:287388/1‑150 (MQ=255)
                                                                  ccTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATa                                                                                     >  1:4515452/1‑150 (MQ=255)
                                                                        aTTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGc                                                                               >  2:2342279/1‑150 (MQ=255)
                                                                           tCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGc                                                                            >  1:2634866/1‑150 (MQ=255)
                                                                           tCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGc                                                                            >  1:3104007/1‑150 (MQ=255)
                                                                              aTGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGt                                                                         <  2:2539400/150‑1 (MQ=255)
                                                                                gATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTc                                                                       <  1:2466599/150‑1 (MQ=255)
                                                                                gATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTc                                                                       <  1:2343583/150‑1 (MQ=255)
                                                                                 aTCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCg                                                                      >  2:370332/1‑150 (MQ=255)
                                                                                     ccGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCAt                                                                  <  1:2298454/150‑1 (MQ=255)
                                                                                     ccGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCAt                                                                  >  1:1021050/1‑150 (MQ=255)
                                                                                       ggCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATcc                                                                <  1:4014089/150‑1 (MQ=255)
                                                                                       ggCGTCAGCAACGAGGTCCCTGCTTGTAATAAGCGCGGCGCACGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATcc                                                                <  1:792168/150‑1 (MQ=255)
                                                                                          gTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAg                                                             >  1:1450607/1‑150 (MQ=255)
                                                                                            cAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGtt                                                           <  1:940409/150‑1 (MQ=255)
                                                                                            cAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGtt                                                           <  2:883795/150‑1 (MQ=255)
                                                                                            cAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGtt                                                           <  2:3862530/150‑1 (MQ=255)
                                                                                             aGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTg                                                          >  1:2192277/1‑150 (MQ=255)
                                                                                             aGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCATCATCCCAGTTg                                                          >  1:1611067/1‑150 (MQ=255)
                                                                                               cAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGcc                                                        <  2:3434941/150‑1 (MQ=255)
                                                                                               cAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGcc                                                        >  1:902949/1‑150 (MQ=255)
                                                                                               cAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGcc                                                        >  1:376213/1‑150 (MQ=255)
                                                                                                aaCGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGccc                                                       <  2:4281638/150‑1 (MQ=255)
                                                                                                aaCGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGccc                                                       >  1:3816299/1‑150 (MQ=255)
                                                                                                  cGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCgg                                                     >  1:2714758/1‑150 (MQ=255)
                                                                                                    aTGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCg                                                   >  1:2166783/1‑150 (MQ=255)
                                                                                                    aTGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCg                                                   >  1:644305/1‑150 (MQ=255)
                                                                                                    aTGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCg                                                   <  1:1251386/150‑1 (MQ=255)
                                                                                                     tGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGt                                                  >  2:2077910/1‑150 (MQ=255)
                                                                                                       tCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTa                                                <  1:3333305/150‑1 (MQ=255)
                                                                                                       tCCCTGCTAGTAATAAGCGCGGCGCAAGCAGGCGTCGTCCGCCCATCGCTTCCAGTAGCTGCCATTCAGAAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCAAGTTGCCCTGCGTTa                                                >  2:1151066/1‑150 (MQ=255)
                                                                                                        cccTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTaa                                               >  1:2553235/1‑150 (MQ=255)
                                                                                                         ccTGCTTTTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCACATGCAGCCAAGCAGTAACCACCTGCTGTCCGGCCTGATCAGAAATCAGCCCGCCAATAAACTGCAGCGGATCGTCATCCCACTTCCCCGGCGTTAAt                                              >  2:295370/1‑150 (MQ=255)
                                                                                                          cTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATc                                             >  2:4283689/1‑150 (MQ=255)
                                                                                                           tGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCg                                            >  1:3215193/1‑150 (MQ=255)
                                                                                                              ttGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctg                                         <  2:267815/150‑1 (MQ=255)
                                                                                                              ttGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctg                                         <  2:4400420/150‑1 (MQ=255)
                                                                                                               tGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctgg                                        <  2:525081/150‑1 (MQ=255)
                                                                                                                  aataaGCGCGGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCgctggctg                                     >  2:3175775/1‑150 (MQ=255)
                                                                                                                        cgcgGCGCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCa                               <  1:3277837/150‑1 (MQ=255)
                                                                                                                           ggcgCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAAcc                            <  1:1114493/150‑1 (MQ=255)
                                                                                                                           ggcgCAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAAcc                            <  2:3932937/150‑1 (MQ=255)
                                                                                                                              gcAAGCAGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGa                         <  2:2464989/150‑1 (MQ=255)
                                                                                                                                   cAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGc                    >  2:1893972/1‑150 (MQ=255)
                                                                                                                                    aGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCa                   >  2:3608459/1‑150 (MQ=255)
                                                                                                                                    aGGGGTCGTCCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTAAATCGCTGGCTGACGGCAACCAGACGAGCa                   >  2:654173/1‑150 (MQ=255)
                                                                                                                                           gtcCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACg            >  2:3832952/1‑150 (MQ=255)
                                                                                                                                           gtcCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACg            >  2:1621307/1‑150 (MQ=255)
                                                                                                                                           gtcCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACg            >  2:2686461/1‑150 (MQ=255)
                                                                                                                                            tcCGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGc           >  2:475868/1‑150 (MQ=255)
                                                                                                                                             ccGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCa          >  2:4080138/1‑150 (MQ=255)
                                                                                                                                              cGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCAt         <  1:3726166/150‑1 (MQ=255)
                                                                                                                                              cGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCAt         <  1:2648103/150‑1 (MQ=255)
                                                                                                                                              cGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCAt         >  1:1925460/1‑150 (MQ=255)
                                                                                                                                              cGCCCATAGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCAt         >  1:4253454/1‑150 (MQ=255)
                                                                                                                                                cccATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCg       >  2:1233656/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGTATCATAACGCCAGACGCGCAATAACGGTCCGAACACCTCTTCATCTGGTACGCCAGCAACGCCTGTCATTTCAATGATCCCCGGCGTCAGCAACGATGTCCCTGCTTGTAATAAGCGCGGCGCAAGCAGGGGTCGTCCGCCCATCGCTTCCAGTTGCTGCCATGCAGTAACCACCTGCTGTGCGGCCTGTTCAGAAATCAGCCCGCCAATAAACGGCTGCGGTTCGTCATCCCAGTTGCCCGGCGTTAATCGCTGGCTGACGGCAACCAGACGAGCAAGAAACGCATCGCCCTG  >  CP009273/1822826‑1823124

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: