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 2,284,519 T→C S181P (TCC→CCC)  narP → nitrate/nitrite response regulator protein NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,284,5190TC98.1% 169.4 / ‑4.1 53S181P (TCC→CCC) narPnitrate/nitrite response regulator protein NarP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (31/21);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

GCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTT  >  NZ_CP009273/2284433‑2284601
                                                                                      |                                                                                  
gCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCg                                                                                 >  1:1002831/1‑90 (MQ=255)
gCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCg                                                                                 >  1:789703/1‑90 (MQ=255)
gCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCg                                                                                 >  1:20911/1‑90 (MQ=255)
  gTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAg                                                                               >  2:1066587/1‑90 (MQ=255)
   tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGc                                                                              >  2:882259/1‑90 (MQ=255)
   tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGc                                                                              <  1:1066587/90‑1 (MQ=255)
   tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGc                                                                              <  1:352426/90‑1 (MQ=255)
     cTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAg                                                                            <  2:1401584/90‑1 (MQ=255)
        aCGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACa                                                                         <  1:1044862/90‑1 (MQ=255)
        aCGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGACCAGACa                                                                         <  2:1041675/90‑1 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    >  2:1012292/1‑90 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    >  1:350510/1‑90 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    <  2:714392/90‑1 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    >  2:734565/1‑90 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    >  1:119436/1‑90 (MQ=255)
                    cTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACat                                                             <  1:5332/90‑1 (MQ=255)
                       gATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATAtt                                                          <  2:548123/90‑1 (MQ=255)
                            tCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCaa                                                     >  2:1545208/1‑90 (MQ=255)
                            tCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCaa                                                     >  1:1447280/1‑90 (MQ=255)
                              tGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATc                                                   >  2:1465242/1‑90 (MQ=255)
                              tGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATc                                                   >  1:986654/1‑90 (MQ=255)
                                cACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATctg                                                 <  1:1415319/90‑1 (MQ=255)
                                  cGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGCTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATctgct                                               <  1:1545208/90‑1 (MQ=255)
                                    aGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACat                                                             >  2:375593/1‑74 (MQ=255)
                                    aGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACat                                                             <  1:375593/74‑1 (MQ=255)
                                    aGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCtgctgc                                             >  1:692588/1‑90 (MQ=255)
                                      cTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGt                                           <  2:892853/90‑1 (MQ=255)
                                        ggCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                         <  2:351952/90‑1 (MQ=255)
                                         gCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaaa                                        >  2:1339626/1‑90 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAg                                                                            >  2:1030542/1‑49 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAg                                                                            <  1:1030542/49‑1 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                         <  2:1293386/84‑1 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                         >  1:1293386/1‑84 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                   >  2:1160533/1‑90 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                   >  2:1200916/1‑90 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                   >  1:49455/1‑90 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                   >  1:1197445/1‑90 (MQ=255)
                                                  cTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTc                               >  2:1539381/1‑90 (MQ=255)
                                                   tGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTcc                              >  2:594295/1‑90 (MQ=255)
                                                       aaaTAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTc                          <  1:707964/90‑1 (MQ=255)
                                                        aaTAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCa                         >  2:1436788/1‑90 (MQ=255)
                                                              cAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTg                   >  2:741816/1‑90 (MQ=255)
                                                              cAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTg                   >  1:113278/1‑90 (MQ=255)
                                                                gATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                 >  1:177299/1‑90 (MQ=255)
                                                                gATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                 >  2:1469225/1‑90 (MQ=255)
                                                                gATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                 >  2:679312/1‑90 (MQ=255)
                                                                 aTTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcg                <  2:78517/90‑1 (MQ=255)
                                                                 aTTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcg                <  1:594295/90‑1 (MQ=255)
                                                                   tGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcggc              <  2:177299/90‑1 (MQ=255)
                                                                    gCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGcc             <  1:259197/90‑1 (MQ=255)
                                                                       tCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacc          <  2:1581328/90‑1 (MQ=255)
                                                                             ttGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTg    <  1:208914/90‑1 (MQ=255)
                                                                               gAATATTCCCGAGCAGACAATAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGtt  >  1:462762/1‑90 (MQ=255)
                                                                                      |                                                                                  
GCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTT  >  NZ_CP009273/2284433‑2284601

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ 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
AAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGG  >  NZ_CP009273/2284421‑2284616
                                                                                                  |                                                                                                 
AAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCC                                                                                                  <  SRR3722078.473435/100‑1 (MQ=60)
  GATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCG                                                                                                >  SRR3722078.1017494/1‑100 (MQ=60)
  GATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCG                                                                                                >  SRR3722078.800031/1‑100 (MQ=60)
  GATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCG                                                                                                >  SRR3722078.21175/1‑100 (MQ=60)
               TGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAA                                                                                   <  SRR3722078.1082674/100‑1 (MQ=60)
               TGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAA                                                                                   >  SRR3722078.120835/1‑100 (MQ=60)
               TGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAA                                                                                   <  SRR3722078.356730/100‑1 (MQ=60)
               TGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAA                                                                                   >  SRR3722078.354805/1‑100 (MQ=60)
                 gtcggcagcgtcagatgtgtataagagacaGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAG                                                                                 <  SRR3722078.1045840/70‑1 (MQ=60)
                    ACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTAC                                                                              <  SRR3722078.1060478/100‑1 (MQ=60)
                              AGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAA                                                                    >  SRR3722078.1470869/1‑100 (MQ=60)
                                agacagGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATC                                                                  <  SRR3722078.380181/94‑1 (MQ=60)
                                CTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATC                                                                  >  SRR3722078.1000942/1‑100 (MQ=60)
                                CTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATC                                                                  <  SRR3722078.5388/100‑1 (MQ=60)
                                      GTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGC                                                            >  SRR3722078.700815/1‑100 (MQ=60)
                                            CACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAAC                                                      <  SRR3722078.1438304/100‑1 (MQ=60)
                                              CGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGCTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTC                                                    <  SRR3722078.1570386/100‑1 (MQ=60)
                                                AGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAA                                                  >  SRR3722078.1216225/1‑100 (MQ=60)
                                                AGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAA                                                  >  SRR3722078.50031/1‑100 (MQ=60)
                                                AGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAA                                                  >  SRR3722078.1313980/1‑100 (MQ=60)
                                                                GTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTG                                  >  SRR3722078.114613/1‑100 (MQ=60)
                                                                  CAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGC                                >  SRR3722078.179374/1‑100 (MQ=60)
                                                                   AAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCG                               <  SRR3722078.716558/100‑1 (MQ=60)
                                                                             ATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTC                     <  SRR3722078.601190/100‑1 (MQ=60)
                                                                                GCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGT                  <  SRR3722078.262366/100‑1 (MQ=60)
                                                                                 CCTCGGTGTTGAATATTCCCGAGCAGACAATAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTT                 >  SRR3722078.468287/1‑100 (MQ=60)
                                                                                      GTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCtgtctcttatacacatctgacgctgccgacgat            >  SRR3722078.752821/1‑67 (MQ=60)
                                                                                      GTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCctgtctctt            >  SRR3722078.1524349/1‑91 (MQ=60)
                                                                                      GTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGC            >  SRR3722078.1308244/1‑100 (MQ=60)
                                                                                        GTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAA          >  SRR3722078.694530/1‑100 (MQ=60)
                                                                                         TTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAAC         <  SRR3722078.211441/100‑1 (MQ=60)
                                                                                                TTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGG  >  SRR3722078.1022694/1‑100 (MQ=60)
                                                                                                TTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGG  >  SRR3722078.516371/1‑100 (MQ=60)
                                                                                                TTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGG  >  SRR3722078.164956/1‑100 (MQ=60)
                                                                                                  |                                                                                                 
AAGATCCCTTCAGCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTTCCTGCAACAACGCGG  >  NZ_CP009273/2284421‑2284616

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