Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 477,267 G→T 100% Q865K (CAG→AAG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273477,2670GT98.1% 515.1 / 0.8 161Q865K (CAG→AAG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (78/80);  total (81/80)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 161 total aligned reads displayed.

CGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTC  >  CP009273/477124‑477411
                                                                                                                                               |                                                                                                                                                
     cGATAACCCCCAGCGGAACGACCAGCATACCGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACAt                                                                                                                                         <  2:1718534/148‑1 (MQ=255)
      gATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATc                                                                                                                                        <  2:525424/148‑1 (MQ=255)
       aTAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCTTAGAAATCGCGTTCCGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATcc                                                                                                                                       <  1:3448053/148‑1 (MQ=255)
        tAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATccc                                                                                                                                      <  2:792472/148‑1 (MQ=255)
        tAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATccc                                                                                                                                      <  2:3385883/148‑1 (MQ=255)
         aaCCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATcccc                                                                                                                                     >  2:1924160/1‑148 (MQ=255)
           cccccAGCGTAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGt                                                                                                                                   <  1:726634/148‑1 (MQ=255)
           cccccAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGt                                                                                                                                   <  1:109858/148‑1 (MQ=255)
             cccAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTcc                                                                                                                                 <  2:3407516/148‑1 (MQ=255)
              ccAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCa                                                                                                                                >  1:542136/1‑148 (MQ=255)
              ccAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCa                                                                                                                                >  1:1646854/1‑148 (MQ=255)
                aGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGt                                                                                                                              <  1:221501/148‑1 (MQ=255)
                  cgcaccGACCAGAATACCGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCa                                                                                                                            <  1:1504080/143‑1 (MQ=255)
                   ggAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCAt                                                                                                                           <  2:2696594/148‑1 (MQ=255)
                    gAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATa                                                                                                                          >  2:3459405/1‑148 (MQ=255)
                    gAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATa                                                                                                                          >  1:2346486/1‑148 (MQ=255)
                    gAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATa                                                                                                                          <  1:3492539/148‑1 (MQ=255)
                    gAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATa                                                                                                                          >  1:1461519/1‑148 (MQ=255)
                    gAACGACCAGCATAACGGAGAACGGAACCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATa                                                                                                                          >  1:1461557/1‑148 (MQ=255)
                          ccAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAac                                                                                                                    >  1:232985/1‑148 (MQ=255)
                            aGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAacac                                                                                                                  <  1:1036293/148‑1 (MQ=255)
                               aTAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACACTCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCgg                                                                                                               <  2:1150539/148‑1 (MQ=255)
                                tAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGt                                                                                                              <  2:2788139/148‑1 (MQ=255)
                                 aaCGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTa                                                                                                             >  1:2997134/1‑148 (MQ=255)
                                 aaCGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTa                                                                                                             >  1:3242401/1‑148 (MQ=255)
                                  aCGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAg                                                                                                            >  2:1136292/1‑148 (MQ=255)
                                  aCGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAg                                                                                                            <  2:1644038/148‑1 (MQ=255)
                                  aCGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAg                                                                                                            >  2:1136330/1‑148 (MQ=255)
                                   cGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAgg                                                                                                           <  2:1567093/148‑1 (MQ=255)
                                        aaCGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGtt                                                                                                      <  1:1090807/148‑1 (MQ=255)
                                        aaCGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGtt                                                                                                      <  1:2060402/148‑1 (MQ=255)
                                        aaCGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGtt                                                                                                      <  2:2645986/148‑1 (MQ=255)
                                           ggAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGGAGGCAGTTTGc                                                                                                   <  1:1789785/148‑1 (MQ=255)
                                            gAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCt                                                                                                  >  1:2873910/1‑148 (MQ=255)
                                            gAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCt                                                                                                  >  1:2874967/1‑148 (MQ=255)
                                            gAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCt                                                                                                  <  1:1277177/148‑1 (MQ=255)
                                             aaTCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTc                                                                                                 >  2:1778496/1‑148 (MQ=255)
                                              aTCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCg                                                                                                <  1:1069678/148‑1 (MQ=255)
                                               tCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGc                                                                                               >  1:2939323/1‑148 (MQ=255)
                                               tCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGc                                                                                               <  2:2973163/148‑1 (MQ=255)
                                                cGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTCGGCAGTTTGCTCGcc                                                                                              >  2:3238994/1‑148 (MQ=255)
                                                  aCCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAg                                                                                            >  1:1466347/1‑148 (MQ=255)
                                                   ccAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgt                                                                                           <  1:2264134/148‑1 (MQ=255)
                                                      gCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttgt                                                                                        <  2:3101640/148‑1 (MQ=255)
                                                       ctctCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttgtt                                                                                       <  1:2218424/148‑1 (MQ=255)
                                                        tctcGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTc                                                                                      <  2:536590/148‑1 (MQ=255)
                                                          tcGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCa                                                                                    >  2:162677/1‑148 (MQ=255)
                                                              aCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAg                                                                                >  1:1380750/1‑148 (MQ=255)
                                                                 gcgcCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTc                                                                             >  2:2022035/1‑148 (MQ=255)
                                                                 gcgcCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTc                                                                             >  2:2760458/1‑148 (MQ=255)
                                                                  cgccgcCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTcc                                                                            <  1:3311545/148‑1 (MQ=255)
                                                                      gccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTg                                                                        <  1:1673819/148‑1 (MQ=255)
                                                                       ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGc                                                                       <  2:1574361/148‑1 (MQ=255)
                                                                       ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGc                                                                       >  1:850651/1‑148 (MQ=255)
                                                                       ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGc                                                                       <  1:672688/148‑1 (MQ=255)
                                                                       ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGc                                                                       >  1:1220119/1‑148 (MQ=255)
                                                                       ccAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGc                                                                       >  1:1220097/1‑148 (MQ=255)
                                                                         aGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCtt                                                                     <  2:3472482/148‑1 (MQ=255)
                                                                             acaGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCAcc                                                                 <  1:2715611/148‑1 (MQ=255)
                                                                              caGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCg                                                                >  1:694388/1‑148 (MQ=255)
                                                                              caGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCg                                                                >  1:414783/1‑148 (MQ=255)
                                                                                ggAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGt                                                              <  1:2076361/148‑1 (MQ=255)
                                                                                ggAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGt                                                              <  2:694366/148‑1 (MQ=255)
                                                                                   acacGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCTGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTACATTGCTTCACCGGTACt                                                           <  2:1288169/148‑1 (MQ=255)
                                                                                    cacGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACtt                                                          >  1:1595024/1‑148 (MQ=255)
                                                                                     acGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCAGTGCTTCACCGGTACtgt                                                         >  2:701351/1‑146 (MQ=255)
                                                                                         cAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTAcc                                                     >  1:93132/1‑148 (MQ=255)
                                                                                         cAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTAcc                                                     >  1:1870561/1‑148 (MQ=255)
                                                                                          aaTCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTAccc                                                    <  2:613236/148‑1 (MQ=255)
                                                                                             cAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGt                                                 >  2:2677129/1‑148 (MQ=255)
                                                                                               aCGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTgc                                               >  2:1282290/1‑148 (MQ=255)
                                                                                               aCGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTgc                                               >  1:396376/1‑148 (MQ=255)
                                                                                               aCGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTgc                                               >  2:3499608/1‑148 (MQ=255)
                                                                                               aCGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTgc                                               >  2:1370895/1‑148 (MQ=255)
                                                                                                 gAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTgccg                                             <  2:805785/148‑1 (MQ=255)
                                                                                                  aaaTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCTCCGGTACTTTTACCCGGTgccgc                                            >  2:2044728/1‑148 (MQ=255)
                                                                                                   aaTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTgccgcc                                           <  1:549758/148‑1 (MQ=255)
                                                                                                   aaTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTgccgcc                                           >  1:3226453/1‑148 (MQ=255)
                                                                                                    aTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCt                                          >  1:736701/1‑148 (MQ=255)
                                                                                                    aTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCt                                          >  2:1033910/1‑148 (MQ=255)
                                                                                                    aTCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCt                                          >  2:1867333/1‑148 (MQ=255)
                                                                                                       gcgTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGc                                       <  1:1924160/148‑1 (MQ=255)
                                                                                                          tACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTa                                    >  1:2252201/1‑148 (MQ=255)
                                                                                                             aGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGa                                 >  2:1739332/1‑148 (MQ=255)
                                                                                                             aGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGa                                 >  2:2277938/1‑148 (MQ=255)
                                                                                                             aGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGa                                 >  2:1788081/1‑148 (MQ=255)
                                                                                                             aGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCAATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGa                                 >  2:2279458/1‑148 (MQ=255)
                                                                                                                 aaGGTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTc                             >  2:1719089/1‑148 (MQ=255)
                                                                                                                   ggTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCa                           <  1:1730672/148‑1 (MQ=255)
                                                                                                                    gTGCCTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCat                          >  1:704410/1‑148 (MQ=255)
                                                                                                                        cTGGTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCatggat                      <  2:2378322/148‑1 (MQ=255)
                                                                                                                          ggTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCatggatgg                    <  1:3497665/148‑1 (MQ=255)
                                                                                                                           gTTGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGc                   >  1:1137429/1‑148 (MQ=255)
                                                                                                                            ttGCCGGAGAGACGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGct                  <  2:902702/148‑2 (MQ=255)
                                                                                                                                  gagagaCGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGt            >  1:1848496/1‑148 (MQ=255)
                                                                                                                                       aCGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTaa       >  2:1870816/1‑148 (MQ=255)
                                                                                                                                        cGTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAAc      >  2:1297216/1‑148 (MQ=255)
                                                                                                                                         gTTCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACg     <  2:1646854/148‑1 (MQ=255)
                                                                                                                                           tCCTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGtt   <  2:1461557/148‑1 (MQ=255)
                                                                                                                                            ccTTATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTc  <  1:1297018/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTC  >  CP009273/477124‑477411

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: