Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,507,551 G→T 6.1% R168S (CGT→AGT)  fecI ← RNA polymerase, sigma 19 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,507,5510GT 6.1% 1137.7 / 29.5 296R168S (CGT→AGT) fecIRNA polymerase, sigma 19 factor
Reads supporting (aligned to +/- strand):  ref base G (146/132);  new base T (9/9);  total (155/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 304 total aligned reads displayed.

CCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCA  >  CP009273/4507402‑4507700
                                                                                                                                                     |                                                                                                                                                     
cctgatcctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACg                                                                                                                                                       <  2:2066012/150‑1 (MQ=255)
  tgatcctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACTGa                                                                                                                                                     >  2:2772152/1‑150 (MQ=255)
  tgatcctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGa                                                                                                                                                     >  2:4364612/1‑150 (MQ=255)
   gatcctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGACTCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGaa                                                                                                                                                    >  2:4484456/1‑150 (MQ=255)
   gatcctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGaa                                                                                                                                                    <  2:666082/150‑1 (MQ=255)
     tcctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAAca                                                                                                                                                  >  1:3913460/1‑150 (MQ=255)
     tcctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAAca                                                                                                                                                  >  2:89125/1‑150 (MQ=255)
      cctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAAcag                                                                                                                                                 >  2:2576627/1‑150 (MQ=255)
       ctgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAAcagc                                                                                                                                                <  2:496054/150‑1 (MQ=255)
        tgTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAAcagca                                                                                                                                               <  2:3957493/150‑1 (MQ=255)
            cATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAgcaggca                                                                                                                                           >  2:2537/1‑150 (MQ=255)
              tACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGt                                                                                                                                         <  2:599216/150‑1 (MQ=255)
                   cTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCg                                                                                                                                    <  1:2576278/150‑1 (MQ=255)
                    tGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGa                                                                                                                                   <  1:2033881/150‑1 (MQ=255)
                           aGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCttt                                                                                                                            >  2:4048322/1‑150 (MQ=255)
                            gcgcgcTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTc                                                                                                                           <  1:4212/150‑1 (MQ=255)
                            gcgcgcTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTc                                                                                                                           <  1:1414790/150‑1 (MQ=255)
                            gcgcgcTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTc                                                                                                                           <  1:4188/150‑1 (MQ=255)
                             cgcgcgTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCg                                                                                                                          <  1:3119203/150‑1 (MQ=255)
                             cgcgcTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCg                                                                                                                          <  1:3898545/150‑1 (MQ=255)
                               cgcTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGcc                                                                                                                        >  1:3760084/1‑150 (MQ=255)
                                 cgtCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCAc                                                                                                                      <  1:3167159/148‑1 (MQ=255)
                                 cgtCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCAc                                                                                                                      <  1:849615/148‑1 (MQ=255)
                                   tCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGt                                                                                                                    <  2:2154149/150‑1 (MQ=255)
                                    ccTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTa                                                                                                                   >  1:4525586/1‑150 (MQ=255)
                                      tgttgtGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACTGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTAtt                                                                                                                 <  2:843562/150‑1 (MQ=255)
                                       gttgtGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTAttt                                                                                                                <  1:2374710/150‑1 (MQ=255)
                                          gtgGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTAtttttt                                                                                                             <  2:1448199/150‑1 (MQ=255)
                                          gtgGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTAtttttt                                                                                                             <  2:632516/150‑1 (MQ=255)
                                           tgGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTc                                                                                                            <  1:4326774/150‑1 (MQ=255)
                                                   cgcgCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCt                                                                                                    <  2:1514106/150‑1 (MQ=255)
                                                   cgcgCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCt                                                                                                    <  2:3174977/150‑1 (MQ=255)
                                                    gcgcTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTg                                                                                                   <  1:1219315/150‑1 (MQ=255)
                                                        tCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATgg                                                                                               <  2:2194526/150‑1 (MQ=255)
                                                         cgccgcTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGa                                                                                              >  2:2309963/1‑150 (MQ=255)
                                                          gcctcTTAGCACGGCATACCAGTGGGAAGCTGCACGCAGCGCCTGACGGCGTGTATCGGTTAACTAAGGATTCATAACCCATACTCAAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCAAGTATTTTTTCACGGAGCTGATGGaa                                                                                             <  1:2215812/150‑1 (MQ=255)
                                                          gccgcTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGaa                                                                                             >  1:1160741/1‑150 (MQ=255)
                                                             gcTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAaca                                                                                          <  2:2470911/150‑1 (MQ=255)
                                                               ttAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACAcc                                                                                        >  1:566923/1‑150 (MQ=255)
                                                                tAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCg                                                                                       >  1:746069/1‑150 (MQ=255)
                                                                    aCGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGtt                                                                                   >  2:1233492/1‑150 (MQ=255)
                                                                    aCGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGtt                                                                                   >  2:1803970/1‑150 (MQ=255)
                                                                      ggCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTtg                                                                                 >  1:1743609/1‑150 (MQ=255)
                                                                         aTACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTgc                                                                              >  2:1989942/1‑150 (MQ=255)
                                                                         aTACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTgc                                                                              >  2:4317629/1‑150 (MQ=255)
                                                                           aCCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTgcgc                                                                            >  1:2486973/1‑150 (MQ=255)
                                                                           aCCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTgcgc                                                                            >  1:1172483/1‑150 (MQ=255)
                                                                            ccAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCa                                                                           <  1:4099392/150‑1 (MQ=255)
                                                                            ccAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCa                                                                           >  1:2798877/1‑150 (MQ=255)
                                                                             cAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCaa                                                                          >  2:3720498/1‑150 (MQ=255)
                                                                              aGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAAt                                                                         >  2:3892537/1‑150 (MQ=255)
                                                                              aGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAAt                                                                         >  1:3354863/1‑150 (MQ=255)
                                                                               gTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAAtc                                                                        >  2:4556601/1‑150 (MQ=255)
                                                                                 gggAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAAtctc                                                                      >  2:1714496/1‑150 (MQ=255)
                                                                                       cTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTa                                                                <  1:1125222/150‑1 (MQ=255)
                                                                                        tGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTAt                                                               <  1:3778104/150‑1 (MQ=255)
                                                                                         gAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATg                                                              >  1:352804/1‑150 (MQ=255)
                                                                                         gAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATg                                                              >  2:2101118/1‑150 (MQ=255)
                                                                                          aaCGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGt                                                             <  2:1578513/150‑1 (MQ=255)
                                                                                          aaCGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGt                                                             <  1:3057563/150‑1 (MQ=255)
                                                                                           aCGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACTGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTc                                                            <  1:2311545/150‑1 (MQ=255)
                                                                                             gCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAg                                                          <  1:1915552/150‑1 (MQ=255)
                                                                                              cAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGa                                                         >  2:3402126/1‑150 (MQ=255)
                                                                                                  gcCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCAt                                                     <  1:90063/150‑1 (MQ=255)
                                                                                                    cTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATcc                                                   >  1:1633740/1‑150 (MQ=255)
                                                                                                     tGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACTGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCa                                                  >  1:3678070/1‑150 (MQ=255)
                                                                                                     tGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGGTTGTGCGCAATCTCGCTGTATGTCAGACCATCCa                                                  >  2:1089588/1‑150 (MQ=255)
                                                                                                       aCGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGt                                                <  2:1631610/150‑1 (MQ=255)
                                                                                                         ggcggGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTg                                              >  2:131461/1‑150 (MQ=255)
                                                                                                         ggcggGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTg                                              >  2:28909/1‑150 (MQ=255)
                                                                                                           cggGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCg                                            <  2:3876737/150‑1 (MQ=255)
                                                                                                             ggAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGaa                                          >  1:540091/1‑150 (MQ=255)
                                                                                                             ggAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGaa                                          >  2:70425/1‑150 (MQ=255)
                                                                                                              gAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGaaa                                         >  2:4169128/1‑150 (MQ=255)
                                                                                                              gAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGaaa                                         <  2:1901906/150‑1 (MQ=255)
                                                                                                                aTCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGc                                       >  2:2330116/1‑150 (MQ=255)
                                                                                                                aTCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGc                                       <  2:3913460/150‑1 (MQ=255)
                                                                                                                aTCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGc                                       <  2:1784418/150‑1 (MQ=255)
                                                                                                                   ggTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGtt                                                                                   >  1:2454105/1‑103 (MQ=255)
                                                                                                                   ggTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGtt                                                                                   <  2:2454105/103‑1 (MQ=255)
                                                                                                                     ttAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGaaa                                  <  2:3738992/150‑1 (MQ=255)
                                                                                                                        aCAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGc                               >  2:484049/1‑150 (MQ=255)
                                                                                                                         cAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCt                              <  1:3001155/150‑1 (MQ=255)
                                                                                                                         cAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCt                              >  2:89896/1‑150 (MQ=255)
                                                                                                                         cAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCt                              >  2:2205892/1‑150 (MQ=255)
                                                                                                                         cAAAGGATTCATAACCCATACCCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCt                              >  1:2239824/1‑150 (MQ=255)
                                                                                                                              gATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACg                         <  2:1062118/150‑1 (MQ=255)
                                                                                                                                ttCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACgtg                       >  2:45982/1‑150 (MQ=255)
                                                                                                                                 tCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACgtgt                      <  2:2708766/150‑1 (MQ=255)
                                                                                                                                    tAACCCATACTCCAGACTGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGtttt                   >  1:291171/1‑150 (MQ=255)
                                                                                                                                    tAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGtttt                   >  1:2851743/1‑150 (MQ=255)
                                                                                                                                    tAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGtttt                   >  1:2679343/1‑150 (MQ=255)
                                                                                                                                     aaCCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTg                  >  2:3772709/1‑150 (MQ=255)
                                                                                                                                     aaCCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTg                  >  2:2509743/1‑150 (MQ=255)
                                                                                                                                       cccATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGcc                <  1:3403285/150‑1 (MQ=255)
                                                                                                                                            aCTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTa           <  2:428429/150‑1 (MQ=255)
                                                                                                                                            aCTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTa           <  2:3002179/150‑1 (MQ=255)
                                                                                                                                            aCTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTa           >  1:1012550/1‑150 (MQ=255)
                                                                                                                                            aCTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGGTTTGCCGTTTa           >  2:4147401/1‑150 (MQ=255)
                                                                                                                                                     ggAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCa  <  1:4039876/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAACAAAGGATTCATAACCCATACTCCAGACGGAACAGCAGGCAGTGCTCGACGGCTTTCGCCACGTATTTTTTCACGGAGCTGATGGAAACACCGAGTTTGTGCGCAATCTCGCTGTATGTCAGACCATCCAGTTGCGAAAGCAGAAACGCTTCACGTGTTTTGCCGTTTAGCCCGTCCA  >  CP009273/4507402‑4507700

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: