Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,732,978 C→T 6.1% intergenic (+91/‑87) malS → / → avtA alpha‑amylase/valine‑pyruvate aminotransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,732,9780CT 6.1% 1042.2 / 35.0 296intergenic (+91/‑87)malS/avtAalpha‑amylase/valine‑pyruvate aminotransferase 1
Reads supporting (aligned to +/- strand):  ref base C (142/136);  new base T (13/5);  total (155/141)
Fisher's exact test for biased strand distribution p-value = 9.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 297 total aligned reads displayed.

GCTTTGTTCGTGAGCATGGCGACGATAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGA  >  CP009273/3732829‑3733126
                                                                                                                                                     |                                                                                                                                                    
 cTTTGTTCGTGAGCATGGCGACGATAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCt                                                                                                                                                     <  1:377721/150‑1 (MQ=255)
   ttGTTCGTGAGCATGGCGACGATAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGTTCa                                                                                                                                                   >  1:1428513/1‑150 (MQ=255)
    tGTTCGTGAGCATGGCGACGATAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCAc                                                                                                                                                  <  2:939775/150‑1 (MQ=255)
         gTGAGCATGGCGACGATAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTc                                                                                                                                             >  1:1258131/1‑150 (MQ=255)
          tGAGCATGGCGACGATAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCa                                                                                                                                            <  1:3440782/150‑1 (MQ=255)
           gAGCATGGCGACGATAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCaa                                                                                                                                           >  2:4115527/1‑150 (MQ=255)
                   cgacgaTAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGa                                                                                                                                   >  1:2723417/1‑150 (MQ=255)
                      cgaTAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGAttt                                                                                                                                >  2:1267042/1‑150 (MQ=255)
                      cgaTAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGAttt                                                                                                                                <  1:3652609/150‑1 (MQ=255)
                        aTAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACACTTCAACAGACCGATTTGc                                                                                                                              >  2:2507140/1‑150 (MQ=255)
                              tGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCa                                                                                                                        <  1:969492/150‑1 (MQ=255)
                              tGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCa                                                                                                                        <  2:1480577/150‑1 (MQ=255)
                               gCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAg                                                                                                                       >  1:2761949/1‑150 (MQ=255)
                               gCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAg                                                                                                                       <  2:2966782/150‑1 (MQ=255)
                               gCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAg                                                                                                                       >  1:1312288/1‑150 (MQ=255)
                               gCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAg                                                                                                                       >  1:3584865/1‑150 (MQ=255)
                                cTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGc                                                                                                                      >  1:2637661/1‑150 (MQ=255)
                                 tGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGAACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCa                                                                                                                     >  1:4444131/1‑150 (MQ=255)
                                    tcgtcTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAAt                                                                                                                  <  2:425547/150‑1 (MQ=255)
                                    tcgtcTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAAt                                                                                                                  <  2:3422236/150‑1 (MQ=255)
                                     cgtcTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATg                                                                                                                 >  2:3463534/1‑150 (MQ=255)
                                           ggcagggcaACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGt                                                                                                           >  1:903622/1‑150 (MQ=255)
                                            gcagggcaACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGtt                                                                                                          <  2:2898828/150‑1 (MQ=255)
                                             cagggcaACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTa                                                                                                         >  2:2118405/1‑150 (MQ=255)
                                              agggcaACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGTTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTAt                                                                                                        >  1:661363/1‑150 (MQ=255)
                                               gggcaACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTAtt                                                                                                       <  1:1141168/150‑1 (MQ=255)
                                               gggcaACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTAtt                                                                                                       <  2:678170/150‑1 (MQ=255)
                                                   aaCAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTa                                                                                                   <  2:1025596/150‑1 (MQ=255)
                                                    aCAGTAACTTTTCCGTCTTCCCGTTCGTCAGTACCTCTGGAAGCCGCCAACCAGGATAAAATGTCAGCCATAATCAGCGTTGCAGGATACAGCAACGCTCACTCTTCAACAGACCGATTTGCACAACAGCAAATGTAGCGTTATTGTTAc                                                                                                  >  2:2623693/1‑150 (MQ=255)
                                                     cAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTAcc                                                                                                 >  1:3258194/1‑150 (MQ=255)
                                                          aCTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTAccttcct                                                                                            >  1:1853912/1‑150 (MQ=255)
                                                          aCTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTAccttcct                                                                                            <  2:3947960/150‑1 (MQ=255)
                                                          aCTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTAccttcct                                                                                            >  1:4109228/1‑150 (MQ=255)
                                                          aCTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTAccttcct                                                                                            <  2:3373058/150‑1 (MQ=255)
                                                              ttCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGAATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCt                                                                                        >  2:3739440/1‑150 (MQ=255)
                                                               tCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTa                                                                                       >  2:3969546/1‑150 (MQ=255)
                                                                 cGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACa                                                                                     >  2:142018/1‑150 (MQ=255)
                                                                 cGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAAACAGCGTTGCAGGATAAAGCTCCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTAATGTTACCTTCCTTGCTACa                                                                                     >  2:3029508/1‑150 (MQ=255)
                                                                      tCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGtt                                                                                <  1:3972780/150‑1 (MQ=255)
                                                                      tCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGtt                                                                                <  1:3592078/150‑1 (MQ=255)
                                                                      tCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGtt                                                                                <  1:3592939/150‑1 (MQ=255)
                                                                       cccGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTc                                                                               >  2:3470318/1‑150 (MQ=255)
                                                                          gTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGAc                                                                            <  2:178646/150‑1 (MQ=255)
                                                                            tCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAg                                                                          <  1:819807/150‑1 (MQ=255)
                                                                            tCGTCAGGACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAg                                                                          <  2:1447044/150‑1 (MQ=255)
                                                                                 aGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGTTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATc                                                                     >  2:309829/1‑150 (MQ=255)
                                                                                 aGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGTTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATc                                                                     >  2:1808171/1‑150 (MQ=255)
                                                                                 aGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATc                                                                     <  1:3623294/150‑1 (MQ=255)
                                                                                  gTACCTCGGGAAGCCGCCCACCATGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATcc                                                                    <  2:2266513/150‑1 (MQ=255)
                                                                                   tACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGTTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATccc                                                                   <  2:4029370/150‑1 (MQ=255)
                                                                                     ccTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGc                                                                 <  1:384240/150‑1 (MQ=255)
                                                                                     ccTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGc                                                                 <  2:667717/150‑1 (MQ=255)
                                                                                      cTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCt                                                                <  1:3659784/150‑1 (MQ=255)
                                                                                        cGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTAGCTACAGAGTTCGACAGATATCCCGCTAt                                                              <  2:327911/150‑1 (MQ=255)
                                                                                         gggAAGCCGCCAACCAGGATAAAATGTCGGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATg                                                             >  1:76417/1‑150 (MQ=255)
                                                                                          ggAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGa                                                            <  1:2394463/150‑1 (MQ=255)
                                                                                            aaGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTTTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACa                                                          >  2:1117261/1‑150 (MQ=255)
                                                                                                cgccAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATtct                                                      >  2:2528521/1‑150 (MQ=255)
                                                                                                 gccAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATtctc                                                     >  1:4368123/1‑150 (MQ=255)
                                                                                                 gccAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATtctc                                                     >  1:4400924/1‑150 (MQ=255)
                                                                                                  ccAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGTTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTcc                                                    >  1:2963363/1‑150 (MQ=255)
                                                                                                     aCCCGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCtt                                                 <  1:461977/150‑1 (MQ=255)
                                                                                                     aCCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCtt                                                 <  1:2907052/150‑1 (MQ=255)
                                                                                                     aCCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCtt                                                 >  2:2776204/1‑150 (MQ=255)
                                                                                                      ccAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCttt                                                >  2:1733060/1‑150 (MQ=255)
                                                                                                         ggATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTg                                             >  1:830826/1‑150 (MQ=255)
                                                                                                          gATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTgg                                            <  1:2881398/150‑1 (MQ=255)
                                                                                                           aTAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGt                                           <  1:1726735/150‑1 (MQ=255)
                                                                                                           aTAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGt                                           <  1:4219278/150‑1 (MQ=255)
                                                                                                            tAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTg                                          >  2:2654266/1‑150 (MQ=255)
                                                                                                             aaaaTGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGa                                         >  1:3101332/1‑150 (MQ=255)
                                                                                                              aaaTGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGAc                                        >  1:414818/1‑150 (MQ=255)
                                                                                                              aaaTGTCAGCCCTAATCAGAGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGAc                                        <  2:1797694/150‑1 (MQ=255)
                                                                                                                 tGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACaaa                                     <  1:1000098/150‑1 (MQ=255)
                                                                                                                   tCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAAtt                                   >  2:3428369/1‑150 (MQ=255)
                                                                                                                     aGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTa                                 >  1:4047327/1‑150 (MQ=255)
                                                                                                                       cccTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTAcc                               >  1:492325/1‑150 (MQ=255)
                                                                                                                        ccTAATCAGCGTTGCAGGATAAAGCACCGTTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTAccc                              >  2:641365/1‑150 (MQ=255)
                                                                                                                        ccTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTAccc                              <  1:2118434/150‑1 (MQ=255)
                                                                                                                         cTAATCAGCGTTGCAGGATAAAGCACCGTTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCg                             <  2:2411022/150‑1 (MQ=255)
                                                                                                                          tAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGc                            >  1:3335154/1‑150 (MQ=255)
                                                                                                                           aaTCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGcc                           >  2:4082430/1‑150 (MQ=255)
                                                                                                                             tCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCAc                         >  2:3772136/1‑150 (MQ=255)
                                                                                                                             tCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCAc                         >  1:1980520/1‑150 (MQ=255)
                                                                                                                             tCAGCGTTGCAGGATAAAGCACCGCACACACTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTAATTGTTACCATCCTAGCACCAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCAc                         >  1:4229943/1‑150 (MQ=255)
                                                                                                                                gCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTcc                      <  2:479742/150‑1 (MQ=255)
                                                                                                                                  gTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCgg                    <  1:4273051/150‑1 (MQ=255)
                                                                                                                                     gCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCAt                 >  2:4288032/1‑150 (MQ=255)
                                                                                                                                      cAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACTGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCAtt                >  1:766485/1‑150 (MQ=255)
                                                                                                                                      cAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCAtt                >  2:200096/1‑150 (MQ=255)
                                                                                                                                        ggATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTAc              <  1:463002/150‑1 (MQ=255)
                                                                                                                                         gATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACg             >  2:300133/1‑150 (MQ=255)
                                                                                                                                          aTAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGc            >  1:2574257/1‑150 (MQ=255)
                                                                                                                                          aTAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGc            >  2:1347451/1‑150 (MQ=255)
                                                                                                                                           tAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCt           >  2:329726/1‑150 (MQ=255)
                                                                                                                                           tAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCt           <  2:2637661/150‑1 (MQ=255)
                                                                                                                                           tAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCt           <  1:3733180/150‑1 (MQ=255)
                                                                                                                                              aGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCACTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGtt        >  2:1906012/1‑150 (MQ=255)
                                                                                                                                                 aCCGTTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGAt     >  2:3493200/1‑150 (MQ=255)
                                                                                                                                                    gCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGa  >  1:1057081/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTTTGTTCGTGAGCATGGCGACGATAAAGTGCTGGTCGTCTGGGCAGGGCAACAGTAACTTTTCCGGCTTCCCGTTCGTCAGTACCTCGGGAAGCCGCCAACCAGGATAAAATGTCAGCCCTAATCAGCGTTGCAGGATAAAGCACCGCTCACTCTTCAACAGACCGATTTGCACCCCAGCAAATGTAGCGTTATTGTTACCTTCCTTGCTACAGAGTTCGACAGATATCCCGCTATGACATTCTCCCTTTTTGGTGACAAATTTACCCGCCACTCCGGCATTACGCTGTTGATGGA  >  CP009273/3732829‑3733126

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: