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,995,567 A→G T635A (ACC→GCC)  xdhA → xanthine dehydrogenase molybdenum‑binding subunit XdhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,995,5670AG96.8% 192.1 / ‑2.0 62T635A (ACC→GCC) xdhAxanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (29/31);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAAT  >  NZ_CP009273/2995480‑2995649
                                                                                       |                                                                                  
gAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGcc                                                                                  >  2:488641/1‑90 (MQ=255)
     cTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCaa                                                                             >  1:1562646/1‑90 (MQ=255)
         aTCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGc                                                                         <  2:1024037/90‑1 (MQ=255)
          tCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCa                                                                        <  2:1225139/90‑1 (MQ=255)
               accaccACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTc                                                                   <  2:769553/90‑1 (MQ=255)
                 caccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCaa                                                                 <  2:1467762/90‑1 (MQ=255)
                 caccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCaa                                                                 <  2:811247/90‑1 (MQ=255)
                 caccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCaa                                                                 <  2:833704/90‑1 (MQ=255)
                   ccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACg                                                               >  2:179330/1‑90 (MQ=255)
                   ccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACg                                                               <  1:422274/90‑1 (MQ=255)
                        aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  2:346046/90‑1 (MQ=255)
                        aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  1:203295/90‑1 (MQ=255)
                        aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  1:910876/90‑1 (MQ=255)
                              ccGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCa                                                    <  1:292793/90‑1 (MQ=255)
                                  cGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGc                                                <  1:179330/90‑1 (MQ=255)
                                   gTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCa                                               <  2:611571/90‑1 (MQ=255)
                                    tttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCat                                              <  1:1376534/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  2:867358/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  1:1005660/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  2:126010/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  1:780144/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  1:404050/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  2:1177242/90‑1 (MQ=255)
                                           gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTa                                       >  2:416739/1‑90 (MQ=255)
                                           gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTa                                       >  2:262343/1‑90 (MQ=255)
                                               cTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATcc                                   <  2:38242/90‑1 (MQ=255)
                                                tttGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCa                                  >  1:146034/1‑90 (MQ=255)
                                                tttGTTGATCTGACGATCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCa                                  >  1:65229/1‑90 (MQ=255)
                                                 ttgttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCAc                                 >  2:1430/1‑90 (MQ=255)
                                                 ttgttgATCTGACGGTCGATATTGCGCTGTGCAAAGCCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCAc                                 >  2:590196/1‑90 (MQ=255)
                                                      gATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCActgctg                            >  2:45721/1‑90 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  2:678386/58‑1 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          >  1:678386/1‑58 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  1:342914/71‑1 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             >  2:342914/1‑71 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGc                          >  2:621256/1‑90 (MQ=255)
                                                           gACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGa                       >  2:1383327/1‑90 (MQ=255)
                                                            aCGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGaa                      <  2:270662/90‑1 (MQ=255)
                                                             cGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAg                     >  2:1023541/1‑90 (MQ=255)
                                                              ggTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAgg                    <  1:874985/90‑1 (MQ=255)
                                                                    tatTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGt              <  1:45721/90‑1 (MQ=255)
                                                                        gcgcTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTacac          >  2:283047/1‑90 (MQ=255)
                                                                           gtgtgCCAAGTCGCCATCGACCGCACCCCCAACGTTCATGATTCAGGCCATCTTCTTAATCCACTGCTGGCAGAAGATCAGGTACAcggc       <  1:821688/89‑1 (MQ=255)
                                                                           cTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggc       <  1:748857/90‑1 (MQ=255)
                                                                           cTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggc       <  2:1405450/90‑1 (MQ=255)
                                                                           cTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggc       <  2:333362/90‑1 (MQ=255)
                                                                             gtgCAAAGTCGCCATCGACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:1418594/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  2:742028/37‑1 (MQ=38)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          >  1:742028/1‑37 (MQ=38)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             >  1:856412/1‑50 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  2:856412/50‑1 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:860500/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:108864/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:480016/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:374905/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:239398/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:1486175/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:590968/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:250761/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTCCTGGCAGAAGGTCAGGTACAcggcgg     >  2:108497/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACCTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:655827/1‑90 (MQ=255)
                                                                                cAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCggaat  <  2:1364884/90‑1 (MQ=255)
                                                                                       |                                                                                  
GAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAAT  >  NZ_CP009273/2995480‑2995649

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
CTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGC  >  NZ_CP009273/2995475‑2995665
                                                                                            |                                                                                                  
CTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAA                                                                                             >  SRR3722078.1588081/1‑100 (MQ=60)
 TGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAAC                                                                                            <  SRR3722078.1318221/100‑1 (MQ=60)
                        CCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTC                                                                     <  SRR3722078.427366/100‑1 (MQ=60)
                             AACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGC                                                                <  SRR3722078.923172/100‑1 (MQ=60)
                             AACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGC                                                                <  SRR3722078.205747/100‑1 (MQ=60)
                                   CCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTC                                                          <  SRR3722078.296361/100‑1 (MQ=60)
                                       CGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAA                                                      <  SRR3722078.181447/100‑1 (MQ=60)
                                         TTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATC                                                    <  SRR3722078.1398731/100‑1 (MQ=60)
                                          aagagacagCCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCC                                                   <  SRR3722078.347108/91‑1 (MQ=60)
                                          TTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCC                                                   <  SRR3722078.1020401/100‑1 (MQ=60)
                                          TTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCC                                                   <  SRR3722078.790273/100‑1 (MQ=60)
                                          TTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCC                                                   <  SRR3722078.408951/100‑1 (MQ=60)
                                           TGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCA                                                  >  SRR3722078.147730/1‑100 (MQ=60)
                                           TGGCTGTACCTTTGTTGATCTGACGATCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCA                                                  >  SRR3722078.65993/1‑100 (MQ=60)
                                                   CCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCctgtctcttatacacatctgac                                          >  SRR3722078.686247/1‑78 (MQ=60)
                                                                 atgtgtataagagacagGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTAC                            <  SRR3722078.294662/83‑1 (MQ=60)
                                                                   GGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAC                          <  SRR3722078.886730/100‑1 (MQ=60)
                                                                        ATATTGCGCTGTGCAAAGTCGCCATCGACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG                     >  SRR3722078.1441633/1‑100 (MQ=60)
                                                                        ATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCctgtctcttatacacatctgacgctgccgacgatggcggtcgt                     >  SRR3722078.751413/1‑57 (MQ=60)
                                                                        ATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCctgtctcttatacacatctgacgctgccga                     >  SRR3722078.867860/1‑70 (MQ=60)
                                                                        ATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG                     >  SRR3722078.597832/1‑100 (MQ=60)
                                                                        ATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG                     >  SRR3722078.253820/1‑100 (MQ=60)
                                                                        ATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG                     >  SRR3722078.110150/1‑100 (MQ=60)
                                                                        ATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACCTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG                     >  SRR3722078.663338/1‑100 (MQ=60)
                                                                         TATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGA                    <  SRR3722078.46260/100‑1 (MQ=60)
                                                                                agGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAA             <  SRR3722078.584586/98‑1 (MQ=60)
                                                                                GTGTGCCAAGTCGCCATCGACCGCACCCCCAACGTTCATGATTCAGGCCATCTTCTTAATCCACTGCTGGCAGAAGATCAGGTACACGGCGGAATGGGAC             <  SRR3722078.832560/100‑1 (MQ=60)
                                                                                CTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAA             <  SRR3722078.758397/100‑1 (MQ=60)
                                                                                  GTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATG           >  SRR3722078.434712/1‑100 (MQ=60)
                                                                                  GTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATG           >  SRR3722078.72082/1‑100 (MQ=60)
                                                                                  GTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATG           >  SRR3722078.879826/1‑100 (MQ=60)
                                                                                  GTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATG           >  SRR3722078.1436556/1‑100 (MQ=60)
                                                                                           CGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGC  <  SRR3722078.1038704/100‑1 (MQ=60)
                                                                                           CGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGC  <  SRR3722078.242317/100‑1 (MQ=60)
                                                                                            |                                                                                                  
CTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGC  >  NZ_CP009273/2995475‑2995665

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