Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 69,546 A→T 16.4% F488Y (TTC→TAC)  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
*CP00927369,5460AT16.4% 785.4 / 132.3 285F488Y (TTC→TAC) thiPthiamine, thiamine pyrophosphate ABC transporter permease; ThiBPQ thiamine, thiamine pyrophosphate importer
Reads supporting (aligned to +/- strand):  ref base A (128/110);  new base T (17/30);  total (145/140)
Fisher's exact test for biased strand distribution p-value = 3.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

TCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACAT  >  CP009273/69398‑69695
                                                                                                                                                    |                                                                                                                                                     
tcagtcagTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGACAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGaa                                                                                                                                                      >  1:531977/1‑150 (MQ=255)
    tcagtTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAatca                                                                                                                                                  <  2:1841558/150‑1 (MQ=255)
      atttttAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATCGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAatcatc                                                                                                                                                <  1:4460789/148‑1 (MQ=255)
      agtTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGTAatcatc                                                                                                                                                <  1:4579382/150‑1 (MQ=255)
      agtTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAatcatc                                                                                                                                                >  2:1344315/1‑150 (MQ=255)
       gtTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCg                                                                                                                                               <  2:49823/150‑1 (MQ=255)
         tttAACATTTCGCCCCGGTTTTTTTACAATCCCGGTAAACAGCAGAAAACCGAGCTGCCGAAGAATTACCGCGGTTACCGCACCGTCCTGGCTGCGATTGGAGCAATTTTGCTGGTTGAGATAACACGGCAGGGTGCGGAAATCATCGtt                                                                                                                                             <  1:1951479/150‑1 (MQ=255)
          ttAACATTTCGCCCCGTTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGTAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTa                                                                                                                                            <  2:108342/150‑1 (MQ=255)
             aCATTTCGCCCTGGTTGTTTTTCAATCACGGTAAACACCATAAAACAGAGCAGCAGCATAATTACCGCGGTGACCGCACCGTCCTGGCTGCTATTGGATAAAATTTGCTGGTATTTATAAAACGGCAGGGTGCGGAAATCATCGTTACCg                                                                                                                                         <  1:3551385/150‑1 (MQ=255)
              cATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGTAATCATCGTTACCGa                                                                                                                                        <  2:4344851/150‑1 (MQ=255)
                 ttCGCCCCGGTTGTTTTTAAATCACGGTAAACAGCAGAAAACAGAGCCGCAGCAGATTTACCGCGGTGACCGCACCGTCCTGGCTGCTTTAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGaaca                                                                                                                                     <  1:3621774/150‑1 (MQ=255)
                  tCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGaacaa                                                                                                                                    >  2:4472921/1‑150 (MQ=255)
                   cGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGTTAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGaacaac                                                                                                                                   <  2:1270383/150‑1 (MQ=255)
                   cGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGaacaac                                                                                                                                   >  2:752260/1‑150 (MQ=255)
                    gCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGCCCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGCGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACg                                                                                                                                  <  2:631085/150‑1 (MQ=255)
                    gCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACg                                                                                                                                  >  1:1009558/1‑150 (MQ=255)
                      cccGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGcc                                                                                                                                >  1:2420176/1‑150 (MQ=255)
                       ccGGTAGATTTTCAATAACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCGTGGCTGCGATAGGAGCCAATTTGCTGGTAGTGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGcca                                                                                                                               <  2:2903778/150‑1 (MQ=255)
                            aGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACg                                                                                                                          >  1:2629781/1‑150 (MQ=255)
                                tttCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCATAATTAACGCGGTGACCGCACCGTCCGGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCaa                                                                                                                      <  1:2443946/150‑1 (MQ=255)
                                  tCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCaaaa                                                                                                                    >  2:873169/1‑150 (MQ=255)
                                   cAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAAt                                                                                                                   >  2:458591/1‑150 (MQ=255)
                                       cACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATATGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCAcc                                                                                                               >  1:3792981/1‑150 (MQ=255)
                                       cACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCAcc                                                                                                               >  2:2662601/1‑150 (MQ=255)
                                       cACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCAcc                                                                                                               >  1:1415299/1‑150 (MQ=255)
                                           gTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATc                                                                                                           >  1:918931/1‑150 (MQ=255)
                                               aCAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACa                                                                                                       <  2:4222758/150‑1 (MQ=255)
                                                cagcagAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAg                                                                                                      >  2:3173148/1‑150 (MQ=255)
                                                   aataaaaCAGAGCAGCAGCATAATTACCGGGGTTACCGCACCGTCCTGGCTGCGATAGGAGCCATTTTGCTGGTAGTTATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCAc                                                                                                   <  1:3688886/147‑1 (MQ=255)
                                                    agAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACg                                                                                                  >  2:3220360/1‑150 (MQ=255)
                                                      aaaaCAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCa                                                                                                >  1:3449065/1‑150 (MQ=255)
                                                      aaaaCAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCa                                                                                                <  1:232418/150‑1 (MQ=255)
                                                      aaaaCAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCa                                                                                                >  1:3942487/1‑150 (MQ=255)
                                                       aaaCTGAGCAGCAGCATAATTAACGCGGTGACCGCACCTTCCTGGCTGCTATTGTAGAAAATTTGCTGGTAGATATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAt                                                                                               <  1:1638331/150‑1 (MQ=255)
                                                       aaaCAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAt                                                                                               >  2:3030104/1‑150 (MQ=255)
                                                       aaaCAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAt                                                                                               <  2:2101349/150‑1 (MQ=255)
                                                       aaaCAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAt                                                                                               >  2:387375/1‑150 (MQ=255)
                                                          cAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCa                                                                                            >  1:41047/1‑150 (MQ=255)
                                                              gcagcagcaGAATTAACGCGGTGACCGCACCGTCCTGTCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCAACAATCGACAGCACGCATGCAAAgg                                                                                        >  1:3351876/1‑150 (MQ=255)
                                                              gcagcagcaGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAgg                                                                                        <  2:1357455/150‑1 (MQ=255)
                                                               cagcagcagAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGc                                                                                       <  2:2035384/150‑1 (MQ=255)
                                                                agcagcagAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGcc                                                                                      <  1:3825052/150‑1 (MQ=255)
                                                                agcagcagAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGcc                                                                                      >  2:2737912/1‑150 (MQ=255)
                                                                  cagcagAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAg                                                                                    >  1:4025459/1‑150 (MQ=255)
                                                                       gAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCt                                                                               <  2:4392515/150‑1 (MQ=255)
                                                                          ttACCGGGGTGACCGCACGGTCCTGGCTGCGTTCGGAGCAAATTGGCGGTTAGATATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTgcg                                                                            <  1:4388574/150‑1 (MQ=255)
                                                                                 ggTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTgg                                                                     >  1:1334878/1‑150 (MQ=255)
                                                                                 ggTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTgg                                                                     >  2:3936261/1‑150 (MQ=255)
                                                                                 ggTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTgg                                                                     >  1:4163284/1‑150 (MQ=255)
                                                                                   tGACCGCACCGTCCTGGCTGCGATAGGAGCCATTTTGCTGGTAGATATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGAc                                                                   <  1:3478053/150‑1 (MQ=255)
                                                                                    gACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACg                                                                  >  1:2744299/1‑150 (MQ=255)
                                                                                     aCCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGt                                                                 >  2:4235869/1‑150 (MQ=255)
                                                                                      ccGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGtt                                                                <  1:2168475/150‑1 (MQ=255)
                                                                                         cACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCa                                                             >  2:3772448/1‑150 (MQ=255)
                                                                                         cACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCa                                                             >  1:791648/1‑150 (MQ=255)
                                                                                            cGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAggg                                                          >  1:3530758/1‑150 (MQ=255)
                                                                                             gTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGgc                                                         >  1:121612/1‑150 (MQ=255)
                                                                                                 tGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGgcgcgc                                                     >  2:648963/1‑150 (MQ=255)
                                                                                                 tGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGgcgcgc                                                     >  1:3415493/1‑150 (MQ=255)
                                                                                                   gCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAg                                                   <  1:1293943/150‑1 (MQ=255)
                                                                                                     tGCGATAGGAGCCAATTTGCTGGTAGAGATAACACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCt                                                 <  1:2924840/150‑1 (MQ=255)
                                                                                                      gCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTc                                                >  2:16582/1‑150 (MQ=255)
                                                                                                        gATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTcca                                              <  2:14134/150‑1 (MQ=255)
                                                                                                           aGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTccacca                                           >  1:3042686/1‑150 (MQ=255)
                                                                                                           aGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTccacca                                           >  2:3372525/1‑150 (MQ=255)
                                                                                                                ccAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACtttt                                      >  1:1610011/1‑150 (MQ=255)
                                                                                                                  aaTTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTaa                                    >  1:1017312/1‑150 (MQ=255)
                                                                                                                  aaTTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTaa                                    <  2:2501798/150‑1 (MQ=255)
                                                                                                                  aaTTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTaa                                    >  2:14903/1‑150 (MQ=255)
                                                                                                                    tttGCTGGTAGAGATAAAACGGCAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGc                                  <  1:2605775/150‑1 (MQ=255)
                                                                                                                     ttGCTGTTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCg                                 <  1:2496054/150‑1 (MQ=255)
                                                                                                                      tGCTGGTAGAGATAAAACGGCAGCGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGt                                <  2:2201872/150‑1 (MQ=255)
                                                                                                                       gCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTg                               <  2:2629781/150‑1 (MQ=255)
                                                                                                                        cTGGTAGAGATAAAACGGCAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGa                              >  2:4107457/1‑150 (MQ=255)
                                                                                                                         tGGTAGAGATAAAACGGCAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGAc                             <  2:2822467/150‑1 (MQ=255)
                                                                                                                         tGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGAc                             >  2:2583140/1‑150 (MQ=255)
                                                                                                                            tAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAg                          >  2:3613272/1‑150 (MQ=255)
                                                                                                                               agaTAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCt                       >  1:1918689/1‑150 (MQ=255)
                                                                                                                                gaTAAAACGGCAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCtt                      <  2:2426802/150‑1 (MQ=255)
                                                                                                                                gaTAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCtt                      >  1:126315/1‑150 (MQ=255)
                                                                                                                                gaTAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCtt                      <  1:1856543/150‑1 (MQ=255)
                                                                                                                                 aTAAAACGGCAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTc                     <  2:741681/150‑1 (MQ=255)
                                                                                                                                  tAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCa                    <  2:1304265/150‑1 (MQ=255)
                                                                                                                                   aaaaCGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCaa                   >  2:3318006/1‑150 (MQ=255)
                                                                                                                                    aaaCGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAAt                  <  2:3543544/150‑1 (MQ=255)
                                                                                                                                       cGGGAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCTCGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATccc               <  2:531977/150‑1 (MQ=255)
                                                                                                                                       cGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATccc               <  1:714760/150‑1 (MQ=255)
                                                                                                                                       cGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATccc               <  2:3161489/150‑1 (MQ=255)
                                                                                                                                       cGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATccc               <  2:4559883/150‑1 (MQ=255)
                                                                                                                                        tgCAGTGTGCGGATATCATATTTTACTATCAACGCCACCACGGCATTATCACAAATCGAGAGCACGAATGCAATGGCCAGCGCCTGCGCCAGTGGACTTTTCAGTGCGAGCAGCTCCACCACTTTTATGAGTGACCATCCTTCAATcccc              <  2:2704503/149‑1 (MQ=255)
                                                                                                                                        ggCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATcccc              >  1:3092361/1‑150 (MQ=255)
                                                                                                                                        ggCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATcccc              <  1:1119469/150‑1 (MQ=255)
                                                                                                                                        ggCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATcccc              <  1:4227531/150‑1 (MQ=255)
                                                                                                                                          cAGGGTGCGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAg            <  2:1915325/150‑1 (MQ=255)
                                                                                                                                          cAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAg            <  2:2393507/150‑1 (MQ=255)
                                                                                                                                            gggTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCg          <  2:1133182/150‑1 (MQ=255)
                                                                                                                                               tGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACt       <  2:161087/150‑1 (MQ=255)
                                                                                                                                                gCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTg      >  2:154633/1‑150 (MQ=255)
                                                                                                                                                 cGGTAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGa     <  2:1123384/150‑1 (MQ=255)
                                                                                                                                                    aaaTCATCGTTACCGAACAACGCCACCACGCCAAAATCACCACTCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACAt  <  1:3254468/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACAT  >  CP009273/69398‑69695

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑