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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,786,991 T→C D18G (GAT→GGT)  yibB ← protein YibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,786,9910TC100.0% 179.6 / NA 50D18G (GAT→GGT) yibBprotein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/32);  total (18/32)

TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786904‑3787073
                                                                                       |                                                                                  
tCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAAccc                                                                                  <  1:613558/90‑1 (MQ=255)
tCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAAccc                                                                                  <  1:211333/90‑1 (MQ=255)
  aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                >  1:392314/1‑90 (MQ=255)
     tGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                             >  1:756588/1‑90 (MQ=255)
     tGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                             >  2:381403/1‑90 (MQ=255)
      gCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTa                                                                            >  2:29659/1‑90 (MQ=255)
         gCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCa                                                                         <  1:684845/90‑1 (MQ=255)
          ccAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCaa                                                                        <  1:167672/90‑1 (MQ=255)
               cGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                   <  2:392314/90‑1 (MQ=255)
               cGCACAAAGTACCTAAAATAAACATCACCCCAACGATCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCa                                                                   <  2:379807/90‑1 (MQ=255)
                    aaaGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAata                                                              <  2:633559/90‑1 (MQ=255)
                      aGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATg                                                            <  1:231774/90‑1 (MQ=255)
                       gTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGc                                                           <  1:652963/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  2:477798/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  2:367728/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  1:582812/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  1:600132/90‑1 (MQ=255)
                          aCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGt                                                        >  1:628353/1‑90 (MQ=255)
                           cTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  2:678886/90‑1 (MQ=255)
                             aaaaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGAt                                                     <  1:29659/90‑1 (MQ=255)
                               aaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATaa                                                   >  1:250940/1‑90 (MQ=255)
                                aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                  <  1:686994/90‑1 (MQ=255)
                                aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                  <  2:893807/90‑1 (MQ=255)
                                aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                  <  2:359110/90‑1 (MQ=255)
                                aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                                  <  1:179730/90‑1 (MQ=255)
                                       tCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtggtg                                           >  2:255141/1‑90 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  1:682987/90‑1 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  1:798534/90‑1 (MQ=255)
                                            tGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatga                                      >  2:345732/1‑90 (MQ=255)
                                             gAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatgat                                     >  1:743570/1‑90 (MQ=255)
                                                cGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTc                                  <  1:381403/90‑1 (MQ=255)
                                                  aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                >  1:664266/1‑90 (MQ=255)
                                                  aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                <  1:658286/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:98857/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  1:345732/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  2:302998/90‑1 (MQ=255)
                                                     tAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAc                             <  2:664266/90‑1 (MQ=255)
                                                          tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                        <  1:225781/90‑1 (MQ=255)
                                                            tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAAt                                                                       >  1:451739/1‑41 (MQ=255)
                                                            tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAAt                                                                       <  2:451739/41‑1 (MQ=255)
                                                            tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCGACCAATGt                                                                     >  1:138165/1‑43 (MQ=255)
                                                            tCACCGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGt                                                                     <  2:138165/43‑1 (MQ=255)
                                                             cACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGa                     <  1:565746/90‑1 (MQ=255)
                                                                              gCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATaaaaaaa    >  1:598952/1‑90 (MQ=255)
                                                                               cAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  2:836013/38‑1 (MQ=39)
                                                                               cAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       >  1:836013/1‑38 (MQ=255)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  1:429141/37‑1 (MQ=39)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       >  2:429141/1‑37 (MQ=39)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa  >  1:442373/1‑90 (MQ=255)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa  >  2:118190/1‑90 (MQ=255)
                                                                                 gTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  2:540518/52‑1 (MQ=255)
                                                                                 gTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       >  1:540518/1‑52 (MQ=255)
                                                                                       |                                                                                  
TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786904‑3787073

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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
TCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATAAATATT  >  NZ_CP009273/3786896‑3787090
                                                                                               |                                                                                                   
TCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCC                                                                                                 >  SRR3722072.399365/1‑100 (MQ=60)
   CATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCT                                                                                              >  SRR3722072.769942/1‑100 (MQ=60)
     TTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAC                                                                                            <  SRR3722072.430857/100‑1 (MQ=60)
      TTTCAAGTGCGGCCAGCCGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGACCAACCCCCTCTACC                                                                                           <  SRR3722072.45233/100‑1 (MQ=60)
      TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722072.905839/100‑1 (MQ=60)
      TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722072.890379/100‑1 (MQ=60)
      TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722072.684165/100‑1 (MQ=60)
      TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                                                           <  SRR3722072.326790/100‑1 (MQ=60)
        TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAA                                                                                         <  SRR3722072.214930/100‑1 (MQ=60)
        TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAA                                                                                         <  SRR3722072.624494/100‑1 (MQ=60)
                 GCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAAT                                                                                <  SRR3722072.697104/100‑1 (MQ=60)
                  CCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATA                                                                               <  SRR3722072.170569/100‑1 (MQ=60)
                        GCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGT                                                                         >  SRR3722072.639572/1‑100 (MQ=60)
                             AAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAA                                                                    >  SRR3722072.255190/1‑100 (MQ=60)
                              AGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAT                                                                   <  SRR3722072.235654/100‑1 (MQ=60)
                               GTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATT                                                                  <  SRR3722072.664616/100‑1 (MQ=60)
                                TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722072.593217/100‑1 (MQ=60)
                                TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                                 <  SRR3722072.610837/100‑1 (MQ=60)
                                   taatcggacgtctcgtgggctcggagatgtgtataagagacagCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGG                                                              <  SRR3722072.436953/57‑1 (MQ=60)
                                     AAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTG                                                            <  SRR3722072.30181/100‑1 (MQ=60)
                                        ATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGAT                                                         <  SRR3722072.699289/100‑1 (MQ=60)
                                        ATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGAT                                                         <  SRR3722072.182842/100‑1 (MQ=60)
                                           AACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGAT                                                      >  SRR3722072.756727/1‑100 (MQ=60)
                                                CAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAT                                                 >  SRR3722072.676142/1‑100 (MQ=60)
                                                   CTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAC                                              <  SRR3722072.812492/100‑1 (MQ=60)
                                                   CTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAC                                              <  SRR3722072.695208/100‑1 (MQ=60)
                                                        CGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATA                                         <  SRR3722072.388220/100‑1 (MQ=60)
                                                          AGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATctgtctcttatacacatctccgagcccacgagacgtccg                                       >  SRR3722072.460007/1‑61 (MQ=60)
                                                          AGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAG                                       <  SRR3722072.670039/100‑1 (MQ=60)
                                                          AGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCGACCAATGTCAAAATAggcctgaatcttatacacatctccgagtccacgggacgtc                                       >  SRR3722072.140627/1‑60 (MQ=60)
                                                            CTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGAC                                     <  SRR3722072.351820/100‑1 (MQ=60)
                                                                  TTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAA                               <  SRR3722072.229599/100‑1 (MQ=60)
                                                                     CACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATAT                            <  SRR3722072.575882/100‑1 (MQ=60)
                                                                            CCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAA                     >  SRR3722072.609635/1‑100 (MQ=60)
                                                                             CCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGctgtctcttatacacatctccgagcccacgagacgtccgatt                    >  SRR3722072.850497/1‑58 (MQ=60)
                                                                              CTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA                   >  SRR3722072.450470/1‑100 (MQ=60)
                                                                               TTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACctgtctcttatacacatctccgagccca                  >  SRR3722072.550252/1‑72 (MQ=60)
                                                                                         GTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATA        >  SRR3722072.679338/1‑100 (MQ=60)
                                                                                               CCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATAAATATT  >  SRR3722072.876061/1‑100 (MQ=60)
                                                                                               |                                                                                                   
TCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATAAATATT  >  NZ_CP009273/3786896‑3787090

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