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 1,047,923 (T)5→4 coding (896/1074 nt) yccM ← 4Fe‑4S binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,047,9190T.98.0% 217.9 / ‑2.0 51coding (900/1074 nt)yccM4Fe‑4S binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (27/23);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCG  >  NZ_CP009273/1047782‑1048048
                                                                                                                                         |                                                                                                                                 
cGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGtt                                                                                                                                  >  2:475313/1‑139 (MQ=255)
               ggCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGt                                                                                                                  <  1:475313/139‑1 (MQ=255)
                           tGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAgg                                                                                                       <  2:234773/138‑1 (MQ=255)
                                       gCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGa                                                                                          >  1:7652/1‑139 (MQ=255)
                                          tACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGActc                                                                                       >  1:88161/1‑139 (MQ=255)
                                             aTTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGt                                                                                                                  >  2:795715/1‑109 (MQ=255)
                                             aTTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGt                                                                                                                  <  1:795715/109‑1 (MQ=255)
                                             aTTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTAc                                                                                    >  2:537445/1‑139 (MQ=255)
                                             aTTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTAc                                                                                    >  2:436360/1‑139 (MQ=255)
                                               tGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTAc                                                                                                   >  2:620321/1‑122 (MQ=255)
                                               tGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTAc                                                                                                   <  1:620321/122‑1 (MQ=255)
                                                gCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGc                                                                                  <  2:87779/138‑1 (MQ=255)
                                                gCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGc                                                                                  <  1:1164410/138‑1 (MQ=255)
                                                 cGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGc                                                                                  >  1:311732/1‑137 (MQ=255)
                                                 cGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGc                                                                                  >  1:878734/1‑137 (MQ=255)
                                                  gTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAg                                                                               >  1:414490/1‑139 (MQ=255)
                                                     aCCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACa                                                                                             <  2:171516/122‑1 (MQ=255)
                                                     aCCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACa                                                                                             >  1:171516/1‑122 (MQ=255)
                                                      ccAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCAt                                                                           <  2:1287984/139‑1 (MQ=255)
                                                      ccAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCAt                                                                           <  1:977757/139‑1 (MQ=255)
                                                         aTCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATa                                                                          >  1:1055844/1‑137 (MQ=255)
                                                         aTCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATa                                                                          >  1:369363/1‑137 (MQ=255)
                                                                cAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTa                                                                 >  1:1126376/1‑139 (MQ=255)
                                                                       aaTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATAGCTCTACGCAAGTCATACAGCCGGTACATTCCa                                                          <  2:797828/139‑1 (MQ=255)
                                                                        aTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTa                                                                                                    <  1:1112546/96‑1 (MQ=255)
                                                                        aTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTa                                                                                                    >  2:1112546/1‑96 (MQ=255)
                                                                        aTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCa                                                          <  2:1126376/138‑1 (MQ=255)
                                                                        aTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCa                                                          <  1:953174/138‑1 (MQ=255)
                                                                        aTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCa                                                          <  2:626295/138‑1 (MQ=255)
                                                                        aTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCAc                                                         >  2:534742/1‑139 (MQ=255)
                                                                             cAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTc                                                    >  2:489837/1‑139 (MQ=255)
                                                                                     gTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTGCCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTg                                             >  1:1014128/1‑138 (MQ=255)
                                                                                     gTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACa                                                                                             <  1:18632/90‑1 (MQ=255)
                                                                                     gTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACa                                                                                             >  2:18632/1‑90 (MQ=255)
                                                                                     gTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTg                                             >  2:1086156/1‑138 (MQ=255)
                                                                                     gTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTg                                             >  1:643842/1‑138 (MQ=255)
                                                                                             aTTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGt                                                                  <  2:1017158/109‑1 (MQ=255)
                                                                                             aTTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGt                                                                  >  1:1017158/1‑109 (MQ=255)
                                                                                             aTTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGaattaatt                                     >  2:153901/1‑138 (MQ=255)
                                                                                                aaCCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTAt                                  <  1:17463/138‑1 (MQ=255)
                                                                                                   cagccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAc                              <  1:1086156/139‑1 (MQ=255)
                                                                                                   cagccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAc                              <  1:903318/139‑1 (MQ=255)
                                                                                                    agccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAcc                             <  2:7652/139‑1 (MQ=255)
                                                                                                    agccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAcc                             <  1:534742/139‑1 (MQ=255)
                                                                                                     gccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAcc                             <  2:252521/139‑1 (MQ=255)
                                                                                                     gccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCg                            >  2:446092/1‑139 (MQ=255)
                                                                                                     gccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCg                            <  2:414490/139‑1 (MQ=255)
                                                                                                      ccagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCgg                           >  2:418766/1‑139 (MQ=255)
                                                                                                              cGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGt                    >  1:399145/1‑138 (MQ=255)
                                                                                                                               ttGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTgcg  >  1:836704/1‑139 (MQ=255)
                                                                                                                                tGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCgg                           <  1:1236207/113‑1 (MQ=255)
                                                                                                                                tGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCgg                           >  2:1236207/1‑113 (MQ=255)
                                                                                                                                         |                                                                                                                                 
CGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCG  >  NZ_CP009273/1047782‑1048048

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATAC  >  NZ_CP009273/1047776‑1048067
                                                                                                                                                   |                                                                                                                                                
AAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAG                                                                                                                                              <  SRR3722129.582187/151‑1 (MQ=60)
                    gagatgtgtataagagacagCTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGG                                                                                                                          <  SRR3722129.807230/131‑1 (MQ=60)
                    CGGCACCGGGCTTTGCCATCCTCCCGCATACATTGCGTAACCAATCACCGCAAACATCATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGG                                                                                                                          <  SRR3722129.463983/151‑1 (MQ=60)
                    CGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGG                                                                                                                          <  SRR3722129.482202/151‑1 (MQ=60)
                                  GCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGAC                                                                                                            >  SRR3722129.7752/1‑151 (MQ=60)
                                   gagacagCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACT                                                                                                           <  SRR3722129.629542/144‑1 (MQ=60)
                                     cgtgggctcggagatgtgtataagagacagCCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGAC                                                                                                            <  SRR3722129.1128077/118‑1 (MQ=60)
                                     ATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCT                                                                                                         >  SRR3722129.89306/1‑151 (MQ=60)
                                         ccttagtctcgtgggctcggagatgtgtataagagacagCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGC                                                                                                     <  SRR3722129.18852/112‑1 (MQ=60)
                                            CGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCA                                                                                                    >  SRR3722129.891264/1‑149 (MQ=60)
                                            CGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCA                                                                                                    >  SRR3722129.316188/1‑149 (MQ=60)
                                             GCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGT                                                                                                 >  SRR3722129.420533/1‑151 (MQ=60)
                                                TACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCctgtctc                                                                                              >  SRR3722129.173878/1‑144 (MQ=60)
                                                    TTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATAC                                                                                            >  SRR3722129.1070692/1‑149 (MQ=60)
                                                    TTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATAC                                                                                            >  SRR3722129.374773/1‑149 (MQ=60)
                                                     TGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAGGTCATACAG                                                                                          <  SRR3722129.1180649/150‑1 (MQ=60)
                                                           ACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTAC                                                                                   <  SRR3722129.991592/151‑1 (MQ=60)
                                                           ACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTAC                                                                                   >  SRR3722129.1142073/1‑151 (MQ=60)
                                                                        AACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCctgtctcttatacacatctccgagcccacgagact                                                                       >  SRR3722129.557816/1‑115 (MQ=60)
                                                                             AATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACC                                                                  <  SRR3722129.966645/150‑1 (MQ=60)
                                                                                CCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTGCCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGA                                                               >  SRR3722129.1028445/1‑150 (MQ=60)
                                                                                CCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGA                                                               >  SRR3722129.653442/1‑150 (MQ=60)
                                                                                        CCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGctgtctcttatacacatctc                                                      >  SRR3722129.1031518/1‑131 (MQ=60)
                                                                                                     TAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATC                                          <  SRR3722129.17670/150‑1 (MQ=60)
                                                                                                        CCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTG                                      <  SRR3722129.1101357/151‑1 (MQ=60)
                                                                                                        CCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTG                                      <  SRR3722129.916125/151‑1 (MQ=60)
                                                                                                         CAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGA                                     <  SRR3722129.542555/151‑1 (MQ=60)
                                                                                                         CAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTG                                      >  SRR3722129.404949/1‑150 (MQ=60)
                                                                                                            gtgtataagagacAGGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATG                                   <  SRR3722129.1253375/137‑1 (MQ=60)
                                                                                                                          GGCTTTTTTATTTGCCGCAGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGC                    >  SRR3722129.848713/1‑151 (MQ=60)
                                                                                                                                            AGGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATAC  <  SRR3722129.156054/151‑1 (MQ=60)
                                                                                                                                              GGTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATAC  >  SRR3722129.1175580/1‑150 (MQ=60)
                                                                                                                                                   |                                                                                                                                                
AAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATAC  >  NZ_CP009273/1047776‑1048067

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