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 F5 I175 R1 328 64.7 2582366 96.2% 2484236 120.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 635,912 G→A G235D (GGC→GAC)  ahpF → alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273635,9120GA91.1% 167.5 / 6.7 56G235D (GGC→GAC) ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (23/28);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

AAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATA  >  NZ_CP009273/635777‑636047
                                                                                                                                       |                                                                                                                                       
aaaTCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACAt                                                                                                                                      <  1:888723/139‑1 (MQ=255)
aaaTCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACAt                                                                                                                                      <  2:790792/139‑1 (MQ=255)
aaaTCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACAt                                                                                                                                      <  2:748659/139‑1 (MQ=255)
     gTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTa                                                                                                                                 <  2:238971/139‑1 (MQ=255)
      ttGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTAc                                                                                                                                >  2:456384/1‑139 (MQ=255)
        gCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTAcc                                                                                                                               <  2:73677/138‑1 (MQ=255)
                    aCTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTAATCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCg                                                                                                                  >  2:342803/1‑139 (MQ=255)
                        gcgcGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACg                                                                                                              >  1:299632/1‑139 (MQ=255)
                               aaaaCGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTtggt                                                                                                       <  2:299632/139‑1 (MQ=255)
                                   cGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTc                                                                                                   <  1:879022/139‑1 (MQ=255)
                                     tGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAg                                                                                                 <  2:443669/139‑1 (MQ=255)
                                     tGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAg                                                                                                 <  1:986221/139‑1 (MQ=255)
                                        ggCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCCATTTAATCCGCACGTAAAGACATCCGTACCGGTGTGATGGGCGAAAGTTTTGGGGGGCAGATc                                                                                              >  2:915895/1‑139 (MQ=255)
                                        ggCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGGGGTCAGATc                                                                                              >  2:915840/1‑139 (MQ=255)
                                          cAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCt                                                                                            <  2:1093621/139‑1 (MQ=255)
                                            gaagaGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCg                                                                                          <  2:155704/139‑1 (MQ=255)
                                                   tGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGt                                                                                   >  1:986970/1‑139 (MQ=255)
                                                            gTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCg                                                                           >  2:1135331/1‑138 (MQ=255)
                                                             tGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaa                                                                         <  2:986970/139‑1 (MQ=255)
                                                             tGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaa                                                                         <  1:700328/139‑1 (MQ=255)
                                                               aTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCCGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGGGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       >  1:211439/1‑139 (MQ=255)
                                                               aTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       >  2:294054/1‑139 (MQ=255)
                                                               aTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       >  1:297765/1‑139 (MQ=255)
                                                               aTGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCAGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaa                                                                        >  2:1198529/1‑138 (MQ=255)
                                                                  cTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTa                                                                    <  1:342803/139‑1 (MQ=255)
                                                                     aTGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACAt                                                                 >  1:568506/1‑139 (MQ=255)
                                                                      tGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACAtt                                                                >  2:1062409/1‑139 (MQ=255)
                                                                         cGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCt                                                             <  2:568506/139‑1 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGTGCCGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATTGAAAACTACATTTCTGTa                                                          >  1:484079/1‑139 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCg                                                                                          >  1:689465/1‑107 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCg                                                                                          <  2:689465/107‑1 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTa                                                          >  1:1230871/1‑139 (MQ=255)
                                                                            aTTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGGTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTa                                                          >  2:918448/1‑139 (MQ=255)
                                                                             ttAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTAc                                                         >  2:342239/1‑139 (MQ=255)
                                                                               aaTCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTAcc                                                        <  1:862130/138‑1 (MQ=255)
                                                                               aaTCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTAcc                                                        <  1:456384/138‑1 (MQ=255)
                                                                                 tcgtcgGTTGCGGCACGGCGGGTGCAGCCGCAGCAAGTGACTCTGCACGTAAAGACATCCGCACCGGTCTGATGGACGAACGTTTTGGTGGTCAGATCCTCGAGACCGTTGATATCGAAAACTTCATTTCTGTACCGaa                                                     <  2:1219159/139‑1 (MQ=255)
                                                                                 tcgtcgGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGaa                                                     <  2:297765/139‑1 (MQ=255)
                                                                                     cgGTTCCGGTCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACt                                                 >  1:882448/1‑139 (MQ=255)
                                                                                     cgGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACt                                                 >  1:1070960/1‑139 (MQ=255)
                                                                                       gTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGa                                               <  1:342239/139‑1 (MQ=255)
                                                                                       gTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGa                                               <  1:1135331/139‑1 (MQ=255)
                                                                                               ccggcGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGa                                       <  1:116698/139‑1 (MQ=255)
                                                                                                  gcggGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGc                                    >  1:1290965/1‑139 (MQ=255)
                                                                                                         cAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAgg                             <  1:294054/139‑1 (MQ=255)
                                                                                                           gCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGgcg                           <  2:665518/139‑1 (MQ=255)
                                                                                                            cGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGgcg                           >  2:1097537/1‑138 (MQ=255)
                                                                                                                 gcaATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       >  1:1135639/1‑89 (MQ=255)
                                                                                                                 gcaATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGaaaa                                                                       <  2:1135639/89‑1 (MQ=255)
                                                                                                                 gcaATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTg                      >  2:31990/1‑138 (MQ=255)
                                                                                                                         cTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCa              <  1:1062409/138‑1 (MQ=255)
                                                                                                                            cGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACg            <  2:484079/137‑1 (MQ=255)
                                                                                                                            cGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGtt          >  2:205593/1‑139 (MQ=255)
                                                                                                                             gCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGtt          >  2:14093/1‑138 (MQ=255)
                                                                                                                             gCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTtg         <  2:291950/139‑1 (MQ=255)
                                                                                                                                    aaGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATa  <  2:105295/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
AAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATA  >  NZ_CP009273/635777‑636047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTG  >  NZ_CP009273/635763‑636055
                                                                                                                                                     |                                                                                                                                               
CATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGAC                                                                                                                                                <  SRR3722129.838809/151‑1 (MQ=60)
         GACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTAC                                                                                                                                        <  SRR3722129.995647/150‑1 (MQ=60)
             GAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTC                                                                                                                                   <  SRR3722129.901365/151‑1 (MQ=60)
                           AATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGT                                                                                                                     >  SRR3722129.303912/1‑151 (MQ=60)
                                gttgtgatattgaaggtttgcctgttcttcctgtgtctttttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGG                                                                                                                <  SRR3722129.1284476/35‑1 (MQ=60)
                                                ACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGAT                                                                                                <  SRR3722129.891553/151‑1 (MQ=60)
                                                  GTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATAC                                                                                              <  SRR3722129.1000165/151‑1 (MQ=60)
                                                      GGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTT                                                                                          >  SRR3722129.1000923/1‑151 (MQ=60)
                                                        ttttttttttttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagCGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTG                                                                                         <  SRR3722129.962646/64‑1 (MQ=60)
                                                                  gaaaaagcataaaGCTTATGACGTATTAATCGTCGGTTCCGGCCAGGCGGGGGAAGCGGCACACATTTACTCCGCACGTAAAGACATCCGTACCGGTCCGATGGGCGAACGTTTTGGTGGTCAGActgtctcttatacacatctccgagcc                                                                              >  SRR3722129.322138/14‑125 (MQ=60)
                                                                  GAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCCGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGGGGTCAGATCCTCGATACCGTTGATATCGAAAAC                                                                              >  SRR3722129.214380/1‑151 (MQ=60)
                                                                  GAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAAC                                                                              >  SRR3722129.302026/1‑151 (MQ=60)
                                                                        GCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATT                                                                        >  SRR3722129.576926/1‑151 (MQ=60)
                                                                          GTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTC                                                                      <  SRR3722129.710654/151‑1 (MQ=60)
                                                                               GCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCCGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATTGAAAACTACATTTCTGTAC                                                                 >  SRR3722129.491077/1‑151 (MQ=60)
                                                                               GCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATctgtctcttatacacatctccg                                                                 >  SRR3722129.699638/1‑129 (MQ=60)
                                                                               GCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTAC                                                                 <  SRR3722129.347749/151‑1 (MQ=60)
                                                                               GCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGGGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTAC                                                                 >  SRR3722129.1247975/1‑151 (MQ=60)
                                                                                      agagacagATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGAC                                                          <  SRR3722129.631268/143‑1 (MQ=60)
                                                                                        GTATTAATCGTCGGTTCCGGTCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTG                                                        >  SRR3722129.895022/1‑151 (MQ=60)
                                                                                        GTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTG                                                        >  SRR3722129.1085970/1‑151 (MQ=60)
                                                                                            TAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAG                                                     <  SRR3722129.462974/150‑1 (MQ=60)
                                                                                            TAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAG                                                     <  SRR3722129.874411/150‑1 (MQ=60)
                                                                                                    GGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGC                                            <  SRR3722129.347172/151‑1 (MQ=60)
                                                                                                    GGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGC                                            <  SRR3722129.1151137/151‑1 (MQ=60)
                                                                                                     GTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCT                                           >  SRR3722129.1308845/1‑151 (MQ=60)
                                                                                                            CCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGC                                    <  SRR3722129.118300/151‑1 (MQ=60)
                                                                                                                    GTGCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGctgtctcttatacacatctccgagcccacgagactaaggc                            >  SRR3722129.1151450/1‑111 (MQ=60)
                                                                                                                      GCAGCGGCAGCAATTTACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAG                          <  SRR3722129.298273/151‑1 (MQ=60)
                                                                                                                        cgtgggctcggagatgtgtataagagacaGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAG                          <  SRR3722129.455544/120‑1 (MQ=60)
                                                                                                                                      ACTCCGCACGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAAT           <  SRR3722129.1077328/150‑1 (MQ=60)
                                                                                                                                              CGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgc  >  SRR3722129.497630/1‑89 (MQ=60)
                                                                                                                                              CGTAAAGACATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTG  >  SRR3722129.879086/1‑151 (MQ=60)
                                                                                                                                                     |                                                                                                                                               
CATGACGTTGACTGAAATCGTTGCCAAAATTGATACTGGCGCGGAAAAACGTGCGGCAGAAGAGCTGAACAAGCGTGATGCTTATGACGTATTAATCGTCGGTTCCGGCCCGGCGGGTGCAGCGGCAGCAATTTACTCCGCACGTAAAGGCATCCGTACCGGTCTGATGGGCGAACGTTTTGGTGGTCAGATCCTCGATACCGTTGATATCGAAAACTACATTTCTGTACCGAAGACTGAAGGGCAGAAGCTGGCAGGCGCACTGAAAGTTCACGTTGATGAATACGACGTTG  >  NZ_CP009273/635763‑636055

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