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,396,099 G→C A413G (GCC→GGC)  nuoC ← NADH‑quinone oxidoreductase subunit C/D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,396,0990GC96.0% 158.7 / ‑0.7 50A413G (GCC→GGC) nuoCNADH‑quinone oxidoreductase subunit C/D
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (27/21);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCA  >  NZ_CP009273/2396013‑2396181
                                                                                      |                                                                                  
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:856268/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:1176102/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:95634/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:483688/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:227604/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:462457/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:580064/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:737050/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTCCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:1457116/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGGAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:693177/1‑90 (MQ=255)
    cTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAAc                                                                             <  1:1121416/90‑1 (MQ=255)
           aCGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTgg                                                                      >  1:404859/1‑90 (MQ=255)
             gTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGa                                                                    >  2:272247/1‑90 (MQ=255)
               cGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAAc                                                                  <  1:462457/90‑1 (MQ=255)
               cGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAAc                                                                  <  1:201322/90‑1 (MQ=255)
               cGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAAc                                                                  <  1:580064/90‑1 (MQ=255)
                  aGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGAc                                                               <  2:1237438/90‑1 (MQ=255)
                     cGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCtt                                                            <  2:1146324/90‑1 (MQ=255)
                      gATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCttt                                                           >  1:458236/1‑90 (MQ=255)
                       aTCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTc                                                          <  2:1360744/90‑1 (MQ=255)
                       aTCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTc                                                          <  2:759090/90‑1 (MQ=255)
                       aTCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTc                                                          <  2:1426144/90‑1 (MQ=255)
                        tCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCa                                                         >  1:1128734/1‑90 (MQ=255)
                        tCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCa                                                         >  2:1014865/1‑90 (MQ=255)
                            ggTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                     <  1:1417044/90‑1 (MQ=255)
                            ggTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                     <  2:762890/90‑1 (MQ=255)
                            ggTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                     <  1:62267/90‑1 (MQ=255)
                              tAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATgg                                                   <  1:834715/90‑1 (MQ=255)
                                  aCGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGtt                                               >  1:21256/1‑90 (MQ=255)
                                     cAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTg                                            >  1:76234/1‑90 (MQ=255)
                                     cAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTg                                            >  2:443696/1‑90 (MQ=255)
                                     cAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTg                                            >  1:1286818/1‑90 (MQ=255)
                                      aGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGc                                           <  1:869145/90‑1 (MQ=255)
                                      aGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGc                                           >  2:443736/1‑90 (MQ=255)
                                          ccGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAgcg                                       >  1:46687/1‑90 (MQ=255)
                                            gcgcCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGcgcc                                     >  1:564342/1‑90 (MQ=255)
                                            gcgcCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGcgcc                                     >  2:1366380/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCttt                                >  2:984500/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCttt                                >  1:1264867/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCttt                                >  1:1518363/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTGTGCAGCGCCGCttt                                >  1:1217955/1‑90 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          >  1:992755/1‑90 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          >  2:87066/1‑90 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          <  1:1095787/90‑1 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          >  1:1118628/1‑90 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          >  1:315899/1‑90 (MQ=255)
                                                        ccccACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTaa                         <  1:242150/90‑1 (MQ=255)
                                                            aCTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACg                     >  2:610369/1‑90 (MQ=255)
                                                            aCTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACg                     >  1:18482/1‑90 (MQ=255)
                                                              tCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGcc                   >  2:1106037/1‑90 (MQ=255)
                                                               ccAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCa                  <  1:1094655/90‑1 (MQ=255)
                                                                cAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAg                 <  2:227604/90‑1 (MQ=255)
                                                                       tctTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTc          >  2:796142/1‑90 (MQ=255)
                                                                       tctTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTc          >  1:777521/1‑90 (MQ=255)
                                                                         tttCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCgg        >  2:384803/1‑90 (MQ=255)
                                                                         tttCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCgg        >  2:1011933/1‑90 (MQ=255)
                                                                             gcgcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATc    >  1:285156/1‑90 (MQ=255)
                                                                              cgcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATcc   <  2:21256/90‑1 (MQ=255)
                                                                              cgcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATcc   <  2:1256882/90‑1 (MQ=255)
                                                                               gcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCa  <  1:507936/90‑1 (MQ=255)
                                                                                      |                                                                                  
GCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCA  >  NZ_CP009273/2396013‑2396181

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTCGAGGAACTCACGC  >  NZ_CP009273/2396003‑2396197
                                                                                                |                                                                                                  
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGG                                                                                                 >  SRR3722078.1070286/1‑100 (MQ=60)
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGG                                                                                                 >  SRR3722078.96754/1‑100 (MQ=60)
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGG                                                                                                 >  SRR3722078.489429/1‑100 (MQ=60)
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGG                                                                                                 >  SRR3722078.230366/1‑100 (MQ=60)
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGG                                                                                                 >  SRR3722078.1308922/1‑100 (MQ=60)
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGGAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGG                                                                                                 >  SRR3722078.701416/1‑100 (MQ=60)
   gagatctacacgaccgccatcgtcggcagcgtcagatgtgtataagagacaGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAA                                                                                              <  SRR3722078.1533374/49‑1 (MQ=60)
           CGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGG                                                                                      >  SRR3722078.409771/1‑100 (MQ=60)
              CTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAA                                                                                   <  SRR3722078.1138533/100‑1 (MQ=60)
                      CGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTT                                                                           >  SRR3722078.463709/1‑100 (MQ=60)
                        TCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCA                                                                         >  SRR3722078.1146005/1‑100 (MQ=60)
                         CGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAG                                                                        <  SRR3722078.467979/100‑1 (MQ=60)
                         CGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAG                                                                        <  SRR3722078.586828/100‑1 (MQ=60)
                         CGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAG                                                                        <  SRR3722078.203745/100‑1 (MQ=60)
                                  TCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTT                                                               >  SRR3722078.21526/1‑100 (MQ=60)
                                     CGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTG                                                            >  SRR3722078.77130/1‑100 (MQ=60)
                                     CGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTG                                                            >  SRR3722078.1307301/1‑100 (MQ=60)
                                      GGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGC                                                           <  SRR3722078.1440060/100‑1 (MQ=60)
                                      GGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGC                                                           <  SRR3722078.62996/100‑1 (MQ=60)
                                        TAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAG                                                         <  SRR3722078.845818/100‑1 (MQ=60)
                                          GCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCG                                                       >  SRR3722078.47239/1‑100 (MQ=60)
                                            ACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCC                                                     >  SRR3722078.570960/1‑100 (MQ=60)
                                                AGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTT                                                 <  SRR3722078.880784/100‑1 (MQ=60)
                                                 GGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTT                                                >  SRR3722078.1543127/1‑100 (MQ=60)
                                                 GGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTT                                                >  SRR3722078.1284936/1‑100 (MQ=60)
                                                 GGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTGTGCAGCGCCGCTTT                                                >  SRR3722078.1237108/1‑100 (MQ=60)
                                                       CGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTA                                          >  SRR3722078.1007181/1‑100 (MQ=60)
                                                       CGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTA                                          >  SRR3722078.319762/1‑100 (MQ=60)
                                                       CGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTA                                          >  SRR3722078.1135708/1‑100 (MQ=60)
                                                            GTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACG                                     >  SRR3722078.18709/1‑100 (MQ=60)
                                                                 GCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGA                                <  SRR3722078.1112448/100‑1 (MQ=60)
                                                                  CCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGAC                               <  SRR3722078.245110/100‑1 (MQ=60)
                                                                       CTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTC                          >  SRR3722078.787601/1‑100 (MQ=60)
                                                                         CCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGG                        <  SRR3722078.1111293/100‑1 (MQ=60)
                                                                             CGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATC                    >  SRR3722078.288633/1‑100 (MQ=60)
                                                                                         GCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTCGAGGAAC        <  SRR3722078.513945/100‑1 (MQ=60)
                                                                                               GCCGGCAACGCCCTGGGACCGACCTTTCAGAATGGTGTTTTGCAGCGTCGCTTTCTCGTAAGCCGCCAGTCGTTTCGGCATCCAGTCGAGCAACTCACGC  <  SRR3722078.1590997/100‑1 (MQ=60)
                                                                                               GCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTCGAGGAACTCACGC  <  SRR3722078.204694/100‑1 (MQ=60)
                                                                                               GCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTCGAGGAACTCACGC  <  SRR3722078.1122909/100‑1 (MQ=60)
                                                                                                |                                                                                                  
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTCGAGGAACTCACGC  >  NZ_CP009273/2396003‑2396197

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