Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,038,316 G→T 88.2% R451S (CGT→AGT)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,038,3160GT88.2% 390.4 / 50.8 161R451S (CGT→AGT) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (7/12);  new base T (70/72);  total (77/84)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 164 total aligned reads displayed.

GATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGC  >  CP009273/1038170‑1038460
                                                                                                                                                  |                                                                                                                                                
gATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCacga                                                                                                                                                 >  1:660459/1‑148 (MQ=255)
   atTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGc                                                                                                                                              >  1:1398648/1‑148 (MQ=255)
      tttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAAc                                                                                                                                           <  1:2343381/148‑1 (MQ=255)
      tttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAAc                                                                                                                                           <  1:2146795/148‑1 (MQ=255)
      tttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAAc                                                                                                                                           <  1:1747133/148‑1 (MQ=255)
        tAGTACGATGGCGCTGCTGATTAGAAAATACATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                         <  2:888770/148‑1 (MQ=255)
        tAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                         >  2:1176688/1‑148 (MQ=255)
        tAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                         <  1:1800973/148‑1 (MQ=255)
        tAGTACGAGGGCGCTGCTGATTAGAAAATAAATTTTTCTTCCGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                         <  2:749566/148‑1 (MQ=255)
           tACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGc                                                                                                                                      >  2:1567711/1‑148 (MQ=255)
            aCGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCa                                                                                                                                     >  1:1026168/1‑148 (MQ=255)
            aCGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCa                                                                                                                                     <  2:681386/148‑1 (MQ=255)
            aCGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCa                                                                                                                                     <  2:1829839/148‑1 (MQ=255)
            aCGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCa                                                                                                                                     <  1:129423/148‑1 (MQ=255)
                tGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCAcgcg                                                                                                                                 <  1:2237253/148‑1 (MQ=255)
                tGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCAcgcg                                                                                                                                 >  2:3059469/1‑148 (MQ=255)
                  gcgcTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGcc                                                                                                                               >  2:1507791/1‑148 (MQ=255)
                     ctgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGc                                                                                                                            >  1:658943/1‑148 (MQ=255)
                     ctgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGc                                                                                                                            >  1:536373/1‑148 (MQ=255)
                     ctgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGc                                                                                                                            >  1:3149006/1‑148 (MQ=255)
                      tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCa                                                                                                                           <  2:773487/148‑1 (MQ=255)
                      tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCa                                                                                                                           <  2:1128311/148‑1 (MQ=255)
                      tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCa                                                                                                                           <  2:812969/148‑1 (MQ=255)
                         tgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACgg                                                                                                                        >  1:2156308/1‑148 (MQ=255)
                          gATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGc                                                                                                                       <  1:3292779/148‑1 (MQ=255)
                           aTTAGAAAATAATTTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCa                                                                                                                      <  1:2740831/148‑1 (MQ=255)
                           aTTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCa                                                                                                                      <  1:2247721/148‑1 (MQ=255)
                              aGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCAccc                                                                                                                   >  2:1079242/1‑148 (MQ=255)
                                aaaaTAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCcac                                                                                                                 <  1:1453859/148‑1 (MQ=255)
                                   ataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGa                                                                                                              >  1:1524446/1‑148 (MQ=255)
                                   ataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGa                                                                                                              >  1:1697110/1‑148 (MQ=255)
                                    taaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaa                                                                                                             <  1:1706531/148‑1 (MQ=255)
                                    taaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaa                                                                                                             <  2:1176842/148‑1 (MQ=255)
                                        tttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaa                                                                                                         <  2:536373/148‑1 (MQ=255)
                                         ttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaaa                                                                                                        >  2:405953/1‑148 (MQ=255)
                                           ttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                      <  2:2913866/148‑1 (MQ=255)
                                           ttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                      <  1:2548326/148‑1 (MQ=255)
                                           ttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                      >  1:1109882/1‑148 (MQ=255)
                                               tACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaa                                                                                                         >  2:457304/1‑141 (MQ=255)
                                               tACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaa                                                                                                         <  1:457304/141‑1 (MQ=255)
                                                    aGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGaa                                                                                             <  1:1954790/148‑1 (MQ=255)
                                                     gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                            <  2:1410787/148‑1 (MQ=255)
                                                      aCGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATa                                                                                           >  1:102033/1‑148 (MQ=255)
                                                       cGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaa                                                                                          <  2:1235620/148‑1 (MQ=255)
                                                       cGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaa                                                                                          >  1:113170/1‑148 (MQ=255)
                                                         ggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaaaa                                                                                        >  1:957843/1‑148 (MQ=255)
                                                              gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                   <  2:2409001/148‑1 (MQ=255)
                                                              gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                   <  2:2409159/148‑1 (MQ=255)
                                                              gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                   <  2:1740139/148‑1 (MQ=255)
                                                               tttATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                                  <  2:1239690/148‑1 (MQ=255)
                                                                ttATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGca                                                                                 >  2:420938/1‑148 (MQ=255)
                                                                ttATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGca                                                                                 >  1:2494617/1‑148 (MQ=255)
                                                                  aTCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcacc                                                                               <  1:377066/148‑1 (MQ=255)
                                                                   tCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcacca                                                                              <  1:191385/148‑1 (MQ=255)
                                                                   tCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcacca                                                                              <  2:621760/148‑1 (MQ=255)
                                                                     cAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACg                                                                            <  2:3143803/148‑1 (MQ=255)
                                                                     cAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACg                                                                            <  2:3065037/148‑1 (MQ=255)
                                                                        ccACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTc                                                                         >  1:2425741/1‑148 (MQ=255)
                                                                         cACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCa                                                                        <  2:2537254/148‑1 (MQ=255)
                                                                         cACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCa                                                                        <  2:2535952/148‑1 (MQ=255)
                                                                           cTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAAc                                                                      >  2:1492142/1‑148 (MQ=255)
                                                                             cGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGc                                                                    >  2:1314854/1‑148 (MQ=255)
                                                                                 cATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttt                                                                >  2:1006483/1‑148 (MQ=255)
                                                                                   ttGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttttt                                                              >  1:113477/1‑148 (MQ=255)
                                                                                   ttGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttttt                                                              >  2:2807866/1‑148 (MQ=255)
                                                                                     gggATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCg                                                            <  1:269903/148‑1 (MQ=255)
                                                                                     gggATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCg                                                            <  1:270104/148‑1 (MQ=255)
                                                                                      ggATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCgg                                                           >  2:2542733/1‑148 (MQ=255)
                                                                                      ggATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCgg                                                           <  2:1667081/148‑1 (MQ=255)
                                                                                       gATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGt                                                          >  2:2288294/1‑148 (MQ=255)
                                                                                              gCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAATCAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTg                                                   <  2:2964934/148‑1 (MQ=255)
                                                                                                attcaCGGTGATGGCGTGGTCTTCCAGTTGTTCGGGGGCTTCTAGACAAGTACGCAACTTCGCACGCGCCAGCACTGCACCCACAGAAATAAACAGTCCAAGAATAAAACCAAGCACCACTTTAAAAGCTTTTTTCGGTTTCACTGGc                                                 <  1:153884/144‑1 (MQ=255)
                                                                                                  aaaCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTg                                               <  1:1945002/148‑1 (MQ=255)
                                                                                                  aaaCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTg                                               >  2:1517668/1‑148 (MQ=255)
                                                                                                    aCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCg                                             >  2:2369710/1‑148 (MQ=255)
                                                                                                     cGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCgg                                            <  2:546291/148‑1 (MQ=255)
                                                                                                        tGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTg                                         <  1:346428/148‑1 (MQ=255)
                                                                                                            gCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTg                                     >  2:69426/1‑148 (MQ=255)
                                                                                                                tGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTg                                 >  2:2050816/1‑148 (MQ=255)
                                                                                                                    ctTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCgg                             >  1:1535776/1‑148 (MQ=255)
                                                                                                                    ctTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCgg                             >  1:1810382/1‑148 (MQ=255)
                                                                                                                     ttCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCggg                            <  1:2933541/148‑1 (MQ=255)
                                                                                                                      tCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGt                           <  2:1109882/148‑1 (MQ=255)
                                                                                                                        cAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCg                         <  2:1903661/148‑1 (MQ=255)
                                                                                                                          gttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGAt                       >  2:266817/1‑148 (MQ=255)
                                                                                                                          gttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGCCCAATAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGTCTGCGGCTGAGTGACTGCCGGGTCGAt                       >  2:260556/1‑148 (MQ=255)
                                                                                                                           ttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATa                      >  2:2311843/1‑148 (MQ=255)
                                                                                                                           ttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATa                      <  1:754400/148‑1 (MQ=255)
                                                                                                                            tgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCGCGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATaa                     <  2:2555846/148‑1 (MQ=255)
                                                                                                                            tgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATaa                     >  1:3061171/1‑148 (MQ=255)
                                                                                                                            tgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATaa                     <  2:2555824/148‑1 (MQ=255)
                                                                                                                             gttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAAt                    <  2:2494617/148‑1 (MQ=255)
                                                                                                                                 cGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGc                <  1:1223366/148‑1 (MQ=255)
                                                                                                                                   gggCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCAc              >  2:574922/1‑148 (MQ=255)
                                                                                                                                    ggCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACg             <  2:1398648/148‑1 (MQ=255)
                                                                                                                                    ggCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACg             >  2:563398/1‑148 (MQ=255)
                                                                                                                                          tACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCa       <  1:2168220/148‑1 (MQ=255)
                                                                                                                                           acacCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCaa      >  1:1212931/1‑148 (MQ=255)
                                                                                                                                             accacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg    <  2:1026168/148‑1 (MQ=255)
                                                                                                                                               cacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGc  <  2:2915476/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGC  >  CP009273/1038170‑1038460

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATG  >  CP009273/1038167‑1038458
                                                                                                                                                     |                                                                                                                                              
GGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACT                                                                                                                                                >  SRR19973070.1402709/1‑150 (MQ=60)
GGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACT                                                                                                                                                >  SRR19973070.1402593/1‑150 (MQ=60)
 GGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGA                                                                                                                                               >  SRR19973070.665161/1‑150 (MQ=60)
    ATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGC                                                                                                                                            >  SRR19973070.1408067/1‑150 (MQ=60)
         TTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAT                                                                                                                                       <  SRR19973070.1759091/150‑1 (MQ=60)
         TTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAT                                                                                                                                       <  SRR19973070.1759247/150‑1 (MQ=60)
         TTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAT                                                                                                                                       <  SRR19973070.2358925/150‑1 (MQ=60)
         TTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAT                                                                                                                                       <  SRR19973070.2161103/150‑1 (MQ=60)
           TAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCG                                                                                                                                     <  SRR19973070.1813412/150‑1 (MQ=60)
            AGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCCTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGC                                                                                                                                    >  SRR19973070.2847162/1‑150 (MQ=60)
             GTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCA                                                                                                                                   >  SRR19973070.1033444/1‑150 (MQ=60)
             GTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCA                                                                                                                                   >  SRR19973070.1034115/1‑150 (MQ=60)
               ACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACG                                                                                                                                 <  SRR19973070.130561/150‑1 (MQ=60)
                  ATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGC                                                                                                                              >  SRR19973070.1769143/1‑150 (MQ=60)
                   TGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCC                                                                                                                             <  SRR19973070.2252014/150‑1 (MQ=60)
                      CGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGC                                                                                                                          >  SRR19973070.3169353/1‑150 (MQ=60)
                      CGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGC                                                                                                                          >  SRR19973070.539970/1‑150 (MQ=60)
                      CGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGC                                                                                                                          >  SRR19973070.663636/1‑150 (MQ=60)
                      CGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGC                                                                                                                          >  SRR19973070.1518984/1‑150 (MQ=60)
                          GCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGG                                                                                                                      >  SRR19973070.2170645/1‑150 (MQ=60)
                          GCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGG                                                                                                                      >  SRR19973070.2170484/1‑150 (MQ=60)
                             GATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCAC                                                                                                                   <  SRR19973070.3313897/150‑1 (MQ=60)
                              ATTAGAAAATAATTTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACC                                                                                                                  <  SRR19973070.2758464/150‑1 (MQ=60)
                              ATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACC                                                                                                                  <  SRR19973070.2262588/150‑1 (MQ=60)
                                   AAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAG                                                                                                             <  SRR19973070.1463984/150‑1 (MQ=60)
                                    AAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGA                                                                                                            >  SRR19973070.1708796/1‑150 (MQ=60)
                                    AAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGA                                                                                                            >  SRR19973070.1535151/1‑150 (MQ=60)
                                      ATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAA                                                                                                          <  SRR19973070.351308/150‑1 (MQ=60)
                                       TAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAAT                                                                                                         <  SRR19973070.1718308/150‑1 (MQ=60)
                                         AATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTNCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAA                                                                                                       >  SRR19973070.1691567/1‑150 (MQ=60)
                                         AATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAA                                                                                                       >  SRR19973070.510177/1‑150 (MQ=60)
                                           gatcTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAAC                                                                                                     <  SRR19973070.460578/146‑1 (MQ=60)
                                           TTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAAC                                                                                                     >  SRR19973070.2060976/1‑150 (MQ=60)
                                           TTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAAC                                                                                                     <  SRR19973070.3118832/150‑1 (MQ=60)
                                            TTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACA                                                                                                    >  SRR19973070.1117477/1‑150 (MQ=60)
                                              TTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGG                                                                                                  <  SRR19973070.2564990/150‑1 (MQ=60)
                                                CTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCC                                                                                                >  SRR19973070.2438986/1‑150 (MQ=60)
                                                       AGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATA                                                                                         >  SRR19973070.1833846/1‑150 (MQ=60)
                                                       AGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATA                                                                                         <  SRR19973070.1967944/150‑1 (MQ=60)
                                                       AGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATA                                                                                         >  SRR19973070.103004/1‑150 (MQ=60)
                                                        GACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAA                                                                                        >  SRR19973070.114219/1‑150 (MQ=60)
                                                        GACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAA                                                                                        >  SRR19973070.114211/1‑150 (MQ=60)
                                                          CGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAA                                                                                      >  SRR19973070.964677/1‑150 (MQ=60)
                                                                 GTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCA                                                                               >  SRR19973070.2510796/1‑150 (MQ=60)
                                                                   TTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACC                                                                             <  SRR19973070.2753584/150‑1 (MQ=60)
                                                                    TATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCA                                                                            >  SRR19973070.783514/1‑150 (MQ=60)
                                                                     ATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCAC                                                                           <  SRR19973070.379828/150‑1 (MQ=60)
                                                                      TCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACG                                                                          <  SRR19973070.192859/150‑1 (MQ=60)
                                                                         AGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTC                                                                       >  SRR19973070.2441684/1‑150 (MQ=60)
                                                                                    CATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTT                                                            >  SRR19973070.114523/1‑150 (MQ=60)
                                                                                        GGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGT                                                        <  SRR19973070.271677/150‑1 (MQ=60)
                                                                                        GGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGT                                                        <  SRR19973070.271879/150‑1 (MQ=60)
                                                                                          GATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTT                                                      >  SRR19973070.2613907/1‑150 (MQ=60)
                                                                                           ATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTC                                                     >  SRR19973070.2522242/1‑150 (MQ=60)
                                                                                             AGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAC                                                   >  SRR19973070.2206815/1‑150 (MQ=60)
                                                                                                   attcacggTGATGGCGTGGTCTTCCAGTTGTTCGGGGGCTTCTAGACAAGTACGCAACTTCGCACGCGCCAGCACTGCACCCACAGAAATAAACAGTCCAAGAATAAAACCAAGCACCACTTTAAAAGCTTTTTTCGGTTTCACTGGCTG                                             <  SRR19973070.155169/142‑1 (MQ=60)
                                                                                                     AAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCG                                           <  SRR19973070.1958064/150‑1 (MQ=60)
                                                                                                     AAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCG                                           <  SRR19973070.1957927/150‑1 (MQ=60)
                                                                                                      AACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGG                                          >  SRR19973070.1283424/1‑150 (MQ=60)
                                                                                                          CTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGA                                      <  SRR19973070.724501/150‑1 (MQ=60)
                                                                                                           TGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAG                                     <  SRR19973070.348950/150‑1 (MQ=60)
                                                                                                                CCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACT                                <  SRR19973070.1539081/150‑1 (MQ=60)
                                                                                                                 CGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTG                               <  SRR19973070.321017/150‑1 (MQ=60)
                                                                                                                     CTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGG                           >  SRR19973070.1546559/1‑150 (MQ=60)
                                                                                                                     CTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGG                           >  SRR19973070.1822862/1‑150 (MQ=60)
                                                                                                                        TTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTC                        <  SRR19973070.2952462/150‑1 (MQ=60)
                                                                                                                           CAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGAT                     >  SRR19973070.610334/1‑150 (MQ=60)
                                                                                                                             GTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAA                   >  SRR19973070.3080766/1‑150 (MQ=60)
                                                                                                                              TTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAAT                  <  SRR19973070.759723/150‑1 (MQ=60)
                                                                                                                                    CGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCAC            <  SRR19973070.1231929/150‑1 (MQ=60)
                                                                                                                                     GGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACG           <  SRR19973070.2262556/150‑1 (MQ=60)
                                                                                                                                       GGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTT         <  SRR19973070.1358887/150‑1 (MQ=60)
                                                                                                                                            CTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAA    >  SRR19973070.1221445/1‑150 (MQ=60)
                                                                                                                                             TACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAT   <  SRR19973070.2182601/150‑1 (MQ=60)
                                                                                                                                              ACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATG  <  SRR19973070.1184941/150‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                              
GGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATG  >  CP009273/1038167‑1038458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑