Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 69,367 C→A 9.1% M16I (ATG→ATT thiQ ← thiamine, thiamine pyrophosphate ABC transporter ATPase; ThiBPQ thiamine, thiamine pyrophosphate importer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927369,3670CA 9.1% 1332.4 / 66.8 359M16I (ATG→ATTthiQthiamine, thiamine pyrophosphate ABC transporter ATPase; ThiBPQ thiamine, thiamine pyrophosphate importer
Reads supporting (aligned to +/- strand):  ref base C (151/175);  new base A (24/9);  total (175/184)
Fisher's exact test for biased strand distribution p-value = 5.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 365 total aligned reads displayed.

ATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTT  >  CP009273/69218‑69516
                                                                                                                                                     |                                                                                                                                                     
aTAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGc                                                                                                                                                       >  1:4882958/1‑150 (MQ=255)
 tAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCa                                                                                                                                                      >  2:2775176/1‑150 (MQ=255)
  aGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCAt                                                                                                                                                     <  1:2408395/150‑1 (MQ=255)
   gTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGTGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTTCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATc                                                                                                                                                    <  1:3777335/150‑1 (MQ=255)
    ttGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCg                                                                                                                                                   <  1:510225/150‑1 (MQ=255)
    ttGTGTGATCTACGCCATCGAGAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCg                                                                                                                                                   <  2:911090/150‑1 (MQ=255)
      gtgtgATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGc                                                                                                                                                 <  2:2442346/150‑1 (MQ=255)
           aTCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAAtg                                                                                                                                            >  1:1290462/1‑150 (MQ=255)
           aTCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAAtg                                                                                                                                            <  1:3197290/150‑1 (MQ=255)
           aTCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACCCATCGGCAAAtg                                                                                                                                            >  1:4785076/1‑150 (MQ=255)
              tACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAAtggtg                                                                                                                                         <  2:1599860/150‑1 (MQ=255)
              tACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAAtggtg                                                                                                                                         >  2:2793996/1‑150 (MQ=255)
                cGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATggtggt                                                                                                                                       <  1:333514/150‑1 (MQ=255)
                cGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATggtggt                                                                                                                                       <  1:881979/150‑1 (MQ=255)
                 gCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTa                                                                                                                                      <  2:4414158/150‑1 (MQ=255)
                   cATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTaaa                                                                                                                                    >  1:358114/1‑150 (MQ=255)
                    aTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAg                                                                                                                                   >  1:289155/1‑150 (MQ=255)
                     tCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGc                                                                                                                                  <  1:4024496/150‑1 (MQ=255)
                     tCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGc                                                                                                                                  >  2:4044812/1‑150 (MQ=255)
                     tCGAGAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGc                                                                                                                                  <  1:3524749/150‑1 (MQ=255)
                      cGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGcc                                                                                                                                 <  1:2054633/150‑1 (MQ=255)
                      cGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGAATCGGCAAATGGTGGTAAAGcc                                                                                                                                 >  2:1909202/1‑150 (MQ=255)
                       gATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCa                                                                                                                                <  2:206192/150‑1 (MQ=255)
                       gATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCa                                                                                                                                <  2:4002028/150‑1 (MQ=255)
                       gATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCa                                                                                                                                >  1:1047192/1‑150 (MQ=255)
                        aTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAg                                                                                                                               >  2:3132724/1‑150 (MQ=255)
                           gTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTg                                                                                                                            <  1:1588111/150‑1 (MQ=255)
                             cAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGat                                                                                                                          <  1:868786/150‑1 (MQ=255)
                              aGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGata                                                                                                                         <  2:1923018/150‑1 (MQ=255)
                               gCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGAATCGGCAAATGGTGGTAAAGCCAGGTGatat                                                                                                                        <  2:2214164/150‑1 (MQ=255)
                                  aaCCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATAtcag                                                                                                                     >  1:1132433/1‑150 (MQ=255)
                                   aCCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATAtcagt                                                                                                                    >  2:639115/1‑150 (MQ=255)
                                   aCCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGAATCGGCAAATGGTGGTAAAGCCAGGTGATAtcagt                                                                                                                    >  2:1090206/1‑150 (MQ=255)
                                    ccGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATAtcagtc                                                                                                                   <  1:3386652/150‑1 (MQ=255)
                                     cGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATAtcagtca                                                                                                                  >  2:4888396/1‑150 (MQ=255)
                                     cGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATAtcagtca                                                                                                                  >  2:3634918/1‑150 (MQ=255)
                                        tggctggcGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGtt                                                                                                               <  2:725568/150‑1 (MQ=255)
                                         ggctggcGTCAGAAACCCGGCGATCAAAGTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGAATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGttt                                                                                                              <  2:226587/150‑1 (MQ=255)
                                             ggcGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAAc                                                                                                          <  2:801287/150‑1 (MQ=255)
                                              gcGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACa                                                                                                         >  2:5073112/1‑150 (MQ=255)
                                              gcGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACa                                                                                                         <  2:4507017/150‑1 (MQ=255)
                                               cGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAt                                                                                                        <  1:696780/150‑1 (MQ=255)
                                                gTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCACATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAtt                                                                                                       >  2:3692179/1‑150 (MQ=255)
                                                  cAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTc                                                                                                     <  1:3913661/150‑1 (MQ=255)
                                                   aGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCg                                                                                                    <  2:3479645/150‑1 (MQ=255)
                                                   aGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCg                                                                                                    >  2:2556003/1‑150 (MQ=255)
                                                   aGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCg                                                                                                    <  2:1245003/150‑1 (MQ=255)
                                                     aaaaCCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGcc                                                                                                  >  1:1347307/1‑150 (MQ=255)
                                                     aaaaCCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGcc                                                                                                  >  2:4216477/1‑150 (MQ=255)
                                                       aaCCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGcccc                                                                                                >  1:4358699/1‑150 (MQ=255)
                                                        aCCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCg                                                                                               >  2:1032379/1‑150 (MQ=255)
                                                        aCCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCg                                                                                               <  2:1397869/150‑1 (MQ=255)
                                                                   aaTTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCa                                                                                    <  1:2530213/150‑1 (MQ=255)
                                                                    aTTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCaa                                                                                   <  1:4038183/150‑1 (MQ=255)
                                                                    aTTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCaa                                                                                   >  1:2012863/1‑150 (MQ=255)
                                                                     ttCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGAATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAAt                                                                                  >  1:75539/1‑150 (MQ=255)
                                                                      tCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATc                                                                                 >  2:2082479/1‑150 (MQ=255)
                                                                       cagcagGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCa                                                                                >  1:1228517/1‑150 (MQ=255)
                                                                        agcagGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCAc                                                                               >  2:3829773/1‑150 (MQ=255)
                                                                         gcagGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACg                                                                              <  1:2332696/150‑1 (MQ=255)
                                                                         gcagGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACg                                                                              <  2:3984932/150‑1 (MQ=255)
                                                                          cagGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACgg                                                                             >  2:4339815/1‑150 (MQ=255)
                                                                            gggTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTa                                                                           <  2:2959454/150‑1 (MQ=255)
                                                                                 cTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAAcag                                                                      >  2:3325800/1‑150 (MQ=255)
                                                                                  ttttACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagc                                                                     >  1:1873517/1‑150 (MQ=255)
                                                                                   tttACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagca                                                                    <  2:1603256/150‑1 (MQ=255)
                                                                                       cccGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGaaa                                                                >  2:4952231/1‑150 (MQ=255)
                                                                                       cccGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGAATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGaaa                                                                <  2:1725664/150‑1 (MQ=255)
                                                                                        ccGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGaaaa                                                               <  2:1290462/150‑1 (MQ=255)
                                                                                                 ttGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagcagc                                                      >  1:1165540/1‑150 (MQ=255)
                                                                                                 ttGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagcagc                                                      >  2:4785462/1‑150 (MQ=255)
                                                                                                   ggCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagcagcag                                                    <  1:3634918/150‑1 (MQ=255)
                                                                                                      ccGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGcagcagcag                                                 >  2:2623596/1‑150 (MQ=255)
                                                                                                              cGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAAcg                                         <  1:649672/150‑1 (MQ=255)
                                                                                                                  aCCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAATTCGCCCAGGTATTTTTTCAATCACGTTAAACAGCAGAAAACAGAGAAGCAGAAGAATTAACGCGGt                                     >  1:2641375/1‑150 (MQ=255)
                                                                                                                   ccTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTg                                    <  1:4004827/150‑1 (MQ=255)
                                                                                                                     tGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGAc                                  <  1:3106193/150‑1 (MQ=255)
                                                                                                                      gCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGAcc                                 >  2:3278060/1‑150 (MQ=255)
                                                                                                                      gCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGAcc                                 <  1:4989290/150‑1 (MQ=255)
                                                                                                                      gCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGAcc                                 <  1:943280/150‑1 (MQ=255)
                                                                                                                        tCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGc                               >  1:3506518/1‑150 (MQ=255)
                                                                                                                          gccgcGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCAc                             >  1:4554949/1‑150 (MQ=255)
                                                                                                                           ccgcGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCAcc                            <  2:3235693/150‑1 (MQ=255)
                                                                                                                            cgcgTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCg                           >  1:2763077/1‑150 (MQ=255)
                                                                                                                               gTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTcc                        <  2:4143195/150‑1 (MQ=255)
                                                                                                                                ttCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCt                       <  2:3884932/150‑1 (MQ=255)
                                                                                                                                   cACCGTTAAGCTAAAACGAATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGc                    >  2:2325981/1‑150 (MQ=255)
                                                                                                                                    aCCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTggat                   <  1:1712084/150‑3 (MQ=255)
                                                                                                                                    aCCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCt                   <  1:4215363/150‑1 (MQ=255)
                                                                                                                                      cGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGc                 <  2:3619832/150‑1 (MQ=255)
                                                                                                                                      cGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGc                 <  2:3424451/150‑1 (MQ=255)
                                                                                                                                       gTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCg                >  1:4309051/1‑150 (MQ=255)
                                                                                                                                        ttAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGa               >  1:2139893/1‑150 (MQ=255)
                                                                                                                                           aGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAg            <  1:1477942/150‑1 (MQ=255)
                                                                                                                                            gCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAgg           >  2:1109940/1‑150 (MQ=255)
                                                                                                                                              tAAAACGAATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAg         >  2:4125808/1‑150 (MQ=255)
                                                                                                                                               aaaaCGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGc        <  1:1370097/150‑1 (MQ=255)
                                                                                                                                                aaaCGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCTGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGcc       >  2:2687829/1‑150 (MQ=255)
                                                                                                                                                aaaCGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGcc       >  2:45961/1‑150 (MQ=255)
                                                                                                                                                     cATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGAAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAAttt  >  1:563427/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTT  >  CP009273/69218‑69516

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: