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 F1 I1 R1 343 268.7 10873516 97.0% 10547310 150.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,537,108 A→G Y105Y (TAT→TAC RPPX_20270 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,537,1080AG100.0% 240.2 / NA 76Y105Y (TAT→TACRPPX_20270hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/38);  total (38/38)

GTGGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACC  >  CP009974/4536960‑4537166
                                                                                                                                                    |                                                          
gtggcgCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTa                                                          <  1:2876947/151‑1 (MQ=255)
gtggcgCTGCTCGGACGTTGAACAGGCAGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTa                                                          >  2:3575391/1‑151 (MQ=255)
 tggcgCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAg                                                         >  1:2879290/1‑151 (MQ=255)
  ggcgCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGt                                                        >  2:740448/1‑151 (MQ=255)
  ggcgCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGt                                                        >  2:4897773/1‑151 (MQ=255)
  ggcgCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGt                                                        <  1:3708871/151‑1 (MQ=255)
   gcgcTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTg                                                       <  1:3287269/151‑1 (MQ=255)
   gcgcTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTg                                                       <  1:2309106/151‑1 (MQ=255)
    cgcTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTgg                                                      >  2:1739895/1‑151 (MQ=255)
     gctgctCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGgc                                                     >  1:1151823/1‑151 (MQ=255)
      ctgctCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGgcg                                                    >  1:5237390/1‑151 (MQ=255)
      ctgctCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGgcg                                                    <  1:1573133/151‑1 (MQ=255)
       tgctCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGgcgc                                                   >  1:2457673/1‑151 (MQ=255)
        gctCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGcc                                                  <  2:1443735/151‑1 (MQ=255)
        gctCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGcc                                                  <  2:1444921/151‑1 (MQ=255)
         ctCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCa                                                 <  2:2745987/151‑1 (MQ=255)
          tCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAg                                                >  1:3141584/1‑151 (MQ=255)
          tCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAg                                                <  1:1298000/151‑1 (MQ=255)
          tCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAg                                                >  1:3105943/1‑151 (MQ=255)
           cGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAgg                                               >  2:133986/1‑151 (MQ=255)
           cGGACGTTGAACAGGCGGCAGAGGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCGCCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAgg                                               <  2:2665340/151‑1 (MQ=255)
             gACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCt                                             <  2:145369/151‑1 (MQ=255)
             gACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCt                                             >  1:1388750/1‑151 (MQ=255)
              aCGTTGAACAGGCGGCAGATGCCATACCAGCAGGGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGACCGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCtt                                            <  1:3221763/151‑1 (MQ=255)
              aCGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCtt                                            <  1:3221850/151‑1 (MQ=255)
               cGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTc                                           <  2:4826592/151‑1 (MQ=255)
                gTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCt                                          <  2:4846666/151‑1 (MQ=255)
                gTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCt                                          >  1:3241325/1‑151 (MQ=255)
                 ttGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTg                                         <  2:971897/151‑1 (MQ=255)
                 ttGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCACGCTTCTg                                         >  1:4817373/1‑151 (MQ=255)
                  tGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGc                                        >  1:2165827/1‑151 (MQ=255)
                  tGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGc                                        <  1:490112/151‑1 (MQ=255)
                  tGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGc                                        <  2:18129/151‑1 (MQ=255)
                   gAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGcc                                       <  1:3315893/151‑1 (MQ=255)
                   gAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGcc                                       >  1:1004745/1‑151 (MQ=255)
                    aaCAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTCGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCg                                      <  2:5252040/151‑1 (MQ=255)
                    aaCAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCg                                      <  1:4938015/151‑1 (MQ=255)
                     aCAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCgg                                     >  1:1159422/1‑151 (MQ=255)
                      cAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCggg                                    <  2:297056/151‑1 (MQ=255)
                       aGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCgggg                                   <  1:3356152/151‑1 (MQ=255)
                       aGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCgggg                                   >  1:2626109/1‑151 (MQ=255)
                       aGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCgggg                                   >  2:2392291/1‑151 (MQ=255)
                        ggcggcAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGt                                  <  1:4621095/151‑1 (MQ=255)
                            gcAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAgg                              >  1:1994011/1‑151 (MQ=255)
                            gcAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAgg                              >  1:578167/1‑151 (MQ=255)
                              aGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCGGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTCCTGCCGGGGGGAGGCa                            >  2:2756843/1‑151 (MQ=255)
                              aGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCa                            >  1:2628271/1‑151 (MQ=255)
                               gATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAg                           >  1:1886982/1‑151 (MQ=255)
                                aTGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGa                          >  2:1534799/1‑151 (MQ=255)
                                 tGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGAt                         >  2:4736345/1‑151 (MQ=255)
                                 tGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGAt                         >  2:290011/1‑151 (MQ=255)
                                  gCCATACCCGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGACCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATg                        <  2:5237917/151‑1 (MQ=255)
                                  gCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATg                        <  2:3286991/151‑1 (MQ=255)
                                   ccATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATgg                       <  1:1402505/151‑1 (MQ=255)
                                   ccATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGGGAGGCAGATgg                       >  2:3259880/1‑151 (MQ=255)
                                     aTACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCg                     >  1:1460668/1‑151 (MQ=25)
                                     aTACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCg                     <  2:408391/151‑1 (MQ=255)
                                      tACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGt                    <  1:868632/151‑1 (MQ=19)
                                      tACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGt                    <  2:363728/151‑1 (MQ=255)
                                      tACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGt                    <  1:3121326/151‑1 (MQ=16)
                                       aCCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTg                   <  2:3736077/151‑1 (MQ=255)
                                       aCCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTg                   >  1:3883433/1‑151 (MQ=21)
                                        gcaggagAGGGCGACGGTCAGGTTGGAGGGGCACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCAGTTCACCGGGCGCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGc                  <  2:4458608/150‑1 (MQ=11)
                                         cagcagAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCa                 <  2:503933/151‑1 (MQ=18)
                                          agcagAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAg                >  1:840767/1‑151 (MQ=18)
                                          agcagAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAg                >  1:1690204/1‑151 (MQ=18)
                                           gcagAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGt               >  2:285286/1‑151 (MQ=18)
                                           gcagAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGt               >  1:1190226/1‑151 (MQ=16)
                                             agagGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTgg             >  2:3504514/1‑151 (MQ=17)
                                              gagGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGc            <  2:1139526/151‑1 (MQ=9)
                                              gagGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGc            <  1:2883896/151‑1 (MQ=12)
                                               agGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCt           >  1:1913058/1‑151 (MQ=9)
                                               agGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCt           >  1:2351626/1‑151 (MQ=9)
                                                gggCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCCGATGGCGTGCAGTGGCtt          >  2:4838059/1‑151 (MQ=9)
                                                gggCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCtt          <  1:432094/151‑1 (MQ=11)
                                                        gtaaGGTTGGAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGccacc  <  2:1538798/148‑1 (MQ=2)
                                                                                                                                                    |                                                          
GTGGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACC  >  CP009974/4536960‑4537166

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGTGGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAGTGG  >  CP009974/4536959‑4537219
                                                                                                                                                     |                                                                                                               
GGTGGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGT                                                                                                                <  SRR10918230.1153180/151‑1 (MQ=60)
GGTGGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGT                                                                                                                >  SRR10918230.2343688/1‑151 (MQ=60)
 GTGGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTA                                                                                                               <  SRR10918230.2889470/151‑1 (MQ=60)
  TGGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAG                                                                                                              >  SRR10918230.2891818/1‑151 (MQ=60)
   GGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGT                                                                                                             <  SRR10918230.3724982/151‑1 (MQ=60)
    GCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTG                                                                                                            <  SRR10918230.3301553/151‑1 (MQ=60)
    GCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTG                                                                                                            <  SRR10918230.2318927/151‑1 (MQ=60)
      GCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGC                                                                                                          >  SRR10918230.1156780/1‑151 (MQ=60)
       CTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCG                                                                                                         <  SRR10918230.1580039/151‑1 (MQ=60)
       CTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCG                                                                                                         >  SRR10918230.5261638/1‑151 (MQ=60)
        TGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGC                                                                                                        >  SRR10918230.2468219/1‑151 (MQ=60)
           TCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAG                                                                                                     >  SRR10918230.3119477/1‑151 (MQ=60)
           TCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAG                                                                                                     <  SRR10918230.1303840/151‑1 (MQ=60)
           TCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAG                                                                                                     >  SRR10918230.3155261/1‑151 (MQ=60)
              GACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCT                                                                                                  >  SRR10918230.1395012/1‑151 (MQ=60)
               ACGTTGAACAGGCGGCAGATGCCATACCAGCAGGGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGACCGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTT                                                                                                 <  SRR10918230.3235771/151‑1 (MQ=60)
               ACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTT                                                                                                 <  SRR10918230.3235858/151‑1 (MQ=60)
                 GTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCT                                                                                               >  SRR10918230.3255409/1‑151 (MQ=60)
                  TTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCACGCTTCTG                                                                                              >  SRR10918230.4839281/1‑151 (MQ=60)
                   TGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGC                                                                                             <  SRR10918230.491815/151‑1 (MQ=60)
                   TGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGC                                                                                             >  SRR10918230.2174967/1‑151 (MQ=60)
                    GAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCC                                                                                            <  SRR10918230.3330288/151‑1 (MQ=60)
                    GAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCC                                                                                            >  SRR10918230.1008857/1‑151 (MQ=60)
                     AACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCG                                                                                           <  SRR10918230.4960459/151‑1 (MQ=59)
                      ACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGG                                                                                          >  SRR10918230.1164412/1‑151 (MQ=57)
                        AGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGG                                                                                        <  SRR10918230.3370689/151‑1 (MQ=52)
                        AGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGG                                                                                        >  SRR10918230.2637537/1‑151 (MQ=52)
                         GGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGT                                                                                       <  SRR10918230.4641998/151‑1 (MQ=50)
                             GCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGG                                                                                   >  SRR10918230.2002437/1‑151 (MQ=41)
                             GCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGG                                                                                   >  SRR10918230.580229/1‑151 (MQ=41)
                               AGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCA                                                                                 >  SRR10918230.2639715/1‑151 (MQ=37)
                                GATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAG                                                                                >  SRR10918230.1895002/1‑151 (MQ=35)
                                    CCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGG                                                                            <  SRR10918230.1408807/151‑1 (MQ=27)
                                      ATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCG                                                                          >  SRR10918230.1467184/1‑151 (MQ=22)
                                       TACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGT                                                                         <  SRR10918230.872052/151‑1 (MQ=20)
                                       TACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGT                                                                         <  SRR10918230.3134917/151‑1 (MQ=20)
                                        ACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTG                                                                        >  SRR10918230.3900402/1‑151 (MQ=20)
                                           AGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAG                                                                     >  SRR10918230.844021/1‑151 (MQ=11)
                                           AGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAG                                                                     >  SRR10918230.1697586/1‑151 (MQ=11)
                                            GCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGT                                                                    >  SRR10918230.1195370/1‑151 (MQ=9)
                                               GAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGC                                                                 <  SRR10918230.2896443/151‑1 (MQ=5)
                                                AGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCT                                                                >  SRR10918230.1921178/1‑151 (MQ=3)
                                                AGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCT                                                                >  SRR10918230.2361656/1‑151 (MQ=3)
                                                 GGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTT                                                               <  SRR10918230.433588/151‑1 (MQ=1)
                                                       CGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGCCA                                                         <  SRR10918230.5306077/151‑1 (MQ=0)
                                                        GGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGCCAC                                                        >  SRR10918230.3097296/1‑151 (MQ=0)
                                                         GTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGCCACC                                                       >  SRR10918230.3058446/1‑151 (MQ=0)
                                                             GGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTT                                                   <  SRR10918230.4878635/151‑1 (MQ=0)
                                                                   AGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGCTGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTC                                             >  SRR10918230.577152/1‑151 (MQ=0)
                                                                   AGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTC                                             >  SRR10918230.2985321/1‑151 (MQ=0)
                                                                    GGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCA                                            <  SRR10918230.1591783/151‑1 (MQ=0)
                                                                     GGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCAC                                           <  SRR10918230.484972/151‑1 (MQ=0)
                                                                        AACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCGTAGTGGCGCCAGGCTTCTGCCGGGGTGAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGG                                        >  SRR10918230.4471198/1‑151 (MQ=0)
                                                                        AACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGG                                        <  SRR10918230.1368756/151‑1 (MQ=0)
                                                                              ACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGG                                  <  SRR10918230.1207396/151‑1 (MQ=0)
                                                                               CGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGG                                 >  SRR10918230.708189/1‑151 (MQ=0)
                                                                                GGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGG                                >  SRR10918230.5312419/1‑151 (MQ=0)
                                                                                 GCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGA                               <  SRR10918230.1304675/151‑1 (MQ=0)
                                                                                 GCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGA                               <  SRR10918230.5265625/151‑1 (MQ=0)
                                                                                   CCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGT                             >  SRR10918230.2415460/1‑151 (MQ=0)
                                                                                      GAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTG                          >  SRR10918230.571654/1‑151 (MQ=0)
                                                                                        ACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTC                        >  SRR10918230.1415070/1‑151 (MQ=0)
                                                                                         CGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCC                       >  SRR10918230.3788546/1‑151 (MQ=0)
                                                                                           GCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGG                     <  SRR10918230.2228612/151‑1 (MQ=0)
                                                                                           GCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGG                     <  SRR10918230.1358751/151‑1 (MQ=0)
                                                                                            CAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGG                    >  SRR10918230.3270502/1‑151 (MQ=0)
                                                                                              ATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCC                  <  SRR10918230.2594528/151‑1 (MQ=0)
                                                                                                TGCCAGTTCTTCTTGCCCTGCCTTTCACCCGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCCCCGGGCGGGGGGAGTTTGTCCGGGCCGT                <  SRR10918230.124435/151‑1 (MQ=0)
                                                                                                 GCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTG               <  SRR10918230.4294735/151‑1 (MQ=0)
                                                                                                   CAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGT             <  SRR10918230.3427437/151‑1 (MQ=0)
                                                                                                    AGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTC            <  SRR10918230.4350265/151‑1 (MQ=0)
                                                                                                       TCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCAT         <  SRR10918230.5381461/151‑1 (MQ=0)
                                                                                                        CTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATG        >  SRR10918230.1054288/1‑151 (MQ=0)
                                                                                                         TTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGA       <  SRR10918230.2218854/151‑1 (MQ=0)
                                                                                                          TTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAA      >  SRR10918230.2739905/1‑151 (MQ=0)
                                                                                                           TTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAG     >  SRR10918230.2704470/1‑151 (MQ=0)
                                                                                                              GCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAGTGG  <  SRR10918230.4743845/151‑1 (MQ=0)
                                                                                                              GCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAGTGG  <  SRR10918230.4742110/151‑1 (MQ=0)
                                                                                                                                                     |                                                                                                               
GGTGGCGCTGCTCGGACGTTGAACAGGCGGCAGATGCCATACCAGCAGAGGGCGACGGTCAGGTTGTAGGGGAACAGCACGGCCCAGAACGGCAATTGCCAGTTCTTTTTGCCCTGCATTTCACCGGGCTCGCCGGTTCGCCAGGGGGCATAGTGGCGCCAGGCTTCTGCCGGGGTCAGGCAGATGGCGTGCAGTGGCTTGTGCCACCAGTTGGGTTCACCGGGCGGGGGGAGTTTGTCCGGGCCGTGGTCCATGAAGTGG  >  CP009974/4536959‑4537219

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