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 F434 I0 R1 1130 215.9 7031662 97.1% 6827743 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 619,048 C→A 100% A307S (GCC→TCC)  fepB ← iron‑enterobactin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273619,0480CA97.2% 440.7 / 4.1 142A307S (GCC→TCC) fepBiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (66/71);  total (68/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 142 total aligned reads displayed.

CGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCC  >  CP009273/618904‑619192
                                                                                                                                                |                                                                                                                                                
cGGTTTTGGTTTGTTGATGATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCt                                                                                                                                               <  2:3425149/148‑1 (MQ=255)
   ttttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgta                                                                                                                                            >  1:1655099/1‑148 (MQ=255)
    tttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgtag                                                                                                                                           >  2:1126565/1‑148 (MQ=255)
    tttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgtag                                                                                                                                           >  2:1126746/1‑148 (MQ=255)
        gTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAAt                                                                                                                                       <  2:1785873/148‑1 (MQ=255)
         tttGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATc                                                                                                                                      >  1:2621844/1‑148 (MQ=255)
         tttGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCATATCCAGCACAAGCATGGAGCTGTAGTAATc                                                                                                                                      >  2:1774502/1‑148 (MQ=255)
         tttGTTGATTATCGGCGTGGTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATc                                                                                                                                      >  1:2513502/1‑148 (MQ=255)
          ttgttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATTCTTAAAAAAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATcc                                                                                                                                     >  1:321767/1‑148 (MQ=255)
          ttgttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATcc                                                                                                                                     >  1:192188/1‑148 (MQ=255)
          ttgttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATcc                                                                                                                                     >  2:1238399/1‑148 (MQ=255)
            gttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAg                                                                                                                                   <  1:3444384/148‑1 (MQ=255)
            gttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAg                                                                                                                                   <  2:1062695/148‑1 (MQ=255)
             ttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGa                                                                                                                                  >  1:532662/1‑148 (MQ=255)
                 ttATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGa                                                                                                                              >  2:2527989/1‑148 (MQ=255)
                    tCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACg                                                                                                                           <  1:1822786/148‑1 (MQ=255)
                     cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          <  1:1811783/148‑1 (MQ=255)
                     cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          <  1:1757699/148‑1 (MQ=255)
                     cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACACCGCCTTAATCCTATCCACCACTTGCAAGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          >  1:1084011/1‑148 (MQ=255)
                     cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGAGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          <  1:2903862/148‑1 (MQ=255)
                     cGGCGTGGTGTTATTTCTGGTGCTGGTGGAGTTGCGACAATTTCGCCAGACGCCGCCGCACGTGACAGCGTCCGACACTTAATGCTTAAAAAAGCGCCATAAGCCAATACACCACATGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          >  1:1083927/1‑148 (MQ=255)
                       gCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGtct                                                                                                                        <  2:797059/148‑1 (MQ=255)
                         gtgtTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCg                                                                                                                      <  1:1313270/148‑1 (MQ=255)
                         gtgtTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCg                                                                                                                      <  1:3058781/148‑1 (MQ=255)
                          tgttgtTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCgg                                                                                                                     <  1:2169233/148‑1 (MQ=255)
                           gttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGt                                                                                                                    >  1:3258900/1‑148 (MQ=255)
                            ttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGtt                                                                                                                   >  2:2639795/1‑148 (MQ=255)
                              gttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTcc                                                                                                                 <  1:1237156/148‑1 (MQ=255)
                              gttATTGCTGGTGCTGGTGGAGGTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTcc                                                                                                                 <  2:851896/148‑1 (MQ=255)
                               ttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTccc                                                                                                                <  1:3093684/148‑1 (MQ=255)
                                tATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                               >  1:26530/1‑148 (MQ=255)
                                tATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                               >  1:1205475/1‑148 (MQ=255)
                                       ggtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATa                                                                                                        >  1:1363067/1‑148 (MQ=255)
                                        gtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaa                                                                                                       >  1:46589/1‑148 (MQ=255)
                                        gtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaa                                                                                                       <  2:263330/148‑1 (MQ=255)
                                         tgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaaa                                                                                                      >  1:768017/1‑148 (MQ=255)
                                         tgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaaa                                                                                                      >  1:8679/1‑148 (MQ=255)
                                                ggAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCtt                                                                                               >  1:1150935/1‑148 (MQ=255)
                                                      gCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTg                                                                                         <  1:2639795/148‑1 (MQ=255)
                                                        gACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTa                                                                                       <  2:3116887/148‑1 (MQ=255)
                                                             tttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgca                                                                                  <  2:2343290/148‑1 (MQ=255)
                                                             tttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgca                                                                                  >  2:1299293/1‑148 (MQ=255)
                                                             tttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgca                                                                                  >  2:1299166/1‑148 (MQ=255)
                                                                    aGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGt                                                                           <  2:2101522/148‑1 (MQ=255)
                                                                    aGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGt                                                                           <  2:2101498/148‑1 (MQ=255)
                                                                     gACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTg                                                                          <  1:992144/148‑1 (MQ=255)
                                                                      aCGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgc                                                                         >  2:2449774/1‑148 (MQ=255)
                                                                      aCGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgc                                                                         >  2:2449961/1‑148 (MQ=255)
                                                                      aCGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgc                                                                         >  1:3478807/1‑148 (MQ=255)
                                                                           gccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGag                                                                    <  1:1890664/148‑1 (MQ=255)
                                                                            ccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagc                                                                   >  1:1375978/1‑148 (MQ=255)
                                                                                 ggTGACAGCTTAAGAAATTTAATGCTTAAAACAGCGCCTTAAGCCTTTCCAGCACTTGCTTGGTGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATATACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCgg                                                              <  2:2633592/148‑1 (MQ=255)
                                                                                  gTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                             <  2:1682159/148‑1 (MQ=255)
                                                                                  gTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                             <  1:1774502/148‑1 (MQ=255)
                                                                                    gACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGAtt                                                           <  2:3271400/148‑1 (MQ=255)
                                                                                     aCAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                          <  1:1689399/148‑1 (MQ=255)
                                                                                     aCAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                          <  2:3096979/148‑1 (MQ=255)
                                                                                     aCAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                          <  1:46949/148‑1 (MQ=255)
                                                                                       aGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGc                                                        <  2:3306024/148‑1 (MQ=255)
                                                                                        gCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCa                                                       <  2:1069928/148‑1 (MQ=255)
                                                                                        gCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCa                                                       <  1:2805652/148‑1 (MQ=255)
                                                                                           tCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGANCGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaa                                                    <  1:1783842/148‑1 (MQ=255)
                                                                                           tCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaa                                                    <  1:1156878/148‑1 (MQ=255)
                                                                                           tCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaa                                                    <  1:2996305/148‑1 (MQ=255)
                                                                                            ccGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaaa                                                   <  2:532662/148‑1 (MQ=255)
                                                                                              gACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATa                                                 >  2:1521046/1‑148 (MQ=255)
                                                                                                cAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                               <  2:509099/148‑1 (MQ=255)
                                                                                                cAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                               >  1:1809862/1‑148 (MQ=255)
                                                                                                cAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                               <  1:1917728/148‑1 (MQ=255)
                                                                                                 aGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCa                                              <  2:1493025/148‑1 (MQ=255)
                                                                                                  gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                             >  2:1885242/1‑148 (MQ=255)
                                                                                                  gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                             <  2:209117/148‑1 (MQ=255)
                                                                                                  gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                             <  1:561126/148‑1 (MQ=255)
                                                                                                  gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                             >  2:765770/1‑148 (MQ=255)
                                                                                                   ttAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATc                                            <  2:239954/148‑1 (MQ=255)
                                                                                                       tGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                        >  2:1872618/1‑148 (MQ=255)
                                                                                                       tGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                        >  2:1808191/1‑148 (MQ=255)
                                                                                                        gCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCAt                                       >  2:1833636/1‑148 (MQ=255)
                                                                                                        gCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCAt                                       >  2:2829240/1‑148 (MQ=255)
                                                                                                         cTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATc                                      <  2:657427/148‑1 (MQ=255)
                                                                                                          ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                     >  1:3042597/1‑148 (MQ=255)
                                                                                                          ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                     >  2:551547/1‑148 (MQ=255)
                                                                                                          ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                     >  1:1295924/1‑148 (MQ=255)
                                                                                                            aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCttt                                   <  2:2825157/148‑1 (MQ=255)
                                                                                                            aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCttt                                   <  2:1072145/148‑1 (MQ=255)
                                                                                                              aaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                 >  1:2486655/1‑148 (MQ=255)
                                                                                                              aaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                 >  1:2012325/1‑148 (MQ=255)
                                                                                                                  gcgcCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATc                             >  1:349908/1‑148 (MQ=255)
                                                                                                                   cgcCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCa                            >  2:2960039/1‑148 (MQ=255)
                                                                                                                       ttAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCgg                        <  1:1238399/148‑1 (MQ=255)
                                                                                                                         aaGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCg                      <  2:2360199/148‑1 (MQ=255)
                                                                                                                           gCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGaa                    <  1:1126565/148‑1 (MQ=255)
                                                                                                                            ccTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAc                   <  2:2493643/148‑1 (MQ=255)
                                                                                                                                 ccAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGaa              <  1:2896702/148‑1 (MQ=255)
                                                                                                                                  cAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAAt             <  1:2154912/148‑1 (MQ=255)
                                                                                                                                  cAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAAt             >  1:753216/1‑148 (MQ=255)
                                                                                                                                     cACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGt          >  1:2858483/1‑148 (MQ=255)
                                                                                                                                            aTGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGc   <  2:1655099/148‑1 (MQ=255)
                                                                                                                                             tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGcc  >  2:3424876/1‑148 (MQ=255)
                                                                                                                                             tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGcc  >  2:1655752/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCC  >  CP009273/618904‑619192

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
AGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTTA  >  CP009273/618902‑619197
                                                                                                                                                  |                                                                                                                                                     
AGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCT                                                                                                                                                    <  SRR19973065.727139/150‑1 (MQ=60)
   GGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTA                                                                                                                                                 >  SRR19973065.1666192/1‑150 (MQ=60)
         GGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATC                                                                                                                                           >  SRR19973065.2638072/1‑150 (MQ=60)
         GGTTTGTTGATTATCGGCGTGGTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATC                                                                                                                                           >  SRR19973065.2529296/1‑150 (MQ=60)
          GTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATTCTTAAAAAAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCC                                                                                                                                          >  SRR19973065.323993/1‑150 (MQ=60)
          GTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCC                                                                                                                                          >  SRR19973065.193689/1‑150 (MQ=60)
           TTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCA                                                                                                                                         <  SRR19973065.754432/150‑1 (MQ=60)
            TTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAG                                                                                                                                        >  SRR19973065.1417084/1‑150 (MQ=60)
             TGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGA                                                                                                                                       >  SRR19973065.536140/1‑150 (MQ=60)
              GTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGAC                                                                                                                                      <  SRR19973065.3465387/150‑1 (MQ=60)
                     ATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACACCGCCTTAATCCTATCCACCACTTGCAAGGAGCTGTAGTAATCCAGACGGAACGT                                                                                                                               >  SRR19973065.1091163/1‑150 (MQ=60)
                     ATCGGCGTGGTGTTATTTCTGGTGCTGGTGGAGTTGCGACAATTTCGCCAGACGCCGCCGCACGTGACAGCGTCCGACACTTAATGCTTAAAAAAGCGCCATAAGCCAATACACCACATGCATGGAGCTGTAGTAATCCAGACGGAACGT                                                                                                                               >  SRR19973065.1091079/1‑150 (MQ=60)
                      TCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTC                                                                                                                              <  SRR19973065.1834789/150‑1 (MQ=60)
                       CGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCT                                                                                                                             <  SRR19973065.1823704/150‑1 (MQ=60)
                       CGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCT                                                                                                                             <  SRR19973065.1769295/150‑1 (MQ=60)
                       CGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCT                                                                                                                             <  SRR19973065.1769251/150‑1 (MQ=60)
                       CGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGAGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCT                                                                                                                             <  SRR19973065.2921568/150‑1 (MQ=60)
                        GGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTC                                                                                                                            >  SRR19973065.3135882/1‑150 (MQ=60)
                           GTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGT                                                                                                                         <  SRR19973065.1321953/150‑1 (MQ=60)
                           GTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGT                                                                                                                         <  SRR19973065.3077477/150‑1 (MQ=60)
                           GTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGT                                                                                                                         >  SRR19973065.3278627/1‑150 (MQ=60)
                            TGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTT                                                                                                                        <  SRR19973065.2183163/150‑1 (MQ=60)
                                GTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCA                                                                                                                    >  SRR19973065.1213227/1‑150 (MQ=60)
                                GTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCA                                                                                                                    >  SRR19973065.26881/1‑150 (MQ=60)
                                GTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCA                                                                                                                    <  SRR19973065.1245272/150‑1 (MQ=60)
                                 TTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAG                                                                                                                   <  SRR19973065.3112573/150‑1 (MQ=60)
                                       CTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATA                                                                                                             >  SRR19973065.1372087/1‑150 (MQ=60)
                                        TGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAA                                                                                                            >  SRR19973065.47128/1‑150 (MQ=60)
                                         GGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAA                                                                                                           >  SRR19973065.8817/1‑150 (MQ=60)
                                         GGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAA                                                                                                           >  SRR19973065.773122/1‑150 (MQ=60)
                                                GTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTT                                                                                                    <  SRR19973065.851155/150‑1 (MQ=60)
                                                GTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTT                                                                                                    >  SRR19973065.1158347/1‑150 (MQ=60)
                                                 TNGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTG                                                                                                   >  SRR19973065.486272/1‑150 (MQ=60)
                                                  GGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGT                                                                                                  <  SRR19973065.1355068/150‑1 (MQ=60)
                                                  GGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGT                                                                                                  >  SRR19973065.2641182/1‑150 (MQ=60)
                                                       TGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGT                                                                                             >  SRR19973065.3453969/1‑150 (MQ=60)
                                                        GCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTA                                                                                            <  SRR19973065.2656094/150‑1 (MQ=60)
                                                             ATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCA                                                                                       >  SRR19973065.1181153/1‑150 (MQ=60)
                                                                    CCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGT                                                                                >  SRR19973065.1893723/1‑150 (MQ=60)
                                                                      AGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGC                                                                              >  SRR19973065.3499974/1‑150 (MQ=60)
                                                                      AGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTGTTGTACTGCAGGCAGGTGC                                                                              >  SRR19973065.3500700/1‑150 (MQ=60)
                                                                       GACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCG                                                                             <  SRR19973065.998734/150‑1 (MQ=60)
                                                                            CGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGC                                                                        >  SRR19973065.1385053/1‑150 (MQ=60)
                                                                             GCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCA                                                                       <  SRR19973065.1903094/150‑1 (MQ=60)
                                                                                    GTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATT                                                                <  SRR19973065.1786110/150‑1 (MQ=60)
                                                                                       ACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGC                                                             <  SRR19973065.1700679/150‑1 (MQ=60)
                                                                                       ACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGC                                                             <  SRR19973065.47490/150‑1 (MQ=60)
                                                                                          GCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATA                                                          <  SRR19973065.2823019/150‑1 (MQ=60)
                                                                                             TCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGANCGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAAT                                                       <  SRR19973065.1795506/150‑1 (MQ=60)
                                                                                             TCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGANCGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAAT                                                       <  SRR19973065.1795569/150‑1 (MQ=60)
                                                                                             TCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAAT                                                       <  SRR19973065.3014658/150‑1 (MQ=60)
                                                                                             TCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAAT                                                       <  SRR19973065.1164303/150‑1 (MQ=60)
                                                                                                GACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGC                                                    >  SRR19973065.1821771/1‑150 (MQ=60)
                                                                                                  CAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAT                                                  <  SRR19973065.1930307/150‑1 (MQ=60)
                                                                                                  CAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAT                                                  <  SRR19973065.1930362/150‑1 (MQ=60)
                                                                                                    GTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCG                                                <  SRR19973065.564698/150‑1 (MQ=60)
                                                                                                     TTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGG                                               >  SRR19973065.2524925/1‑150 (MQ=60)
                                                                                                          GCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCT                                          >  SRR19973065.1304462/1‑150 (MQ=60)
                                                                                                          GCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCT                                          >  SRR19973065.3061218/1‑150 (MQ=60)
                                                                                                              AAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCT                                      >  SRR19973065.2025302/1‑150 (MQ=60)
                                                                                                              AAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCT                                      >  SRR19973065.2502299/1‑150 (MQ=60)
                                                                                                                AACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGA                                    <  SRR19973065.2770280/150‑1 (MQ=60)
                                                                                                                  CAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATC                                  >  SRR19973065.3031024/1‑150 (MQ=60)
                                                                                                                  CAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATC                                  >  SRR19973065.352357/1‑150 (MQ=60)
                                                                                                                  CAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATC                                  >  SRR19973065.2132626/1‑150 (MQ=60)
                                                                                                                         TTAAGCCTATCCAGCACGTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCG                           <  SRR19973065.1246526/150‑1 (MQ=60)
                                                                                                                             GCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACA                       <  SRR19973065.1288889/150‑1 (MQ=60)
                                                                                                                             GCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACA                       <  SRR19973065.1134075/150‑1 (MQ=60)
                                                                                                                             GCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACA                       <  SRR19973065.1133894/150‑1 (MQ=60)
                                                                                                                                  TCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAAT                  >  SRR19973065.758241/1‑150 (MQ=60)
                                                                                                                                   CCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATA                 <  SRR19973065.2914442/150‑1 (MQ=60)
                                                                                                                                   CCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATA                 <  SRR19973065.2914386/150‑1 (MQ=60)
                                                                                                                                    CAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAG                <  SRR19973065.2168771/150‑1 (MQ=60)
                                                                                                                                     AGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGT               >  SRR19973065.2876041/1‑150 (MQ=60)
                                                                                                                                               TGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAT     >  SRR19973065.2993642/1‑150 (MQ=60)
                                                                                                                                               TGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAT     >  SRR19973065.2993664/1‑150 (MQ=60)
                                                                                                                                                 GAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTT   <  SRR19973065.658857/150‑1 (MQ=60)
                                                                                                                                                 GAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTT   <  SRR19973065.809979/150‑1 (MQ=60)
                                                                                                                                                  AGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTTA  >  SRR19973065.1492526/1‑150 (MQ=60)
                                                                                                                                                  AGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTTA  >  SRR19973065.3212214/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                     
AGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTTA  >  CP009273/618902‑619197

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: