Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,173,145 A→T 100% Q658L (CAG→CTG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,173,1450AT98.4% 1301.7 / 5.7 381Q658L (CAG→CTG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base T (205/170);  total (209/172)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 393 total aligned reads displayed.

CGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGT  >  CP009273/4172999‑4173290
                                                                                                                                                  |                                                                                                                                                 
 gaagaaGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCAACCctgca                                                                                                                                                <  1:961520/149‑1 (MQ=255)
  aagaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgcag                                                                                                                                               <  2:4175525/149‑1 (MQ=255)
  aagaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgcag                                                                                                                                               >  2:3118816/1‑149 (MQ=255)
   agaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAgg                                                                                                                                              >  1:1664773/1‑149 (MQ=255)
    gaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggt                                                                                                                                             >  1:1122606/1‑149 (MQ=255)
    gaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggt                                                                                                                                             >  1:1757529/1‑149 (MQ=255)
    gaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggt                                                                                                                                             >  1:2908088/1‑149 (MQ=255)
      agGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtgg                                                                                                                                           >  1:621482/1‑149 (MQ=255)
      agGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtgg                                                                                                                                           >  2:2528954/1‑149 (MQ=255)
       ggCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtggt                                                                                                                                          >  2:4683655/1‑149 (MQ=255)
       ggCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtggt                                                                                                                                          <  1:2499995/149‑1 (MQ=255)
       ggCAACTACGTTATCGCCCAGGCGAACTCCAACGTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtggt                                                                                                                                          >  2:1683478/1‑149 (MQ=255)
          aaCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATc                                                                                                                                       >  2:712131/1‑149 (MQ=255)
          aaCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATc                                                                                                                                       <  2:3939397/149‑1 (MQ=255)
           aCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                      <  2:3298726/149‑1 (MQ=255)
           aCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                      <  1:4478977/149‑1 (MQ=255)
           aCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                      >  2:2080994/1‑149 (MQ=255)
           aCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                      <  1:2080994/149‑1 (MQ=255)
           aCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                      >  2:4478977/1‑149 (MQ=255)
             tACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGt                                                                                                                                    >  1:2980568/1‑149 (MQ=255)
              aCGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTc                                                                                                                                   >  1:1239867/1‑149 (MQ=255)
              aCGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTc                                                                                                                                   <  1:4268858/149‑1 (MQ=255)
                gTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCgg                                                                                                                                 <  2:4750965/149‑1 (MQ=255)
                   aTCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGc                                                                                                                              <  2:857161/149‑1 (MQ=255)
                    tCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCg                                                                                                                             <  2:3443201/149‑1 (MQ=255)
                     cGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGt                                                                                                                            <  2:19230/149‑1 (MQ=255)
                      gCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTc                                                                                                                           >  1:4296272/1‑149 (MQ=255)
                       cccAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTcc                                                                                                                          >  1:2689484/1‑149 (MQ=255)
                        ccAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTccc                                                                                                                         <  1:4344068/149‑1 (MQ=255)
                          aGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTg                                                                                                                       >  1:1577255/1‑149 (MQ=255)
                          aGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTg                                                                                                                       >  1:1353136/1‑149 (MQ=255)
                          aGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTg                                                                                                                       >  2:4314360/1‑149 (MQ=255)
                           ggCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGa                                                                                                                      >  2:3108452/1‑149 (MQ=255)
                           ggCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGa                                                                                                                      >  2:1712475/1‑149 (MQ=255)
                             cGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATc                                                                                                                    <  2:1277381/149‑1 (MQ=255)
                              gAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATcc                                                                                                                   <  2:3040749/149‑1 (MQ=255)
                               aaCTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATccc                                                                                                                  >  2:2885691/1‑149 (MQ=255)
                               aaCTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATccc                                                                                                                  >  1:580713/1‑149 (MQ=255)
                                aCTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATccct                                                                                                                 >  1:3238806/1‑148 (MQ=255)
                                    cAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTcc                                                                                                             <  2:4258223/149‑1 (MQ=255)
                                     aaCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCt                                                                                                            >  2:2304346/1‑149 (MQ=255)
                                      aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTg                                                                                                           <  1:4246590/149‑1 (MQ=255)
                                          ggATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAac                                                                                                       <  1:2360453/149‑1 (MQ=255)
                                          ggATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAac                                                                                                       >  2:1396062/1‑149 (MQ=255)
                                           gATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAaca                                                                                                      >  2:4066795/1‑149 (MQ=255)
                                            aTGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAacac                                                                                                     >  1:753678/1‑149 (MQ=255)
                                              gaagaaGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGa                                                                                                   <  2:3891304/149‑1 (MQ=255)
                                                  aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGAc                                                                                               <  2:3227641/149‑1 (MQ=255)
                                                  aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGAc                                                                                               >  1:3476757/1‑149 (MQ=255)
                                                  aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGAc                                                                                               >  2:2146742/1‑149 (MQ=255)
                                                   agGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACg                                                                                              >  2:3848181/1‑149 (MQ=255)
                                                   agGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACg                                                                                              >  1:1447000/1‑149 (MQ=255)
                                                      ccACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCa                                                                                           >  1:681084/1‑149 (MQ=255)
                                                         cTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAAcc                                                                                        <  2:1366152/149‑1 (MQ=255)
                                                         cTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATTACGCCAAcc                                                                                        >  2:1106525/1‑149 (MQ=255)
                                                             gTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGc                                                                                    <  2:4612480/149‑1 (MQ=255)
                                                                  agaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGa                                                                               >  2:2978743/1‑149 (MQ=255)
                                                                   gaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGAACAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGAt                                                                              <  1:1716268/149‑1 (MQ=255)
                                                                      cTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATggg                                                                           <  2:585436/149‑1 (MQ=255)
                                                                      cTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATggg                                                                           >  1:4043792/1‑149 (MQ=255)
                                                                       tGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGt                                                                          <  2:1522405/149‑1 (MQ=255)
                                                                       tGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGt                                                                          >  1:162922/1‑149 (MQ=255)
                                                                        ggTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTg                                                                         <  2:4393791/149‑1 (MQ=255)
                                                                         gTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGc                                                                        <  2:4022634/149‑1 (MQ=255)
                                                                          tAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCg                                                                       >  2:4725647/1‑149 (MQ=255)
                                                                            aCTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGaa                                                                     <  1:4675479/149‑1 (MQ=255)
                                                                             cTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGacc                                                                    >  2:3776317/1‑147 (MQ=255)
                                                                               tGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACAt                                                                  <  1:3938362/149‑1 (MQ=255)
                                                                                         aaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATGGATGGGTGCGAACATGCAACGTCAg                                                        >  1:226160/1‑149 (MQ=255)
                                                                                          aGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                       >  2:855921/1‑149 (MQ=255)
                                                                                          aGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                       <  1:4786031/149‑1 (MQ=255)
                                                                                            gCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGcc                                                     >  1:684572/1‑149 (MQ=255)
                                                                                               aaTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGtt                                                  >  2:1163547/1‑149 (MQ=255)
                                                                                               aaTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGtt                                                  >  1:1040952/1‑149 (MQ=255)
                                                                                                 tCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTcc                                                >  2:2461061/1‑149 (MQ=255)
                                                                                                    aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAc                                             >  2:2792168/1‑149 (MQ=255)
                                                                                                    aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAc                                             <  2:1324151/149‑1 (MQ=255)
                                                                                                    aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAc                                             <  1:1949498/149‑1 (MQ=255)
                                                                                                      cTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGActc                                           <  2:3877440/149‑1 (MQ=255)
                                                                                                         gTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgc                                        <  2:4141931/149‑1 (MQ=255)
                                                                                                             aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                    >  2:1892721/1‑149 (MQ=255)
                                                                                                                cgcgACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGa                                 <  1:3329680/149‑1 (MQ=255)
                                                                                                                 gcgACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGAt                                <  1:4174968/149‑1 (MQ=255)
                                                                                                                    aCCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                             <  1:4601725/149‑1 (MQ=255)
                                                                                                                    aCCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                             <  1:4076963/149‑1 (MQ=255)
                                                                                                                     ccAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGc                            <  2:2387441/149‑1 (MQ=255)
                                                                                                                      cAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGcc                           <  1:3065392/149‑1 (MQ=255)
                                                                                                                       aGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                          >  2:2936943/1‑149 (MQ=255)
                                                                                                                        ggTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                         >  1:2187374/1‑149 (MQ=255)
                                                                                                                        ggTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                         >  1:2769164/1‑149 (MQ=255)
                                                                                                                        ggTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTGTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                         <  1:3880402/149‑1 (MQ=255)
                                                                                                                         gTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCt                        >  2:354626/1‑149 (MQ=255)
                                                                                                                            gACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                     >  2:4352299/1‑149 (MQ=255)
                                                                                                                            gACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                     >  2:1546521/1‑149 (MQ=255)
                                                                                                                                 cATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                >  1:703863/1‑149 (MQ=255)
                                                                                                                                   tGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt              <  1:833896/149‑1 (MQ=255)
                                                                                                                                     gACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg            <  2:4101002/149‑1 (MQ=255)
                                                                                                                                        gTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTAt         <  1:3658354/149‑1 (MQ=255)
                                                                                                                                          agccACCCTGCCGGTGGGATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCACCGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATgg       <  2:2647837/147‑1 (MQ=255)
                                                                                                                                              aCCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg   >  1:3552471/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGT  >  CP009273/4172999‑4173290

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: ‑