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 547,953 T→C D33G (GAC→GGC)  purK ← 5‑(carboxyamino)imidazole ribonucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273547,9530TC100.0% 132.7 / NA 38D33G (GAC→GGC) purK5‑(carboxyamino)imidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/21);  total (17/21)

GAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTTCGT  >  NZ_CP009273/547829‑548087
                                                                                                                            |                                                                                                                                      
gAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGc                                                                                                                          >  1:47259/1‑139 (MQ=255)
   cacaTCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGc                                                                                                                          >  2:681828/1‑136 (MQ=255)
      aTCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGa                                                                                                                      >  1:185875/1‑137 (MQ=255)
      aTCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGa                                                                                                                      >  2:739919/1‑137 (MQ=255)
       tCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAg                                                                                                                   <  1:191390/139‑1 (MQ=255)
       tCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAg                                                                                                                   <  1:739919/139‑1 (MQ=255)
         gcgGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCa                                                                                                                 <  2:883550/139‑1 (MQ=255)
         gcgGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCa                                                                                                                 <  2:477774/139‑1 (MQ=255)
               cacaAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGcc                                                                                                            >  2:1026196/1‑138 (MQ=255)
                  aaaGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAAc                                                                                                        >  1:540232/1‑139 (MQ=255)
                    aGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACgg                                                                                                      >  2:331472/1‑139 (MQ=255)
                       ccGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTc                                                                                                   >  1:1233693/1‑139 (MQ=255)
                               cgcgcCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCgc                                                                                                 >  1:686757/1‑133 (MQ=255)
                               cgcgcCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCgc                                                                                                 <  2:686757/133‑1 (MQ=255)
                                     aGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACgca                                                                                     <  1:22899/139‑1 (MQ=255)
                                          gcgGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATa                                                                                <  2:696513/139‑1 (MQ=255)
                                          gcgGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATa                                                                                <  2:185875/139‑1 (MQ=255)
                                          gcgGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGATGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATa                                                                                <  1:965201/139‑1 (MQ=255)
                                                                     ttCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACg                                                     <  1:766788/139‑1 (MQ=255)
                                                                             cAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTg                                             >  1:27138/1‑139 (MQ=255)
                                                                                cGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCa                                                                                                                 <  1:837525/68‑1 (MQ=255)
                                                                                cGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCa                                                                                                                 >  2:837525/1‑68 (MQ=255)
                                                                                cGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGtt                                           <  1:377798/138‑1 (MQ=255)
                                                                                cGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGtt                                           >  2:377798/1‑138 (MQ=255)
                                                                                    aaTCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATg                                      >  2:1090447/1‑139 (MQ=255)
                                                                                     aTCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGc                                     >  2:179408/1‑139 (MQ=255)
                                                                                         cGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCa                                 <  2:1195461/139‑1 (MQ=255)
                                                                                              tgttgAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGc                                     <  1:1146720/130‑1 (MQ=255)
                                                                                              tgttgAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGc                                     >  2:1146720/1‑130 (MQ=255)
                                                                                                        gCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCggg                  <  2:27138/139‑1 (MQ=255)
                                                                                                        gCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCggg                  <  2:1043278/139‑1 (MQ=255)
                                                                                                          aCCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAg                                                                                                                   >  1:76089/1‑40 (MQ=255)
                                                                                                          aCCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAg                                                                                                                   <  2:76089/40‑1 (MQ=255)
                                                                                                                      tCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTtc    >  2:119425/1‑139 (MQ=255)
                                                                                                                        aGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTtcgt  <  1:369124/139‑1 (MQ=255)
                                                                                                                        aGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTtcg   <  2:540232/138‑1 (MQ=255)
                                                                                                                        aGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTtcg   <  1:636015/138‑1 (MQ=255)
                                                                                                                        aGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTtcg   <  1:1074947/138‑1 (MQ=255)
                                                                                                                            |                                                                                                                                      
GAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTTCGT  >  NZ_CP009273/547829‑548087

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ 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
AGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTTCGTCGGTCTGGGCT  >  NZ_CP009273/547809‑548098
                                                                                                                                                |                                                                                                                                                 
AGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCC                                                                                                                                             <  SRR3722129.685960/151‑1 (MQ=60)
 GACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCC                                                                                                                                            <  SRR3722129.4962/151‑1 (MQ=60)
         GCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCC                                                                                                                                    >  SRR3722129.47893/1‑151 (MQ=60)
               ATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGAC                                                                                                                                >  SRR3722129.188493/1‑149 (MQ=60)
                          ATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAAC                                                                                                                   <  SRR3722129.194067/151‑1 (MQ=60)
                          ATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAAC                                                                                                                   <  SRR3722129.750747/151‑1 (MQ=60)
                           TCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACG                                                                                                                  >  SRR3722129.548139/1‑151 (MQ=60)
                            agaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacaGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGG                                                                                                                 <  SRR3722129.849544/91‑1 (MQ=60)
                                GTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCG                                                                                                             >  SRR3722129.1250833/1‑151 (MQ=60)
                                        AGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTG                                                                                                     >  SRR3722129.696878/1‑151 (MQ=60)
                                                       tttttttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacaGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGC                                                                                       <  SRR3722129.1028777/70‑1 (MQ=60)
                                                        CAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCT                                                                                     <  SRR3722129.23187/151‑1 (MQ=60)
                                                             CGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGATGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTG                                                                                <  SRR3722129.978872/151‑1 (MQ=60)
                                                                                      GCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGT                                                       >  SRR3722129.27484/1‑151 (MQ=60)
                                                                                        GTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTT                                                     <  SRR3722129.777948/151‑1 (MQ=60)
                                                                                                  AGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCAC                                           <  SRR3722129.383316/151‑1 (MQ=60)
                                                                                                        AATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCG                                     <  SRR3722129.1162658/151‑1 (MQ=60)
                                                                                                                   GTTGAAAAGGCACCGCCGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaagaggctatact                          >  SRR3722129.77076/1‑62 (MQ=60)
                                                                                                                                   CGCCGGTTCAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTTCGTCGG          >  SRR3722129.807090/1‑151 (MQ=60)
                                                                                                                                           CAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTTCGTCGGTCTGGGCT  <  SRR3722129.374530/151‑1 (MQ=60)
                                                                                                                                           CAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTTCGTCGGTCTGGGC   <  SRR3722129.645468/150‑1 (MQ=60)
                                                                                                                                           CAGCGCCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTTCGTCGGTCTGGGC   <  SRR3722129.1090008/150‑1 (MQ=60)
                                                                                                                                                |                                                                                                                                                 
AGACGGTCAGCAATAATCGGGAACACATCGCGGTTCACAAAGGCCGGATGGCGCGCCAGCTCGCGGGTTAATGCGGTTTCCGGCCAGCGTTCTATCTCAGCGGTAATCACGCTTTGTTGAAAAGGCACCGCCGCCGGTTCAGCGTCCAGCCCGACTGGCCAGACAGCAATGCCTAACGGTTCGCCTGCCTGACGCAGCATACGGCCTAACTGCCCGTTACCGAGGACGCAAACCTGTTTCATGCCGCACCTCGCGGGTCCGGGTTTTCCAGCACTTCGTCGGTCTGGGCT  >  NZ_CP009273/547809‑548098

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