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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 674,442 T→G T141P (ACT→CCT)  ybeT ← tetratricopeptide repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273674,4420TG100.0% 174.9 / NA 50T141P (ACT→CCT) ybeTtetratricopeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/28);  total (22/28)

GTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAG  >  NZ_CP009273/674355‑674523
                                                                                        |                                                                                 
gttgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGgat                                                                                 >  1:1193561/1‑90 (MQ=255)
gttgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGgat                                                                                 >  1:842275/1‑90 (MQ=255)
  tgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGatat                                                                               <  2:369964/90‑1 (MQ=255)
  tgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGatat                                                                               <  2:1043415/90‑1 (MQ=255)
  tgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGatat                                                                               >  1:926508/1‑90 (MQ=255)
  tgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGatat                                                                               <  1:655316/90‑1 (MQ=255)
   gttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATAtt                                                                              <  1:369677/90‑1 (MQ=255)
     ttGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtga                                                                            >  2:1227738/1‑90 (MQ=255)
        gCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagca                                                                         <  2:558710/90‑1 (MQ=255)
        gCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagca                                                                         <  2:511106/90‑1 (MQ=255)
        gCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagca                                                                         <  2:275106/90‑1 (MQ=255)
         cGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                        <  2:1202851/90‑1 (MQ=255)
         cGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                        <  2:1159743/90‑1 (MQ=255)
         cGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                        <  2:722419/90‑1 (MQ=255)
          ggTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcatg                                                                       <  2:181709/90‑1 (MQ=255)
          ggTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcatg                                                                       <  2:296029/90‑1 (MQ=255)
          ggTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcatg                                                                       >  1:1199073/1‑90 (MQ=255)
             cAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagc                                                                          >  2:1222426/1‑84 (MQ=255)
             cAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagc                                                                          <  1:1222426/84‑1 (MQ=255)
                    aaCCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCt                                                             <  1:1011825/90‑1 (MQ=255)
                         tttttGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGGTCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTc                                                        >  2:1132814/1‑90 (MQ=255)
                         tttttGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTc                                                        >  2:338255/1‑90 (MQ=255)
                               agaTTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgc                                                  <  2:34150/90‑1 (MQ=255)
                               aagatCAAAAAGGTTAAAAATTAATCCCTGGGCCTTTTTTTAAGATAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgc                                                  <  1:241432/86‑1 (MQ=255)
                                gaTTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgct                                                 <  1:911830/90‑1 (MQ=255)
                                           gTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCa                                      >  1:249982/1‑90 (MQ=255)
                                             tAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCaga                                    >  1:30267/1‑90 (MQ=255)
                                             tAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCaga                                    >  2:839478/1‑90 (MQ=255)
                                                cATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGaa                                 <  2:586397/90‑1 (MQ=255)
                                                cATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGaa                                 >  1:315314/1‑90 (MQ=255)
                                                  aTTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAt                                <  1:728991/90‑1 (MQ=255)
                                                  aTTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAt                                <  1:338255/90‑1 (MQ=255)
                                                   ttAATACCTGAGCCATTTTTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAtt                               <  2:738831/90‑1 (MQ=255)
                                                     aaTACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGc                             >  1:239458/1‑90 (MQ=255)
                                                     aaTACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGc                             >  1:1034018/1‑90 (MQ=255)
                                                     aaTACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATCGc                             <  2:147585/90‑1 (MQ=255)
                                                          cTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTAc                        <  2:842275/90‑1 (MQ=255)
                                                           tGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                        <  1:1027343/41‑1 (MQ=255)
                                                           tGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                        >  2:1027343/1‑41 (MQ=255)
                                                           tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                       >  2:10576/1‑90 (MQ=255)
                                                           tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                       >  1:272572/1‑90 (MQ=255)
                                                           tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                       >  2:1001276/1‑90 (MQ=255)
                                                           tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                       >  1:491259/1‑90 (MQ=255)
                                                            gAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACgg                      >  1:331192/1‑90 (MQ=255)
                                                             aGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACggg                     >  1:1036011/1‑90 (MQ=255)
                                                                 aTTTCTGTAGATCAAACCTAAAGGATATTgagcatgagcat                                                                  >  2:944624/1‑41 (MQ=255)
                                                                 aTTTCTGTAGATCAAACCTAAAGGATATTgagcatgagcat                                                                  <  1:944624/41‑1 (MQ=255)
                                                                         aGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGa         <  1:839478/90‑1 (MQ=255)
                                                                         aGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGa         <  1:1001276/90‑1 (MQ=255)
                                                                               aaCCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCa   >  1:794966/1‑90 (MQ=255)
                                                                               aaCCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCa   >  1:499363/1‑90 (MQ=255)
                                                                                aCCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAg  <  2:331192/90‑1 (MQ=255)
                                                                                        |                                                                                 
GTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAG  >  NZ_CP009273/674355‑674523

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGGTAGATCC  >  NZ_CP009273/674345‑674540
                                                                                                 |                                                                                                  
TTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGAT                                                                                                  >  SRR3722089.854219/1‑100 (MQ=60)
TTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGAT                                                                                                  >  SRR3722089.1211706/1‑100 (MQ=60)
  TTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATAT                                                                                                >  SRR3722089.940029/1‑100 (MQ=60)
       AATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGC                                                                                           <  SRR3722089.53960/100‑1 (MQ=60)
          GTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATG                                                                                        >  SRR3722089.1217303/1‑100 (MQ=60)
            TGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAG                                                                                      <  SRR3722089.663748/100‑1 (MQ=60)
             GTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGC                                                                                     <  SRR3722089.374220/100‑1 (MQ=60)
                 GGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATAC                                                                                 <  SRR3722089.1240999/100‑1 (MQ=60)
                   agatggtcgtcggcagcgtcagatgtgtataagagacagCATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCC                                                                               <  SRR3722089.1042673/61‑1 (MQ=60)
                         agatggtcgtcggcagcgtcagatgtgtataagagacagTACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTC                                                                         <  SRR3722089.958472/61‑1 (MQ=60)
                              AACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTG                                                                    <  SRR3722089.1026896/100‑1 (MQ=60)
                                         aagatcaaaaaggttaaaaattaatccctgggcctttttttaagatAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATAC                                                         <  SRR3722089.244504/54‑1 (MQ=60)
                                          GATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACC                                                        <  SRR3722089.925102/100‑1 (MQ=60)
                                           ATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCA                                                       >  SRR3722089.253171/1‑100 (MQ=60)
                                             TCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGA                                                     >  SRR3722089.30636/1‑100 (MQ=60)
                                                TAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAA                                                  >  SRR3722089.319311/1‑100 (MQ=60)
                                                    GGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGC                                              >  SRR3722089.1049474/1‑100 (MQ=60)
                                                    GGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGC                                              >  SRR3722089.242498/1‑100 (MQ=60)
                                                          CATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACG                                        >  SRR3722089.276055/1‑100 (MQ=60)
                                                          CATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACG                                        >  SRR3722089.497078/1‑100 (MQ=60)
                                                           ATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGG                                       <  SRR3722089.738475/100‑1 (MQ=60)
                                                           ATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGG                                       <  SRR3722089.342474/100‑1 (MQ=60)
                                                           ATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGG                                       >  SRR3722089.335335/1‑100 (MQ=60)
                                                            TTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGG                                      >  SRR3722089.1051504/1‑100 (MQ=60)
                                                                             TCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCctgtctcttatacacatctgacgctgccgacgaccatcttagtgtagatctc                     >  SRR3722089.479565/1‑48 (MQ=60)
                                                                              CTGTGGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCA                    >  SRR3722089.805920/1‑100 (MQ=60)
                                                                              CTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCA                    >  SRR3722089.505271/1‑100 (MQ=60)
                                                                                  AGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACC                <  SRR3722089.1016162/100‑1 (MQ=60)
                                                                                  AGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACC                <  SRR3722089.851369/100‑1 (MQ=60)
                                                                                        AACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAG          >  SRR3722089.9740/1‑100 (MQ=60)
                                                                                         ACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGG         >  SRR3722089.1096400/1‑100 (MQ=60)
                                                                                         ACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGG         >  SRR3722089.218727/1‑100 (MQ=60)
                                                                                          CCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGGT        >  SRR3722089.552592/1‑100 (MQ=60)
                                                                                               AGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGGctgtct   >  SRR3722089.740452/1‑94 (MQ=60)
                                                                                               AGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGGTAGATC   >  SRR3722089.533811/1‑100 (MQ=60)
                                                                                                GGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGGTAGATCC  >  SRR3722089.675270/1‑100 (MQ=60)
                                                                                                GGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGGTAGATCC  <  SRR3722089.1141098/100‑1 (MQ=60)
                                                                                                GGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGGTAGATCC  >  SRR3722089.502267/1‑100 (MQ=60)
                                                                                                 |                                                                                                  
TTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGACCCCCAGACCTTTGAGGTAGATCC  >  NZ_CP009273/674345‑674540

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