Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 70,249 A→C 8.4% C254G (TGC→GGC)  thiP ← thiamine, thiamine pyrophosphate ABC transporter permease; ThiBPQ thiamine, thiamine pyrophosphate importer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927370,2490AC 8.4% 645.1 / 51.2 227C254G (TGC→GGC) thiPthiamine, thiamine pyrophosphate ABC transporter permease; ThiBPQ thiamine, thiamine pyrophosphate importer
Reads supporting (aligned to +/- strand):  ref base A (113/95);  new base C (11/8);  total (124/103)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 239 total aligned reads displayed.

CAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAA  >  CP009273/70100‑70398
                                                                                                                                                     |                                                                                                                                                     
   caACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGgcagca                                                                                                                                                    >  1:3008565/1‑150 (MQ=255)
     aCAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGcaca                                                                                                                                                  >  2:2424061/1‑150 (MQ=255)
       agcagcagcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACAcc                                                                                                                                                >  2:185946/1‑150 (MQ=255)
         cagcagcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCCGCACACCAt                                                                                                                                              <  1:3030221/150‑1 (MQ=255)
         cagcagcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCAt                                                                                                                                              <  1:1654141/150‑1 (MQ=255)
            cagcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTg                                                                                                                                           <  2:3450195/150‑1 (MQ=255)
             agcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTgg                                                                                                                                          <  2:1491031/150‑1 (MQ=255)
               cgcCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCCGCACACCATCTGGag                                                                                                                                        <  2:2952348/150‑1 (MQ=255)
                 ccAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGACGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGagca                                                                                                                                      >  2:1175621/1‑150 (MQ=255)
                   aGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGagcagc                                                                                                                                    >  1:1107613/1‑150 (MQ=255)
                   aGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGagcagc                                                                                                                                    <  1:379481/150‑1 (MQ=255)
                   aGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGagcagc                                                                                                                                    <  2:3567762/150‑1 (MQ=255)
                    gCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAgcg                                                                                                                                   >  2:2006212/1‑150 (MQ=255)
                    gCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAgcg                                                                                                                                   >  2:1719368/1‑150 (MQ=255)
                    gCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAgcg                                                                                                                                   >  2:1719299/1‑150 (MQ=255)
                      acaATTAACACCGTGTCGCTAATGCGGCTATGCAGACGATCGTGCGGGTCGAGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGAAATGGCATTACTCAATCTCTGACTCAACTGCACCAGCCCGTGGCAGCACTCCATCTGGAGCAGCGcc                                                                                                                                 <  2:4167938/150‑1 (MQ=255)
                        aaTTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAg                                                                                                                               <  2:3467821/150‑1 (MQ=255)
                            aaCACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCAtt                                                                                                                           <  1:2625882/150‑1 (MQ=255)
                                ccGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCcg                                                                                                                       >  1:2815816/1‑150 (MQ=255)
                                 cGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCCGCACACCATCTGGAGCAGCGCCAGCATTGCcgc                                                                                                                      >  1:4539511/1‑150 (MQ=255)
                                 cGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCcgc                                                                                                                      >  2:4115118/1‑150 (MQ=255)
                                  gtgtCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCcgcg                                                                                                                     >  2:3765355/1‑150 (MQ=255)
                                      cGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCggg                                                                                                                 <  2:1155029/150‑1 (MQ=255)
                                      cGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCggg                                                                                                                 >  2:1214022/1‑150 (MQ=255)
                                        cAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCa                                                                                                               >  2:689994/1‑150 (MQ=255)
                                         aaaTGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAg                                                                                                              >  2:3311700/1‑150 (MQ=255)
                                         aaaTGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAg                                                                                                              >  2:3518165/1‑150 (MQ=255)
                                         aaaTGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAg                                                                                                              >  2:3518021/1‑150 (MQ=255)
                                              cGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCg                                                                                                         <  2:3365297/150‑1 (MQ=255)
                                               ggCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCCTCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGAGCGt                                                                                                        >  2:2881893/1‑150 (MQ=255)
                                               ggCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGt                                                                                                        >  1:4093970/1‑150 (MQ=255)
                                                 cTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAg                                                                                                      <  2:1306949/150‑1 (MQ=255)
                                                 cTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAg                                                                                                      >  1:974551/1‑150 (MQ=255)
                                                 cTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAg                                                                                                      >  2:2243701/1‑150 (MQ=255)
                                                  tATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCCGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGt                                                                                                     >  2:2625592/1‑150 (MQ=255)
                                                   aTGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTc                                                                                                    <  2:2010068/150‑1 (MQ=255)
                                                      cAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTa                                                                                                 <  1:3620571/150‑1 (MQ=255)
                                                      cAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTa                                                                                                 <  2:3683715/150‑1 (MQ=255)
                                                      cAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTa                                                                                                 >  1:1366126/1‑150 (MQ=255)
                                                       aGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTaa                                                                                                <  1:1527304/150‑1 (MQ=255)
                                                       aGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTaa                                                                                                >  2:2603092/1‑150 (MQ=255)
                                                         aCGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACt                                                                                              >  2:72549/1‑150 (MQ=255)
                                                           gATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCa                                                                                            >  2:2271097/1‑150 (MQ=255)
                                                            aTCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAg                                                                                           >  1:3268577/1‑150 (MQ=255)
                                                            aTCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAg                                                                                           >  1:991840/1‑150 (MQ=255)
                                                             tCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAgc                                                                                          >  1:1100101/1‑150 (MQ=255)
                                                              cGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAgcg                                                                                         >  2:1391034/1‑150 (MQ=255)
                                                                  cGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCNg                                                                                     <  1:420746/150‑1 (MQ=255)
                                                                    ggTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGAt                                                                                   <  1:695961/150‑1 (MQ=255)
                                                                     gTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATa                                                                                  <  2:3121217/150‑1 (MQ=255)
                                                                     gTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATa                                                                                  <  2:3595870/150‑1 (MQ=255)
                                                                       cgcgCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGa                                                                                >  1:2883378/1‑150 (MQ=255)
                                                                       cgcgCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGa                                                                                >  1:642623/1‑150 (MQ=255)
                                                                         cgcCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGAtt                                                                              >  1:2511981/1‑150 (MQ=255)
                                                                             aGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCa                                                                          <  2:1051443/150‑1 (MQ=255)
                                                                             aGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCa                                                                          >  2:739021/1‑150 (MQ=255)
                                                                                cTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTAATCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTATTCGTAACTCAGCGCCTGATAGATTGCCAGcg                                                                       >  2:4408738/1‑149 (MQ=255)
                                                                                    cagcagCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGAt                                                                   >  2:1262984/1‑150 (MQ=255)
                                                                                     agcagcGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATa                                                                  <  2:1365924/150‑1 (MQ=255)
                                                                                        agcGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCCTCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTg                                                               <  1:3972355/150‑1 (MQ=255)
                                                                                          cGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGt                                                             >  2:2485633/1‑150 (MQ=255)
                                                                                          cGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGt                                                             <  1:187180/150‑1 (MQ=255)
                                                                                           gtggtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTc                                                            >  1:3197577/1‑150 (MQ=255)
                                                                                            tggtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCg                                                           >  2:305823/1‑150 (MQ=255)
                                                                                            tggtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCg                                                           >  1:1148634/1‑150 (MQ=255)
                                                                                             ggtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCCGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGc                                                          >  2:434733/1‑150 (MQ=255)
                                                                                              gtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCCGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGcc                                                         >  1:754201/1‑150 (MQ=255)
                                                                                              gtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGcc                                                         >  2:1966110/1‑150 (MQ=255)
                                                                                              gtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGcc                                                         >  2:1966923/1‑150 (MQ=255)
                                                                                               tgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCt                                                        <  1:116573/150‑1 (MQ=255)
                                                                                               tgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCt                                                        >  2:4009107/1‑150 (MQ=255)
                                                                                               tgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCt                                                        >  1:4383535/1‑150 (MQ=255)
                                                                                                     gcgcAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGAc                                                  <  1:2791627/150‑1 (MQ=255)
                                                                                                            ggCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGccccc                                           >  1:206733/1‑150 (MQ=255)
                                                                                                             gCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGcccccc                                          <  2:364732/150‑1 (MQ=255)
                                                                                                             gCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGcccccc                                          <  2:359888/150‑1 (MQ=255)
                                                                                                             gCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCAGTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGcccccc                                          <  2:3281276/150‑1 (MQ=255)
                                                                                                              ccTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                         <  1:2787710/150‑1 (MQ=255)
                                                                                                                ttACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGc                                       <  1:4457109/150‑1 (MQ=255)
                                                                                                                ttACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGc                                       >  1:1486558/1‑150 (MQ=255)
                                                                                                                    tCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATa                                   >  2:1525444/1‑150 (MQ=255)
                                                                                                                         cGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACg                              <  2:936902/150‑1 (MQ=255)
                                                                                                                           cTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGt                            <  1:3263498/150‑1 (MQ=255)
                                                                                                                            tGACTCAACAGCACCAGCCCGAGGCCGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTg                           >  1:2025432/1‑150 (MQ=255)
                                                                                                                             gACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTgg                          >  2:2398024/1‑150 (MQ=255)
                                                                                                                             gACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTgg                          >  1:2242442/1‑150 (MQ=255)
                                                                                                                               cTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCg                        <  2:87546/150‑1 (MQ=255)
                                                                                                                               cTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCg                        <  1:2424061/150‑1 (MQ=255)
                                                                                                                               cTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATAGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCg                        <  1:3765355/150‑1 (MQ=255)
                                                                                                                                       cACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGc                <  1:2446382/150‑1 (MQ=255)
                                                                                                                                           aGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGaaa            >  1:3809452/1‑150 (MQ=255)
                                                                                                                                             cccGAGGCCGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACa          <  1:2344031/150‑1 (MQ=255)
                                                                                                                                             cccGAGGCCGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACa          <  1:39809/150‑1 (MQ=255)
                                                                                                                                               cGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACaga        >  2:2737377/1‑150 (MQ=255)
                                                                                                                                                gAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACagag       <  1:2485633/150‑1 (MQ=255)
                                                                                                                                                gAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACagag       <  1:2722158/150‑1 (MQ=255)
                                                                                                                                                 aGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGc      <  1:1339398/150‑1 (MQ=255)
                                                                                                                                                 aGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGc      >  1:2861336/1‑150 (MQ=255)
                                                                                                                                                  ggCCGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCa     <  2:2313926/150‑1 (MQ=255)
                                                                                                                                                    cagcaCACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATa   <  1:185946/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAA  >  CP009273/70100‑70398

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: