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,057,664 T→C G95G (GGA→GGG cobT ← nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,057,6640TC98.3% 184.6 / ‑3.8 58G95G (GGA→GGGcobTnicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (29/28);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCAC  >  NZ_CP009273/2057578‑2057749
                                                                                      |                                                                                     
gggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCa                                                                                    >  2:1382332/1‑90 (MQ=255)
gggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCa                                                                                    >  1:646030/1‑90 (MQ=255)
  gATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACg                                                                                  <  2:333772/90‑1 (MQ=255)
      ggCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTc                                                                              <  2:332732/90‑1 (MQ=255)
        cTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCat                                                                            <  1:1397089/90‑1 (MQ=255)
          cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCatat                                                                          >  1:784986/1‑90 (MQ=255)
          cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCatat                                                                          >  2:553749/1‑90 (MQ=255)
          cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCatat                                                                          >  1:929195/1‑90 (MQ=255)
          cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCatat                                                                          >  2:1076459/1‑90 (MQ=255)
           agcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATAtt                                                                         >  1:685811/1‑90 (MQ=255)
              agTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTc                                                                      <  1:814234/90‑1 (MQ=255)
              agTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTc                                                                      <  1:85953/90‑1 (MQ=255)
               gTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCa                                                                     >  2:768609/1‑90 (MQ=255)
               gTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCa                                                                     >  2:1096926/1‑90 (MQ=255)
               gTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCa                                                                     >  2:1315261/1‑90 (MQ=255)
                 aTCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGc                                                                   >  2:1283245/1‑90 (MQ=255)
                 aTCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGc                                                                   >  1:20549/1‑90 (MQ=255)
                  tCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                  <  2:1314014/90‑1 (MQ=255)
                  tCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                  <  2:1361506/90‑1 (MQ=255)
                    aaTACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTg                                                                >  2:1101930/1‑90 (MQ=255)
                    aaTACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTg                                                                >  1:136032/1‑90 (MQ=255)
                      tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                              >  1:900962/1‑90 (MQ=255)
                      tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                              <  2:20549/90‑1 (MQ=255)
                         cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                           >  1:938779/1‑90 (MQ=255)
                          aaCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGc                                                          >  1:268347/1‑90 (MQ=255)
                          aaCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGc                                                          >  1:43803/1‑90 (MQ=255)
                             aTCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGt                                                       <  2:238502/90‑1 (MQ=255)
                               cAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCa                                                     >  2:838672/1‑90 (MQ=255)
                               cAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCa                                                     >  1:685791/1‑90 (MQ=255)
                                 aTTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAct                                                   <  2:199075/90‑1 (MQ=255)
                                      gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACGCACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                              >  1:653727/1‑90 (MQ=255)
                                      gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                              >  2:1075874/1‑90 (MQ=255)
                                      gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                              >  1:1112268/1‑90 (MQ=255)
                                       tGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                             <  1:1076459/90‑1 (MQ=255)
                                              tAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGaaa                                      >  1:636936/1‑90 (MQ=255)
                                              tAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGaaa                                      >  1:160226/1‑90 (MQ=255)
                                                     cGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAc                               >  1:67662/1‑90 (MQ=255)
                                                     cGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAc                               <  1:47004/90‑1 (MQ=255)
                                                     cGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAc                               <  1:1101930/90‑1 (MQ=255)
                                                     cGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAc                               <  1:1194148/90‑1 (MQ=255)
                                                      gCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                              <  2:537037/90‑1 (MQ=255)
                                                       cTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAccc                             >  1:1433031/1‑90 (MQ=255)
                                                          gTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                          <  2:1231738/90‑1 (MQ=255)
                                                          gTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                          <  2:636936/90‑1 (MQ=255)
                                                            tCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTtcc                        <  1:1075874/90‑1 (MQ=255)
                                                                   gCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACCTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGa                 >  1:1303284/1‑90 (MQ=255)
                                                                      cacacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGc              <  2:1084235/90‑1 (MQ=255)
                                                                       acacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGcc             >  1:15882/1‑90 (MQ=255)
                                                                        cacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCg            <  2:15882/90‑1 (MQ=255)
                                                                        cacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCg            <  1:287566/90‑1 (MQ=255)
                                                                             gCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCCCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATc       <  2:616191/90‑1 (MQ=255)
                                                                             gCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATc       <  1:1096926/90‑1 (MQ=255)
                                                                                ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCCCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  1:998097/90‑1 (MQ=255)
                                                                                ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  2:1428066/90‑1 (MQ=255)
                                                                                ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  2:136032/90‑1 (MQ=255)
                                                                                ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    >  1:185611/1‑90 (MQ=255)
                                                                                ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  1:353843/90‑1 (MQ=255)
                                                                                ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  1:524662/90‑1 (MQ=255)
                                                                                  ttgtcCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGcac  <  1:305410/90‑1 (MQ=255)
                                                                                  ttgtcCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGcac  <  2:383244/90‑1 (MQ=255)
                                                                                      |                                                                                     
GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCAC  >  NZ_CP009273/2057578‑2057749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAG  >  NZ_CP009273/2057565‑2057759
                                                                                                   |                                                                                               
gatgtgtataagagacagAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTC                                                                                                 <  SRR3722078.1529117/82‑1 (MQ=60)
   GATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCA                                                                                              >  SRR3722078.653458/1‑100 (MQ=60)
   GATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCA                                                                                              <  SRR3722078.1474423/100‑1 (MQ=60)
      AAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGG                                                                                           <  SRR3722078.1094086/100‑1 (MQ=60)
             GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATAT                                                                                    >  SRR3722078.941981/1‑100 (MQ=60)
             GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATAT                                                                                    >  SRR3722078.795218/1‑100 (MQ=60)
              GGATGGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATT                                                                                   >  SRR3722078.693875/1‑100 (MQ=60)
                    GCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGC                                                                             >  SRR3722078.20809/1‑100 (MQ=60)
                     CTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCC                                                                            <  SRR3722078.1419697/100‑1 (MQ=60)
                       CAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTG                                                                          >  SRR3722078.137607/1‑100 (MQ=60)
                         GCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTA                                                                        >  SRR3722078.913056/1‑100 (MQ=60)
                           AGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATG                                                                      <  SRR3722078.86964/100‑1 (MQ=60)
                           AGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATG                                                                      <  SRR3722078.825006/100‑1 (MQ=60)
                            GTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGG                                                                     >  SRR3722078.951807/1‑100 (MQ=60)
                             TNTCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGC                                                                    >  SRR3722078.44318/1‑100 (MQ=60)
                             TATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGC                                                                    >  SRR3722078.271625/1‑100 (MQ=60)
                                  ATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCA                                                               >  SRR3722078.693855/1‑100 (MQ=60)
                                         CATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACGCACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTT                                                        >  SRR3722078.661216/1‑100 (MQ=60)
                                         CATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTT                                                        >  SRR3722078.1129248/1‑100 (MQ=60)
                                                 ACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAA                                                >  SRR3722078.162098/1‑100 (MQ=60)
                                                 ACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAA                                                >  SRR3722078.644256/1‑100 (MQ=60)
                                                    TGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAG                                             <  SRR3722078.1092727/100‑1 (MQ=60)
                                                        CGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAC                                         >  SRR3722078.68453/1‑100 (MQ=60)
                                                          TTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCC                                       >  SRR3722078.1456337/1‑100 (MQ=60)
                                                                  CGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCC                               <  SRR3722078.47558/100‑1 (MQ=60)
                                                                  CGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCC                               <  SRR3722078.1212854/100‑1 (MQ=60)
                                                                  CGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCC                               <  SRR3722078.1118723/100‑1 (MQ=60)
                                                                      TGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACCTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGA                           >  SRR3722078.1324064/1‑100 (MQ=60)
                                                                         TCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGC                        <  SRR3722078.1092129/100‑1 (MQ=60)
                                                                          CTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCC                       >  SRR3722078.16080/1‑100 (MQ=60)
                                                                                   CACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGC              >  SRR3722078.187819/1‑100 (MQ=60)
                                                                                     CACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCAC            <  SRR3722078.291073/100‑1 (MQ=60)
                                                                                          GCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATA       <  SRR3722078.1113615/100‑1 (MQ=60)
                                                                                             GGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCCCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACC    <  SRR3722078.1012646/100‑1 (MQ=60)
                                                                                             GGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACC    <  SRR3722078.530866/100‑1 (MQ=60)
                                                                                             GGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACC    <  SRR3722078.358158/100‑1 (MQ=60)
                                                                                               TTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAG  <  SRR3722078.309139/100‑1 (MQ=60)
                                                                                                   |                                                                                               
GTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAG  >  NZ_CP009273/2057565‑2057759

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: