Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,978,089 A→C 24.7% M288R (ATG→AGG)  yqeF ← short chain acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,978,0890AC24.7% 150.1 / 18.6 89M288R (ATG→AGG) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (19/47);  new base C (17/6);  total (36/53)
Fisher's exact test for biased strand distribution p-value = 2.07e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACC  >  CP009273/2977941‑2978237
                                                                                                                                                    |                                                                                                                                                     
gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCa                                                                                                                                                       <  2:1444310/149‑1 (MQ=255)
 acCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCAt                                                                                                                                                      <  2:1172746/149‑1 (MQ=255)
    cGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAAt                                                                                                                                                  <  2:760302/150‑1 (MQ=255)
     gACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATg                                                                                                                                                 <  1:616968/150‑1 (MQ=255)
       cGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGcc                                                                                                                                               <  1:647747/150‑1 (MQ=255)
       cGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCCGACACCCCCGCCCAATTTCCATCaaatcc                                                                                                                                               >  1:186956/1‑146 (MQ=255)
             tCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCt                                                                                                                                         <  2:366604/150‑1 (MQ=255)
             tCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCt                                                                                                                                         <  1:586715/150‑1 (MQ=255)
                  cTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCCGGCAACGGCGGGGCGCAAAAACCCGCGCAAATCCCCTCCAAGCCCGTTCTacac                                                                                                                                     >  2:1339845/1‑149 (MQ=255)
                    cAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCAATTCCCCACAAAGCCCGGTCTACACCaa                                                                                                                                  >  1:848260/1‑150 (MQ=255)
                     aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCa                                                                                                                                   <  1:901518/148‑1 (MQ=255)
                     aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAAc                                                                                                                                 <  1:825025/150‑1 (MQ=255)
                       gCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGAACTACACCAACGc                                                                                                                               >  1:394149/1‑150 (MQ=255)
                         aTCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTg                                                                                                                             <  1:405809/150‑1 (MQ=255)
                               ccAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCCCAAACACCGGCCCACATTTCATCAATACCGGAACTACACAAAAACAGGCaaac                                                                                                                       >  1:614357/1‑149 (MQ=255)
                                 aaCCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCACTTCCCATCACTTCCGGCTCTACACCCACACCGGAAAATgc                                                                                                                     >  2:930747/1‑150 (MQ=255)
                                      aaaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCAATTCCCATCAATGCCCGATCTACACCAACCCTGGCAAATGcgcgg                                                                                                                 >  2:1210507/1‑148 (MQ=255)
                                       aaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATATTCCGGCACAATACCGAGCACTGCGAGATCTACACCAACGCTGGCAAATGCGCGAAt                                                                                                               <  2:1370667/150‑1 (MQ=255)
                                       aaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCTGAGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAAt                                                                                                               <  2:1073668/150‑1 (MQ=255)
                                        aGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATg                                                                                                              <  1:33827/150‑1 (MQ=255)
                                        aGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCAAAACCCGCGCCAATCCCCCCCATGCCGGGTCTTCACCCACCCCGGCCAAAGCCCGGATg                                                                                                              >  2:1276659/1‑150 (MQ=255)
                                           gcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCACCTACCCGCGCTATTTCCCTTAATGCCGGGTCTTCACCAACACCCGCAAAAGCACGaaatgtc                                                                                                           >  1:76339/1‑145 (MQ=255)
                                           gcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCCATTTCCATCCATGCCGGGTCTACACCAACGCCGGCCAATGCGCGaa                                                                                                                >  2:1437515/1‑145 (MQ=255)
                                            ccTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAAt                                                                                                               <  1:1437515/145‑1 (MQ=255)
                                                 gcagcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCAATGCCCGATCTACACCCACGCTGGCAAAAACGCGCATGCGGGCCCGc                                                                                                     >  2:225658/1‑150 (MQ=255)
                                                  cagcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCACCGTAGCAATACCGCGCCATCCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCa                                                                                                    <  2:242904/150‑1 (MQ=255)
                                                  cagcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGc                                                                                                     <  1:767221/149‑1 (MQ=255)
                                                   agcaAGCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCAc                                                                                                  <  1:610950/150‑1 (MQ=255)
                                                      aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCAATTCCCCTCAATTGCGGGTCTACACCAAACCTGGCCAAAGCGCGAAAGCGGGCCAGCACggg                                                                                               >  1:1164736/1‑150 (MQ=255)
                                                        aCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTa                                                                                             <  1:111056/150‑1 (MQ=255)
                                                        aCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTa                                                                                             <  1:977455/150‑1 (MQ=255)
                                                          gCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTaaa                                                                                           <  1:406008/150‑1 (MQ=255)
                                                           cTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAAt                                                                                          <  2:395393/150‑1 (MQ=255)
                                                              cATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCGCAATTTCCATCAATGCCGGAGCTACCCCAACACTGGCCAAAGCGCGAAAACGGGCCACCACCGGGAAAAtt                                                                                         >  1:470325/1‑149 (MQ=255)
                                                               aTTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGAGCGATTCGCCGCAATGCGGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCa                                                                                       <  2:473823/149‑1 (MQ=255)
                                                                ttAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTTCCATCAATGCCGGAACTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAAATCAAc                                                                                     >  1:1418545/1‑150 (MQ=255)
                                                                 tAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTTGCCTAAACCCGCGCCAATCCCCTCCATGCCCGGTCTCCACCACCCCTGGGAAATGCGCGGATGTGGGGGCCCACGGGGAAATTTaaag                                                                                    >  2:1016561/1‑148 (MQ=255)
                                                                 tAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCACATTCCACCAAAGCCGGGTCTACCCCAACCCTGGGAAAAGCGCGAAAACCGGGCAGCACGGGGAAATTCAACg                                                                                    >  1:857336/1‑150 (MQ=255)
                                                                   gCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCt                                                                                  <  2:304000/150‑1 (MQ=255)
                                                                     cTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCGCACTTCCCCTCAATGCCGGATCTACACCAACGCCGGCAAAAGCGCGAATGCGGGCCAGCACCGGGAAAATCAACGCTcc                                                                                >  1:593953/1‑149 (MQ=255)
                                                                        gATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCGCAGCAATGCGGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTg                                                                              <  2:327116/149‑1 (MQ=255)
                                                                        gATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGt                                                                               <  1:496962/148‑1 (MQ=255)
                                                                        gATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGt                                                                               <  2:1153323/148‑1 (MQ=255)
                                                                           aaGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCGCAATTCCCCTCAAAGCCCGATCTACACCAACGCTGGCCAATACGCGAAAGCGGGCCAGCACGGGGAAAATCAAAGCCCGTGCttt                                                                          >  1:710003/1‑150 (MQ=255)
                                                                            aGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCttt                                                                          <  1:1103643/149‑1 (MQ=255)
                                                                             gATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTgg                                                                        <  1:225658/150‑1 (MQ=255)
                                                                                 gACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCTTCTGAATGATACGGCGCCAACCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTc                                                                    <  2:717060/150‑1 (MQ=255)
                                                                                 gACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCCATTTCCCCCCATGCCGGATCTTCACCCAAGCTGGCAAATGCGCGGAAGGGGGCCCGCACCGGGAAAATTAAAGCTCGTGCTTTTGCttt                                                                    >  1:758579/1‑149 (MQ=255)
                                                                                  aCTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCCATTCCCCTCCATTCCGGATCCACACCAACGCCGGGAAAATCGCGCAATCGGGCCAGCCCGGGGAAATTCAACGCTCGTGCTTTTGCtttg                                                                   >  2:874777/1‑148 (MQ=255)
                                                                                       aGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTTCCCTCAATGCCGGGTCTACACCAACCCTGGCAAATGCGCGAATGCGGGGCCGCACGGGGAAAATCAACGCCCGTGCTTTTGCTTCGCtcat                                                              >  1:389360/1‑150 (MQ=255)
                                                                                           aaCTGCCAGCCTACACGCTCCCGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcatc                                                          <  1:1339845/150‑1 (MQ=255)
                                                                                             cTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACCCCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTcatcatcat                                                        <  1:1315994/150‑1 (MQ=255)
                                                                                              tGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCATCCATTCCGGGTCTACACCAACACTGGGAAAAGCGCGAAAACGGGCCCGCACGGGTGAAATTAACGCTCGCGCTTTGGCCTCTCTCATCATCAtt                                                       >  2:609140/1‑150 (MQ=255)
                                                                                                ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTTTCTTTAATGATTCGGCGCCCACCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTa                                                      <  2:1257325/149‑1 (MQ=255)
                                                                                                ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTCCGGAGCAATTCCCCGCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCAtt                                                       <  2:1314425/148‑1 (MQ=255)
                                                                                                ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCAtt                                                       <  1:1210507/148‑1 (MQ=255)
                                                                                                ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACAACCGCGCAATTTCCATCAATGCCGGATCTACACCAACGCCGGCAAAAGCGCGAAAGCGGGCCAGCACGGGGAAAATTAAAGCTCGTGCTTTTGCTTCGCTCATCATCATTAc                                                     >  2:783288/1‑150 (MQ=255)
                                                                                                ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACAACCGCGCAATTCCCCTCCATGCCGGAGCTACACCAACGCCGGCAAAAGCGCGAATGCGGGGCAGCACGGGGAAATTCAACACCCCTGCTTTTGTTTCGCTCCTCATCAtttc                                                     >  2:157841/1‑148 (MQ=255)
                                                                                                 cAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTAcc                                                    <  1:1016561/150‑1 (MQ=255)
                                                                                                 cAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTAcc                                                    <  1:1245298/150‑1 (MQ=255)
                                                                                                    ccTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGAGCAATACGCAGCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCa                                                 <  2:186956/150‑1 (MQ=255)
                                                                                                      tACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCCTCAATGCCTGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGGCAGCACGGGGAAATTCAACCCGCGTGCTTTGGCTTCGCTCATCATCATTACCGCaac                                               >  1:906951/1‑148 (MQ=255)
                                                                                                      tACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAAATCCCATCAATGCCGGATCTACACCAACACCGGGAAAAACGCGAATGCGGGCCCGCACGGGGAAAATCAAAGCTCGTGCTTTTGCTTCGCCCCCCATCATTACCGCaac                                               >  2:1345892/1‑148 (MQ=255)
                                                                                                       acacGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCCATGCCCGGTCTACACCAACCCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAAATTAACGCCCGCGCTTTGGCTTTGCTCATCATCCTTTACGCAGcc                                              >  2:543749/1‑149 (MQ=255)
                                                                                                        cacGCTCCAGGCAACGGCGGGTCGCATACTCCGGCACGATTCCCAGCCAAGCGGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTg                                             <  2:1437148/150‑1 (MQ=255)
                                                                                                        cacGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTg                                             <  1:583800/150‑1 (MQ=255)
                                                                                                         acGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgc                                            <  2:985999/150‑1 (MQ=255)
                                                                                                         acGCTCCAGGCAACGGCGGGTCGCATACAACCGCGCAAATCCCCTCAATCCCGGAGCTACACCAACGCCGGCAAAAGCGCGGATGCGGGCCCGCACGGGGAAATTTAACACCCGTGCTTTTGCTTCGCTCATCATCATTACCGCAGCTgc                                            >  2:636905/1‑150 (MQ=255)
                                                                                                            cTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCGCATCAATGCCAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGcc                                         <  2:1215359/150‑1 (MQ=255)
                                                                                                               cAGGCAACGGCGGGTCGCATACACCCGCGCAATTTCCCTCAATGCCGGATCTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCCGCACGGGGAAATTCAACGCCCGTGCTTTTGCTTCGCTCATCATCATTACCGCCGCTGCGCCCTc                                      >  2:114505/1‑150 (MQ=255)
                                                                                                                aGGCAACGGCGGGGCGCATACACCCGCGCAAATCCCATCAATGCCGGATCTACACCCACGCCGGCAAAAGCGCGAATGCGGGCCAGCCCGGGGAAATTCAAACCCCGTGCTTTTGCTTCGCTCCTCATCATTACCCCCGCTGCGcccccg                                     >  2:1218615/1‑146 (MQ=255)
                                                                                                                     aCGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTAt                                <  1:157841/150‑1 (MQ=255)
                                                                                                                       ggcggGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATgg                              <  2:20297/150‑1 (MQ=255)
                                                                                                                        gcggGTCGCATACTCCGGCGCAATTCCCAGCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGga                             <  2:65177/150‑1 (MQ=255)
                                                                                                                          gggTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatg                           <  1:754488/150‑1 (MQ=255)
                                                                                                                            gTCGCATTCACCGGCGCGATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatgat                         <  2:614357/150‑1 (MQ=255)
                                                                                                                               gCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCa                      >  1:3151/1‑150 (MQ=255)
                                                                                                                                cATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAt                     <  2:710003/150‑1 (MQ=255)
                                                                                                                                 aTACACCGGCGCAATTCCCCTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt                    <  2:796469/150‑1 (MQ=255)
                                                                                                                                 aTACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAt                     >  1:151432/1‑149 (MQ=255)
                                                                                                                                  tACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTa                   >  1:574215/1‑150 (MQ=255)
                                                                                                                                      ccGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCAt                                       >  2:1017494/1‑126 (MQ=255)
                                                                                                                                      ccGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCg               <  1:543749/150‑1 (MQ=255)
                                                                                                                                       cGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATc                                      <  1:1017494/126‑1 (MQ=255)
                                                                                                                                       cGGAGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCg               <  2:288746/149‑1 (MQ=255)
                                                                                                                                        ggCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt             <  2:848260/150‑1 (MQ=255)
                                                                                                                                          cgcAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGt           <  1:52292/150‑1 (MQ=255)
                                                                                                                                             aaTTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTca         >  2:1102739/1‑149 (MQ=255)
                                                                                                                                                 cccATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGaa    >  2:687311/1‑150 (MQ=255)
                                                                                                                                                 cccATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGa     <  1:1206072/149‑1 (MQ=255)
                                                                                                                                                    aTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAcc  >  2:549482/1‑149 (MQ=255)
                                                                                                                                                    aTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAcc  >  1:968123/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACC  >  CP009273/2977941‑2978237

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: