Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 479,244 C→A 100% A206S (GCC→TCC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273479,2440CA95.2% 856.6 / 28.5 272A206S (GCC→TCC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (5/8);  new base A (123/136);  total (128/144)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 272 total aligned reads displayed.

GATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACG  >  CP009273/479097‑479389
                                                                                                                                                   |                                                                                                                                                 
 aTTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGt                                                                                                                                                >  1:435791/1‑150 (MQ=255)
   tCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCCGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTaa                                                                                                                                              >  1:1976949/1‑150 (MQ=255)
   tCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTaa                                                                                                                                              >  1:1980146/1‑150 (MQ=255)
    cACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAAt                                                                                                                                             >  1:564922/1‑150 (MQ=255)
       tttCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGAc                                                                                                                                          <  2:277674/150‑1 (MQ=255)
        ttCAGCAGGATTTTGCCTAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACa                                                                                                                                         >  2:2532261/1‑150 (MQ=255)
        ttCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACa                                                                                                                                         >  1:3543701/1‑150 (MQ=255)
        ttCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACa                                                                                                                                         <  2:2179449/150‑1 (MQ=255)
          cagcagGATTTTGCCTAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATc                                                                                                                                       >  1:692186/1‑150 (MQ=255)
          cagcagGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATc                                                                                                                                       >  1:1797930/1‑150 (MQ=255)
                  ttttGCCTAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCg                                                                                                                               <  1:3000888/150‑1 (MQ=255)
                   tttGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGt                                                                                                                              <  1:565717/150‑1 (MQ=255)
                    ttGCCTAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGtt                                                                                                                             >  2:517628/1‑150 (MQ=255)
                     tGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTa                                                                                                                            <  1:1711739/150‑1 (MQ=255)
                     tGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTa                                                                                                                            >  2:2946907/1‑150 (MQ=255)
                      gCCTAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAg                                                                                                                           <  2:125074/150‑1 (MQ=255)
                      gCCGAACTCTTCAGTAGAGGTCAGCCGCGTCAGAGCCATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAg                                                                                                                           <  2:3632627/150‑1 (MQ=255)
                       cctaaCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGc                                                                                                                          <  1:1244303/147‑1 (MQ=255)
                       ccGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGc                                                                                                                          >  2:3798808/1‑150 (MQ=255)
                            ctctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGaa                                                                                                                     <  1:2974019/150‑1 (MQ=255)
                            ctctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGaa                                                                                                                     <  1:2975143/150‑1 (MQ=255)
                             tctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAAt                                                                                                                    >  1:2798445/1‑150 (MQ=255)
                              ctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAAtt                                                                                                                   <  1:2630773/150‑1 (MQ=255)
                              ctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAAtt                                                                                                                   <  2:646752/150‑1 (MQ=255)
                                tCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTg                                                                                                                 <  1:899093/150‑1 (MQ=255)
                                tCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTg                                                                                                                 >  1:805345/1‑150 (MQ=255)
                                tCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTg                                                                                                                 <  1:2347276/150‑1 (MQ=255)
                                 cAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt                                                                                                                <  2:2223529/150‑1 (MQ=255)
                                 cAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt                                                                                                                >  1:3382462/1‑150 (MQ=255)
                                   gTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTc                                                                                                              <  2:511707/150‑1 (MQ=255)
                                     agagGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAg                                                                                                            <  1:226753/150‑1 (MQ=255)
                                       agGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCt                                                                                                          <  1:2621891/150‑1 (MQ=255)
                                        ggTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCTACCGGCGTTAGCTGGAATTTGTTCAGCTc                                                                                                         <  1:1215578/150‑1 (MQ=255)
                                           cAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCAtt                                                                                                      >  1:3985832/1‑150 (MQ=255)
                                            aGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTc                                                                                                     >  1:3526174/1‑150 (MQ=255)
                                             gACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCg                                                                                                    <  2:376284/150‑1 (MQ=255)
                                              aCGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCgg                                                                                                   >  2:2598291/1‑150 (MQ=255)
                                              aCGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCgg                                                                                                   >  2:679770/1‑150 (MQ=255)
                                              aCGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCgg                                                                                                   >  2:581923/1‑150 (MQ=255)
                                               cgcgTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCggg                                                                                                  <  2:426680/150‑1 (MQ=255)
                                                 cgTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGtt                                                                                                <  1:3965326/150‑1 (MQ=255)
                                                 cgTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGtt                                                                                                <  1:2720656/150‑1 (MQ=255)
                                                  gTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTc                                                                                               <  1:4060808/150‑1 (MQ=255)
                                                  gTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTc                                                                                               <  2:4349321/150‑1 (MQ=255)
                                                     tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATc                                                                                            >  2:255629/1‑150 (MQ=255)
                                                     tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATc                                                                                            >  2:861590/1‑150 (MQ=255)
                                                      gAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATcc                                                                                           <  2:1734866/150‑1 (MQ=255)
                                                       aGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCa                                                                                          >  1:121474/1‑150 (MQ=255)
                                                            taataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATAc                                                                                     <  1:1376369/150‑1 (MQ=255)
                                                               taGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCa                                                                                  >  1:941563/1‑150 (MQ=255)
                                                                   ggCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATcgc                                                                              >  2:501098/1‑150 (MQ=255)
                                                                     cGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGt                                                                            >  2:3434999/1‑150 (MQ=255)
                                                                     cGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGt                                                                            >  2:1944409/1‑150 (MQ=255)
                                                                      gTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTa                                                                           >  2:2303556/1‑150 (MQ=255)
                                                                        tAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACt                                                                         <  2:4480850/150‑1 (MQ=255)
                                                                        tAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACt                                                                         <  2:3400391/150‑1 (MQ=255)
                                                                           gCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACtgtg                                                                      >  1:3349662/1‑150 (MQ=255)
                                                                             tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGaa                                                                    <  2:239569/150‑1 (MQ=255)
                                                                              gTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAAc                                                                   <  2:692186/150‑1 (MQ=255)
                                                                                tGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCg                                                                 <  2:1708020/150‑1 (MQ=255)
                                                                                tGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCg                                                                 >  2:754621/1‑150 (MQ=255)
                                                                                 ggCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGa                                                                >  2:3992791/1‑150 (MQ=255)
                                                                                      ttCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaacaac                                                           <  1:4288853/150‑1 (MQ=255)
                                                                                            ggcggcGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACa                                                     >  2:2806556/1‑150 (MQ=255)
                                                                                             gcggcgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACAt                                                    <  2:4484574/150‑1 (MQ=255)
                                                                                             gcggcgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACAt                                                    <  2:1977049/150‑1 (MQ=255)
                                                                                              cggcgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATc                                                   >  2:1780577/1‑150 (MQ=255)
                                                                                               ggcgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCa                                                  >  2:3261442/1‑150 (MQ=255)
                                                                                                 cgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCAcc                                                >  2:3170716/1‑150 (MQ=255)
                                                                                                 cgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCAcc                                                <  1:2532261/150‑1 (MQ=255)
                                                                                                  gTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCAccc                                               >  1:3969633/1‑150 (MQ=255)
                                                                                                    accaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCAc                                             <  1:561828/150‑1 (MQ=255)
                                                                                                         cGAGCAGACCCGCCGCAACCGGCGCGTTCAGCGCTTTGATGGAGGTAATGACAACACCCGGCGTTAGCTGGAATATGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCg                                        <  2:4236065/150‑1 (MQ=255)
                                                                                                          gAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGa                                       <  2:1864825/150‑1 (MQ=255)
                                                                                                           aGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGAc                                      >  2:1585482/1‑150 (MQ=255)
                                                                                                            gCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACg                                     >  1:4393088/1‑150 (MQ=255)
                                                                                                             cTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGt                                    <  1:3434999/150‑1 (MQ=255)
                                                                                                               gACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTAc                                  <  1:643692/150‑1 (MQ=255)
                                                                                                               gACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTAc                                  <  1:1624104/150‑1 (MQ=255)
                                                                                                                 cccGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACgg                                <  2:504058/150‑1 (MQ=255)
                                                                                                                   cgccgcAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCt                              >  1:3616135/1‑150 (MQ=255)
                                                                                                                   cgccgcAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCt                              <  1:1887430/150‑1 (MQ=255)
                                                                                                                    gccgcAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTg                             >  1:2406520/1‑150 (MQ=255)
                                                                                                                      cgcAACCTGGGCGTTCTGCGCTTTGATTGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGAt                           <  2:4027015/150‑1 (MQ=255)
                                                                                                                       gcAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATg                          >  1:2881400/1‑150 (MQ=255)
                                                                                                                       gcAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATg                          <  1:3798814/150‑1 (MQ=255)
                                                                                                                       gcAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATg                          <  2:2120091/150‑1 (MQ=255)
                                                                                                                        cAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTAGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATgg                         >  1:1985049/1‑150 (MQ=255)
                                                                                                                         aaCCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGc                        <  2:1965514/150‑1 (MQ=255)
                                                                                                                         aaCCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGc                        >  2:3228026/1‑150 (MQ=255)
                                                                                                                          aCCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGTAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCa                       >  1:4036671/1‑150 (MQ=255)
                                                                                                                             tGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCt                    >  1:3910528/1‑150 (MQ=255)
                                                                                                                              gggCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCtt                   >  1:896058/1‑150 (MQ=255)
                                                                                                                                  gTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCat               >  1:147657/1‑150 (MQ=255)
                                                                                                                                     cTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATAtt            <  2:1575701/150‑1 (MQ=255)
                                                                                                                                       gcgcTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTcg          >  2:2322246/1‑150 (MQ=255)
                                                                                                                                         gcTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTcgcc        <  2:147657/150‑1 (MQ=255)
                                                                                                                                            ttGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCgccgcc     <  2:3813474/150‑1 (MQ=255)
                                                                                                                                             tGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCa    <  2:2933192/150‑1 (MQ=255)
                                                                                                                                              gATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCAc   <  2:3223165/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACG  >  CP009273/479097‑479389

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: