Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,736,722 C→A 100% G210G (GGC→GGA yiaK → 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,736,7220CA98.5% 866.3 / 1.4 280G210G (GGC→GGAyiaK2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (113/161);  total (117/163)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 282 total aligned reads displayed.

CGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTC  >  CP009273/3736575‑3736871
                                                                                                                                                   |                                                                                                                                                     
    aTTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATc                                                                                                                                                 <  1:2898454/150‑1 (MQ=255)
        cttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCgaga                                                                                                                                             <  2:3603349/150‑1 (MQ=255)
          tccACGCCGATCACCCTGGTCGATATGTCGATGTCGATGTTCTCTGACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGCCCAAAGAACCTGGAGTTATCGagaag                                                                                                                                           <  2:667916/150‑1 (MQ=255)
          tccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGagaag                                                                                                                                           <  2:2463389/150‑1 (MQ=255)
             aCGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAAt                                                                                                                                        <  2:1806295/150‑1 (MQ=255)
                 cGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTTCGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATcgcc                                                                                                                                    <  2:2337760/150‑1 (MQ=255)
                 cGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGTTGAGGGCAATTTGACCAAAGAACCAGGAGTTATCGAGAAGAATcgcc                                                                                                                                    >  1:1823775/1‑150 (MQ=255)
                    tCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCCATTTGACCAAAGAACCGGGAGTTATCGAGAAGAATCGCCGCa                                                                                                                                 >  1:4093227/1‑150 (MQ=255)
                    tCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCa                                                                                                                                 >  1:1214495/1‑150 (MQ=255)
                     cACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCGGGAGTTATCGAGAAGAATCGCCGCAt                                                                                                                                >  2:4916641/1‑150 (MQ=255)
                      aCCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAtt                                                                                                                               <  1:5146816/150‑1 (MQ=255)
                      aCCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAtt                                                                                                                               <  1:1460391/150‑1 (MQ=255)
                      aCCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAtt                                                                                                                               <  2:1030118/150‑1 (MQ=255)
                       ccATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCAttt                                                                                                                              <  2:3214050/150‑1 (MQ=255)
                        cctGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGACGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTGGAGGATGAGGGCCAATTTACCCAACAACCCGGCGTTATCGAGAAGAATCGCCGCAtttt                                                                                                                             >  1:1892667/3‑150 (MQ=255)
                           ggTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGcc                                                                                                                          <  1:4139201/150‑1 (MQ=255)
                                    gtcgatgtcgatGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTa                                                                                                                 <  1:3782200/150‑1 (MQ=255)
                                      cgatgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACt                                                                                                               <  1:2987533/150‑1 (MQ=255)
                                      cgatgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACt                                                                                                               >  1:141137/1‑150 (MQ=255)
                                       gatgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTg                                                                                                              >  2:1345933/1‑150 (MQ=255)
                                        atgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTgg                                                                                                             >  1:2023846/1‑150 (MQ=255)
                                         tgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGa                                                                                                            >  2:2817062/1‑150 (MQ=255)
                                          gtcgatgtTCTCTTACGGCATGTTAGAAGTTACCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGaa                                                                                                           <  1:1837794/150‑1 (MQ=255)
                                            cgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAg                                                                                                         >  2:353165/1‑150 (MQ=255)
                                            cgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAg                                                                                                         <  2:3302831/150‑1 (MQ=255)
                                             gatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAgg                                                                                                        >  2:2984337/1‑150 (MQ=255)
                                              atgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGt                                                                                                       <  1:1378092/150‑1 (MQ=255)
                                              atgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAAATTGACCAAAGAACCTGGAGTTATCGAGTAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGt                                                                                                       >  1:4487749/1‑150 (MQ=255)
                                               tgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGTCAGGTCGTCAGCTCCCGTTCGATGGTGGCTTTTATTATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGtt                                                                                                      <  1:5222443/150‑1 (MQ=255)
                                               tgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGtt                                                                                                      <  1:2204613/150‑1 (MQ=255)
                                                gtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTc                                                                                                     >  2:1973209/1‑150 (MQ=255)
                                                 ttCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCt                                                                                                    <  2:1773375/150‑1 (MQ=255)
                                                 ttCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCt                                                                                                    >  2:1244412/1‑150 (MQ=255)
                                                    tctTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGc                                                                                                 <  2:2634349/150‑1 (MQ=255)
                                                     ctTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCa                                                                                                >  1:3106728/1‑150 (MQ=255)
                                                      ttACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCAt                                                                                               >  2:5058046/1‑150 (MQ=255)
                                                        aCGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGt                                                                                             >  2:5298267/1‑150 (MQ=255)
                                                        aCGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGt                                                                                             <  2:2749577/150‑1 (MQ=255)
                                                         cGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTc                                                                                            <  1:5308522/150‑1 (MQ=255)
                                                          ggCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCg                                                                                           >  2:4586965/1‑150 (MQ=255)
                                                           gCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGa                                                                                          <  1:806214/150‑1 (MQ=255)
                                                            cATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAt                                                                                         <  2:1141197/150‑1 (MQ=255)
                                                             aTGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAtt                                                                                        <  1:1245738/150‑1 (MQ=255)
                                                             aTGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAtt                                                                                        >  1:3059489/1‑150 (MQ=255)
                                                               gTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtgt                                                                                      <  2:2687178/150‑1 (MQ=255)
                                                                ttAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtgtg                                                                                     <  2:5053647/150‑1 (MQ=255)
                                                                   gAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTg                                                                                  <  2:4733078/150‑1 (MQ=255)
                                                                    aaGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCGTTGATGATGAGGGCAAATTGACCAAAGAACCGGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCGACTGGAAAGGGTCTGGCCTGTCGAATGTGCTgg                                                                                 >  1:4036269/1‑150 (MQ=255)
                                                                       ttAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGata                                                                              <  1:3946092/150‑1 (MQ=255)
                                                                         aaCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                            <  1:2092559/150‑1 (MQ=255)
                                                                         aaCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                            <  2:2531990/150‑1 (MQ=255)
                                                                         aaCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                            >  1:3161842/1‑150 (MQ=255)
                                                                         aaCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAAATTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATg                                                                            <  1:1722670/150‑1 (MQ=255)
                                                                          aCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGa                                                                           >  1:913623/1‑150 (MQ=255)
                                                                           ccGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGAt                                                                          <  1:1679529/150‑1 (MQ=255)
                                                                            cGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATc                                                                         <  2:4332578/150‑1 (MQ=255)
                                                                            cGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATc                                                                         <  1:4792251/150‑1 (MQ=255)
                                                                                tGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTa                                                                     >  1:3890988/1‑150 (MQ=255)
                                                                                    aGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTActct                                                                 <  2:1208141/150‑1 (MQ=255)
                                                                                      gtcgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTcc                                                               <  2:3632371/150‑1 (MQ=255)
                                                                                       tcgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCt                                                              <  2:56120/150‑1 (MQ=255)
                                                                                       tcgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCt                                                              <  2:562833/150‑1 (MQ=255)
                                                                                       tcgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCt                                                              <  1:1133524/150‑1 (MQ=255)
                                                                                        cgtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCtt                                                             >  1:2605209/1‑150 (MQ=255)
                                                                                         gtcAGCTCCCGGTCGATGGTGGCTTTGATGTTGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCttt                                                            >  1:2817129/1‑150 (MQ=255)
                                                                                         gtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCttt                                                            <  2:4780492/150‑1 (MQ=255)
                                                                                         gtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCttt                                                            <  1:1560798/150‑1 (MQ=255)
                                                                                           cAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAATAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCtttac                                                          >  1:2412520/1‑148 (MQ=255)
                                                                                           cAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTcc                                                          <  2:974443/150‑1 (MQ=255)
                                                                                            aGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCg                                                         >  1:1177808/1‑150 (MQ=255)
                                                                                               tCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACg                                                      >  2:2874279/1‑150 (MQ=255)
                                                                                               tCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACg                                                      <  2:479552/150‑1 (MQ=255)
                                                                                                         tggtggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGt                                            >  2:88469/1‑150 (MQ=255)
                                                                                                             ggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGcc                                        <  2:3380795/150‑1 (MQ=255)
                                                                                                             ggCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGcc                                        <  2:4097481/150‑1 (MQ=255)
                                                                                                              gCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCg                                       <  1:4247415/150‑1 (MQ=255)
                                                                                                                tttGATGATGAGGGCAATTTGACCAAAGAACCTGGTGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGaa                                     >  2:1841496/1‑150 (MQ=255)
                                                                                                                  tgatgatgaGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGt                                   <  1:1442723/150‑1 (MQ=255)
                                                                                                                     tgatgaGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCAc                                <  2:3173449/150‑1 (MQ=255)
                                                                                                                       atgaGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCAccc                              >  2:4772360/1‑150 (MQ=255)
                                                                                                                        tgaGGGCAATGTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCa                             <  2:4666546/150‑1 (MQ=255)
                                                                                                                         gaGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAg                            <  1:4840269/150‑1 (MQ=255)
                                                                                                                          aGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAgg                           >  2:3955059/1‑150 (MQ=255)
                                                                                                                           gggCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGa                          >  2:2409691/1‑150 (MQ=255)
                                                                                                                             gCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGaca                        <  1:963459/150‑1 (MQ=255)
                                                                                                                             gCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGaca                        <  1:4337418/150‑1 (MQ=255)
                                                                                                                              cAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaa                       <  1:2005783/150‑1 (MQ=255)
                                                                                                                               aaTTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaac                      <  1:4998005/150‑1 (MQ=255)
                                                                                                                               aaTTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTAGCCGAAGTCACCCAGGacaac                      <  2:681070/150‑1 (MQ=255)
                                                                                                                                aTTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaaca                     <  2:2266264/150‑1 (MQ=255)
                                                                                                                                aTTTGACCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGacaaca                     <  2:1568185/150‑1 (MQ=255)
                                                                                                                                     aCCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAATGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcgac                <  1:3204186/150‑1 (MQ=255)
                                                                                                                                     aCCAAAGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGcgac                <  2:3842009/150‑1 (MQ=255)
                                                                                                                                        aaaGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGaa             >  2:3756227/1‑150 (MQ=255)
                                                                                                                                          aGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATa           <  1:3407763/150‑1 (MQ=255)
                                                                                                                                          aGAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATa           >  2:4514406/1‑150 (MQ=255)
                                                                                                                                           gAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATAc          <  2:1177808/150‑1 (MQ=255)
                                                                                                                                           gAACCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATAc          <  2:4063361/150‑1 (MQ=255)
                                                                                                                                            aaCCTGGAGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACg         <  1:1345933/150‑1 (MQ=255)
                                                                                                                                                   cGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTc  >  2:2018002/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTCCGACGGCGCATCCGTTGCCGAAGTCACCCAGGACAACAGCGACGAATACGGCATTTC  >  CP009273/3736575‑3736871

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: