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% 871.4 / 1.6 265Q658L (CAG→CTG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (144/116);  total (147/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 276 total aligned reads displayed.

GAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACG  >  CP009273/4173000‑4173289
                                                                                                                                                 |                                                                                                                                                
gaagaaGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgc                                                                                                                                                >  2:2295741/1‑148 (MQ=255)
  agaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgcag                                                                                                                                              >  1:2918312/1‑148 (MQ=255)
  agaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgcag                                                                                                                                              >  2:2655044/1‑148 (MQ=255)
  agaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCctgcag                                                                                                                                              <  2:2931710/148‑1 (MQ=255)
   gaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAgg                                                                                                                                             >  2:3263386/1‑148 (MQ=255)
   gaagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAgg                                                                                                                                             >  1:986028/1‑148 (MQ=255)
    aagGCAACTACGTTATCTCCCATGCGAACTCCAACTTGGATGAATAAGGCCACTTCGTAGAAGACCTTGTAACTTGACTTAGCAAAGGCGAATACAGCTAGTTCAGCAGCGACCAGGTTGACTACATTGACGTATCCACCCATCAggt                                                                                                                                            >  1:1884881/1‑148 (MQ=255)
    aagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggt                                                                                                                                            >  2:2542183/1‑148 (MQ=255)
    aagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggt                                                                                                                                            >  2:3296916/1‑148 (MQ=255)
    aagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggt                                                                                                                                            >  2:742116/1‑148 (MQ=255)
     agGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAggtg                                                                                                                                           <  2:162886/148‑1 (MQ=255)
       gCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtggt                                                                                                                                         <  1:2712443/148‑1 (MQ=255)
       gCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAggtggt                                                                                                                                         <  1:2252928/148‑1 (MQ=255)
        cAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTa                                                                                                                                        <  2:2154881/148‑1 (MQ=255)
         aaCTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTAt                                                                                                                                       <  2:2697656/148‑1 (MQ=255)
           cTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                     >  2:363764/1‑148 (MQ=255)
           cTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATcc                                                                                                                                     >  2:3122880/1‑148 (MQ=255)
             aCGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGt                                                                                                                                   <  1:218976/148‑1 (MQ=255)
             aCGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGt                                                                                                                                   <  2:2926380/148‑1 (MQ=255)
              cGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTc                                                                                                                                  <  2:539325/148‑1 (MQ=255)
               gTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCg                                                                                                                                 <  2:2705524/148‑1 (MQ=255)
                ttATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCgg                                                                                                                                >  1:1232584/1‑148 (MQ=255)
                ttATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCgg                                                                                                                                <  1:252582/148‑1 (MQ=255)
                 tATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGt                                                                                                                               >  2:2380600/1‑148 (MQ=255)
                   tCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGc                                                                                                                             >  1:1407624/1‑148 (MQ=255)
                   tCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGc                                                                                                                             >  2:678681/1‑148 (MQ=255)
                       ccAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTcc                                                                                                                         <  1:3169152/148‑1 (MQ=255)
                        cAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTccc                                                                                                                        <  2:1741470/148‑1 (MQ=255)
                         aGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCt                                                                                                                       >  1:230795/1‑148 (MQ=255)
                         aGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCt                                                                                                                       >  1:229977/1‑148 (MQ=255)
                               aCTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCTCCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATccc                                                                                                                 >  2:2122482/1‑148 (MQ=255)
                               aCTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATccc                                                                                                                 <  1:1993075/148‑1 (MQ=255)
                               aCTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATccc                                                                                                                 >  2:882932/1‑148 (MQ=255)
                                   cAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTc                                                                                                             <  1:2357954/148‑1 (MQ=255)
                                       ttGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTgg                                                                                                         >  2:1394965/1‑148 (MQ=255)
                                         ggATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGaa                                                                                                       <  1:725751/148‑1 (MQ=255)
                                               agaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGAt                                                                                                 <  2:1693911/148‑1 (MQ=255)
                                                gaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATg                                                                                                >  1:2782637/1‑148 (MQ=255)
                                                 aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGa                                                                                               <  1:2851579/148‑1 (MQ=255)
                                                 aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGa                                                                                               >  2:3176017/1‑148 (MQ=255)
                                                    gCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGc                                                                                            <  1:453065/148‑1 (MQ=255)
                                                    gCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGc                                                                                            >  2:1846249/1‑148 (MQ=255)
                                                      cACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCa                                                                                          >  2:1409861/1‑148 (MQ=255)
                                                       aCTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCaa                                                                                         >  2:3023547/1‑148 (MQ=255)
                                                        cTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAAc                                                                                        >  2:1318178/1‑148 (MQ=255)
                                                        cTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTGTCCTGGAACACGATGACGCCAAc                                                                                        >  2:827479/1‑148 (MQ=255)
                                                              agaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCa                                                                                  >  2:1303149/1‑148 (MQ=255)
                                                               gaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAt                                                                                 <  1:2580270/148‑1 (MQ=255)
                                                                 agaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTg                                                                               >  2:4710/1‑148 (MQ=255)
                                                                    ccTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATg                                                                            <  2:655336/148‑1 (MQ=255)
                                                                      tGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATggg                                                                          <  1:2118390/148‑1 (MQ=255)
                                                                       ggTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTAGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGt                                                                         <  1:3343737/148‑1 (MQ=255)
                                                                         tAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGc                                                                       >  2:1978961/1‑148 (MQ=255)
                                                                          aaCTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCg                                                                      >  2:2328277/1‑148 (MQ=255)
                                                                            cTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGaa                                                                    >  1:927629/1‑148 (MQ=255)
                                                                             ttGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                                   >  1:3067662/1‑148 (MQ=255)
                                                                             ttGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                                   <  1:2939138/148‑1 (MQ=255)
                                                                              tGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACa                                                                  >  1:679578/1‑148 (MQ=255)
                                                                                ccGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATg                                                                <  1:750899/148‑1 (MQ=255)
                                                                                  gTAGCAAATTCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGTACGTATCCACCCTGAAGGTGGTATCCGACGGTGCGTCCCTGATCCCTTTCATGGAACACGAAGACGCCAACCGAGCATTCATGGGTGCGAACATGCa                                                              >  2:3344251/1‑148 (MQ=255)
                                                                                  gTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCa                                                              <  1:950830/148‑1 (MQ=255)
                                                                                   tAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCaa                                                             <  1:1846249/148‑1 (MQ=255)
                                                                                       aaaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTc                                                         <  2:1186271/148‑1 (MQ=255)
                                                                                        aaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                        >  1:3442761/1‑148 (MQ=255)
                                                                                          ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                      >  2:2410274/1‑148 (MQ=255)
                                                                                          ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                      <  2:1213891/148‑1 (MQ=255)
                                                                                           gCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                                     >  1:542646/1‑148 (MQ=255)
                                                                                              aaTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGt                                                  <  2:2960328/148‑1 (MQ=255)
                                                                                                tCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTc                                                >  1:2231507/1‑148 (MQ=255)
                                                                                                 ccAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTcc                                               <  1:1326399/148‑1 (MQ=255)
                                                                                                   aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                             >  1:1503911/1‑148 (MQ=255)
                                                                                                     cTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAct                                           >  1:1678234/1‑148 (MQ=255)
                                                                                                            aGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcg                                    <  1:1451082/148‑1 (MQ=255)
                                                                                                             gCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                   >  1:2061808/1‑148 (MQ=255)
                                                                                                             gCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                   >  2:1118552/1‑148 (MQ=255)
                                                                                                             gCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                   <  2:1492572/148‑1 (MQ=255)
                                                                                                             gCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                   >  1:2549880/1‑148 (MQ=255)
                                                                                                                 cgACCAGGTTGACTACATGGACGTATCCACCCTGCTGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGAt                               >  2:2893223/1‑148 (MQ=255)
                                                                                                                    ccAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                            >  2:3482103/1‑148 (MQ=255)
                                                                                                                     gaGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGc                           <  2:2066691/147‑1 (MQ=255)
                                                                                                                       ggTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                         <  2:2254123/148‑1 (MQ=255)
                                                                                                                       ggTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                         >  2:1998988/1‑148 (MQ=255)
                                                                                                                       ggTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                         <  1:1751592/148‑1 (MQ=255)
                                                                                                                        gTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                        >  2:408361/1‑148 (MQ=255)
                                                                                                                        gTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGc                        >  2:2524597/1‑148 (MQ=255)
                                                                                                                        gTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCGGATAAGCCGc                        >  1:1300566/1‑148 (MQ=255)
                                                                                                                           gACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtgg                     >  1:1073219/1‑148 (MQ=255)
                                                                                                                             cTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                   <  1:305052/148‑1 (MQ=255)
                                                                                                                                cATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt                <  1:633463/148‑1 (MQ=255)
                                                                                                                                  tGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTAc              >  1:1040821/1‑148 (MQ=255)
                                                                                                                                    gACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg            >  1:1638725/1‑148 (MQ=255)
                                                                                                                                    gACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg            <  1:3488266/148‑1 (MQ=255)
                                                                                                                                     aCGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTgg           >  2:603633/1‑148 (MQ=255)
                                                                                                                                      cGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGt          <  1:3296815/148‑1 (MQ=255)
                                                                                                                                         aTCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATg       >  2:2832424/1‑148 (MQ=255)
                                                                                                                                          tCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATgg      >  1:2475553/1‑148 (MQ=255)
                                                                                                                                           ccACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGa     >  2:2738514/1‑148 (MQ=255)
                                                                                                                                           ccACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGa     <  2:1560185/148‑1 (MQ=255)
                                                                                                                                              cccTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg  >  1:1775383/1‑148 (MQ=255)
                                                                                                                                              cccTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg  >  2:1159120/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACG  >  CP009273/4173000‑4173289

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: