Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 70,465 G→T 43.4% H182N (CAT→AAT)  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,4650GT43.4% 203.4 / 277.7 228H182N (CAT→AAT) thiPthiamine, thiamine pyrophosphate ABC transporter permease; ThiBPQ thiamine, thiamine pyrophosphate importer
Reads supporting (aligned to +/- strand):  ref base G (63/66);  new base T (50/49);  total (113/115)
Fisher's exact test for biased strand distribution p-value = 8.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 230 total aligned reads displayed.

TAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCC  >  CP009273/70317‑70614
                                                                                                                                                    |                                                                                                                                                     
tAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGc                                                                                                                                                      >  2:3649915/1‑150 (MQ=255)
 aGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGcc                                                                                                                                                     <  1:2643460/150‑1 (MQ=255)
 aGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGcc                                                                                                                                                     >  2:1590146/1‑150 (MQ=255)
 aGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGcc                                                                                                                                                     <  2:2333710/150‑1 (MQ=255)
  gATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCa                                                                                                                                                    <  1:819744/150‑1 (MQ=255)
  gATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCa                                                                                                                                                    <  1:434733/150‑1 (MQ=255)
  gATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCa                                                                                                                                                    <  1:433285/150‑1 (MQ=255)
  gATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCa                                                                                                                                                    <  1:1263002/150‑1 (MQ=255)
   aTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAg                                                                                                                                                   >  2:1046760/1‑150 (MQ=255)
    ttGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCCGCGATAGCACGGGGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGc                                                                                                                                                  <  2:2496073/150‑1 (MQ=255)
    ttGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGc                                                                                                                                                  >  2:4417769/1‑150 (MQ=255)
      gCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTa                                                                                                                                                <  1:2243701/150‑1 (MQ=255)
         aGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGc                                                                                                                                             <  2:220485/150‑1 (MQ=255)
           cTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCAt                                                                                                                                           >  1:738054/1‑150 (MQ=255)
           cTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCAt                                                                                                                                           <  1:689994/150‑1 (MQ=255)
           cTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCAt                                                                                                                                           <  1:1391034/150‑1 (MQ=255)
             cGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATcc                                                                                                                                         >  1:4509103/1‑150 (MQ=255)
             cGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATcc                                                                                                                                         <  2:1107613/150‑1 (MQ=255)
              gATAGTGGTCGCCTGCGGACCGCCCCACAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATccc                                                                                                                                        <  1:3656814/150‑1 (MQ=255)
               aTAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCa                                                                                                                                       <  1:164123/150‑1 (MQ=255)
                 aGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAg                                                                                                                                     >  2:1256691/1‑150 (MQ=255)
                      tCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTggg                                                                                                                                >  2:605994/1‑150 (MQ=255)
                      tCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTggg                                                                                                                                >  2:233995/1‑150 (MQ=255)
                         ccTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGggcgg                                                                                                                             >  1:2388399/1‑150 (MQ=255)
                         ccTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGggcgg                                                                                                                             >  1:83680/1‑150 (MQ=255)
                               gACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCACCCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGtt                                                                                                                       >  2:2732335/1‑150 (MQ=255)
                               gACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGtt                                                                                                                       <  2:3268686/150‑1 (MQ=255)
                               gACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGtt                                                                                                                       <  2:3197875/150‑1 (MQ=255)
                                aCCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTg                                                                                                                      <  1:1214022/150‑1 (MQ=255)
                                        ccAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGtt                                                                                                              >  1:602067/1‑150 (MQ=255)
                                         cAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTc                                                                                                             <  2:3008565/150‑1 (MQ=255)
                                            cGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGcc                                                                                                          >  2:2303020/1‑150 (MQ=255)
                                             gATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCg                                                                                                         <  1:4138157/150‑1 (MQ=255)
                                                aGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCgggg                                                                                                      <  2:1760164/150‑1 (MQ=255)
                                                     ggTGGCTAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGtt                                                                                                 >  1:3275154/1‑150 (MQ=255)
                                                      gTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGttt                                                                                                <  1:2972780/150‑1 (MQ=255)
                                                      gTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGttt                                                                                                <  1:4009107/150‑1 (MQ=255)
                                                      gTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGttt                                                                                                <  1:2603087/150‑1 (MQ=255)
                                                       tGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGtttt                                                                                               <  2:3186982/150‑1 (MQ=255)
                                                         gCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTcc                                                                                             >  2:3884851/1‑150 (MQ=255)
                                                           gAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAg                                                                                           <  2:2389761/150‑1 (MQ=255)
                                                            aaGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGt                                                                                          <  2:2428450/150‑1 (MQ=255)
                                                            aaGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGt                                                                                          <  2:1317837/150‑1 (MQ=255)
                                                             aGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTg                                                                                         <  1:3859691/150‑1 (MQ=255)
                                                               cTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGcc                                                                                       >  2:4507535/1‑150 (MQ=255)
                                                               cTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGcc                                                                                       >  2:2271584/1‑150 (MQ=255)
                                                                   cGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGa                                                                                   >  1:221740/1‑150 (MQ=255)
                                                                     aaaCAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGt                                                                                 >  2:88095/1‑150 (MQ=255)
                                                                     aaaCAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGt                                                                                 <  1:3445638/150‑1 (MQ=255)
                                                                     aaaCAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGt                                                                                 >  2:1982984/1‑150 (MQ=255)
                                                                        cAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGtaat                                                                              <  2:792905/150‑1 (MQ=255)
                                                                         agagCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGtaata                                                                             <  2:3711156/150‑1 (MQ=255)
                                                                         agagCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGtaata                                                                             <  1:1733903/150‑1 (MQ=255)
                                                                         agagCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGtaata                                                                             <  2:754201/150‑1 (MQ=255)
                                                                          gagCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGtaataa                                                                            >  2:2771499/1‑150 (MQ=255)
                                                                           agCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAg                                                                           >  2:3470687/1‑150 (MQ=255)
                                                                           agCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAg                                                                           >  2:3470426/1‑150 (MQ=255)
                                                                           agCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAg                                                                           >  1:2820506/1‑150 (MQ=255)
                                                                            gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGc                                                                          >  1:2285711/1‑150 (MQ=255)
                                                                             cATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCg                                                                         >  1:4521781/1‑150 (MQ=255)
                                                                                 aaGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTc                                                                     <  2:2883909/150‑1 (MQ=255)
                                                                                 aaGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTc                                                                     >  2:3582299/1‑150 (MQ=255)
                                                                                 aaGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTc                                                                     >  1:2061993/1‑150 (MQ=255)
                                                                                  aGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCg                                                                    <  1:2006212/150‑1 (MQ=255)
                                                                                  aGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCg                                                                    <  1:1719300/150‑1 (MQ=255)
                                                                                     tAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCa                                                                 <  2:3756962/150‑1 (MQ=255)
                                                                                       aGCGCAGCAACCGGCGTGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAATTTGACGCTGTTCGCCGGGGTTGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATc                                                               <  1:1175621/150‑1 (MQ=255)
                                                                                        gcgcAGGAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCg                                                              <  2:984754/150‑1 (MQ=255)
                                                                                        gcgcAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCg                                                              <  1:1966141/150‑1 (MQ=255)
                                                                                          gcagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGc                                                            >  1:2506662/1‑150 (MQ=255)
                                                                                            agcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAg                                                          <  2:1486638/150‑1 (MQ=255)
                                                                                             gcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGa                                                         <  2:2815816/150‑1 (MQ=255)
                                                                                              caACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGAt                                                        <  2:4172770/150‑1 (MQ=255)
                                                                                               aaCCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGAtt                                                       <  1:2625592/150‑1 (MQ=255)
                                                                                                    gcggGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTaaaaa                                                  <  1:1880544/150‑1 (MQ=255)
                                                                                                     cggGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTaaaaaa                                                 <  1:478529/150‑1 (MQ=255)
                                                                                                      gggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAac                                                <  2:2402694/150‑1 (MQ=255)
                                                                                                      gggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAac                                                >  2:428828/1‑150 (MQ=255)
                                                                                                         aTTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                             >  2:2461385/1‑150 (MQ=255)
                                                                                                         aTTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAAcaca                                             <  2:3979238/150‑1 (MQ=255)
                                                                                                               cGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGcc                                       >  1:1350387/1‑150 (MQ=255)
                                                                                                                gCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCa                                      <  1:2881893/150‑1 (MQ=255)
                                                                                                                gCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCa                                      <  2:4093970/150‑1 (MQ=255)
                                                                                                                  cGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGCTGTTTTCCAGTGCCTGGAGTAATAATCGGCTCACCATCGGCAGATTAAAAAACACATGGGCCAga                                    >  2:631139/1‑149 (MQ=255)
                                                                                                                  cGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGc                                    >  2:631162/1‑150 (MQ=255)
                                                                                                                     aaCCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCaaa                                 >  2:69952/1‑150 (MQ=255)
                                                                                                                         aCGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATAc                             >  1:1743940/1‑150 (MQ=255)
                                                                                                                           ggCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCt                           >  2:910896/1‑150 (MQ=255)
                                                                                                                            gCCATTCGACGAAGCGGAAAAATTGCCAGATACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCtt                          <  1:2142296/150‑1 (MQ=255)
                                                                                                                             ccATTCGACGAAGCGGAAATAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTg                         <  2:2860964/150‑1 (MQ=255)
                                                                                                                             ccATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTg                         <  1:1044842/150‑1 (MQ=255)
                                                                                                                             ccATTCGAAGAAGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGACGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTg                         <  2:2861349/150‑1 (MQ=255)
                                                                                                                                ttCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAg                      >  2:1029152/1‑150 (MQ=255)
                                                                                                                                       aaGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAg               >  2:1906750/1‑150 (MQ=255)
                                                                                                                                        aGCGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAgg              <  1:2270947/150‑1 (MQ=255)
                                                                                                                                          cGGAAAAAATTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGc            >  1:1701705/1‑150 (MQ=255)
                                                                                                                                                aaaTTCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAg      >  1:1176299/1‑150 (MQ=255)
                                                                                                                                                aaaTGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAg      >  2:1456981/1‑150 (MQ=255)
                                                                                                                                                   tGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTc   <  1:3065994/150‑1 (MQ=255)
                                                                                                                                                    gCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTcc  <  2:1147947/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCC  >  CP009273/70317‑70614

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: