Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F2 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,034,572 T→G Q74H (CAA→CAC mobB ← molybdopterin‑guanine dinucleotide biosynthesis protein MobB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,034,5720TG86.0% 141.9 / 12.7 50Q74H (CAA→CACmobBmolybdopterin‑guanine dinucleotide biosynthesis protein MobB
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (17/26);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGG  >  NZ_CP009273/4034486‑4034659
                                                                                      |                                                                                       
tCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgc                                                                                      <  2:444619/90‑1 (MQ=255)
  cAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctg                                                                                    >  1:532398/1‑90 (MQ=255)
    gCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTgctggc                                                                                  >  2:631887/1‑90 (MQ=255)
        tGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCtggctggc                                                                              <  2:658603/90‑1 (MQ=255)
        tGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCtggctggc                                                                              >  1:588271/1‑90 (MQ=255)
         gAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCa                                                                             <  2:417089/90‑1 (MQ=255)
             tATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGa                                                                         <  2:754391/90‑1 (MQ=255)
                  aTTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGtt                                                                    <  1:698947/90‑1 (MQ=255)
                      gACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTtgtg                                                                >  2:450008/1‑90 (MQ=255)
                       aCTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGc                                                               <  2:23056/90‑1 (MQ=255)
                        cTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGcc                                                              <  2:588271/90‑1 (MQ=255)
                          tGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCcgc                                                            >  2:124876/1‑90 (MQ=255)
                             gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                         <  2:4720/90‑1 (MQ=255)
                                 aaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGcc                                                     >  2:265507/1‑90 (MQ=255)
                                     gTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgc                                                 >  1:708625/1‑90 (MQ=255)
                                     gTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgc                                                 >  1:339638/1‑90 (MQ=255)
                                       aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                               <  1:527925/90‑1 (MQ=255)
                                       aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                               <  1:499172/90‑1 (MQ=255)
                                       aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                               <  1:639665/90‑1 (MQ=255)
                                             aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCa                                         >  2:350186/1‑90 (MQ=255)
                                             aGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCa                                         >  1:684816/1‑90 (MQ=255)
                                                 cttcttCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                     <  2:370825/90‑1 (MQ=255)
                                                 cttcttCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                     <  1:734663/90‑1 (MQ=255)
                                                 cttcttCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                     <  1:381693/90‑1 (MQ=255)
                                                   tcttcGTCTGGTGTTTCTGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTa                                   >  1:119408/1‑90 (MQ=255)
                                                     ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATc                                 <  1:373533/90‑1 (MQ=255)
                                                               gtTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc                       >  1:676202/1‑90 (MQ=255)
                                                                tttCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                      <  1:820628/90‑1 (MQ=255)
                                                                tttCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                      <  2:221484/90‑1 (MQ=255)
                                                                 ttCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCtt                     <  2:505270/90‑1 (MQ=255)
                                                                 ttCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCtt                     >  2:120495/1‑90 (MQ=255)
                                                                 ttCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCtt                     >  1:142376/1‑90 (MQ=255)
                                                                       catcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc               >  2:621622/1‑90 (MQ=255)
                                                                        atcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCGGCTCATAGCTATCTTTGCCTGGCTTATCAACa              <  2:228690/90‑1 (MQ=255)
                                                                          caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc            <  2:390646/90‑1 (MQ=255)
                                                                          caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc            <  1:291500/90‑1 (MQ=255)
                                                                          caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc            <  1:388500/90‑1 (MQ=255)
                                                                          caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc            <  1:528247/90‑1 (MQ=255)
                                                                            aGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCa          >  2:750591/1‑90 (MQ=255)
                                                                            aGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCa          <  2:730634/90‑1 (MQ=255)
                                                                            aGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCa          <  2:838293/90‑1 (MQ=255)
                                                                             ggCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCat         <  2:666341/90‑1 (MQ=255)
                                                                                 cATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATc                                 <  2:220427/62‑1 (MQ=255)
                                                                                 cATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATc                                 >  1:220427/1‑62 (MQ=255)
                                                                                  aTCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                      >  2:549800/1‑72 (MQ=255)
                                                                                  aTCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                      <  1:549800/72‑1 (MQ=255)
                                                                                  aTCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCAt    >  2:628738/1‑90 (MQ=255)
                                                                                   tCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg   <  1:523444/90‑1 (MQ=255)
                                                                                    cGTTGCTGGCTGGCAACGATGGTTTGTTCCGCGCCAGCCTTGCGCGGCTCATAGCTCTCTTTGCCTGGCTTATCAACATCCATATCATgg  <  2:220664/90‑1 (MQ=255)
                                                                                    cGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATgg  <  1:120495/90‑1 (MQ=255)
                                                                                      |                                                                                       
TCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGG  >  NZ_CP009273/4034486‑4034659

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  NZ_CP009273/4034476‑4034671
                                                                                                |                                                                                                   
CAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGC                                                                                                  <  SRR3722079.766037/100‑1 (MQ=60)
CAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGC                                                                                                  <  SRR3722079.429441/100‑1 (MQ=60)
  GAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTG                                                                                                >  SRR3722079.537366/1‑100 (MQ=60)
        AATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGC                                                                                          >  SRR3722079.594131/1‑100 (MQ=60)
                            ATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGC                                                                      <  SRR3722079.706335/100‑1 (MQ=60)
                                     GCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGC                                                             >  SRR3722079.716148/1‑100 (MQ=60)
                                     GCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGC                                                             >  SRR3722079.342317/1‑100 (MQ=60)
                                             CTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCA                                                     >  SRR3722079.692015/1‑100 (MQ=60)
                                                 AGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGC                                                 <  SRR3722079.532837/100‑1 (MQ=60)
                                                 AGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGC                                                 <  SRR3722079.646254/100‑1 (MQ=60)
                                                 AGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGC                                                 <  SRR3722079.503640/100‑1 (MQ=60)
                                                   ATCCAGCTCTTCTTCGTCTGGTGTTTCTGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTA                                               >  SRR3722079.120390/1‑100 (MQ=60)
                                                           CTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCC                                       <  SRR3722079.384800/100‑1 (MQ=60)
                                                           CTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCC                                       <  SRR3722079.742520/100‑1 (MQ=60)
                                                               TTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGC                                   <  SRR3722079.376555/100‑1 (MQ=60)
                                                               TTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGC                                   >  SRR3722079.683297/1‑100 (MQ=60)
                                                                 CGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTT                                 >  SRR3722079.143587/1‑100 (MQ=60)
                                                                          agagacagATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATC                        <  SRR3722079.555066/92‑1 (MQ=60)
                                                                          TTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATC                        <  SRR3722079.829419/100‑1 (MQ=60)
                                                                                 CATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCctgtctcttatacacatc                 >  SRR3722079.222304/1‑82 (MQ=60)
                                                                                    CAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATNG              <  SRR3722079.435556/100‑1 (MQ=60)
                                                                                    CAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGG              <  SRR3722079.391707/100‑1 (MQ=60)
                                                                                    CAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGG              <  SRR3722079.533161/100‑1 (MQ=60)
                                                                                    CAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGG              <  SRR3722079.293901/100‑1 (MQ=60)
                                                                                       GGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGC           >  SRR3722079.136381/1‑100 (MQ=60)
                                                                                             TCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGC     <  SRR3722079.528291/100‑1 (MQ=60)
                                                                                              CGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCT    <  SRR3722079.121484/100‑1 (MQ=60)
                                                                                              CGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCT    <  SRR3722079.352944/100‑1 (MQ=60)
                                                                                                GTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  SRR3722079.610292/1‑100 (MQ=60)
                                                                                                GTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  SRR3722079.8444/1‑100 (MQ=60)
                                                                                                |                                                                                                   
CAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTA  >  NZ_CP009273/4034476‑4034671

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