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
A3 F28 I0 R1 51 55.9 1320120 86.3% 1139263 218.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,132,787 A→C 100% I204S (ATC→AGC)  wcaA ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC98.1% 134.3 / ‑5.3 48I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (26/21);  minor base T (1/0);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-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.

TCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGAC  >  NC_000913/2132577‑2133010
                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
tCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCAGGAGAATTTTTTCGGAGACGAGGGGGTCTGCATTTCACCGTGATTGATATACAGGCCCTGCGGCGCCtcttct                                                                                                                                                                                                                >  2:159321/1‑228 (MQ=255)
                  ttaccgcccGAGAGGAGGGTCAGAAAAGTGCGCCAGGTCAGACGTTTAGTGCGGGATTGGGAGAGGGTAAACATGTGATATTTTTTGCTGGCGCGGTCGAATTTGTTTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                                                                                                              <  2:464653/221‑1 (MQ=255)
                               gcagcGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAATccaccacca                                                                                                                                                                                 >  1:98696/1‑228 (MQ=255)
                                    gTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGCATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGGCTCTGCTTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                                                                                                            >  1:498661/1‑228 (MQ=255)
                                            aGTGCGCCCGGTCATACGTTTATTGCGGTTCTGGTATAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                                                                                                     <  1:159321/227‑1 (MQ=255)
                                                   cAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTAGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGGGATCTGCATTTCAGCGAGATTGATATGAAGGTCCTGCGTCGCCTCTTTTACTTTCCACGGTTCGCCGTCCTCACCCACCATCTCCAGGAAAATGTGGtca                                                                                                                                                             >  2:135292/1‑226 (MQ=255)
                                                          tACGTTTATTGCGGATCTGGTAGAGGGTACACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAGGAAGATGTGGTAATCCtgtg                                                                                                                                                      >  2:657895/1‑228 (MQ=255)
                                                          tACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCCCAGGAAGATGTCGGAATCCttgg                                                                                                                                                      >  2:23741/1‑225 (MQ=255)
                                                                  ttGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                                                                                                               <  1:23741/227‑1 (MQ=255)
                                                                      ggagttGGTAGTGGGTAAACAGTTGATTTTTTTTGCTGGCGCGGTCGAATTTGGCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                                                                                                          <  2:236227/223‑1 (MQ=255)
                                                                      ggATCTGGTGGAGGGTAAAAAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTGTCTCGGTGACGAGGTGATCTGCATTGCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                                                                                                          <  1:252541/228‑1 (MQ=255)
                                                                        aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTAACCGTGATTGATGTGCAGGCTCTGGGTCGCCTGTTATACTTTCCACGGTTCCTCGTACTCCACCACCACCCTCAGTCAGATGACTGAATCCTGTGCCGCTTTCATTTCg                                                                                                                                        >  2:20851/1‑228 (MQ=255)
                                                                                  gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCCTTCAGTCCggaacccacca                                                                                                                              >  2:146507/1‑219 (MQ=255)
                                                                                  gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTTGCCCGCTTTCAGTTCGGCATCGaaca                                                                                                                              >  1:158080/1‑228 (MQ=255)
                                                                                     tAAACCGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                                                                                                            <  1:20851/227‑1 (MQ=255)
                                                                                     tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCATGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGAGGTCGCAATCCTGTGCCGCTTTCAGTTCGGTTTCGaacaaa                                                                                                                            >  2:358017/1‑227 (MQ=255)
                                                                                     tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt                                                                                                                                                >  2:440559/1‑207 (MQ=255)
                                                                                      aaaCATGTGATTTTTTGTGCTGGCGCGGTCTAATTTGTTTTTGTGCTTGCGGTAAAAGGGGAAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                                                                                                          <  1:254517/228‑1 (MQ=255)
                                                                                      aaaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                                                                                                               <  1:440559/207‑1 (MQ=255)
                                                                                       aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctct                                                                                                                                                                                                                   >  2:473803/1‑138 (MQ=255)
                                                                                        aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                                                                                                                  <  1:473803/138‑1 (MQ=255)
                                                                                             tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCCCCGGAAGCTTTGCTAATCCTGTGCCGCTTTCGTTTCGGTTTCGAACAAAAACTtttt                                                                                                                    >  2:552073/1‑227 (MQ=255)
                                                                                               atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGAACAAACACCCTTTTaa                                                                                                                 >  2:553512/1‑228 (MQ=255)
                                                                                                     tttGCTGCCGCGGTCTAATTTGTTTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTGTCGGTGACGCGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCa                                                                                                           <  2:174443/228‑1 (MQ=255)
                                                                                                        ggtgtcgcgGTCGGATTTGTCTTTGGTCGTGCGGTAAAATTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                                                                                                            <  1:135292/219‑1 (MQ=255)
                                                                                                            gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCAGCTTTCAGTTCGGTATCGAGCAACAACTCTTTGAAACGCCAGGc                                                                                                        >  1:431221/1‑224 (MQ=255)
                                                                                                                     agtttGTCGTTGTGATTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGGCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc                                                                                            <  1:587177/225‑1 (MQ=255)
                                                                                                                         tGTCTTTGTGCTTGCGCTAAAAGTGGAAGGACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAGCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg                                                                                        <  1:595104/227‑1 (MQ=255)
                                                                                                                           ttttttggCTTGGGGGAAAAGGGGAAGTCCCGGGATAATTTTTGCGGGGGCGAGGTGGTCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGtt                                                                                      <  2:609435/220‑1 (MQ=255)
                                                                                                                              ttGGGCTTGCTGTAAAAGTGTGACTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTTCATTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCTACGGTTCGCCGGTCTCCACCACCATCCTCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCa                                                                                  <  1:611152/228‑1 (MQ=255)
                                                                                                                                     tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGTTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACCCCCTTTGAAAAGCCATCCCCCGGAAAAGCCATCGTTGCCAAtaaaat                                                                           >  2:527136/1‑223 (MQ=255)
                                                                                                                                        ggTAAAAGTTGAATTTCCCGGGAAATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATAGCCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGc                                                                        <  2:117214/228‑1 (MQ=255)
                                                                                                                                        ggTAAAAGTTGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGc                                                                        <  1:553512/228‑1 (MQ=255)
                                                                                                                                                    aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAACCGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACa                                                            >  2:500589/1‑228 (MQ=255)
                                                                                                                                                    aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACa                                                            >  1:172692/1‑228 (MQ=255)
                                                                                                                                                        cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCATTGATATTGCGTTTGTAGACCAggcg                                                        >  1:262420/1‑228 (MQ=255)
                                                                                                                                                                                ggTTATCTGCATTTCCCCTTGGTTGATATGCGGGCTCTGCGTCGCCTTTTCTACTTTCCACGTTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATa                                <  2:366939/228‑1 (MQ=255)
                                                                                                                                                                                    aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGACGTGAATACGGCGATTTGGGATACAGc                            >  1:171226/1‑228 (MQ=255)
                                                                                                                                                                                           tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGc                     >  1:136688/1‑228 (MQ=255)
                                                                                                                                                                                            ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGc                                                                                                        >  2:24784/1‑144 (MQ=255)
                                                                                                                                                                                            ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGGAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGGGAATACGGGGCTTTGGGATACAGCGGCaaaca                    >  2:219658/1‑225 (MQ=255)
                                                                                                                                                                                             tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGcc                                                                                                       <  1:24784/144‑1 (MQ=255)
                                                                                                                                                                                                 cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGGGAAACCGGCGATTTGGGAAACAGCGGAAAACCTGCCg               >  2:385665/1‑228 (MQ=255)
                                                                                                                                                                                                   tgattgatATGCAGGCTCTGGGGCGCCACTTCTTGTTTCCACGGTTCGCCGTACCCCCCCCCCATCCGCAGGAAGATGTCGTAAACCGGGCCCGCTTTCAGTTCGGTATCGAACAAACACGCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGt             <  2:54921/228‑1 (MQ=255)
                                                                                                                                                                                                   tgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACACCGGCAAGCTTGCCGGt             >  1:494828/1‑228 (MQ=255)
                                                                                                                                                                                                   tgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAACCACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGTTACACCGGCAAGCTTGCCg               >  1:322977/1‑226 (MQ=255)
                                                                                                                                                                                                   gttttgatATCCAGGCTCTGGGGCGCCTTTTCTACTTTCCCCGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGAGTTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGt             <  2:158080/225‑1 (MQ=255)
                                                                                                                                                                                                              cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGAc  >  1:72751/1‑228 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                               
TCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGAC  >  NC_000913/2132577‑2133010

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTA  >  NC_000913/2132544‑2133035
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
attcTCCAGCCTCTTTAGCGCCCCCGAATTCCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGGGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTTGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGT                                                                                                                                                                                                                                                   <  M01186:129:000000000‑A5M8V:1:1109:7190:14834/247‑1 (MQ=60)
                                         tcttttgccgttactcaatgacgggagcgtaagcaaagggcgcaaggtctttcggttatgGCGGATCTGGTAGAGGGTAACCAGCTGATATTTTTTGCTGGCGCGGTCGAATTGGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACC                                                                                                                                                                                                           <  M01186:129:000000000‑A5M8V:1:1114:13553:5949/190‑1 (MQ=60)
                                                     ACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAATCCACCACCATCCGCAGGAAGA                                                                                                                                                                                              >  M01186:129:000000000‑A5M8V:1:1104:5818:17799/1‑251 (MQ=60)
                                                          CCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGCATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGGCTCTGCTTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAGGATGTCA                                                                                                                                                                                         >  M01186:129:000000000‑A5M8V:1:2107:27432:16463/1‑251 (MQ=60)
                                                                 acgtGTCAGGAAAGTGCGCCCGGTCATACGTTTATTGCGGTTCTGGTATAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCC                                                                                                                                                                                   <  M01186:129:000000000‑A5M8V:1:1106:25840:24847/246‑1 (MQ=60)
                                                                                       GTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAT                                                                                                                                                             <  M01186:129:000000000‑A5M8V:1:1101:18642:22287/250‑1 (MQ=60)
                                                                                           TACGTTGATTGGGGATCTGGTGGAGGGTAAAAAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTGTCTCGGTGACGAGGTGATCTGCATTGCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAAC                                                                                                                                                        <  M01186:129:000000000‑A5M8V:1:1110:26282:22331/251‑1 (MQ=60)
                                                                                                        GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTTGCCCGCTTTCAGTTCGGCATCGAACAACCACCCTTTGA                                                                                                                                           >  M01186:129:000000000‑A5M8V:1:1106:17143:23453/1‑251 (MQ=60)
                                                                                                          TCTGGTAGAGGGTAAACCGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAA                                                                                                                                          <  M01186:129:000000000‑A5M8V:1:1101:26184:19768/250‑1 (MQ=60)
                                                                                                           TTGGTAGAGGGTAAACATGTGATTTTTTGTGCTGGCGCGGTCTAATTTGTTTTTGTGCTTGCGGTAAAAGGGGAAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAAC                                                                                                                                        <  M01186:129:000000000‑A5M8V:1:1110:18123:24743/251‑1 (MQ=60)
                                                                                                           CTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAT                                                                                                                                                             <  M01186:129:000000000‑A5M8V:1:2105:15928:12211/230‑1 (MQ=60)
                                                                                                             GGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAC                                                                                                                                                                                                                                <  M01186:129:000000000‑A5M8V:1:2106:26019:16959/161‑1 (MQ=60)
                                                                                                                             gtgatatttgtgggtgtCGCGGTCGGATTTGTCTTTGGTCGTGCGGTAAAATTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTA                                                                                                                          <  M01186:129:000000000‑A5M8V:1:1105:13400:28098/230‑1 (MQ=60)
                                                                                                                                  ATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCAGCTTTCAGTTCGGTATCGAGCAACAACTCTTTGAAACGCCAGGCCCAGGTAAAGAC                                                                                                                     >  M01186:129:000000000‑A5M8V:1:2105:17955:2448/1‑247 (MQ=60)
                                                                                                                                          CTGGCTCGGTCGAGTTTGTCGTTGTGATTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGGCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAA                                                                                                          <  M01186:129:000000000‑A5M8V:1:2111:21489:9050/250‑1 (MQ=60)
                                                                                                                                              CGTGGTCGAATTTGTCTTTGTGCTTGCGCTAAAAGTGGAAGGACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAGCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGAT                                                                                                      <  M01186:129:000000000‑A5M8V:1:2111:4898:15990/250‑1 (MQ=60)
                                                                                                                                                   TCTAATTTGTCTTTGGGCTTGCTGTAAAAGTGTGACTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTTCATTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCTACGGTTCGCCGGTCTCCACCACCATCCTCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCG                                                                                                <  M01186:129:000000000‑A5M8V:1:2112:5556:10881/251‑1 (MQ=60)
                                                                                                                                                             TTTTGTGCTTGCGGTAAAAGTTGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAAC                                                                                      <  M01186:129:000000000‑A5M8V:1:2109:28533:20609/251‑1 (MQ=60)
                                                                                                                                                                          GTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAAA                                                                         >  M01186:129:000000000‑A5M8V:1:1107:18933:14560/1‑251 (MQ=60)
                                                                                                                                                                              AAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCATTGATATTGCGTTTGTAGACCAGGCGGCGTGAAAACGG                                                                     >  M01186:129:000000000‑A5M8V:1:1111:7902:10933/1‑251 (MQ=60)
                                                                                                                                                                                            ctcggagatgtttataagcgacagGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATAC                                                        <  M01186:129:000000000‑A5M8V:1:2106:21937:14058/226‑1 (MQ=60)
                                                                                                                                                                                                          GTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGACGTGAATACGGCGATTTGGGATACAGCGGCAAcctcgcg                                         >  M01186:129:000000000‑A5M8V:1:1107:18152:13371/1‑244 (MQ=60)
                                                                                                                                                                                                                 GGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGG                                  >  M01186:129:000000000‑A5M8V:1:1106:17383:4481/1‑251 (MQ=60)
                                                                                                                                                                                                                 GGTGATCTGCATTTCACCGTGATTGATATGCAGGCT                                                                                                                                                                                                                                                         <  M01186:129:000000000‑A5M8V:1:1102:15834:19879/36‑1 (MQ=60)
                                                                                                                                                                                                                  GTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAC                                                                                                                     <  M01186:129:000000000‑A5M8V:1:1101:18031:23375/167‑1 (MQ=60)
                                                                                                                                                                                                                         GCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACACCGGCAAGCTTGCCGGTTGGGAATAaagt                          >  M01186:129:000000000‑A5M8V:1:2107:15959:13401/1‑247 (MQ=60)
                                                                                                                                                                                                                         GCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAACCACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGTTACACCGGCAAGCTTGCCGGTTGGGAATAGA                            >  M01186:129:000000000‑A5M8V:1:1114:16447:21937/1‑249 (MQ=60)
                                                                                                                                                                                                                                    TGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCA               >  M01186:129:000000000‑A5M8V:1:1103:21649:18450/1‑251 (MQ=60)
                                                                                                                                                                                                                                                 GGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGGATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGATGGGAATAGACTTTGCCCTGGCAAACGgctgtctca  >  M01186:129:000000000‑A5M8V:1:1113:19980:13106/1‑242 (MQ=60)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        
ATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTA  >  NC_000913/2132544‑2133035

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