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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,972,691 A→G intergenic (‑237/‑543) flhD ← / → uspC flagellar transcriptional regulator FlhD/universal stress protein UspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,972,6910AG98.0% 160.9 / ‑3.7 52intergenic (‑237/‑543)flhD/uspCflagellar transcriptional regulator FlhD/universal stress protein UspC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (27/23);  total (28/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 = 4.26e-01

TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACT  >  NZ_CP009273/1972605‑1972772
                                                                                       |                                                                                 
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                >  2:1327177/1‑90 (MQ=255)
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                >  1:122503/1‑90 (MQ=255)
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                >  2:1000456/1‑90 (MQ=255)
      tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTg                                                                          >  1:650257/1‑90 (MQ=255)
       aGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGc                                                                         <  1:901577/90‑1 (MQ=255)
       aGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGc                                                                         <  1:1560786/90‑1 (MQ=255)
              ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACTCGCACATACAACGGGGGGGGGCTGCGATTTTc                                                                  <  2:445342/90‑1 (MQ=255)
              ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTc                                                                  >  2:1038045/1‑90 (MQ=255)
              ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTc                                                                  <  1:228505/90‑1 (MQ=255)
                  tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaata                                                              >  1:587368/1‑90 (MQ=255)
                   ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                             >  1:110897/1‑90 (MQ=255)
                   ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaataa                                                             >  1:125028/1‑90 (MQ=255)
                     acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATg                                                           <  1:180216/90‑1 (MQ=255)
                        cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                        >  2:528673/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCCGCACATACAACGGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                        <  1:1099852/90‑1 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                        <  2:1551161/89‑1 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGt                                                        >  1:1551161/1‑89 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                       >  2:172178/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                       >  2:1051709/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                       >  1:1376083/1‑90 (MQ=255)
                          tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGa                                                      >  2:777503/1‑90 (MQ=255)
                           aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGAt                                                     >  1:1509431/1‑90 (MQ=255)
                                   ttttttCTAAATCGACGCAACTGTTCTCGTCACCACCCCGCCATACCACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                             <  1:1449381/90‑1 (MQ=255)
                                     tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATcaca                                           <  2:1023057/90‑1 (MQ=255)
                                             aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAacac                                   <  1:68579/90‑1 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAt                              <  1:159688/90‑1 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAt                              <  2:1239950/90‑1 (MQ=255)
                                                   gCAACTGTTCTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATTCAGATCACACAAAACACTCAAtt                             <  2:1373297/90‑1 (MQ=255)
                                                         gTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTaa                       >  2:76730/1‑90 (MQ=255)
                                                           aCTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                     >  1:1499088‑M1/1‑88 (MQ=255)
                                                           aCTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                     >  1:839213‑M1/1‑88 (MQ=255)
                                                            cTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                    <  2:1456717‑M1/90‑4 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACaa                                        >  1:1356699/1‑68 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACaa                                        <  2:1356699/68‑1 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                  >  2:875710‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                  >  1:722095‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaagg                 >  2:506573‑M1/1‑84 (MQ=255)
                                                               gTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAGTGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaagg                 <  2:371067‑M1/90‑7 (MQ=255)
                                                               gTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaagg                 <  1:510021‑M1/90‑7 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                             <  1:1116453/61‑1 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATca                                             >  2:1116453/1‑61 (MQ=255)
                                                                tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                >  1:430875‑M1/1‑83 (MQ=255)
                                                                   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg             >  2:840751‑M1/1‑80 (MQ=255)
                                                                        cGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa        >  1:983682‑M1/1‑75 (MQ=255)
                                                                        cGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa        >  2:1571058‑M1/1‑75 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   >  2:1264749‑M1/1‑70 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   >  2:1099858‑M1/1‑70 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   <  2:722095‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   <  1:1038045‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   >  1:557014‑M1/1‑70 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc   <  2:86041‑M1/90‑21 (MQ=255)
                                                                              tACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccct  <  1:342933‑M1/90‑22 (MQ=255)
                                                                              tACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccct  <  2:994423‑M1/90‑22 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa        >  2:311011‑M1/1‑65 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa        <  1:311011‑M1/80‑16 (MQ=255)
                                                                                       |                                                                                 
TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACT  >  NZ_CP009273/1972605‑1972772

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 32 ≤ 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
TAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGG  >  NZ_CP009273/1972594‑1972782
                                                                                                  |                                                                                           
cggcgaccaccgagatctacacgaccgccatcgtcggcagcgtcagatgtgtataagagacaGCAACTGTACTCGTCACTACACGCACATACAACGGGGG                                                                                           <  SRR3722078.1008246/38‑1 (MQ=60)
 AATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGG                                                                                          >  SRR3722078.123930/1‑100 (MQ=60)
     TGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGC                                                                                      <  SRR3722078.907765/100‑1 (MQ=60)
       TCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTG                                                                                    >  SRR3722078.657719/1‑100 (MQ=60)
                  AGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAAT                                                                         <  SRR3722078.913675/100‑1 (MQ=60)
                  AGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAAT                                                                         <  SRR3722078.1586188/100‑1 (MQ=60)
                   GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGGCTGCGATTTTCAAT                                                                         <  SRR3722078.545034/100‑1 (MQ=60)
                   GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATA                                                                        >  SRR3722078.594204/1‑100 (MQ=60)
                    AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAA                                                                       >  SRR3722078.112209/1‑100 (MQ=60)
                    AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAA                                                                       >  SRR3722078.126475/1‑100 (MQ=60)
                         GGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGT                                                                  <  SRR3722078.231276/100‑1 (MQ=60)
                          GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTG                                                                 >  SRR3722078.1398270/1‑100 (MQ=60)
                          GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTG                                                                 >  SRR3722078.1576428/1‑100 (MQ=60)
                            CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGAT                                                               >  SRR3722078.1534047/1‑100 (MQ=60)
                                ACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAG                                                           <  SRR3722078.182339/100‑1 (MQ=60)
                                    TTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCCGCACATACAACGGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATC                                                        <  SRR3722078.1116603/100‑1 (MQ=60)
                                              agatgtgtataagagacagACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAC                                             <  SRR3722078.1133499/81‑1 (MQ=60)
                                              TTTTTTCTAAATCGACGCAACTGTTCTCGTCACCACCCCGCCATACCACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAC                                             <  SRR3722078.1473018/100‑1 (MQ=60)
                                                        ATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTT                                   <  SRR3722078.69381/100‑1 (MQ=60)
                                                            ACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                               >  SRR3722078.1523537/1‑98 (MQ=60)
                                                            ACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                               >  SRR3722078.850385/1‑98 (MQ=60)
                                                             CGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                              <  SRR3722078.161555/100‑4 (MQ=60)
                                                               CAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATctgtctcttata                            >  SRR3722078.1378480/1‑88 (MQ=60)
                                                               CAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                            >  SRR3722078.731055/1‑95 (MQ=60)
                                                                 ACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                          >  SRR3722078.436060/1‑93 (MQ=60)
                                                                  atntannnnnnnntaaangaACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtg                         <  SRR3722078.1152316/80‑9 (MQ=60)
                                                                         CGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                  >  SRR3722078.997903/1‑85 (MQ=60)
                                                                          GTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaag                 <  SRR3722078.516061/100‑17 (MQ=60)
                                                                              CTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc             >  SRR3722078.563577/1‑80 (MQ=60)
                                                                                   agagacagCAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata        <  SRR3722078.445776/92‑26 (MQ=60)
                                                                                   CGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata        <  SRR3722078.314820/100‑26 (MQ=60)
                                                                                     CACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg      >  SRR3722078.953539/1‑73 (MQ=60)
                                                                                        ATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga   <  SRR3722078.1053512/100‑31 (MQ=60)
                                                                                         TACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagat  <  SRR3722078.347127/100‑32 (MQ=60)
                                                                                                  |                                                                                           
TAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGG  >  NZ_CP009273/1972594‑1972782

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