Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,197,734 C→T 6.3% R219K (AGA→AAA)  cohE ← e14 prophage; repressor protein phage e14

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,197,7340CT 6.3% 1030.1 / 37.5 286R219K (AGA→AAA) cohEe14 prophage; repressor protein phage e14
Reads supporting (aligned to +/- strand):  ref base C (149/119);  new base T (11/7);  total (160/126)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 309 total aligned reads displayed.

TCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGT  >  CP009273/1197586‑1197882
                                                                                                                                                    |                                                                                                                                                    
tCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtc                                                                                                                                                      <  1:4161486/149‑1 (MQ=255)
tCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtc                                                                                                                                                      <  1:1961174/149‑1 (MQ=255)
  gtacgtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtctt                                                                                                                                                    <  2:3853714/149‑1 (MQ=255)
   tacgtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtcttc                                                                                                                                                   >  1:4434197/1‑149 (MQ=255)
   tacgtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtcttc                                                                                                                                                   <  1:645828/149‑1 (MQ=255)
   tacgtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtcttc                                                                                                                                                   >  1:2136105/1‑149 (MQ=255)
       tacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTgg                                                                                                                                               <  2:1646539/149‑1 (MQ=255)
          cgggcgggCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGttt                                                                                                                                            <  1:925354/149‑1 (MQ=255)
           gggcgggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGtttt                                                                                                                                           >  2:1501401/1‑149 (MQ=255)
               gggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTg                                                                                                                                       >  2:2904255/1‑149 (MQ=255)
                ggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTga                                                                                                                                      >  2:4940247/1‑149 (MQ=255)
                  cTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTgaga                                                                                                                                    >  2:180837/1‑149 (MQ=255)
                                  ctgcctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAAt                                                                                                                    >  1:1731048/1‑149 (MQ=255)
                                   tgcctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAAtt                                                                                                                   <  2:3350852/149‑1 (MQ=255)
                                     cctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATtat                                                                                                                 >  1:1156452/1‑149 (MQ=255)
                                              tgttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAAtt                                                                                                        >  1:3828510/1‑149 (MQ=255)
                                              tgttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAAtt                                                                                                        <  2:683735/149‑1 (MQ=255)
                                              tgttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAAtt                                                                                                        >  2:3758637/1‑149 (MQ=255)
                                              tgttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAAtt                                                                                                        >  2:4333530/1‑149 (MQ=255)
                                              tgttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAAtt                                                                                                        >  1:4786168/1‑149 (MQ=255)
                                               gttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTa                                                                                                       >  2:1011145/1‑149 (MQ=255)
                                                 ttCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTAcc                                                                                                     >  2:1416572/1‑149 (MQ=255)
                                                  tctcTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCg                                                                                                    >  2:134895/1‑149 (MQ=255)
                                                    tctTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGtt                                                                                                  >  1:2429651/1‑149 (MQ=255)
                                                      ttGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGa                                                                                                <  2:3073713/149‑1 (MQ=255)
                                                      ttGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGa                                                                                                <  1:3104276/149‑1 (MQ=255)
                                                            aCATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTaa                                                                                          >  1:328869/1‑149 (MQ=255)
                                                            aCATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTaa                                                                                          >  1:4000869/1‑149 (MQ=255)
                                                              aTTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATg                                                                                        <  2:3452847/149‑1 (MQ=255)
                                                               ttATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGt                                                                                       >  2:4457462/1‑149 (MQ=255)
                                                                tATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTa                                                                                      >  2:1143324/1‑149 (MQ=255)
                                                                  tGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAg                                                                                    >  2:3555821/1‑149 (MQ=255)
                                                                         cAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAgg                                                                             >  2:5099225/1‑149 (MQ=255)
                                                                          aaaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGc                                                                            <  1:1007751/149‑1 (MQ=255)
                                                                           aaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATATAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGcc                                                                           >  2:502110/1‑149 (MQ=255)
                                                                           aaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGcc                                                                           <  2:4389379/149‑1 (MQ=255)
                                                                           aaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGcc                                                                           >  1:3386788/1‑149 (MQ=255)
                                                                           aaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGcc                                                                           >  1:4902360/1‑149 (MQ=255)
                                                                            aGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCa                                                                          >  1:3273241/1‑149 (MQ=255)
                                                                              gTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGt                                                                        >  1:345871/1‑149 (MQ=255)
                                                                              gTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGt                                                                        >  1:1031585/1‑149 (MQ=255)
                                                                              gTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGt                                                                        >  1:4473932/1‑149 (MQ=255)
                                                                              gTACAGTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGCTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGt                                                                        <  2:2822932/149‑1 (MQ=255)
                                                                                aCATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGttt                                                                      <  1:4705504/149‑1 (MQ=255)
                                                                                 cATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTg                                                                     <  1:373398/149‑1 (MQ=255)
                                                                                 cATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTg                                                                     <  1:2335889/149‑1 (MQ=255)
                                                                                   ttGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTggg                                                                   <  1:2412255/149‑1 (MQ=255)
                                                                                   ttGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTggg                                                                   <  2:3912659/149‑1 (MQ=255)
                                                                                    tGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGACAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATTTAAGGCTCCGGCCAGTTTGGGt                                                                  >  2:2405570/1‑149 (MQ=255)
                                                                                       cAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTa                                                               >  2:1878447/1‑149 (MQ=255)
                                                                                         aGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAAc                                                             <  2:613472/149‑1 (MQ=255)
                                                                                              aaaaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCttt                                                        <  1:3156719/149‑1 (MQ=255)
                                                                                               aaaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTg                                                       >  1:3523808/1‑149 (MQ=255)
                                                                                                 aaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTgag                                                     >  2:4017764/1‑149 (MQ=255)
                                                                                                  aaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTgaga                                                    >  1:3724971/1‑149 (MQ=255)
                                                                                                   aaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGAt                                                   >  2:4806917/1‑149 (MQ=255)
                                                                                                   aaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGAt                                                   >  2:211806/1‑149 (MQ=255)
                                                                                                    aCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATa                                                  >  2:2234931/1‑149 (MQ=255)
                                                                                                     ccagCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATaa                                                 <  1:614455/146‑1 (MQ=255)
                                                                                                      cTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAAc                                                >  1:1915956/1‑149 (MQ=255)
                                                                                                        gCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGc                                              <  1:1661041/149‑1 (MQ=255)
                                                                                                          cGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCtg                                            >  1:3916876/1‑149 (MQ=255)
                                                                                                           gAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCtgt                                           >  2:4161824/1‑149 (MQ=255)
                                                                                                            aaGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCtgtg                                          >  1:1223624/1‑149 (MQ=255)
                                                                                                             aGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTgtgt                                         >  1:2708052/1‑149 (MQ=255)
                                                                                                             aGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTgtgt                                         >  1:5025170/1‑149 (MQ=255)
                                                                                                             aGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTgtgt                                         <  1:4361945/149‑1 (MQ=255)
                                                                                                              gCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTc                                        >  1:976390/1‑149 (MQ=255)
                                                                                                                aGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTccc                                      <  1:5054169/149‑1 (MQ=255)
                                                                                                                 ggTTCATAAACATTGATTAGGCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCa                                     >  1:294749/1‑149 (MQ=255)
                                                                                                                    tCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtct                                  <  1:2461422/149‑1 (MQ=255)
                                                                                                                      aTAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtcttc                                <  2:4431001/149‑1 (MQ=255)
                                                                                                                      aTAAACATTGATTAGGCTTTGATTTTGTATCTTCGTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtcttc                                <  2:2592106/149‑1 (MQ=255)
                                                                                                                       tAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtcttct                               >  2:4116907/1‑149 (MQ=255)
                                                                                                                        aaaCATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTa                              <  1:1850945/149‑1 (MQ=255)
                                                                                                                        aaaCATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTa                              <  1:277755/149‑1 (MQ=255)
                                                                                                                         aaCATTGATTAGGCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTAt                             >  2:4615009/1‑149 (MQ=255)
                                                                                                                         aaCATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTAt                             >  2:5146905/1‑149 (MQ=255)
                                                                                                                         aaCATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTAt                             >  1:2809874/1‑149 (MQ=255)
                                                                                                                          aCATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATc                            >  2:4704255/1‑149 (MQ=255)
                                                                                                                          aCATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGGTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATc                            >  2:1815207/1‑149 (MQ=255)
                                                                                                                             ttGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAAc                         <  2:2451631/149‑1 (MQ=255)
                                                                                                                             ttGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAAc                         >  1:818108/1‑149 (MQ=255)
                                                                                                                             ttGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAAc                         >  1:759694/1‑149 (MQ=255)
                                                                                                                             ttGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAAc                         <  1:2075769/149‑1 (MQ=255)
                                                                                                                               gATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCt                       <  2:802996/149‑1 (MQ=255)
                                                                                                                                 ttAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCttt                     >  1:1346305/1‑149 (MQ=255)
                                                                                                                                  tAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCtttt                    <  1:5099287/149‑1 (MQ=255)
                                                                                                                                    ggCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGa                  <  2:153927/149‑1 (MQ=255)
                                                                                                                                     gCTTTGATTTTGTATTTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGaa                 <  1:885330/149‑1 (MQ=255)
                                                                                                                                      cTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAg                >  2:193092/1‑149 (MQ=255)
                                                                                                                                        ttGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggt              <  2:1239538/149‑1 (MQ=255)
                                                                                                                                          gATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggtgg            <  2:2958015/149‑1 (MQ=255)
                                                                                                                                           aTTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggtggt           >  1:1429412/1‑149 (MQ=255)
                                                                                                                                                 tATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCAcc     <  2:3681703/149‑1 (MQ=255)
                                                                                                                                                 tATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCAcc     >  1:2959012/1‑149 (MQ=255)
                                                                                                                                                 tATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCAcc     >  1:1887109/1‑149 (MQ=255)
                                                                                                                                                   tcttctTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTg   >  2:866145/1‑149 (MQ=255)
                                                                                                                                                    cttcttGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGt  <  2:2743040/149‑1 (MQ=255)
                                                                                                                                                    cttcttGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGt  <  1:3758701/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTTTCACCTGT  >  CP009273/1197586‑1197882

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: