Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,118,851 G→T 8.8% T84T (ACG→ACT metB → cystathionine gamma‑synthase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,118,8510GT 8.8% 1157.3 / 59.0 313T84T (ACG→ACTmetBcystathionine gamma‑synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (140/145);  new base T (13/15);  total (153/160)
Fisher's exact test for biased strand distribution p-value = 8.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 319 total aligned reads displayed.

TTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCATTACGGGCAGCGCT  >  CP009273/4118702‑4119000
                                                                                                                                                     |                                                                                                                                                     
tttACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACg                                                                                                                                                       >  2:3970783/1‑150 (MQ=255)
  tACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGAc                                                                                                                                                     >  2:2810617/1‑150 (MQ=255)
  tACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGAc                                                                                                                                                     <  1:1918553/150‑1 (MQ=255)
   aCCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGAcc                                                                                                                                                    >  1:1508592/1‑150 (MQ=255)
    ccGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCg                                                                                                                                                   >  2:100944/1‑150 (MQ=255)
     cGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGt                                                                                                                                                  >  2:2592506/1‑150 (MQ=255)
       gATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCt                                                                                                                                                >  1:536419/1‑150 (MQ=255)
         tttAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCttt                                                                                                                                              <  2:1996912/150‑1 (MQ=255)
         tttAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCttt                                                                                                                                              <  1:2831119/150‑1 (MQ=255)
             aTGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGa                                                                                                                                          >  1:3771286/1‑150 (MQ=255)
                 aCCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAAcc                                                                                                                                      <  2:3468323/150‑1 (MQ=255)
                  ccGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCt                                                                                                                                     <  2:1183970/150‑1 (MQ=255)
                  ccGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCt                                                                                                                                     >  1:67832/1‑150 (MQ=255)
                    gcgcgcgcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTgg                                                                                                                                   <  1:2238873/150‑1 (MQ=255)
                         cgcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATc                                                                                                                              <  1:1821906/150‑1 (MQ=255)
                              gATGACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATctgctg                                                                                                                         <  2:1635131/150‑1 (MQ=255)
                               aTTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTgg                                                                                                                        >  2:3367246/1‑150 (MQ=255)
                                 tACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGtt                                                                                                                      <  2:3263629/150‑1 (MQ=255)
                                  aCTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTg                                                                                                                     <  1:213804/150‑1 (MQ=255)
                                  aCTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTg                                                                                                                     >  1:3447552/1‑150 (MQ=255)
                                   cTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTgc                                                                                                                    <  1:2268547/150‑1 (MQ=255)
                                      gcgtcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgcc                                                                                                                 <  1:1808292/150‑1 (MQ=255)
                                          cgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCAc                                                                                                             <  2:1493992/150‑1 (MQ=255)
                                          cgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCAc                                                                                                             >  2:3591141/1‑150 (MQ=255)
                                           gcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACg                                                                                                            <  2:3801579/150‑1 (MQ=255)
                                            cgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGa                                                                                                           >  1:3568282/1‑150 (MQ=255)
                                             ggCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGAc                                                                                                          >  2:3704930/1‑150 (MQ=255)
                                                aaCCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGc                                                                                                       >  2:685316/1‑150 (MQ=255)
                                                 aCCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCt                                                                                                      <  1:4037427/150‑1 (MQ=255)
                                                  cccAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTa                                                                                                     <  2:1059943/150‑1 (MQ=255)
                                                       cgcgcgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcggcg                                                                                                <  1:892420/150‑1 (MQ=255)
                                                       cgcgcgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcggcg                                                                                                <  2:2729344/150‑1 (MQ=255)
                                                              tgtgGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTa                                                                                         <  1:3128928/150‑1 (MQ=255)
                                                               gtgGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTAt                                                                                        <  1:4444206/150‑1 (MQ=255)
                                                               gtgGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTAt                                                                                        >  1:2278108/1‑150 (MQ=255)
                                                                  gTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGc                                                                                     <  2:3190073/150‑1 (MQ=255)
                                                                   ttCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGcc                                                                                    <  2:744602/150‑1 (MQ=255)
                                                                    tCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCt                                                                                   >  1:4065678/1‑150 (MQ=255)
                                                                      aGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGt                                                                                 >  2:1570186/1‑150 (MQ=255)
                                                                       gCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGtt                                                                                <  2:4520429/150‑1 (MQ=255)
                                                                        cGTGCGCTGGCAGACCTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCGGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTc                                                                               <  1:3088762/150‑1 (MQ=255)
                                                                        cGTGCGCTGGCAGAACTGGTAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTc                                                                               <  1:3681082/150‑1 (MQ=255)
                                                                         gTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCg                                                                              <  2:4275818/150‑1 (MQ=255)
                                                                              cTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATATGCTGGTTTCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGt                                                                         >  2:540252/1‑150 (MQ=255)
                                                                               tGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTc                                                                        >  1:2480683/1‑150 (MQ=255)
                                                                                ggCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCt                                                                       <  2:123378/150‑1 (MQ=255)
                                                                                ggCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCt                                                                       >  1:4111275/1‑150 (MQ=255)
                                                                                  cAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTgg                                                                     >  1:4531573/1‑150 (MQ=255)
                                                                                  cAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTgg                                                                     <  2:1781229/150‑1 (MQ=255)
                                                                                   aGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGc                                                                    <  2:4304298/150‑1 (MQ=255)
                                                                                         ggAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGAAGCTATCGCCTGTTCGACAGTCTGGCGAAAcg                                                              >  2:2161534/1‑150 (MQ=255)
                                                                                          gAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAAcgc                                                             >  1:1180553/1‑150 (MQ=255)
                                                                                            aGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTCTCGCCTGTTCGACAGTCTGGCGAAACGCgg                                                           >  2:2886355/1‑150 (MQ=255)
                                                                                             ggtggtGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGt                                                          <  1:113741/150‑1 (MQ=255)
                                                                                              gtggtgCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGtt                                                         <  1:2885275/150‑1 (MQ=255)
                                                                                               tggtgctggtgcAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTg                                                        >  2:2292855/1‑150 (MQ=255)
                                                                                                 gtgctggtgcAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCt                                                      <  1:1026235/150‑1 (MQ=255)
                                                                                                 gtgctggtgcAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCt                                                      <  2:562518/150‑1 (MQ=255)
                                                                                                    ctggtgcAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATc                                                   <  2:2064864/150‑1 (MQ=255)
                                                                                                    ctggtgcAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATc                                                   >  1:1348298/1‑150 (MQ=255)
                                                                                                      ggtgcAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATcgc                                                 >  2:48124/1‑150 (MQ=255)
                                                                                                       gtgcAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATcgcg                                                <  1:1614999/150‑1 (MQ=255)
                                                                                                        tgcAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCgt                                               >  2:720027/1‑150 (MQ=255)
                                                                                                          cAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGAGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCgtgt                                             >  2:4243607/1‑150 (MQ=255)
                                                                                                           aGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGtgtt                                            <  1:3267669/150‑1 (MQ=255)
                                                                                                            gTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGtgttg                                           <  2:1729031/150‑1 (MQ=255)
                                                                                                            gTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGtgttg                                           <  1:427435/150‑1 (MQ=255)
                                                                                                              acttactAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTgttgtt                                         <  1:713943/150‑1 (MQ=255)
                                                                                                               cttactAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGttgttt                                        <  1:2438795/150‑1 (MQ=255)
                                                                                                               cttactAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGttgttt                                        >  2:4478181/1‑150 (MQ=255)
                                                                                                                ttactAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGttgtttg                                       <  2:303675/150‑1 (MQ=255)
                                                                                                                   ctAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTttgttg                                    >  2:221716/1‑150 (MQ=255)
                                                                                                                    tAATACCGGCATGTCCGCGATTCACCTGGTAACTACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGa                                   <  1:1978121/150‑1 (MQ=255)
                                                                                                                      aTACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATc                                 <  1:253990/150‑1 (MQ=255)
                                                                                                                         ccGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTGTGTTGATCAAg                              >  2:2717923/1‑150 (MQ=255)
                                                                                                                           gggaTGTCGGCGTTTCACCTGGTAACGACGGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGTTCAAGGc                            <  2:187283/147‑1 (MQ=255)
                                                                                                                             cATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGa                          >  1:4481069/1‑150 (MQ=255)
                                                                                                                             cATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGa                          >  1:3417129/1‑150 (MQ=255)
                                                                                                                              aTGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGAt                         <  1:606380/150‑1 (MQ=255)
                                                                                                                              aTGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGAt                         >  2:1383938/1‑150 (MQ=255)
                                                                                                                               tGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATg                        >  2:881852/1‑150 (MQ=255)
                                                                                                                               tGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATg                        <  1:1427525/150‑1 (MQ=255)
                                                                                                                               tGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATg                        <  1:2044598/150‑1 (MQ=255)
                                                                                                                               tGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATg                        >  2:3315062/1‑150 (MQ=255)
                                                                                                                               tGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATg                        <  1:2987126/150‑1 (MQ=255)
                                                                                                                                   cgcgATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGTAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGTCGATGAACa                    <  2:2836823/150‑1 (MQ=255)
                                                                                                                                   cgcgATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACa                    <  2:2369035/150‑1 (MQ=255)
                                                                                                                                    gcgATTCACCTGGTAACGACCGTCTTTTTGAAACCTGTCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGACTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAg                   <  2:1600328/150‑1 (MQ=255)
                                                                                                                                     cgATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAgg                  >  2:3774172/1‑150 (MQ=255)
                                                                                                                                       aTTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCa                >  2:1442012/1‑150 (MQ=255)
                                                                                                                                       aTTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCa                >  2:618822/1‑150 (MQ=255)
                                                                                                                                        ttCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCAt               <  1:2666197/150‑1 (MQ=255)
                                                                                                                                        ttCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCAt               <  1:1776530/150‑1 (MQ=255)
                                                                                                                                        ttCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCAt               <  1:656548/150‑1 (MQ=255)
                                                                                                                                        gtCACCTGGTCACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCAt               <  1:1554221/149‑1 (MQ=255)
                                                                                                                                           aCCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCATTAc            <  2:2334721/150‑1 (MQ=255)
                                                                                                                                           aCCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCATTAc            <  2:3570226/150‑1 (MQ=255)
                                                                                                                                            ccTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCATTACg           <  2:4359864/150‑1 (MQ=255)
                                                                                                                                                  aaCGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCATTACGGGCAgc     >  2:1098644/1‑150 (MQ=255)
                                                                                                                                                   aCGACCGTCTTTTTGAAACCTGGCTATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCATTACGGGCAgcg    <  1:3873994/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACAGGCATTACGGGCAGCGCT  >  CP009273/4118702‑4119000

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: