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 3,934,995 +C intergenic (‑308/‑173) yieP ← / → BW25113_RS19520 FadR/GntR family transcriptional regulator/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,934,9911.C93.2% 123.2 / 3.3 44intergenic (‑304/‑177)yieP/BW25113_RS19520FadR/GntR family transcriptional regulator/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (24/17);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 8.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTT  >  NZ_CP009273/3934912‑3935079
                                                                                  |                                                                                        
aGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGcgcc                                                                                 >  1:1438301/1‑90 (MQ=255)
      gTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTgag                                                                           >  1:1175514/1‑90 (MQ=255)
                  aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGTGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                               >  1:263560/1‑90 (MQ=11)
                  aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                               <  1:513220/90‑1 (MQ=255)
                  aCTCCGGCAGAGAAAGCAAAAATAAATGCNTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGACAAaaagcaaagc                                                               <  2:546192/90‑1 (MQ=14)
                   cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                              >  2:362498/1‑90 (MQ=255)
                   cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                              >  1:769/1‑90 (MQ=255)
                     ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                            >  2:1116623/1‑90 (MQ=255)
                       ggCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt                                                         <  1:1519299/90‑1 (MQ=39)
                          agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGc                                                       <  1:1083780/90‑1 (MQ=255)
                             gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGctct                                                    >  1:742302/1‑90 (MQ=255)
                               aaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCttt                                                  >  1:1117749/1‑90 (MQ=255)
                               aaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCttt                                                  >  1:622408/1‑90 (MQ=255)
                                 gCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTaa                                                >  1:229772/1‑90 (MQ=255)
                                    aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa                                             <  2:1405707/90‑1 (MQ=255)
                                    aaaaTAAATGCTAGACTCTGTAGCGGGAAATGGTATTATGCACACCCCCGTGCCGCTGAGAAAAAGCAAAGCGGCGCTGCTCTTTAACaa                                             <  2:804465/90‑1 (MQ=11)
                                          aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATc                                       <  1:407886/90‑1 (MQ=255)
                                          aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATc                                       <  1:1314090/90‑1 (MQ=255)
                                            tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                     <  1:1229032/90‑1 (MQ=255)
                                            tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                     <  2:224869/90‑1 (MQ=255)
                                               ttGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACa                                  >  2:537319/1‑90 (MQ=255)
                                                tGACTCTGTAGCGGGAAAGGGTATTATGCACCCCCCCGCGCCGCCGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa                                 <  2:1320846/90‑1 (MQ=14)
                                                tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa                                 >  1:599780/1‑90 (MQ=255)
                                                     ctGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgt                            >  2:115736/1‑90 (MQ=255)
                                                           cGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGc                      <  2:1422293/90‑1 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                     >  1:138793/1‑74 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                     <  2:138793/74‑1 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                               <  2:1213696/80‑1 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                               >  1:1213696/1‑80 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                     >  2:1419379/1‑90 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                     >  2:802260/1‑90 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGCGTGGGCa                     >  1:839649/1‑90 (MQ=11)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  1:931688/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  1:1315625/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  2:1558873/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  1:654303/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  2:357853/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  <  1:537319/90‑1 (MQ=255)
                                                                 aGCGTATTATGCACACCCCCGCGCCGCTGAGaaaaa                                                                      <  2:292856/36‑1 (MQ=11)
                                                                 aGCGTATTATGCACACCCCCGCGCCGCTGAGaaaaa                                                                      >  1:292856/1‑36 (MQ=11)
                                                                  gCGTATTATGCAACCCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGac               <  1:293550/90‑2 (MQ=11)
                                                                        tatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg         <  1:1419379/90‑1 (MQ=255)
                                                                          tGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGa       <  1:237035/90‑1 (MQ=255)
                                                                          tGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGa       <  2:1035197/90‑1 (MQ=255)
                                                                                cacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCGTTAACAATTTATCAGACAATCTATGTGGGCACTCGAAGATACGGATTCtt  >  2:1421758/1‑90 (MQ=255)
                                                                                  |                                                                                        
AGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTT  >  NZ_CP009273/3934912‑3935079

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ 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
CCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  NZ_CP009273/3934902‑3935090
                                                                                               |                                                                                               
CCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCC                                                                                            >  SRR3722078.1461710/1‑100 (MQ=7)
      GGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAG                                                                                      >  SRR3722078.1193838/1‑100 (MQ=7)
                  TCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGTGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGC                                                                          >  SRR3722078.266777/1‑100 (MQ=23)
                   CTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCG                                                                         >  SRR3722078.776/1‑100 (MQ=24)
                      tccacgaccgccatcgtcggcagcgtcagatgtgtataagagacagGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAC                                                                     <  SRR3722078.999417/54‑1 (MQ=34)
                            nnnnnnnnnnnnnnnnnnnnAATAAATGCTTGACTATGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGAAAAGCGGCACTGCTC                                                                <  SRR3722078.1517910/80‑1 (MQ=23)
                            ACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTC                                                                <  SRR3722078.519295/100‑1 (MQ=24)
                             CTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCT                                                               >  SRR3722078.751693/1‑100 (MQ=24)
                               CCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTT                                                             >  SRR3722078.1134816/1‑100 (MQ=24)
                               CCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTT                                                             >  SRR3722078.629580/1‑100 (MQ=24)
                                 GGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAA                                                           >  SRR3722078.232558/1‑100 (MQ=24)
                                 GGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAAC                                                          <  SRR3722078.1544074/100‑1 (MQ=26)
                                    AGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAA                                                        <  SRR3722078.1100182/100‑1 (MQ=24)
                                                AATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAA                                            >  SRR3722078.606720/1‑100 (MQ=24)
                                                    AATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTG                                        <  SRR3722078.412821/100‑1 (MQ=24)
                                                    AATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTG                                        <  SRR3722078.1335084/100‑1 (MQ=24)
                                                     tcgtcggcagcgtcagatgtgtataagagacaggcacacCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTG                                        <  SRR3722078.1402986/61‑1 (MQ=10)
                                                      TGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTG                                      <  SRR3722078.1248394/100‑1 (MQ=24)
                                                            ACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGctgtct                                >  SRR3722078.140406/1‑94 (MQ=25)
                                                            ACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCA                                >  SRR3722078.1232775/1‑100 (MQ=24)
                                                            ACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGCGTGGGCA                                >  SRR3722078.850834/1‑100 (MQ=23)
                                                               CTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTC                             >  SRR3722078.944540/1‑100 (MQ=24)
                                                               CTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTC                             >  SRR3722078.1336650/1‑100 (MQ=24)
                                                               CTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTC                             >  SRR3722078.661800/1‑100 (MQ=24)
                                                                 GTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCctgtctcttatacacatctgacgctgccgacgatggcggtcgtg                           >  SRR3722078.296425/1‑56 (MQ=30)
                                                                     ataagagacagATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGAT                       <  SRR3722078.1017135/89‑1 (MQ=7)
                                                                         AAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGG                   <  SRR3722078.543670/100‑1 (MQ=19)
                                                                            GCGTATTATGCAACCCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGACGATACGGATT                <  SRR3722078.297126/100‑1 (MQ=5)
                                                                            GCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATT                >  SRR3722078.85585/1‑100 (MQ=6)
                                                                              GTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGnnnnnnnnnnnnnnnnnnnnn              >  SRR3722078.160201/1‑79 (MQ=7)
                                                                              GTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCT              >  SRR3722078.654753/1‑100 (MQ=6)
                                                                              GTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCT              >  SRR3722078.1337013/1‑100 (MQ=6)
                                                                                ATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTA            >  SRR3722078.580972/1‑100 (MQ=6)
                                                                                 TTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAA           >  SRR3722078.644735/1‑100 (MQ=6)
                                                                                  TATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAAC          <  SRR3722078.1442425/100‑1 (MQ=6)
                                                                                   ATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACctgtctcttatacacatctgacgctgccgacgatggc         >  SRR3722078.685837/1‑63 (MQ=8)
                                                                                   ATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACG         >  SRR3722078.1316023/1‑100 (MQ=6)
                                                                                    TGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGT        <  SRR3722078.239924/100‑1 (MQ=6)
                                                                                    TGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTC       >  SRR3722078.45572/1‑100 (MQ=7)
                                                                                          CCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  SRR3722078.1426774/1‑100 (MQ=7)
                                                                                          CCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  SRR3722078.1452765/1‑100 (MQ=7)
                                                                                               |                                                                                               
CCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAG  >  NZ_CP009273/3934902‑3935090

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