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,526,050 A→G T95A (ACC→GCC)  cysK → cysteine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,526,0500AG100.0% 225.1 / NA 67T95A (ACC→GCC) cysKcysteine synthase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/33);  total (34/33)

GTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGG  >  NZ_CP009273/2525963‑2526132
                                                                                       |                                                                                  
gTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGcc                                                                                  <  1:976031/90‑1 (MQ=255)
   gAACCGACCAGCGGTACTACCGGGATTGCCCTGGCTTATGTGGCTGCCGCTCGCGGTTACAAACTCCCCCTGACCCTGCCAGAAGCCATg                                                                               <  1:1209763/90‑1 (MQ=255)
   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  1:187281/90‑1 (MQ=255)
   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  1:536628/90‑1 (MQ=255)
   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  1:484097/90‑1 (MQ=255)
    aaCCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGa                                                                              >  1:781819/1‑90 (MQ=255)
       cGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTa                                                                           >  2:336007/1‑90 (MQ=255)
       cGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTa                                                                           >  1:653326/1‑90 (MQ=255)
         aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTAtt                                                                         >  1:1073146/1‑90 (MQ=255)
                  aaTACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAcgccgc                                                                >  1:330147/1‑90 (MQ=255)
                    tACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCaa                                                              <  1:292282/90‑1 (MQ=255)
                       cGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct                                                           >  1:1115278/1‑90 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               >  1:800793/1‑69 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               >  1:1500934/1‑69 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  2:1500934/69‑1 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  2:800793/69‑1 (MQ=255)
                        gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctg                                                          >  2:997937/1‑90 (MQ=255)
                              gCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGaaa                                                    >  2:796371/1‑90 (MQ=255)
                              gCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGaaa                                                    >  1:1132061/1‑90 (MQ=255)
                                   ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt                                               >  1:1430739/1‑90 (MQ=255)
                                   ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt                                               <  2:295863/90‑1 (MQ=255)
                                   ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt                                               <  2:330147/90‑1 (MQ=255)
                                   ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGtt                                               >  1:1082801/1‑90 (MQ=255)
                                    gCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              <  1:421973/90‑1 (MQ=255)
                                    gCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              <  2:1215641/90‑1 (MQ=255)
                                        aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAc                                                                     <  2:1217532/63‑1 (MQ=255)
                                        aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAc                                                                     >  1:1217532/1‑63 (MQ=255)
                                        aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg                                          >  1:585696/1‑90 (MQ=255)
                                         tGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGc                                         <  2:275228/90‑1 (MQ=255)
                                          gTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCa                                        <  1:389855/90‑1 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               >  1:47717/1‑48 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  2:47717/48‑1 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  1:330924/1‑90 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  2:1322550/1‑90 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGGGCAAAc                                     >  1:949025/1‑90 (MQ=255)
                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGCATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  1:448033/1‑90 (MQ=255)
                                              cTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAcc                                    >  2:1021179/1‑90 (MQ=255)
                                                gCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTg                                  >  1:1031122/1‑90 (MQ=255)
                                                  cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              >  2:133066/1‑76 (MQ=255)
                                                  cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              <  1:133066/76‑1 (MQ=255)
                                                  cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt                                >  1:569076/1‑90 (MQ=255)
                                                    cTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTTGGTGCAAACCTGGTGc                              <  1:1417869/90‑1 (MQ=255)
                                                    cTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGc                              <  1:62347/90‑1 (MQ=255)
                                                    cTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGc                              <  1:49595/90‑1 (MQ=255)
                                                     tCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCt                             >  2:940718/1‑90 (MQ=255)
                                                            tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa                      <  1:966480/90‑1 (MQ=255)
                                                            tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa                      <  2:177842/90‑1 (MQ=255)
                                                            tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAACCTGGTGCTGACGGAAg                     <  1:924061/90‑1 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              >  2:214313/1‑90 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              >  2:1573147/1‑90 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              <  1:790567/90‑1 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              <  1:760341/90‑1 (MQ=255)
                                                                    cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa              <  2:585696/90‑1 (MQ=255)
                                                                      cccTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              >  1:385453/1‑56 (MQ=255)
                                                                      cccTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                              <  2:385453/56‑1 (MQ=255)
                                                                         tGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCa         <  2:1430739/90‑1 (MQ=255)
                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct                                                           >  2:1158468/1‑39 (MQ=39)
                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct                                                           <  1:1158468/39‑1 (MQ=39)
                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAt        >  1:1227910/1‑90 (MQ=255)
                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAt        >  1:944681/1‑90 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    >  1:1484356/1‑90 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    >  2:980956/1‑90 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    <  1:1005324/90‑1 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    <  2:1006090/90‑1 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    <  1:980956/90‑1 (MQ=255)
                                                                              aTGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaaa    <  1:917768/90‑1 (MQ=255)
                                                                                gCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAgg  >  1:396778/1‑90 (MQ=255)
                                                                                       |                                                                                  
GTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGG  >  NZ_CP009273/2525963‑2526132

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
GGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCAAAAAGCAG  >  NZ_CP009273/2525951‑2526149
                                                                                                   |                                                                                                   
GGCGTTGAATTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAG                                                                                                     <  SRR3722078.607997/100‑1 (MQ=60)
GGCGTTGAACTGGTTGAACCGACCAGCGGTAACACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAG                                                                                                     <  SRR3722078.1184981/100‑1 (MQ=60)
 agatctccacgcccgccatcgtcggcagcgtcagatgtgtataagagacagATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGC                                                                                                    <  SRR3722078.992816/49‑1 (MQ=60)
 agatctacccgcccgccatcgtcggcagcgtcagatgtgtataagagacagATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGC                                                                                                    <  SRR3722078.577498/49‑1 (MQ=60)
 GCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGC                                                                                                    <  SRR3722078.1453326/100‑1 (MQ=60)
      GAACTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGA                                                                                               >  SRR3722078.791994/1‑100 (MQ=60)
        ACTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGT                                                                                             <  SRR3722078.877461/100‑1 (MQ=60)
         CTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTA                                                                                            >  SRR3722078.660810/1‑100 (MQ=60)
           GGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATT                                                                                          >  SRR3722078.1089347/1‑100 (MQ=60)
            GTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTG                                                                                         <  SRR3722078.990058/100‑1 (MQ=60)
               gcagcgtcagatgtgtataagagacagGCACTGGCCTACTTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAC                                                                                      <  SRR3722078.1425938/73‑1 (MQ=60)
               GAACCGACCAGCGGTACTACCGGGATTGCCCTGGCTTATGTGGCTGCCGCTCGCGGTTACAAACTCCCCCTGACCCTGCCAGAAGCCATGAGTATTGAAC                                                                                      <  SRR3722078.1228772/100‑1 (MQ=60)
               GAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAC                                                                                      <  SRR3722078.542972/100‑1 (MQ=60)
               GAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAC                                                                                      <  SRR3722078.489847/100‑1 (MQ=60)
               GAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAC                                                                                      <  SRR3722078.189507/100‑1 (MQ=60)
                    GACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGC                                                                                 >  SRR3722078.334186/1‑100 (MQ=60)
                         GCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCT                                                                            >  SRR3722078.1132299/1‑100 (MQ=60)
                          CGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACctgtctcttat                                                                           >  SRR3722078.811332/1‑89 (MQ=60)
                          CGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACctgtctcttat                                                                           >  SRR3722078.1525411/1‑89 (MQ=60)
                                tctacacgaccgccatcgtcggcagcgtcagatgtgtataagagacagCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAA                                                                     <  SRR3722078.486201/52‑1 (MQ=60)
                                TACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAA                                                                     >  SRR3722078.1149396/1‑100 (MQ=60)
                                TACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAA                                                                     <  SRR3722078.295847/100‑1 (MQ=60)
                                   gcccgccatcgtcggcagcgtcagatgtgtataagagacagAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCG                                                                  <  SRR3722078.1176398/59‑1 (MQ=60)
                                     GGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTT                                                                >  SRR3722078.1454000/1‑100 (MQ=60)
                                     GGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTT                                                                >  SRR3722078.1099182/1‑100 (MQ=60)
                                          GCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTctgtctcttatacacat                                                           >  SRR3722078.1236678/1‑83 (MQ=60)
                                          GCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTG                                                           >  SRR3722078.592518/1‑100 (MQ=60)
                                               GGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACctgtctcttatacacatctgacgctgccgacg                                                      >  SRR3722078.48274/1‑68 (MQ=60)
                                               GGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAC                                                      >  SRR3722078.334970/1‑100 (MQ=60)
                                               GGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGGGCAAAC                                                      >  SRR3722078.962331/1‑100 (MQ=60)
                                               GGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGCATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAC                                                      >  SRR3722078.453392/1‑100 (MQ=60)
                                                acagATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACC                                                     <  SRR3722078.134602/96‑1 (MQ=60)
                                                GCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACC                                                     <  SRR3722078.427060/100‑1 (MQ=60)
                                                  CTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTG                                                   >  SRR3722078.1046428/1‑100 (MQ=60)
                                                    ATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGT                                                 >  SRR3722078.575745/1‑100 (MQ=60)
                                                      GTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGC                                               <  SRR3722078.394596/100‑1 (MQ=60)
                                                                CTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTTGGTGCAAACCTGGTGCTGACGGAAGG                                     <  SRR3722078.1440900/100‑1 (MQ=60)
                                                                CTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGG                                     <  SRR3722078.50173/100‑1 (MQ=60)
                                                                CTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGG                                     <  SRR3722078.63076/100‑1 (MQ=60)
                                                                        TACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCctgtctcttatacacatctgacgc                             >  SRR3722078.390150/1‑76 (MQ=60)
                                                                        TACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAG                             <  SRR3722078.980251/100‑1 (MQ=60)
                                                                        TACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAACCTGGTGCTGACGGAAGGTGCTAAAGG                            <  SRR3722078.936719/100‑1 (MQ=60)
                                                                            AACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAT                         >  SRR3722078.1247253/1‑100 (MQ=60)
                                                                            AACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAT                         >  SRR3722078.957872/1‑100 (MQ=60)
                                                                                CACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAA                     >  SRR3722078.1508556/1‑100 (MQ=60)
                                                                                CACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAA                     <  SRR3722078.770108/100‑1 (MQ=60)
                                                                                CACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAA                     <  SRR3722078.800913/100‑1 (MQ=60)
                                                                                  CCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGG                   >  SRR3722078.401591/1‑100 (MQ=60)
                                                                                          ATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCC           <  SRR3722078.930263/100‑1 (MQ=60)
                                                                                          ATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCC           <  SRR3722078.1020060/100‑1 (MQ=60)
                                                                                          ATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCC           <  SRR3722078.995098/100‑1 (MQ=60)
                                                                                           TGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCA          >  SRR3722078.1425527/1‑100 (MQ=60)
                                                                                            GCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCAA         >  SRR3722078.1499794/1‑100 (MQ=60)
                                                                                                  AGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCAAAAAGCA   >  SRR3722078.210218/1‑100 (MQ=60)
                                                                                                   GCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGctgtctcttatacacatctgacgctgc  >  SRR3722078.775699/1‑73 (MQ=60)
                                                                                                   GCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGctgtctcttatacacatctgacgctgc  >  SRR3722078.112757/1‑73 (MQ=60)
                                                                                                   |                                                                                                   
GGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGCAATCCAAAAAGCAG  >  NZ_CP009273/2525951‑2526149

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