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
A2 F1 I1 R1 404 265.1 10478538 96.8% 10143224 150.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,632,096 G→C intergenic (+658/‑172) emrB → / → RPPX_20645 multidrug resistance protein B/oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,632,0960GC99.1% 360.0 / ‑4.8 105intergenic (+658/‑172)emrB/RPPX_20645multidrug resistance protein B/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (59/45);  minor base A (0/1);  total (59/46)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 105 total aligned reads displayed.

CACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCC  >  CP009974/4632007‑4632243
                                                                                         |                                                                                                                                                   
cACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGa                                                                                        <  2:2071464/151‑1 (MQ=12)
cACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGa                                                                                        <  2:2069710/151‑1 (MQ=12)
cACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGa                                                                                        >  1:2210295/1‑151 (MQ=9)
 aCTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGAt                                                                                       >  2:148870/1‑151 (MQ=9)
 aCTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGAt                                                                                       <  1:1795435/151‑1 (MQ=12)
  cTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAATGCCTGCGTTCGGCCTCCTTAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATc                                                                                      >  1:4109753/1‑151 (MQ=255)
  cTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATc                                                                                      <  1:2088306/151‑1 (MQ=18)
  cTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATc                                                                                      <  1:2092973/151‑1 (MQ=16)
    ccGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCaa                                                                                    >  1:3160507/1‑151 (MQ=18)
       gCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAg                                                                                 <  1:5066247/151‑1 (MQ=255)
       gCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTCGCCGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGCCCGGCCTCCTGCAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAg                                                                                 >  2:4938357/1‑151 (MQ=11)
        ccTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGc                                                                                >  2:1559745/1‑151 (MQ=255)
         cTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGcc                                                                               >  1:3588148/1‑151 (MQ=255)
          ttGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCca                                                                              >  1:2813537/1‑151 (MQ=255)
            gCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGGAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCcaca                                                                            >  2:254758/1‑151 (MQ=255)
            gCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCcaca                                                                            >  2:3888103/1‑151 (MQ=255)
               cccGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGATCGGCCTCCTGAAGTCGCGAAGATCAAGATAAAGATCAAAAGCCAGATCAACAGCCacaaca                                                                         >  2:1876946/1‑151 (MQ=255)
                 cGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAc                                                                       >  1:140166/1‑151 (MQ=255)
                  ggCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCACGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAca                                                                      >  1:1048955/1‑151 (MQ=255)
                  ggCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCGTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAca                                                                      <  2:5046367/151‑1 (MQ=255)
                    cAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagc                                                                    <  1:2472645/151‑1 (MQ=255)
                     aGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagca                                                                   >  1:5100455/1‑151 (MQ=255)
                     aGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagca                                                                   >  1:525651/1‑151 (MQ=255)
                      ggACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagcag                                                                  >  2:4468581/1‑151 (MQ=255)
                       gACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCaga                                                                 <  1:3322712/151‑1 (MQ=255)
                       gACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCaga                                                                 >  1:1508102/1‑151 (MQ=255)
                          cgcgCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAggg                                                              >  1:4007662/1‑151 (MQ=255)
                           gcgcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGt                                                             >  2:4172342/1‑151 (MQ=255)
                             gcCGAAGGCCCCATCCTGGGGCGTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           >  2:4653278/1‑151 (MQ=255)
                             gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           <  2:1944261/151‑1 (MQ=255)
                             gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           >  1:3962324/1‑151 (MQ=255)
                              ccGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTgg                                                          <  1:3563638/151‑1 (MQ=255)
                               cGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGc                                                         <  1:3129743/151‑1 (MQ=255)
                                gAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCt                                                        >  1:5118537/1‑151 (MQ=255)
                                gAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCt                                                        >  1:3526191/1‑151 (MQ=255)
                                 aaGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtg                                                       <  2:4367345/151‑1 (MQ=255)
                                 aaGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtg                                                       <  2:317184/151‑1 (MQ=255)
                                  aGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtgt                                                      <  2:1807477/151‑1 (MQ=255)
                                  aGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtgt                                                      >  1:4155232/1‑151 (MQ=255)
                                  aGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtgt                                                      <  1:2675570/151‑1 (MQ=255)
                                    gCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGc                                                    >  2:1999599/1‑151 (MQ=255)
                                    gCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGc                                                    <  2:2780505/151‑1 (MQ=255)
                                    gCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGACTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGc                                                    <  2:4085522/151‑1 (MQ=255)
                                     ccccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCt                                                   <  1:3881590/151‑1 (MQ=255)
                                     ccccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCt                                                   <  1:2860121/151‑1 (MQ=255)
                                      cccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTg                                                  >  1:4549029/1‑151 (MQ=255)
                                      cccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGCGTTGGCTGTGCTg                                                  >  1:5043365/1‑151 (MQ=255)
                                       ccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGa                                                 >  2:1043490/1‑151 (MQ=255)
                                         aTCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATc                                               <  1:2471276/151‑1 (MQ=255)
                                         aTCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATc                                               <  1:1380664/151‑1 (MQ=255)
                                          tCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCt                                              >  2:4639639/1‑151 (MQ=255)
                                           ccTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCtt                                             <  1:3758634/151‑1 (MQ=255)
                                           ccTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCtt                                             <  2:3262749/151‑1 (MQ=255)
                                            cTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTc                                            <  1:3851457/151‑1 (MQ=255)
                                            cTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTc                                            >  1:3759767/1‑151 (MQ=255)
                                             gggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCg                                           <  2:1218351/150‑1 (MQ=255)
                                              ggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGt                                          <  2:2963515/151‑1 (MQ=255)
                                              ggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAAAAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGt                                          <  1:566569/151‑1 (MQ=255)
                                               gggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTc                                         <  2:1237357/151‑1 (MQ=255)
                                                ggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGCTCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCACAGGGTTGGCTGTGCTGATCTTCGTCa                                        >  2:5099876/1‑151 (MQ=255)
                                                ggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCa                                        >  2:4171906/1‑151 (MQ=255)
                                                 gCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCAt                                       <  2:1918192/151‑1 (MQ=255)
                                                  ccTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCAtt                                      >  2:3018026/1‑151 (MQ=255)
                                                   cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTc                                     >  1:4913475/1‑151 (MQ=255)
                                                   cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTc                                     >  2:1647948/1‑151 (MQ=255)
                                                     cAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCa                                   <  2:3966038/151‑1 (MQ=255)
                                                     cAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCa                                   <  1:2583878/151‑1 (MQ=255)
                                                      ggcgcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCaa                                  <  2:214580/150‑1 (MQ=255)
                                                      aGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCaa                                  <  1:3849624/151‑1 (MQ=255)
                                                      aGCGCTTGACGGGCATCCATGCCCGCCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCACGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCAGTCCaa                                  >  2:3564696/1‑151 (MQ=18)
                                                       gcgcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAAc                                 >  1:2675628/1‑151 (MQ=255)
                                                        cgcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGCGTTGGCTGTGCTGATCTTCGTCATTCCTACg                                >  2:4340517/1‑151 (MQ=255)
                                                          cTTGACGGGTATCCATGCCCGTCACCTGCCTACGCAATGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCACGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCa                              <  2:4092000/151‑1 (MQ=25)
                                                          cTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCa                              <  1:2401209/151‑1 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCac                             >  2:501157/1‑150 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             >  1:1183061/1‑151 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             <  1:3209819/151‑1 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             >  1:477439/1‑151 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             >  1:1443025/1‑151 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             >  1:5061382/1‑151 (MQ=255)
                                                            tGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgc                            >  1:1554659/1‑151 (MQ=255)
                                                             gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcg                           >  1:2456041/1‑151 (MQ=255)
                                                             gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcg                           >  1:3096590/1‑151 (MQ=255)
                                                             gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcg                           <  1:4146292/151‑1 (MQ=255)
                                                              aCGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcgc                          >  1:1125540/1‑151 (MQ=255)
                                                               cGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGcc                         >  1:231518/1‑151 (MQ=255)
                                                                gggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCa                        >  2:3907847/1‑151 (MQ=255)
                                                                gggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCa                        >  2:819589/1‑151 (MQ=255)
                                                                  gcatggaTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCAtt                      >  1:2007793/7‑151 (MQ=255)
                                                                   cATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCAttt                     >  1:1924779/1‑151 (MQ=255)
                                                                    aTCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTc                    >  2:3309056/1‑151 (MQ=255)
                                                                       gaTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGcc                 >  2:5072630/2‑151 (MQ=255)
                                                                          gCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTAt              <  2:525651/151‑1 (MQ=255)
                                                                           cccGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGCCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCACCGCCATTTCGCCTAtt             >  2:2471888/1‑151 (MQ=255)
                                                                             cGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGt           <  2:610928/151‑1 (MQ=255)
                                                                                cACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGa        <  2:217213/151‑1 (MQ=255)
                                                                                 aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       <  1:4468581/151‑1 (MQ=255)
                                                                                 aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       >  2:4685912/1‑151 (MQ=255)
                                                                                 aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       <  1:1912666/151‑1 (MQ=255)
                                                                                      ccTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTcc  <  1:231418/151‑1 (MQ=255)
                                                                                         |                                                                                                                                                   
CACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCC  >  CP009974/4632007‑4632243

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
ACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCCC  >  CP009974/4631950‑4632244
                                                                                                                                                  |                                                                                                                                                    
ACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCA                                                                                                                                                  <  SRR10918231.374754/151‑1 (MQ=0)
    AAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGC                                                                                                                                              <  SRR10918231.3270818/151‑1 (MQ=0)
              CGCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGG                                                                                                                                    <  SRR10918231.2205614/151‑1 (MQ=0)
               GCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGC                                                                                                                                   <  SRR10918231.4859097/151‑1 (MQ=0)
                 TCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCT                                                                                                                                 >  SRR10918231.3292040/1‑151 (MQ=0)
                  CCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTC                                                                                                                                >  SRR10918231.1324239/1‑151 (MQ=0)
                            GGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGC                                                                                                                      <  SRR10918231.1553793/151‑1 (MQ=0)
                                                         CACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGA                                                                                         >  SRR10918231.2220205/1‑151 (MQ=0)
                                                          ACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGAT                                                                                        <  SRR10918231.1803608/151‑1 (MQ=1)
                                                           CTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAATGCCTGCGTTCGGCCTCCTTAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATC                                                                                       >  SRR10918231.4129152/1‑151 (MQ=3)
                                                           CTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATC                                                                                       <  SRR10918231.2097669/151‑1 (MQ=3)
                                                           CTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATC                                                                                       <  SRR10918231.2102368/151‑1 (MQ=3)
                                                             CCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAA                                                                                     >  SRR10918231.3174930/1‑151 (MQ=7)
                                                                GCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAG                                                                                  <  SRR10918231.5090943/151‑1 (MQ=18)
                                                                  CTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCC                                                                                >  SRR10918231.3604519/1‑151 (MQ=22)
                                                                   TTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCA                                                                               >  SRR10918231.2826431/1‑151 (MQ=24)
                                                                    TGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCAC                                                                              <  SRR10918231.3989036/151‑1 (MQ=27)
                                                                          CGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAC                                                                        >  SRR10918231.140675/1‑151 (MQ=39)
                                                                           GGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCACGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACA                                                                       >  SRR10918231.1053671/1‑151 (MQ=40)
                                                                             CAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGC                                                                     <  SRR10918231.2483954/151‑1 (MQ=45)
                                                                              AGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCA                                                                    >  SRR10918231.5125375/1‑151 (MQ=47)
                                                                              AGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCA                                                                    >  SRR10918231.527524/1‑151 (MQ=47)
                                                                                GACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGA                                                                  >  SRR10918231.1515098/1‑151 (MQ=51)
                                                                                GACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGA                                                                  <  SRR10918231.3337777/151‑1 (MQ=51)
                                                                                   CGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGG                                                               >  SRR10918231.4026459/1‑151 (MQ=57)
                                                                                    aagcgctgggccccaggaTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCTCTGCGCTAGGCCGGCGGTCGGCCTGCTGAAGTCGCcagctgcggaggttgcagtatagatcggaagacaacacgtctgaacgccagtcggc                                                              >  SRR10918231.1707151/19‑95 (MQ=3)
                                                                                      GCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTG                                                            >  SRR10918231.3980872/1‑151 (MQ=60)
                                                                                       CCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGG                                                           <  SRR10918231.3579899/151‑1 (MQ=60)
                                                                                        CGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGC                                                          <  SRR10918231.3144033/151‑1 (MQ=60)
                                                                                         GAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCT                                                         >  SRR10918231.5143604/1‑151 (MQ=60)
                                                                                         GAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCT                                                         >  SRR10918231.3542270/1‑151 (MQ=60)
                                                                                           AGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGT                                                       >  SRR10918231.4174864/1‑151 (MQ=60)
                                                                                           AGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGT                                                       <  SRR10918231.2687888/151‑1 (MQ=60)
                                                                                              CCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCT                                                    <  SRR10918231.2873204/151‑1 (MQ=60)
                                                                                              CCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCT                                                    <  SRR10918231.3899630/151‑1 (MQ=60)
                                                                                               CCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTG                                                   >  SRR10918231.4570702/1‑151 (MQ=60)
                                                                                               CCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGCGTTGGCTGTGCTG                                                   >  SRR10918231.5067913/1‑151 (MQ=60)
                                                                                                  ATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATC                                                <  SRR10918231.2482660/151‑1 (MQ=60)
                                                                                                  ATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATC                                                <  SRR10918231.2482582/151‑1 (MQ=60)
                                                                                                  ATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATC                                                <  SRR10918231.1387212/151‑1 (MQ=60)
                                                                                                    CCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTT                                              <  SRR10918231.3775906/151‑1 (MQ=60)
                                                                                                     CTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTC                                             <  SRR10918231.3869297/151‑1 (MQ=60)
                                                                                                     CTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTC                                             >  SRR10918231.3777050/1‑151 (MQ=60)
                                                                                                       GGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAAAAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGT                                           <  SRR10918231.568637/151‑1 (MQ=60)
                                                                                                            CTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTC                                      >  SRR10918231.4937176/1‑151 (MQ=60)
                                                                                                             aagggaggTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCC                                     >  SRR10918231.217993/9‑151 (MQ=60)
                                                                                                              CAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCA                                    <  SRR10918231.2595772/151‑1 (MQ=60)
                                                                                                               AGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAA                                   <  SRR10918231.3867446/151‑1 (MQ=60)
                                                                                                                gcctttgcCGGGCATCCATGCCCGGCACCTGCCTCCGCAAGGCCTGCGTTCAGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAAC                                  <  SRR10918231.629118/143‑1 (MQ=60)
                                                                                                                GCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAAC                                  >  SRR10918231.2687948/1‑151 (MQ=60)
                                                                                                                   CTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCA                               <  SRR10918231.2412170/151‑1 (MQ=60)
                                                                                                                    TTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAG                              >  SRR10918231.1449785/1‑151 (MQ=60)
                                                                                                                    TTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAG                              >  SRR10918231.1188601/1‑151 (MQ=60)
                                                                                                                    TTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAG                              >  SRR10918231.5086037/1‑151 (MQ=60)
                                                                                                                    TTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAG                              <  SRR10918231.3224417/151‑1 (MQ=60)
                                                                                                                    TTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAG                              >  SRR10918231.479157/1‑151 (MQ=60)
                                                                                                                     TGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGC                             >  SRR10918231.1561817/1‑151 (MQ=60)
                                                                                                                      GACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCG                            >  SRR10918231.2467274/1‑151 (MQ=60)
                                                                                                                      GACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCG                            <  SRR10918231.4165874/151‑1 (MQ=60)
                                                                                                                      GACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCG                            >  SRR10918231.3110731/1‑151 (MQ=60)
                                                                                                                       ACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGC                           >  SRR10918231.1130726/1‑151 (MQ=60)
                                                                                                                        CGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCC                          >  SRR10918231.232340/1‑151 (MQ=60)
                                                                                                                           GCATGGATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATT                       >  SRR10918231.2016814/1‑151 (MQ=60)
                                                                                                                            CATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTT                      >  SRR10918231.1933431/1‑151 (MQ=60)
                                                                                                                                          ACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAT        <  SRR10918231.4489816/151‑1 (MQ=60)
                                                                                                                                          ACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAT        <  SRR10918231.1921271/151‑1 (MQ=60)
                                                                                                                                           CCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATG       <  SRR10918231.3358466/151‑1 (MQ=60)
                                                                                                                                               CCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCC   <  SRR10918231.232239/151‑1 (MQ=60)
                                                                                                                                                CTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCCC  <  SRR10918231.5218166/151‑1 (MQ=60)
                                                                                                                                                CTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCCC  >  SRR10918231.858791/1‑151 (MQ=60)
                                                                                                                                                  |                                                                                                                                                    
ACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCCC  >  CP009974/4631950‑4632244

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