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 3,535,370 C→A 87.8% S535R (AGC→AGA)  feoB → fused ferrous iron transporter, protein B: GTP‑binding protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,535,3700CA87.8% 251.0 / 25.6 128S535R (AGC→AGA) feoBfused ferrous iron transporter, protein B: GTP‑binding protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (10/6);  new base A (67/37);  total (77/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

TTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTA  >  CP009273/3535230‑3535517
                                                                                                                                            |                                                                                                                                                   
ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGaaaa                                                                                                                                              >  1:2385706/1‑148 (MQ=255)
 tGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGACAAt                                                                                                                                             >  2:620401/1‑148 (MQ=255)
    cGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAAtcgt                                                                                                                                          <  2:1584920/148‑1 (MQ=255)
    cGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAAtcgt                                                                                                                                          <  2:2287614/148‑1 (MQ=255)
      gTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAAtcgtcg                                                                                                                                        <  2:353451/148‑1 (MQ=255)
      gTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAAtcgtcg                                                                                                                                        >  1:2274811/1‑148 (MQ=255)
       tCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGa                                                                                                                                       <  2:236594/148‑1 (MQ=255)
          aTCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATac                                                                                                                                    >  2:1730887/1‑147 (MQ=255)
          aTCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATaa                                                                                                                                    <  2:803311/148‑1 (MQ=255)
          aTCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATaa                                                                                                                                    <  2:2480220/148‑1 (MQ=255)
            cATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACa                                                                                                                                  >  1:746066/1‑148 (MQ=255)
               gTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCa                                                                                                                               <  1:136465/148‑1 (MQ=255)
               gTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCa                                                                                                                               <  1:2789881/148‑1 (MQ=255)
                   cacaCGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGa                                                                                                                           >  1:727412/1‑148 (MQ=255)
                        tttAAAAGCCTGATTATCCCGACCTGGCAGCGTCTGAACGGCTTGGTTGTGGGTGCTGGTAATGTGATCATCTTGGTCAGAATTTTCCTGAGCGCTTTCAACCGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCgg                                                                                                                      <  2:1051806/147‑1 (MQ=255)
                        gTTAAAAGCCTGATTATCCAGACCGGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGGTGGTAAAGTGATAATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCgg                                                                                                                      <  2:489595/148‑1 (MQ=255)
                         ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGgc                                                                                                                     >  2:813376/1‑148 (MQ=255)
                              aGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTgg                                                                                                                >  2:1230968/1‑148 (MQ=255)
                              aGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTgg                                                                                                                >  2:2258674/1‑148 (MQ=255)
                                ccTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg                                                                                                              >  1:512494/1‑148 (MQ=255)
                                ccTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg                                                                                                              >  1:512458/1‑148 (MQ=255)
                                ccTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg                                                                                                              >  1:2901989/1‑148 (MQ=255)
                                  tGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtc                                                                                                            >  2:1841986/1‑148 (MQ=255)
                                  tGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtc                                                                                                            >  1:1697527/1‑148 (MQ=255)
                                  tGATTACCCACACCTGGCCGCGTCTGAACGGTTTCGTTCTGGGTGCTGTTAAAGTGATCATCATCGTCCGCATTTTCCTGAGCGCTTACAACCGCTTCTCGCTGAGTGGCAAAATCGTCGATAACAACAACGACTCGGCGCTGGcgtc                                                                                                            <  2:3061593/148‑1 (MQ=255)
                                    aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccg                                                                                                          >  2:234057/1‑148 (MQ=255)
                                    aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccg                                                                                                          >  2:234048/1‑148 (MQ=255)
                                     ttATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGACAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgt                                                                                                         >  1:2301713/1‑148 (MQ=255)
                                         ccAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGc                                                                                                     >  2:2307473/1‑148 (MQ=255)
                                             cccgggCCGCGTCTCAACGGGTTCGTTCTGCGTGCTGGTAACGTGATCATCATCGTCAGCTTTTTCCTGCGCGCTTACAACAGCTACGCGCTGAGCGGCAAAATGGTCTATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCggg                                                                                                 <  2:565961/144‑1 (MQ=255)
                                               cTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTg                                                                                               >  2:499147/1‑148 (MQ=255)
                                                tGGCAGCGTCTGACAGGCGTCGTACTGCGTGCTGCTAAAGTTAGCATCATCGTCACAATTTTCCTGAGCGATTTAAACAGCTTCTCGCTGAGAGGGAAAATCGTCCATAACATCACCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGa                                                                                              <  2:180275/148‑1 (MQ=255)
                                                tGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGa                                                                                              >  2:663409/1‑148 (MQ=255)
                                                         cttaaaGGTTTCTTTCTGCGTGCTGTTAATGTTATCAAAATCGTCAGATTTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGTGGTAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCgg                                                                                     <  1:263174/145‑1 (MQ=255)
                                                         cTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCgg                                                                                     >  2:3247399/1‑148 (MQ=255)
                                                          ttaacggCTTCGTTCTGGGTGCTGTTAAAGTGATCATCATCGTCAGCTTTTTCCTGCGCGCTTACAACCGCTTCTCGCTGAGAGGGAAAATCGTCGATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGt                                                                                    <  1:2965320/143‑1 (MQ=255)
                                                           gAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTc                                                                                   >  2:2083764/1‑148 (MQ=255)
                                                            aaaGGCTTCGTTCTGCGTGCTGGTAAATTGATCATCATCGTCAGCATTTTCCTGAGCGCTTACAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCt                                                                                  <  2:2676895/148‑1 (MQ=255)
                                                            aaaGGCTTCGTACTGCGTGCTGGTAAAGTGATCATCATCGGCAGCATTTTCCTGAGCGCTTCCAACCGCTTCTCGCTGAGAGGGAAAATCGTCGATAACACCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCt                                                                                  <  2:2927102/148‑1 (MQ=255)
                                                             aaGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCtt                                                                                 >  1:233320/1‑148 (MQ=255)
                                                             aaGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCtt                                                                                 >  1:233312/1‑148 (MQ=255)
                                                               ggCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCa                                                                               >  2:1885871/1‑148 (MQ=255)
                                                                  ttCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGc                                                                            >  1:920508/1‑148 (MQ=255)
                                                                  ttCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGc                                                                            >  2:2318699/1‑148 (MQ=255)
                                                                  ttCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGc                                                                            >  2:2318850/1‑148 (MQ=255)
                                                                     gTTCTGCTTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCGGGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCaa                                                                         <  2:977989/148‑1 (MQ=255)
                                                                      ttCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAAt                                                                        >  2:1225399/1‑148 (MQ=255)
                                                                        cTGGGTGCTGTTAAAGTGATCATCATCTTCAGATTTTTCCTGAGGGCTTCCAACAGCTTCTCGCTGAGTGGTAAAATCGTCGATAACAACAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTg                                                                      <  1:1072919/148‑1 (MQ=255)
                                                                          gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGc                                                                    >  1:402731/1‑148 (MQ=255)
                                                                          gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGc                                                                    >  1:402772/1‑148 (MQ=255)
                                                                                   tAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGaa                                                           >  2:215892/1‑148 (MQ=255)
                                                                                    aaaGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAg                                                          >  1:177173/1‑148 (MQ=255)
                                                                                    aaaGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAg                                                          >  1:177176/1‑148 (MQ=255)
                                                                                    aaaGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAg                                                          <  2:3248446/148‑1 (MQ=255)
                                                                                    aaaGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAg                                                          >  2:101302/1‑148 (MQ=255)
                                                                                      aGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGAt                                                        >  1:555436/1‑148 (MQ=255)
                                                                                       gTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATa                                                       >  1:2392425/1‑148 (MQ=255)
                                                                                        ttatcatcatcGTCAGCATTTTCCTGAGCGATTTCAACAGCTTCTCGCTGAGAGGTAAATTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATaa                                                      <  2:1546269/146‑1 (MQ=255)
                                                                                        tGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATaa                                                      >  2:3185704/1‑148 (MQ=255)
                                                                                        tGATCATCAACGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGTGGGAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATaa                                                      <  1:2737514/148‑1 (MQ=255)
                                                                                           tcatcatcGTCAGCTTTTTCGTGAGCGCTTTCAACTGCTTCTCGCTTTGAGGCAAAATCGTCTATAACAGCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTg                                                   <  1:2106733/148‑1 (MQ=255)
                                                                                           tcatcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTg                                                   >  1:749755/1‑148 (MQ=255)
                                                                                            catcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCCGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATCACTgg                                                  >  2:948562/1‑148 (MQ=255)
                                                                                            catcatcGTCAGCATTGTCCTGCCCGCTTACACCCGCTTCGCGCTGCGCGGGAAAATCGTGTATAACATCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTgg                                                  <  1:471449/148‑1 (MQ=255)
                                                                                             atcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggc                                                 >  1:2051458/1‑148 (MQ=255)
                                                                                                 tcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggcaggc                                             >  2:2641989/1‑148 (MQ=255)
                                                                                                  cGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggcaggca                                            >  2:812977/1‑148 (MQ=255)
                                                                                                   ttcatatttttCCTGAGCGCTTTCAACCGCTTCTCGGTGAGAGGCAAAATCGTCGATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCaa                                           <  1:2522434/141‑1 (MQ=255)
                                                                                                   gTCAGCTTTTTCCTGAGCGCTTTCAACAGGTTCTCGCTGAGGGGTAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCaa                                           <  1:2175078/148‑1 (MQ=255)
                                                                                                    tCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAAc                                          >  2:1799326/1‑148 (MQ=255)
                                                                                                    tCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAAc                                          >  2:1799273/1‑148 (MQ=255)
                                                                                                    tCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAAc                                          >  2:126461/1‑148 (MQ=255)
                                                                                                    tCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAAc                                          >  2:2018145/1‑148 (MQ=255)
                                                                                                       gctttttCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGt                                       <  2:2808954/145‑1 (MQ=255)
                                                                                                       gCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGt                                       >  1:2982997/1‑148 (MQ=255)
                                                                                                        cagtttCCTGAGCGCTTTCAACAGCTTCTCGCTGAGTGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGtt                                      <  2:1865961/145‑1 (MQ=255)
                                                                                                          ttttCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTgg                                    >  1:2869540/1‑148 (MQ=255)
                                                                                                           tttCCTGAGCGCTTTCAACAGGTTCTCGCTGAGAGGGAAAATCGTCGATAACATCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                                   <  1:2141764/148‑1 (MQ=255)
                                                                                                           gttCCTGAGCGCTTACAACCGCTTCTCGCTGTGTGGTAAAATCGTCTATAACCGCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                                   <  1:1870875/147‑1 (MQ=255)
                                                                                                               cTGAGCGCTTTCCACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACACCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGt                               <  2:3236347/148‑1 (MQ=255)
                                                                                                               cTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGt                               >  1:2220293/1‑148 (MQ=255)
                                                                                                                 gAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                             >  2:1403486/1‑148 (MQ=255)
                                                                                                                 gAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                             >  2:1403515/1‑148 (MQ=255)
                                                                                                                    cgcTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                          >  1:2765330/1‑148 (MQ=255)
                                                                                                                     gcTTACAACGGTTTCGCGGTGACAGTGAAAATCGTTCAAAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAg                         <  1:1422747/148‑1 (MQ=255)
                                                                                                                      cTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                        >  2:2261190/1‑148 (MQ=255)
                                                                                                                      cTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                        >  2:153072/1‑148 (MQ=255)
                                                                                                                      cTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                        >  2:2262653/1‑148 (MQ=255)
                                                                                                                      cTTACAACCGCTTCGCGCTGAGAGGTAAAATCGTCGATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                        <  1:1655071/148‑1 (MQ=255)
                                                                                                                           aaCAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCa                   >  2:943086/1‑148 (MQ=255)
                                                                                                                           aaCAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCa                   >  2:2629636/1‑148 (MQ=255)
                                                                                                                             ccgCTTCTCGCTGAGCGGGAAAATCGTCGATAACAACAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATg                 <  2:14327/146‑1 (MQ=255)
                                                                                                                             cAGCTTCTCGGTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATg                 <  1:1658298/148‑1 (MQ=255)
                                                                                                                             cAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATg                 >  2:1580576/1‑148 (MQ=255)
                                                                                                                                   ctcGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa           <  2:813132/148‑1 (MQ=255)
                                                                                                                                       cTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaag       >  2:270431/1‑148 (MQ=255)
                                                                                                                                        tGAGAGGGAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGt      <  2:1321799/148‑1 (MQ=255)
                                                                                                                                         gagaGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCCTGAAGATAACTGGCAGGCAACGGTTGGCCGGTGTACAGGTGCCATGGCGAAAGACGTg     >  2:1547689/1‑148 (MQ=255)
                                                                                                                                         gagaGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTg     >  2:1267820/1‑148 (MQ=255)
                                                                                                                                            cGGTAAAATCGTCGAAAACAACCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTa  <  1:1536072/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTA  >  CP009273/3535230‑3535517

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ 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
TCATGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGT  >  CP009273/3535222‑3535519
                                                                                                                                                    |                                                                                                                                                     
TCATGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAG                                                                                                                                                      <  SRR19973070.2138526/150‑1 (MQ=60)
TCATGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAG                                                                                                                                                      <  SRR19973070.656896/150‑1 (MQ=60)
      AGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAA                                                                                                                                                >  SRR19973070.2401460/1‑150 (MQ=60)
           CCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTC                                                                                                                                           <  SRR19973070.1783642/150‑1 (MQ=60)
            CGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCG                                                                                                                                          >  SRR19973070.2289924/1‑150 (MQ=60)
                  ATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACA                                                                                                                                    >  SRR19973070.751336/1‑150 (MQ=60)
                       GTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAAC                                                                                                                               <  SRR19973070.137654/150‑1 (MQ=60)
                       GTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAAC                                                                                                                               <  SRR19973070.2731103/150‑1 (MQ=60)
                       GTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAAC                                                                                                                               <  SRR19973070.2807856/150‑1 (MQ=60)
                         ACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGA                                                                                                                             >  SRR19973070.732571/1‑150 (MQ=60)
                         ACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGA                                                                                                                             >  SRR19973070.2343397/1‑150 (MQ=60)
                                      AGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                                                                                                                >  SRR19973070.515977/1‑150 (MQ=60)
                                      AGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                                                                                                                >  SRR19973070.516013/1‑150 (MQ=60)
                                      AGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCG                                                                                                                >  SRR19973070.2920760/1‑150 (MQ=60)
                                        CCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTC                                                                                                              >  SRR19973070.1709214/1‑150 (MQ=60)
                                           GATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGACAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGT                                                                                                           >  SRR19973070.2317023/1‑150 (MQ=60)
                                             TTACCCAGACCTGGCAGCGTCTGAACGGCTTCGTTGGGGGTGCTGTTAATGTTAACATCATCTACAGATTTTTCCTGAGCGCTTAAAACAGCTTCGCGCTGAGTGGCAAAATCGTCTATAACAACAACGACTCGGCGCTGGCGTCCGTCA                                                                                                         <  SRR19973070.1900886/150‑1 (MQ=39)
                                                                 CTTAAAGGTTTCTTTCTGCGTGCTGTTAATGTTATCAAAATCGTCAGATTTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGTGGTAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTC                                                                                     <  SRR19973070.264924/150‑1 (MQ=60)
                                                                  ttaacGGCTTCGTTCTGGGTGCTGTTAAAGTGATCATCATCGTCAGCTTTTTCCTGCGCGCTTACAACCGCTTCTCGCTGAGAGGGAAAATCGTCGATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCT                                                                                    <  SRR19973070.2984384/145‑1 (MQ=60)
                                                                   taacgttttggttctgcgtgcttttaatgttataATCATCTTCAGAATTTTCCTGCGCGTTTACAACAGCTTCTCGCTGATTGTTAAAATCTTCTATAACAACCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTT                                                                                   <  SRR19973070.3167918/116‑1 (MQ=48)
                                                                   gatcggcttcgttctgggtggtggtaatgttataatcctctacagtttttgCCTGAGCGGTTTCAACCGCTTCGCGCTGAGAGGTAAAATCGTCTATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTT                                                                                   <  SRR19973070.2272244/99‑1 (MQ=60)
                                                                   GAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTT                                                                                   >  SRR19973070.234960/1‑150 (MQ=60)
                                                                   GAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTT                                                                                   >  SRR19973070.234952/1‑150 (MQ=60)
                                                                        GCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGC                                                                              >  SRR19973070.927069/1‑150 (MQ=60)
                                                                                CTGGGTGCTGTTAAAGTGATCATCATCTTCAGATTTTTCCTGAGGGCTTCCAACAGCTTCTCGCTGAGTGGTAAAATCGTCGATAACAACAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGC                                                                      <  SRR19973070.1080407/150‑1 (MQ=60)
                                                                                CTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGC                                                                      >  SRR19973070.405690/1‑150 (MQ=60)
                                                                                CTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGC                                                                      >  SRR19973070.405731/1‑150 (MQ=60)
                                                                                          GTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAG                                                            >  SRR19973070.178581/1‑150 (MQ=60)
                                                                                          GTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAG                                                            >  SRR19973070.178578/1‑150 (MQ=60)
                                                                                            AAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGAT                                                          >  SRR19973070.559101/1‑150 (MQ=60)
                                                                                             AAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATA                                                         >  SRR19973070.2407750/1‑150 (MQ=60)
                                                                                             AAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATA                                                         >  SRR19973070.2408220/1‑150 (MQ=60)
                                                                                               GTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAAC                                                       >  SRR19973070.202299/1‑150 (MQ=60)
                                                                                                TGATCATCAACGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGTGGGAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACT                                                      <  SRR19973070.2755138/150‑1 (MQ=60)
                                                                                                 GATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTG                                                     >  SRR19973070.755049/1‑150 (MQ=60)
                                                                                                   TCATCATCGTCAGCTTTTTCGTGAGCGCTTTCAACTGCTTCTCGCTTTGAGGCAAAATCGTCTATAACAGCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGC                                                   <  SRR19973070.2120861/150‑1 (MQ=60)
                                                                                                   TCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGC                                                   >  SRR19973070.2065332/1‑150 (MQ=60)
                                                                                                   TCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGC                                                   >  SRR19973070.389401/1‑150 (MQ=60)
                                                                                                    CATCATCGTCTGCATTTTCCTGCGGGGTTACAACCGCTTCGCGGTGAGCGGTAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCA                                                  <  SRR19973070.474862/150‑1 (MQ=60)
                                                                                                    CATCATCGTCAGCATTGTCCTGCCCGCTTACACCCGCTTCGCGCTGCGCGGGAAAATCGTGTATAACATCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCA                                                  <  SRR19973070.474787/150‑1 (MQ=60)
                                                                                                           ttcatatTTTTCCTGAGCGCTTTCAACCGCTTCTCGGTGAGAGGCAAAATCGTCGATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACG                                           <  SRR19973070.2538889/143‑1 (MQ=60)
                                                                                                           GTCAGCTTTTTCCTGAGCGCTTTCAACAGGTTCTCGCTGAGGGGTAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACG                                           <  SRR19973070.2189489/150‑1 (MQ=60)
                                                                                                             CAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGT                                         >  SRR19973070.3002135/1‑150 (MQ=60)
                                                                                                                CATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGG                                      >  SRR19973070.2888127/1‑150 (MQ=60)
                                                                                                                 ATTTTCCTTAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGC                                     >  SRR19973070.2176124/1‑150 (MQ=60)
                                                                                                                 ATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGC                                     >  SRR19973070.636570/1‑150 (MQ=60)
                                                                                                                   tttcgtgagggatagcaaccGCTTCTCGCTGAGAGGTAAATTCGTCGATAACAACCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCT                                   <  SRR19973070.2050181/130‑1 (MQ=60)
                                                                                                                   TTTCCTGCGCGCTTTCAACAGCTTCTCGCTTAGAGGTAAAATCGTCTATAACAACAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCT                                   <  SRR19973070.1591748/150‑1 (MQ=60)
                                                                                                                   TTTCCTGAGCGCTTTCAACAGGTTCTCGCTGAGAGGGAAAATCGTCGATAACATCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCT                                   <  SRR19973070.2156043/150‑1 (MQ=60)
                                                                                                                   GTTCCTGAGCGCTTACAACCGCTTCTCGCTGTGTGGTAAAATCGTCTATAACCGCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCT                                   <  SRR19973070.1883638/150‑1 (MQ=60)
                                                                                                                     TCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGT                                 >  SRR19973070.2234870/1‑150 (MQ=60)
                                                                                                                          AGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACA                            >  SRR19973070.2783042/1‑150 (MQ=60)
                                                                                                                              CTTACAACCGCTTCGCGCTGAGAGGTAAAATCGTCGATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTG                        <  SRR19973070.1666480/150‑1 (MQ=60)
                                                                                                                                   AACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATG                   >  SRR19973070.2812633/1‑150 (MQ=60)
                                                                                                                                    cccGCTTCTCGGTGAGTGGGAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGG                  <  SRR19973070.526783/147‑1 (MQ=60)
                                                                                                                                     CAGCTTCTCGGTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGC                 <  SRR19973070.1669721/150‑1 (MQ=60)
                                                                                                                                                TGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGG      >  SRR19973070.30334/1‑150 (MQ=60)
                                                                                                                                                  AGTGGTAAAATCGTCTATAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTA    <  SRR19973070.822220/150‑1 (MQ=60)
                                                                                                                                                    CGGTAAAATCGTCGAAAACAACCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGT  <  SRR19973070.1546859/150‑1 (MQ=60)
                                                                                                                                                    AGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGT  >  SRR19973070.1804417/1‑150 (MQ=60)
                                                                                                                                                    |                                                                                                                                                     
TCATGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGT  >  CP009273/3535222‑3535519

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