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,439,101 T→C K255E (AAG→GAG)  mepA ← penicillin‑insensitive murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,439,1010TC97.7% 133.3 / ‑3.6 44K255E (AAG→GAG) mepApenicillin‑insensitive murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (31/12);  total (32/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAA  >  NZ_CP009273/2439019‑2439183
                                                                                  |                                                                                  
cAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGc                                                                             <  1:288837/90‑1 (MQ=255)
 aGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCt                                                                            >  1:1086045/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCtt                                                                           >  2:483589/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  2:836676/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  2:821261/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  2:815648/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  2:188718/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  2:1210730/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  2:1041526/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  1:627225/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  1:1361582/1‑90 (MQ=255)
  ggCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCtt                                                                           >  1:1158293/1‑90 (MQ=255)
   gCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCttt                                                                          >  2:260027/1‑90 (MQ=255)
            cGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTccc                                                                 <  2:734597/90‑1 (MQ=255)
              tttCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCgg                                                               >  1:245030/1‑90 (MQ=255)
              tttCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCgg                                                               >  1:1125008/1‑90 (MQ=255)
              tttCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCgg                                                               >  2:30979/1‑90 (MQ=255)
                tCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGtt                                                             >  1:1302561/1‑90 (MQ=255)
                tCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGtt                                                             >  1:1158634/1‑90 (MQ=255)
                tCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGtt                                                             >  2:1543410/1‑90 (MQ=255)
                           aCGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                                  >  1:98876/1‑90 (MQ=255)
                           aCGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                                  >  1:403733/1‑90 (MQ=255)
                           aCGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                                  >  2:1335501/1‑90 (MQ=255)
                           aCGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                                  >  1:64201/1‑90 (MQ=255)
                           aCGTGCTCATCCAGTAGCGCCTGCCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCa                                                  >  2:1330619/1‑90 (MQ=255)
                            cGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCaa                                                 >  2:1113570/1‑90 (MQ=255)
                            cGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCaa                                                 >  1:178356/1‑90 (MQ=255)
                                  cATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGc                                           <  1:1337654/90‑1 (MQ=255)
                                   aTCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCt                                          >  1:1284730/1‑90 (MQ=255)
                                             gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCt                                >  1:1064680/1‑90 (MQ=255)
                                             gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCt                                >  2:777011/1‑90 (MQ=255)
                                             gcCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCt                                >  1:571626/1‑90 (MQ=255)
                                                tGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGcc                             >  2:1537557/1‑90 (MQ=255)
                                                    aGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGc                         <  1:1330619/90‑1 (MQ=255)
                                                     ggAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCa                        >  1:384621/1‑90 (MQ=255)
                                                          gcggcAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCAt                   <  2:827425/90‑1 (MQ=255)
                                                           cggcAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATc                  <  2:41317/90‑1 (MQ=255)
                                                           cggcAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATc                  <  1:324093/90‑1 (MQ=255)
                                                                   gAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATg          <  1:459681/90‑1 (MQ=255)
                                                                    aGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATgg         <  2:104470/90‑1 (MQ=255)
                                                                    aGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATgg         <  1:836676/90‑1 (MQ=255)
                                                                     ggCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGc        <  2:1086045/90‑1 (MQ=255)
                                                                       cGGTGTCTTCTCCTCAGGCTTTTTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCgg      <  1:926038/90‑1 (MQ=255)
                                                                           gtCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTaaa  >  2:131382/1‑90 (MQ=255)
                                                                                  |                                                                                  
CAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAA  >  NZ_CP009273/2439019‑2439183

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCGGGGAAACCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACAC  >  NZ_CP009273/2439003‑2439198
                                                                                                  |                                                                                                 
tcgtcggcagcgtcagatgtgtataagagacagCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCC                                                                                                  <  SRR3722078.881076/67‑1 (MQ=60)
       AACCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCT                                                                                           >  SRR3722078.1102501/1‑100 (MQ=60)
       AACCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCT                                                                                           <  SRR3722078.767641/100‑1 (MQ=60)
        ACCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTT                                                                                          >  SRR3722078.634442/1‑100 (MQ=60)
        ACCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTT                                                                                          >  SRR3722078.1176221/1‑100 (MQ=60)
        ACCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTT                                                                                          >  SRR3722078.1383466/1‑100 (MQ=60)
                CAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTC                                                                                  <  SRR3722078.292360/100‑1 (MQ=60)
                 AGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCctgtctcttat                                                                                 >  SRR3722078.535191/1‑89 (MQ=60)
                    CTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGG                                                                              >  SRR3722078.248025/1‑100 (MQ=60)
                    CTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGG                                                                              >  SRR3722078.1142211/1‑100 (MQ=60)
                      ATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTT                                                                            >  SRR3722078.1323332/1‑100 (MQ=60)
                      ATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTT                                                                            >  SRR3722078.1176568/1‑100 (MQ=60)
                       ccgagatctacacgaccgccatcgtcggcagcgtcagatgtgtataagagacaGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTT                                                                           <  SRR3722078.243627/47‑1 (MQ=60)
                                 CCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCA                                                                 >  SRR3722078.64956/1‑100 (MQ=60)
                                 CCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCA                                                                 >  SRR3722078.100036/1‑100 (MQ=60)
                                 CCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCA                                                                 >  SRR3722078.408627/1‑100 (MQ=60)
                                  CATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAA                                                                >  SRR3722078.180443/1‑100 (MQ=60)
                                        ATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTctgtctcttatacacatctgacgctgccgacga                                                          >  SRR3722078.1217485/1‑67 (MQ=60)
                                         TCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCT                                                         >  SRR3722078.1305166/1‑100 (MQ=60)
                                                  CATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTC                                                <  SRR3722078.1359098/100‑1 (MQ=60)
                                                   ATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCT                                               >  SRR3722078.578315/1‑100 (MQ=60)
                                                   ATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCT                                               >  SRR3722078.1080720/1‑100 (MQ=60)
                                                           GCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCA                                       >  SRR3722078.389308/1‑100 (MQ=60)
                                                                    AGGAAGGCGGCCACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCC                              <  SRR3722078.1351919/100‑1 (MQ=60)
                                                                           CGGCAACGGAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGC                       <  SRR3722078.328067/100‑1 (MQ=60)
                                                                                   GAGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGG               <  SRR3722078.465175/100‑1 (MQ=60)
                                                                                    AGGCGGTGTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGT              <  SRR3722078.847809/100‑1 (MQ=60)
                                                                                       CGGTGTCTTCTCCTCAGGCTTTTTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGA           <  SRR3722078.938758/100‑1 (MQ=60)
                                                                                           GTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCctgtctcttatacacatctgacg       >  SRR3722078.1520246/1‑77 (MQ=60)
                                                                                           GTCTTCTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTT       >  SRR3722078.405691/1‑100 (MQ=60)
                                                                                                CTCCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACAG  <  SRR3722078.672340/100‑1 (MQ=60)
                                                                                                  |                                                                                                 
GCGGGGAAACCATAAACAGGCTATTAAACGTTTCCATTAGATCACGTGCTCATCCAGTAGCGCCTGGCAGGAAGGCGGCAACGGAGGCGGTGTCTTCTTCTCAGGCTTTGTTGTTCCCGGTTTTGGAGGTTCAAACCAGCTTTGCAGTTCTGCCCCGCAACCATCGCCTGATGGCGGTAAAGGTTGATCTTCACAC  >  NZ_CP009273/2439003‑2439198

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