Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 70,744 G→T 7.0% R89R (CGG→AGG)  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,7440GT 7.0% 1009.9 / 31.7 264R89R (CGG→AGG) thiPthiamine, thiamine pyrophosphate ABC transporter permease; ThiBPQ thiamine, thiamine pyrophosphate importer
Reads supporting (aligned to +/- strand):  ref base G (119/126);  new base T (11/8);  total (130/134)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

TAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCG  >  CP009273/70599‑70893
                                                                                                                                                 |                                                                                                                                                     
tAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGAccc                                                                                                                                                   >  2:3218373/1‑150 (MQ=255)
tAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGAccc                                                                                                                                                   >  2:3278665/1‑150 (MQ=255)
 aGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAGGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATGGCACACAGACGCAACAGCGCCAGCCTACCCg                                                                                                                                                  >  1:1595571/1‑150 (MQ=255)
 aGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCTCACAGACGCAACAGCGCCAGCCGACCCg                                                                                                                                                  <  1:4444527/150‑1 (MQ=255)
  gggCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCgg                                                                                                                                                 >  2:2510868/1‑150 (MQ=255)
      gAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAg                                                                                                                                             <  1:4370599/150‑1 (MQ=255)
      gAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAg                                                                                                                                             >  1:759918/1‑150 (MQ=255)
      gAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAg                                                                                                                                             >  1:277387/1‑150 (MQ=255)
      gAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAg                                                                                                                                             >  1:2392508/1‑150 (MQ=255)
      gAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAg                                                                                                                                             >  2:3469710/1‑150 (MQ=255)
           ggTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCt                                                                                                                                        <  2:1537662/150‑1 (MQ=255)
             tCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGc                                                                                                                                      <  2:2451705/150‑1 (MQ=255)
                aCTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCTACCCGGAAAGCGCCTGCGAt                                                                                                                                   >  1:1575402/1‑150 (MQ=255)
                aCTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGAc                                                                                                                                                     >  2:1444091/1‑132 (MQ=255)
                aCTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGAt                                                                                                                                   >  2:1341676/1‑150 (MQ=255)
                aCTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGAt                                                                                                                                   <  2:3930163/150‑1 (MQ=255)
                  tCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCTACCCGGAAAGCGCCTGCGATag                                                                                                                                 >  2:1030147/1‑150 (MQ=255)
                     aGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAgc                                                                                                                              >  1:2611560/1‑150 (MQ=255)
                          gAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGcgcgcg                                                                                                                         <  1:3269658/150‑1 (MQ=255)
                           aGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCgg                                                                                                                        <  2:3271833/150‑1 (MQ=255)
                           aGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCgg                                                                                                                        >  1:1043070/1‑150 (MQ=255)
                            gCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCggg                                                                                                                       <  1:1384633/150‑1 (MQ=255)
                              gATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCg                                                                                                                     <  2:1110924/150‑1 (MQ=255)
                                 tGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAgg                                                                                                                  >  1:132056/1‑150 (MQ=255)
                                    cAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAat                                                                                                               <  2:334172/150‑1 (MQ=255)
                                     agagTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAACACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCTACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAata                                                                                                              <  2:313474/150‑1 (MQ=255)
                                     agagTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAata                                                                                                              <  1:2954937/150‑1 (MQ=255)
                                       agTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATc                                                                                                            <  1:3407775/150‑1 (MQ=255)
                                       agTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATc                                                                                                            <  2:3885696/150‑1 (MQ=255)
                                       agTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATc                                                                                                            <  1:4623721/150‑1 (MQ=255)
                                        gtgtTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATcg                                                                                                           >  1:2741519/1‑150 (MQ=255)
                                         tgtTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATcgc                                                                                                          <  2:141510/150‑1 (MQ=255)
                                          gtTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATcgcg                                                                                                         <  1:961730/150‑1 (MQ=255)
                                          gtTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATcgcg                                                                                                         <  2:118864/150‑1 (MQ=255)
                                            tGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCggg                                                                                                       <  1:1996816/150‑1 (MQ=255)
                                             gccagccaGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGt                                                                                                      <  1:604314/150‑1 (MQ=255)
                                              ccagccagCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTa                                                                                                     <  1:866847/150‑1 (MQ=255)
                                                 gccagccCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGa                                                                                                  <  2:1914318/150‑1 (MQ=255)
                                                     gccCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACaga                                                                                              >  2:3362731/1‑150 (MQ=255)
                                                         tGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGc                                                                                          <  1:439135/150‑1 (MQ=255)
                                                         tGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCCGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGAGCGGGCGAGGACAATCGCGGGTATGACAGAGAGc                                                                                          >  2:3324230/1‑150 (MQ=255)
                                                            cGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGt                                                                                       <  2:16803/150‑1 (MQ=255)
                                                             gACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTg                                                                                      <  2:1788929/150‑1 (MQ=255)
                                                             gACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTg                                                                                      <  1:1636770/150‑1 (MQ=255)
                                                             gACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTg                                                                                      <  2:2335658/150‑1 (MQ=255)
                                                              aCCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGc                                                                                     >  2:2098278/1‑150 (MQ=255)
                                                               ccATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGcc                                                                                    >  1:738941/1‑150 (MQ=255)
                                                                 aTAGACGCTAAGAATGCCGAAAACAGCGACCAACAGCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGAGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCga                                                                                  <  2:1070226/150‑1 (MQ=255)
                                                                 aTAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCga                                                                                  <  1:3105151/150‑1 (MQ=255)
                                                                     aCGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCgagaga                                                                              <  1:3156825/150‑1 (MQ=255)
                                                                      cGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCTACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGaa                                                                             >  2:1615653/1‑150 (MQ=255)
                                                                        cTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAAc                                                                           <  2:3655036/150‑1 (MQ=255)
                                                                        cTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAAc                                                                           <  1:2676763/150‑1 (MQ=255)
                                                                        cTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAAc                                                                           <  1:2420524/150‑1 (MQ=255)
                                                                         tAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACg                                                                          >  2:4067299/1‑150 (MQ=255)
                                                                           aGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGcc                                                                        >  1:3183187/1‑150 (MQ=255)
                                                                             aaTGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTg                                                                      <  2:2758915/150‑1 (MQ=255)
                                                                             aaTGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTg                                                                      >  1:2393473/1‑150 (MQ=255)
                                                                             aaTGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTg                                                                      >  2:3291568/1‑150 (MQ=255)
                                                                             aaTGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTg                                                                      >  2:3616406/1‑150 (MQ=255)
                                                                                    aaaaCAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAg                                                               >  1:4134130/1‑150 (MQ=255)
                                                                                     aaaCAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAgg                                                              >  2:460604/1‑150 (MQ=255)
                                                                                     aaaCAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAgg                                                              >  2:242112/1‑150 (MQ=255)
                                                                                          gCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGACAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAg                                                         <  2:1725517/150‑1 (MQ=255)
                                                                                            gACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAgcg                                                       >  2:3931552/1‑150 (MQ=255)
                                                                                             aCCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAgcgc                                                      <  2:2062017/150‑1 (MQ=255)
                                                                                              ccAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGca                                                     >  2:943186/1‑150 (MQ=255)
                                                                                               cAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGcac                                                    >  1:2514734/1‑150 (MQ=255)
                                                                                                aaCACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGcacc                                                   >  1:4105589/1‑150 (MQ=255)
                                                                                                  cacCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCTACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGcaccac                                                 <  1:3684782/150‑1 (MQ=255)
                                                                                                   acCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACcaca                                                <  1:3937290/150‑1 (MQ=255)
                                                                                                       ggaggaTCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGcc                                            <  1:1886595/150‑1 (MQ=255)
                                                                                                       ggaggaTCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGcc                                            <  1:3513352/150‑1 (MQ=255)
                                                                                                           gaTCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAg                                        <  2:4224979/150‑1 (MQ=255)
                                                                                                           gaTCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAg                                        <  2:2907321/150‑1 (MQ=255)
                                                                                                            aTCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGa                                       >  1:677825/1‑150 (MQ=255)
                                                                                                             tCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGAt                                      <  2:2983125/150‑1 (MQ=255)
                                                                                                               aaGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAg                                    >  1:137243/1‑150 (MQ=255)
                                                                                                                aGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGACGCGCACCACATGCCACAGATAGc                                   <  1:2179317/150‑1 (MQ=255)
                                                                                                                                aCGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGc                   >  1:3003871/1‑150 (MQ=255)
                                                                                                                                 cGCTACAGCTCCAGGCGTCCGGTAAAGCGCTTTCGATAGAGCGCGAGGGCGAGGAATATCGCGGGTATTACAGAGAGAAGTGCCGAGAGAAAAGCCTGCCAGAAGGAGATGCGCACCTAATGACACAGTTAGCTGTCGTGCCAGACTGcc                  <  2:1308580/150‑1 (MQ=255)
                                                                                                                                   cAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCAc                <  1:499083/150‑1 (MQ=255)
                                                                                                                                    aaCAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCAcc               >  2:2496593/1‑150 (MQ=255)
                                                                                                                                     aCAGCGCCAGCCTACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCAccc              >  1:4255518/1‑150 (MQ=255)
                                                                                                                                       aGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAg            >  2:310527/1‑150 (MQ=255)
                                                                                                                                       aGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAg            >  2:592738/1‑150 (MQ=255)
                                                                                                                                        gcgcCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGt           <  1:3708233/150‑1 (MQ=255)
                                                                                                                                          gcCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCa         >  1:3610368/1‑150 (MQ=255)
                                                                                                                                          gcCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCa         >  1:1031570/1‑150 (MQ=255)
                                                                                                                                            cAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATc       >  1:857308/1‑150 (MQ=255)
                                                                                                                                            cAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATc       >  2:2706770/1‑150 (MQ=255)
                                                                                                                                             aGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATcc      >  1:3823228/1‑150 (MQ=255)
                                                                                                                                             aGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATcc      >  1:3455821/1‑150 (MQ=255)
                                                                                                                                             aGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATcc      >  1:3867059/1‑150 (MQ=255)
                                                                                                                                              gCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCt     <  1:2034938/150‑1 (MQ=255)
                                                                                                                                              gCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCt     <  2:4348615/150‑1 (MQ=255)
                                                                                                                                              gCCGACCCGGAAAGCGCCTGCGAGAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCt     <  2:2021862/150‑1 (MQ=255)
                                                                                                                                               ccGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTg    <  1:1546779/150‑1 (MQ=255)
                                                                                                                                                 gACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTgcg  <  2:4487678/150‑1 (MQ=255)
                                                                                                                                                 gACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTgcg  >  2:1194598/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCGCGCGGGCGAGGAATATCGCGGGTATGACAGAGAGCAGTGCCGAGAGAAACGCCTGCCAGAAGGAGAAGCGCACCACATGCCACAGATAGCTGTCCTGCCAGACTGCCACCCAGTCATCCTGCG  >  CP009273/70599‑70893

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: