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 F5 I97 R1 275 44.0 1733530 96.5% 1672856 122.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 668,149 A→G N697N (AAT→AAC leuS ← leucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273668,1490AG95.9% 170.4 / ‑0.8 49N697N (AAT→AACleuSleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (23/24);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CATCGGTTGGTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGT  >  NZ_CP009273/668013‑668277
                                                                                                                                        |                                                                                                                                
cATCGGTTGGTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGtt                                                                                                                                >  1:373330/1‑139 (MQ=255)
     gTTGGTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTtcag                                                                                                                           >  2:785/1‑139 (MQ=255)
      ttGGTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTtcagt                                                                                                                          <  2:710322/139‑1 (MQ=255)
            cTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAgcg                                                                                                                     <  2:518173/138‑1 (MQ=255)
               tCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCa                                                                                                                   <  1:307699/137‑1 (MQ=255)
                 gCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGTGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCa                                                                                                                <  1:577319/138‑1 (MQ=255)
                 gCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCa                                                                                                                <  1:541104/138‑1 (MQ=255)
                    aGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCa                                                                                                                <  2:513313/135‑1 (MQ=255)
                    aGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCa                                                                                                                >  1:513313/1‑135 (MQ=255)
                         gTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAg                                                                                                        <  1:233489/138‑1 (MQ=255)
                             atcaGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCg                                                                                                   <  1:38654/139‑1 (MQ=255)
                               caGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGc                                                                                                  >  1:687879/1‑138 (MQ=255)
                                  cTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACa                                                                                              <  2:352370/139‑1 (MQ=255)
                                  cTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACa                                                                                              >  2:686673/1‑139 (MQ=255)
                                  cTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACa                                                                                              >  1:736294/1‑139 (MQ=255)
                                     cATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATc                                                                                            <  1:552931/138‑1 (MQ=255)
                                             ccgcAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGc                                                                                   >  2:471310/1‑139 (MQ=255)
                                               gcAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCtg                                                                                 >  1:652335/1‑139 (MQ=255)
                                                              aaGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAAcc                                                                   <  1:438868/138‑1 (MQ=255)
                                                               aGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAg                                                                 <  2:652335/139‑1 (MQ=255)
                                                                     gacgacGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTc                                                             >  2:621654/1‑137 (MQ=255)
                                                                     gacgacGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTc                                                             >  2:10286/1‑137 (MQ=255)
                                                                           ggCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACg                                                     >  1:391810/1‑139 (MQ=255)
                                                                            gCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCtgtg                                                                               <  1:606677/112‑1 (MQ=255)
                                                                            gCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCtgtg                                                                               >  2:606677/1‑112 (MQ=255)
                                                                            gCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGt                                                    >  2:304572/1‑139 (MQ=255)
                                                                                atatCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGttt                                                  >  2:179311/1‑137 (MQ=255)
                                                                                   tcatcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGa                                             <  2:51728/139‑1 (MQ=255)
                                                                                     atcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGaa                                            <  1:686673/138‑1 (MQ=255)
                                                                                     atcGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGaa                                            <  1:302928/138‑1 (MQ=255)
                                                                                            aCTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAg                                    >  1:26274/1‑139 (MQ=255)
                                                                                             cTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGc                                   <  2:550447/139‑1 (MQ=255)
                                                                                             cTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGc                                   <  1:561493/139‑1 (MQ=255)
                                                                                                        gTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCacac                        >  1:348423/1‑139 (MQ=255)
                                                                                                        gTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCa                                                                                                                <  2:680320/51‑1 (MQ=255)
                                                                                                        gTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCa                                                                                                                >  1:680320/1‑51 (MQ=255)
                                                                                                                   aTCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTg              <  2:427848/138‑1 (MQ=255)
                                                                                                                   aTCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTg              >  1:427848/1‑138 (MQ=255)
                                                                                                                   aTCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGc             >  2:182714/1‑139 (MQ=255)
                                                                                                                    tCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGc                                   >  2:93787/1‑116 (MQ=255)
                                                                                                                    tCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGc                                   <  1:93787/116‑1 (MQ=255)
                                                                                                                     cgcgACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCa           >  2:228761/1‑139 (MQ=255)
                                                                                                                     cgcgACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCa           >  1:521320/1‑139 (MQ=255)
                                                                                                                      gcgACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCTCCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCa           <  1:228761/138‑1 (MQ=255)
                                                                                                                        gACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTc        <  2:553428/139‑1 (MQ=255)
                                                                                                                         aCGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCg       <  1:75769/139‑1 (MQ=255)
                                                                                                                               cgcTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTcc                                                            >  2:686427/1‑80 (MQ=255)
                                                                                                                               cgcTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTcc                                                            <  1:686427/80‑1 (MQ=255)
                                                                                                                               cgcTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGt  >  2:499872/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                
CATCGGTTGGTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGT  >  NZ_CP009273/668013‑668277

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATC  >  NZ_CP009273/668002‑668283
                                                                                                                                                   |                                                                                                                                      
ATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTT                                                                                                                                     >  SRR3722120.378841/1‑151 (MQ=60)
                    GTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTT                                                                                                                 >  SRR3722120.520633/1‑151 (MQ=60)
                         TTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAG                                                                                                              <  SRR3722120.312161/149‑1 (MQ=60)
                           CGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGTGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGC                                                                                                           <  SRR3722120.585354/150‑1 (MQ=60)
                           CGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGC                                                                                                           <  SRR3722120.548741/150‑1 (MQ=60)
                               AGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCA                                                                                                       >  SRR3722120.697298/1‑150 (MQ=60)
                                  TTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACAT                                                                                                   >  SRR3722120.746254/1‑151 (MQ=60)
                                   TGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACAT                                                                                                   <  SRR3722120.236795/150‑1 (MQ=60)
                                       CATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCT                                                                                              <  SRR3722120.39126/151‑1 (MQ=60)
                                               CCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGT                                                                                      >  SRR3722120.661316/1‑151 (MQ=60)
                                               CCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTG                                                                                       <  SRR3722120.560704/150‑1 (MQ=60)
                                                           atgtgtataagagacAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAAC                                                                          <  SRR3722120.615098/136‑1 (MQ=60)
                                                                        GAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGAC                                                              <  SRR3722120.445286/150‑1 (MQ=60)
                                                                           GGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGT                                                          >  SRR3722120.397549/1‑151 (MQ=60)
                                                                              acgagattcgccttagtctcgtgggctcggagatgtgtataagagacagTCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTC                                                       <  SRR3722120.695826/102‑1 (MQ=60)
                                                                                            TATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGC                                         >  SRR3722120.26571/1‑151 (MQ=60)
                                                                                               CATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCC                                       <  SRR3722120.696073/150‑1 (MQ=60)
                                                                                               CATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCC                                       <  SRR3722120.307268/150‑1 (MQ=60)
                                                                                                       gtataagagacaGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACAC                              <  SRR3722120.95039/139‑1 (MQ=60)
                                                                                                       ACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACAC                              <  SRR3722120.569354/151‑1 (MQ=60)
                                                                                                        CTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACC                             >  SRR3722120.353463/1‑151 (MQ=60)
                                                                                                        CTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaa                              >  SRR3722120.689659/1‑74 (MQ=60)
                                                                                                                   GTTTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCC                  >  SRR3722120.434094/1‑151 (MQ=60)
                                                                                                                     TTTATGCACATCGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCAT                >  SRR3722120.528741/1‑151 (MQ=60)
                                                                                                                                CGCGACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCTCCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCA      <  SRR3722120.231977/150‑1 (MQ=60)
                                                                                                                                   GACGCAGCGCTTTCTGGTTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATC  <  SRR3722120.76763/151‑1 (MQ=60)
                                                                                                                                                   |                                                                                                                                      
ATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCATCAGCTCCATAATCGCCGCAATTGCGGTGTTGAAGGTCTGACGACGGCCGATATCATCGGTCACTTTAGCGATCGTTTTATGCACATCGCGACGCAGCGCTTTCTGATTTTCAGTCAGCGCATCAACGTTCAGTGCCGCAACATCACCTTTTGCTGTGTGCTCGTAAACCAGTTTCCAGACACGTTTCAGGAAGCGGTTAGCCCCTTCCACACCGGATTCCTGCCATTCGAGAGTCATATC  >  NZ_CP009273/668002‑668283

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