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 F14 I0 R1 248 61.6 2601574 98.7% 2567753 146.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 1,842,575 Δ2 bp 100% coding (2218‑2219/2754 nt) gacS ← sensor protein GacS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC1,842,5750T.97.7% 193.5 / ‑1.5 44coding (2219/2754 nt)gacSsensor protein GacS
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (26/17);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_002947_CJ‑RC1,842,5760G.97.7% 190.1 / ‑2.0 44coding (2218/2754 nt)gacSsensor protein GacS
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (26/17);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CACTTCATCACCACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGC  >  NC_002947_CJ‑RC/1842434‑1842708
                                                                                                                                             ||                                                                                                                                    
cACTTCATCACCACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTtct                                                                                                                              >  6:111295/1‑149 (MQ=255)
  cTTCATCACCACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTcc                                                                                                                            <  6:159686/149‑1 (MQ=255)
      atcaCCACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCca                                                                                                                          >  8:42450/1‑147 (MQ=255)
         accaccTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGc                                                                                                                     <  7:233818/149‑1 (MQ=255)
            accTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGa                                                                                                                  >  6:282066/1‑149 (MQ=255)
                  gCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGtt                                                                                                            <  4:58217/149‑1 (MQ=255)
                    cAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCg                                                                                                          <  6:315024/149‑1 (MQ=255)
                      gTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCgg                                                                                                         >  5:33151/1‑148 (MQ=255)
                       ttGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCCGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTa                                                                                                       <  4:89183/149‑1 (MQ=255)
                         gCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTa                                                                                                       >  1:40974/1‑147 (MQ=255)
                             ttCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGgcg                                                                                                 >  6:127184/1‑149 (MQ=255)
                                  tGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGcc                                                                                            >  4:201755/1‑149 (MQ=255)
                                           tGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGcc                                                                                   >  4:100922/1‑149 (MQ=255)
                                             gTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCgg                                                                                 >  5:14754/1‑149 (MQ=255)
                                              tCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCg                                                                                  <  4:186627/149‑1 (MQ=255)
                                                  gTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATc                                                                            >  8:117895/1‑149 (MQ=255)
                                                     gtcgtcCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTg                                                                          >  4:202639/1‑148 (MQ=255)
                                                     gtcgtcCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGc                                                                         <  5:137116/149‑1 (MQ=255)
                                                      tcgtcCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGca                                                                        <  5:220811/149‑1 (MQ=255)
                                                         tcCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCTGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCAcat                                                                     <  2:140797/149‑1 (MQ=255)
                                                             tGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCAcatccat                                                                 <  3:100922/149‑1 (MQ=255)
                                                              gCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCAGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCAcatccat                                                                 >  3:247330/1‑148 (MQ=255)
                                                                  cTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCa                                                            >  4:138391/1‑149 (MQ=255)
                                                                     tgCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCa                                                         <  5:214143/149‑1 (MQ=255)
                                                                       caacaaCGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAgg                                                       >  7:139533/1‑149 (MQ=255)
                                                                        aacaacGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGt                                                      >  3:276531/1‑149 (MQ=255)
                                                                        aacaacGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGt                                                      >  3:233057/1‑149 (MQ=255)
                                                                              gcgcgTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGTGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAAc                                                 >  5:80014/1‑148 (MQ=255)
                                                                              gcgcgTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAAc                                                 >  7:332931/1‑148 (MQ=255)
                                                                                gcgTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAAt                                              >  4:66195/1‑149 (MQ=255)
                                                                                 cgTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAAtt                                             <  7:42450/149‑1 (MQ=255)
                                                                                   tttttCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAAttct                                           <  4:233057/149‑1 (MQ=255)
                                                                                     tttCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAAttcttc                                         >  4:124210/1‑149 (MQ=255)
                                                                                      ttCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTcc                                        >  6:38878/1‑149 (MQ=255)
                                                                                      ttCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTcc                                        >  4:70746/1‑149 (MQ=255)
                                                                                       tCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTcc                                        <  4:183588/148‑1 (MQ=255)
                                                                                         gTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGc                                     <  7:18120/149‑1 (MQ=255)
                                                                                                atggcatggGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCt                              >  1:114317/1‑149 (MQ=255)
                                                                                                atggcatggGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCt                              <  5:111295/149‑1 (MQ=255)
                                                                                                     atggGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTAcc                         >  7:292265/1‑149 (MQ=255)
                                                                                                                gggCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAAc               >  7:216865/1‑148 (MQ=255)
                                                                                                                     acgaTCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCAttgt         <  5:282066/149‑1 (MQ=255)
                                                                                                                     acgaTCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCAttgt         <  8:139533/149‑1 (MQ=255)
                                                                                                                             cAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGc  >  1:270996/1‑148 (MQ=255)
                                                                                                                                             ||                                                                                                                                    
CACTTCATCACCACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGC  >  NC_002947_CJ‑RC/1842434‑1842708

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CACCACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCC  >  NC_002947_CJ‑RC/1842442‑1842709
                                                                                                                                     |                                                                                                                                      
CACCACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCG                                                                                                                     <  NB501016:282:H2V7GAFX2:4:21502:4852:11374/151‑1 (MQ=60)
             AGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGT                                                                                                         >  NB501016:282:H2V7GAFX2:3:11408:15586:9092/1‑150 (MQ=60)
                TGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTAC                                                                                                       >  NB501016:282:H2V7GAFX2:1:11109:11516:7429/1‑149 (MQ=60)
                                    GGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGG                                                                                 >  NB501016:282:H2V7GAFX2:3:11404:19273:6184/1‑151 (MQ=60)
                                            AGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCA                                                                         <  NB501016:282:H2V7GAFX2:3:11606:5869:12748/151‑1 (MQ=60)
                                             GTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCAC                                                                        <  NB501016:282:H2V7GAFX2:3:21501:13810:12179/151‑1 (MQ=60)
                                                    ATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATC                                                                 <  NB501016:282:H2V7GAFX2:2:11211:19237:13048/151‑1 (MQ=60)
                                                     TGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCAGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATC                                                                 >  NB501016:282:H2V7GAFX2:2:21210:26078:6462/1‑150 (MQ=60)
                                                            GTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAG                                                         <  NB501016:282:H2V7GAFX2:3:21412:3512:2405/151‑1 (MQ=60)
                                                              GCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGT                                                       >  NB501016:282:H2V7GAFX2:4:11606:8135:12607/1‑151 (MQ=60)
                                                               CAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTC                                                      >  NB501016:282:H2V7GAFX2:2:21207:11077:3319/1‑151 (MQ=60)
                                                               CAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTC                                                      >  NB501016:282:H2V7GAFX2:2:21305:9035:9316/1‑151 (MQ=60)
                                                                     CGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGTGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACA                                                 >  NB501016:282:H2V7GAFX2:3:11506:12841:13706/1‑150 (MQ=60)
                                                                     CGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACA                                                 >  NB501016:282:H2V7GAFX2:4:21612:2968:6827/1‑150 (MQ=60)
                                                                        GCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTC                                             <  NB501016:282:H2V7GAFX2:4:11409:6958:16290/151‑1 (MQ=60)
                                                                                CGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCA                                     <  NB501016:282:H2V7GAFX2:4:11404:17483:16780/151‑1 (MQ=60)
                                                                                       CATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTG                              <  NB501016:282:H2V7GAFX2:3:11601:16269:7783/151‑1 (MQ=60)
                                                                                       CATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTG                              >  NB501016:282:H2V7GAFX2:1:11301:21662:1542/1‑151 (MQ=60)
                                                                                            CATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCG                         >  NB501016:282:H2V7GAFX2:4:21603:12302:3593/1‑151 (MQ=60)
                                                                                                       AGGGCGACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACC               >  NB501016:282:H2V7GAFX2:4:21411:25075:1465/1‑150 (MQ=60)
                                                                                                            GACGATCGGCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTT         <  NB501016:282:H2V7GAFX2:3:21603:6309:4189/151‑1 (MQ=60)
                                                                                                                    GCAGCGGGTTGCCGCTC‑‑GGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCC  >  NB501016:282:H2V7GAFX2:1:21211:23220:13004/1‑150 (MQ=60)
                                                                                                                                     |                                                                                                                                      
CACCACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCC  >  NC_002947_CJ‑RC/1842442‑1842709

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