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 F36 I1 R1 10 110.6 3657964 98.0% 3584804 143.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,188,953 G→A E1202K (GAA→AAA)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,188,9530GA98.4% 419.1 / ‑3.2 124E1202K (GAA→AAA) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (66/56);  total (67/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGA  >  NC_000913/4188808‑4189098
                                                                                                                                                 |                                                                                                                                                 
ggCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAAc                                                                                                                                                >  3:151286/1‑149 (MQ=255)
ggCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAAc                                                                                                                                                >  7:263729/1‑149 (MQ=255)
   aaTCGTGGCTGAAATGAGCGGTTTCGTTTCCTTCGGTAAAGAAACCTAAGGTATAGGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTAGGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACgtg                                                                                                                                             <  1:185181/149‑1 (MQ=255)
   aaTCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGAGTCCGAAATGGCGGCCGCTCAACGGGGGCGAAGGTAACCgtg                                                                                                                                             >  6:254343/1‑149 (MQ=255)
    aTCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACgtgt                                                                                                                                            <  5:218445/149‑1 (MQ=255)
    aTCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACgtgt                                                                                                                                            <  5:358525/149‑1 (MQ=255)
    aTCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACgtg                                                                                                                                             >  6:42895/1‑148 (MQ=255)
     tCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTa                                                                                                                                           <  6:97092/149‑1 (MQ=255)
       cTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAg                                                                                                                                          >  3:158202/1‑148 (MQ=255)
        tGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGaa                                                                                                                                        >  5:197098/1‑149 (MQ=255)
        tGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGaa                                                                                                                                        >  4:129360/1‑149 (MQ=255)
            tGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCCCCCCGGTAGACAGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtg                                                                                                                                    >  1:363696/1‑149 (MQ=255)
            tGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtg                                                                                                                                    >  2:324270/1‑149 (MQ=255)
             gAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGGCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTGCGAAGGTAAACGTGTAGAACgtgg                                                                                                                                   >  6:279204/1‑149 (MQ=255)
              aaaTCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtggt                                                                                                                                  >  8:241257/1‑149 (MQ=255)
              aaaTCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtg                                                                                                                                    <  4:151286/147‑1 (MQ=255)
              aaaTCAGCGGTAATGTTTCCTTCGGTAAATAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATACGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACgtgg                                                                                                                                   <  2:78366/148‑1 (MQ=255)
                 tCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGATATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGAc                                                                                                                               >  1:359738/1‑149 (MQ=255)
                      ggaaTCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCACGGTAGACCGTAGCGATCCGTGCGAAGAGATGATTCCGAAATGGCTTCACATCAACGTGTTCTAAGGTAAAAGTCTAGAACGTGCTGAAGTaaa                                                                                                                          >  8:207196/4‑148 (MQ=255)
                      ggTATCGTTTCCGTCGGTAAAGAAACCAAAGGGTCACGTCGGCGGGTTATGTCCCCGGGAGCCGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTCTCGTAAt                                                                                                                          >  3:382285/1‑149 (MQ=255)
                          tCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTcc                                                                                                                      <  1:324270/149‑1 (MQ=255)
                          tCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTcc                                                                                                                      >  7:297858/1‑149 (MQ=255)
                                  ttCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCGCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTcc                                                                                                              >  6:315364/1‑149 (MQ=255)
                                  ttCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTcc                                                                                                              >  7:338647/1‑149 (MQ=255)
                                    cGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTcc                                                                                                              <  4:368982/147‑1 (MQ=255)
                                       tAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAg                                                                                                         >  4:284875/1‑149 (MQ=255)
                                       tAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAg                                                                                                         >  8:190240/1‑149 (MQ=255)
                                              aCCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTTGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCa                                                                                                  <  6:92688/149‑1 (MQ=255)
                                               ccAAAGTTAAACGTCGTCTGGTTATCAACCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGCTGACGTAATTTCCGACGGTCAGGAAGCTCAGCAc                                                                                                 >  7:456435/1‑149 (MQ=255)
                                                cAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCAc                                                                                                 >  7:344933/1‑148 (MQ=255)
                                                cAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCAc                                                                                                 >  5:279514/1‑148 (MQ=255)
                                                  aaGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACATGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCg                                                                                                             >  7:325694/1‑134 (MQ=255)
                                                   aGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACa                                                                                             <  3:129360/149‑1 (MQ=255)
                                                     gTAAACGTCGTCTGGTTATCACCCGGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAtt                                                                                           <  5:267620/149‑1 (MQ=255)
                                                     gTAAACGTCGTCTGGTTATCACCCCGGTCGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAtt                                                                                           <  3:8955/149‑1 (MQ=255)
                                                     gTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAtt                                                                                           <  3:2727/149‑1 (MQ=255)
                                                       aaaCGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACTAGGTGAACTAATTTACAAATATACGAAAGCGCCGCACGACATtct                                                                                         >  6:166976/1‑149 (MQ=255)
                                                       aaaCGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtct                                                                                         >  4:97686/1‑149 (MQ=255)
                                                         aCGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgc                                                                                       <  4:103019/149‑1 (MQ=255)
                                                           gtcgtcTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcgt                                                                                     >  8:244037/1‑149 (MQ=255)
                                                              gtcTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcgtct                                                                                   <  4:138884/148‑1 (MQ=255)
                                                               tcTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcgtctgc                                                                                 >  7:34903/1‑149 (MQ=255)
                                                               tcTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcgtctgc                                                                                 >  2:425768/1‑149 (MQ=255)
                                                                 tGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTctgcgtctgcgt                                                                               >  7:154419/1‑149 (MQ=255)
                                                                 tGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTctgcgtctgcgt                                                                               >  1:40224/1‑149 (MQ=255)
                                                                   gTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCgtgg                                                                             >  5:296672/1‑149 (MQ=255)
                                                                   gTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCgtg                                                                              <  1:28266/148‑1 (MQ=255)
                                                                    ttATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCgtggt                                                                            >  6:62299/1‑149 (MQ=255)
                                                                    ttATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTACGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCAGGAAGCGCCGAACGACAGTCTGCGTCTGCgtggt                                                                            >  6:419615/1‑149 (MQ=255)
                                                                       tCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGgtgt                                                                          <  3:97686/148‑1 (MQ=255)
                                                                       tCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGgtgt                                                                          <  6:125973/148‑1 (MQ=255)
                                                                         aCCCGGGTAGACGGTAGCGATCGTTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCa                                                                       <  4:382285/149‑1 (MQ=255)
                                                                         aCCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCa                                                                       <  5:42895/149‑1 (MQ=255)
                                                                           cccGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATg                                                                     >  8:396624/1‑149 (MQ=255)
                                                                           cccGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATg                                                                     >  2:398198/1‑149 (MQ=255)
                                                                               gTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGt                                                                 <  4:158202/149‑1 (MQ=255)
                                                                                    cGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTc                                                            <  7:241257/149‑1 (MQ=255)
                                                                                      gTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGt                                                          <  5:257662/149‑1 (MQ=255)
                                                                                       tAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGtt                                                         <  6:197098/149‑1 (MQ=255)
                                                                                         gCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTa                                                        <  8:297858/148‑1 (MQ=255)
                                                                                         gCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTa                                                        >  2:115871/1‑148 (MQ=255)
                                                                                           gATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGATCTCGGTGATCTCCTTAACATAAAAAAAAAAAAGCCGCACGACATTATTCGTCTGCGTGGGGTTCATGCTGTTACTCGTTACAt                                                     >  4:435200/1‑149 (MQ=255)
                                                                                             tCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCg                                                   >  6:282133/1‑149 (MQ=255)
                                                                                               cGTACGAAGAGATGATTCCGAAATGGCGTCCGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGtt                                                 <  7:92005/149‑1 (MQ=255)
                                                                                               cGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGTGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGt                                                  <  4:4593/148‑1 (MQ=255)
                                                                                                 tACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTaa                                               <  7:190240/149‑1 (MQ=255)
                                                                                                  aCGAAGAGATGATTCCGAAATGGCGTCAGCTCACCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGAAGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGGTAAc                                              >  3:222005/1‑149 (MQ=255)
                                                                                                   cgacgaGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAAc                                              <  1:100319/148‑1 (MQ=255)
                                                                                                   agAAGAGATGATTCCGAAATGGCGGCAGCTCACCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAAc                                              <  5:279204/147‑1 (MQ=255)
                                                                                                    gAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACTACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACg                                             >  4:274574/1‑148 (MQ=255)
                                                                                                    gAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGt                                                                 <  7:437398/128‑1 (MQ=255)
                                                                                                     aaGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACg                                             >  5:129388/1‑147 (MQ=255)
                                                                                                          atgatTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACa                                      <  5:315364/149‑1 (MQ=255)
                                                                                                             atTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGa                                   <  2:23285/149‑1 (MQ=255)
                                                                                                              ttCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGAc                                  <  2:363696/149‑1 (MQ=255)
                                                                                                                ccGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACg                                 >  3:148418/1‑148 (MQ=255)
                                                                                                                 cGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTctgcgtctgcgt                                                                               <  1:184508/101‑1 (MQ=255)
                                                                                                                 cGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGAAGGTCCGGAAGCGCCGCACGACATtctgagtctgctt                                                                               >  2:184508/1‑99 (MQ=255)
                                                                                                                    aaTGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATAc                            <  2:401336/149‑1 (MQ=255)
                                                                                                                     aTGGCGTCAGCTCAACGTGGTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATAcc                           <  7:207196/149‑1 (MQ=255)
                                                                                                                        gCGTCAGCTCAAGGTGTTCGAATGTAAACGTGTAGAACGTTGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTTGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGt                         <  8:344933/148‑1 (MQ=255)
                                                                                                                        gCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGt                         <  5:245629/148‑1 (MQ=255)
                                                                                                                         cGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCt                       <  1:445810/149‑1 (MQ=255)
                                                                                                                           tCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGCCGGTCCGGAAGCGCCGCACGCCATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                      >  5:364108/1‑148 (MQ=255)
                                                                                                                           tCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                      >  1:404764/1‑148 (MQ=255)
                                                                                                                              gCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAgg                  <  5:62299/149‑1 (MQ=255)
                                                                                                                               cTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                 >  1:106871/1‑149 (MQ=255)
                                                                                                                               cTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                 <  2:29169/149‑1 (MQ=255)
                                                                                                                                  aaCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt              >  8:389538/1‑149 (MQ=255)
                                                                                                                                   aCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt             <  8:34903/149‑1 (MQ=255)
                                                                                                                                   aCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt              >  6:247137/1‑148 (MQ=255)
                                                                                                                                      tgtTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCTCAAGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGAAGGGCGtt             >  3:57895/1‑146 (MQ=255)
                                                                                                                                       gtTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAg                                     <  3:116240/121‑1 (MQ=255)
                                                                                                                                        ttCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGAt        >  1:299038/1‑149 (MQ=255)
                                                                                                                                           gAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTaa                                       <  2:52706/115‑2 (MQ=255)
                                                                                                                                           gAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTaa                                       >  1:52706/1‑114 (MQ=255)
                                                                                                                                           gAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGAGTCTTCGGGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa     <  3:274574/149‑1 (MQ=255)
                                                                                                                                           gAAGGTAAACGTGTAGAACGTGGTGAAGTAATTTCCGAAAGAAAGTAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa     >  4:97648/1‑149 (MQ=255)
                                                                                                                                             aGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTc                                                            <  6:263151/92‑1 (MQ=255)
                                                                                                                                              ggTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTCCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGa  <  2:450431/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGA  >  NC_000913/4188808‑4189098

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ 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
CGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAA  >  NC_000913/4188807‑4189101
                                                                                                                                                  |                                                                                                                                                    
CGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACG                                                                                                                                                  >  NB501016:230:HFVJLAFXY:2:11302:4993:10181/1‑151 (MQ=60)
CGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACG                                                                                                                                                  >  NB501016:230:HFVJLAFXY:4:21404:20755:17529/1‑151 (MQ=60)
   CAATCGTGGCTGAAATGAGCGGTTTCGTTTCCTTCGGTAAAGAAACCTAAGGTATAGGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTAGGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGT                                                                                                                                               <  NB501016:230:HFVJLAFXY:1:11305:3805:4639/151‑1 (MQ=60)
    AATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTA                                                                                                                                              <  NB501016:230:HFVJLAFXY:3:21510:22785:9727/151‑1 (MQ=60)
    AATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTA                                                                                                                                              <  NB501016:230:HFVJLAFXY:3:11611:6589:3793/151‑1 (MQ=60)
       CCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGA                                                                                                                                            >  NB501016:230:HFVJLAFXY:2:11303:21865:12212/1‑150 (MQ=60)
        CTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAAC                                                                                                                                          >  NB501016:230:HFVJLAFXY:3:11607:16547:18861/1‑151 (MQ=60)
            CTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCCCCCCGGTAGACAGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGG                                                                                                                                      >  NB501016:230:HFVJLAFXY:1:21208:17886:11315/1‑151 (MQ=60)
              GAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTG                                                                                                                                    >  NB501016:230:HFVJLAFXY:2:11206:15836:9414/1‑151 (MQ=60)
                 ATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGATATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACG                                                                                                                                 >  NB501016:230:HFVJLAFXY:1:21207:2498:17688/1‑151 (MQ=60)
                   CAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGGCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGCTTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGATCGTGGTGACGT                                                                                                                                >  NB501016:230:HFVJLAFXY:2:21306:1153:4341/1‑150 (MQ=60)
                      CGGTATCGTTTCCGTCGGTAAAGAAACCAAAGGGTCACGTCGGCGGGTTATGTCCCCGGGAGCCGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTCTCGTAATT                                                                                                                            >  NB501016:230:HFVJLAFXY:2:21304:16109:8019/1‑151 (MQ=60)
                          ATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCG                                                                                                                        >  NB501016:230:HFVJLAFXY:4:21409:23795:19431/1‑151 (MQ=60)
                          ATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCG                                                                                                                        <  NB501016:230:HFVJLAFXY:1:21202:16153:3601/151‑1 (MQ=60)
                                  CTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCG                                                                                                                >  NB501016:230:HFVJLAFXY:4:21504:25570:14060/1‑151 (MQ=60)
                                  CTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCG                                                                                                                >  NB501016:230:HFVJLAFXY:4:21410:16301:5124/1‑151 (MQ=60)
                                      GGTAAAGAAACAAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTACCGATCCGTACGAAGAGATGATTCCGAAAAGGCGTCAGCTCAACGAGTTCGAAGGTAAACGAGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGA                                                                                                              >  NB501016:230:HFVJLAFXY:3:21610:8277:12995/1‑149 (MQ=60)
                                               ACCAAAGTTAAACGTCGTCTGGTTATCAACCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGCTGACGTAATTTCCGACGGTCAGGAAGCtcagcact                                                                                                   >  NB501016:230:HFVJLAFXY:4:21611:2328:6609/1‑143 (MQ=60)
                                                CCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACG                                                                                                   >  NB501016:230:HFVJLAFXY:3:21408:14040:19017/1‑150 (MQ=60)
                                                CCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACG                                                                                                   >  NB501016:230:HFVJLAFXY:4:21505:11499:14065/1‑150 (MQ=60)
                                                 CAAAGGTAAACGTCGTCTGGTTATCACGCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGTCGTCAGCTCAACGTGTTCGAAGGTAAACGTTTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGAC                                                                                                 <  NB501016:230:HFVJLAFXY:2:11309:23577:2264/151‑1 (MQ=60)
                                                  AAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACATGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGG                                                                                                               >  NB501016:230:HFVJLAFXY:4:21502:10527:11269/1‑136 (MQ=60)
                                                   AAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAT                                                                                               <  NB501016:230:HFVJLAFXY:2:11210:4780:17436/151‑1 (MQ=60)
                                                     GGTAAACGTCGTCTGGTTATCACCCGGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTC                                                                                             <  NB501016:230:HFVJLAFXY:3:21406:18137:19719/151‑1 (MQ=60)
                                                     GGTAAACGTCGTCTGGTTATCACCCCGGTCGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTC                                                                                             <  NB501016:230:HFVJLAFXY:2:11102:26759:9586/151‑1 (MQ=60)
                                                     GGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTC                                                                                             <  NB501016:230:HFVJLAFXY:2:11101:1281:9559/151‑1 (MQ=60)
                                                               GTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCG                                                                                   >  NB501016:230:HFVJLAFXY:4:11406:9058:5275/1‑151 (MQ=60)
                                                                 CTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTG                                                                                 >  NB501016:230:HFVJLAFXY:4:11511:10812:13207/1‑151 (MQ=60)
                                                                 CTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTG                                                                                 >  NB501016:230:HFVJLAFXY:1:11106:2913:8939/1‑151 (MQ=60)
                                                                   GGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGT                                                                               >  NB501016:230:HFVJLAFXY:3:21411:15070:12791/1‑151 (MQ=60)
                                                                   GGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGT                                                                               >  NB501016:230:HFVJLAFXY:1:11104:24632:17022/1‑151 (MQ=60)
                                                                   GGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGG                                                                                <  NB501016:230:HFVJLAFXY:1:11104:15267:14628/150‑1 (MQ=60)
                                                                       ATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTT                                                                            <  NB501016:230:HFVJLAFXY:2:11205:6573:11394/150‑1 (MQ=60)
                                                                         CACCCCGGTAGACGGTAGCGATCCGTACGAAGAGAGGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCAT                                                                         <  NB501016:230:HFVJLAFXY:3:11408:9497:1566/151‑1 (MQ=60)
                                                                                  AGACGGTAGCGATCCGTACGAAGAGATGATTCCGAACTGGCGTCAGCTCAACGTGTTCGAAGGTAAACGGGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTAC                                                                 <  NB501016:230:HFVJLAFXY:3:21611:3698:12371/150‑1 (MQ=60)
                                                                                    ACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCG                                                              <  NB501016:230:HFVJLAFXY:1:11109:10034:7541/151‑1 (MQ=60)
                                                                                    ACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCG                                                              <  NB501016:230:HFVJLAFXY:4:21401:25734:6712/151‑1 (MQ=60)
                                                                                     CGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGT                                                             <  NB501016:230:HFVJLAFXY:4:21601:6324:15863/151‑1 (MQ=60)
                                                                                      GGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTT                                                            <  NB501016:230:HFVJLAFXY:3:21405:17314:10032/151‑1 (MQ=60)
                                                                                       GTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTT                                                            >  NB501016:230:HFVJLAFXY:1:21302:1783:13087/1‑150 (MQ=60)
                                                                                               CCGTACGAAGAGATGATTCCGAAATGGCGTCCGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTA                                                   <  NB501016:230:HFVJLAFXY:4:11502:3943:3169/151‑1 (MQ=60)
                                                                                                 GTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAAC                                                 <  NB501016:230:HFVJLAFXY:4:11605:9181:5250/151‑1 (MQ=60)
                                                                                                  TACGAAGAGATGATTCCGAAATGGCGTCAGCTCACCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGAAGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGGTAACG                                                >  NB501016:230:HFVJLAFXY:2:21101:2803:3133/1‑151 (MQ=60)
                                                                                                   caGAAGAGATGATTCCGAAATGGCGGCAGCTCACCGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACG                                                <  NB501016:230:HFVJLAFXY:3:21408:1223:17930/148‑1 (MQ=60)
                                                                                                   ACGACGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACG                                                <  NB501016:230:HFVJLAFXY:1:11203:17558:11003/150‑1 (MQ=60)
                                                                                                    CGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTATTTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTT                                                                   <  NB501016:230:HFVJLAFXY:4:21608:7471:5253/130‑1 (MQ=60)
                                                                                                    CGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGA                                               >  NB501016:230:HFVJLAFXY:3:11509:6717:2304/1‑150 (MQ=60)
                                                                                                     GAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGA                                               >  NB501016:230:HFVJLAFXY:3:11509:12346:18349/1‑149 (MQ=60)
                                                                                                          GATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAG                                        <  NB501016:230:HFVJLAFXY:3:21503:13711:12174/151‑1 (MQ=60)
                                                                                                               TTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAAAGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGT                                   >  NB501016:230:HFVJLAFXY:3:21412:20804:4110/1‑151 (MQ=60)
                                                                                                                TCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGT                                   >  NB501016:230:HFVJLAFXY:2:11302:4928:1328/1‑150 (MQ=60)
                                                                                                                 CCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTG                                                                                 <  NB501016:230:HFVJLAFXY:1:11305:25224:2534/103‑1 (MQ=60)
                                                                                                                  CGAAATGGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTA                                                           >  NB501016:230:HFVJLAFXY:4:11605:16232:10374/1‑124 (MQ=60)
                                                                                                                     AATGGCGTCAGCTCAACGTGGTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCG                             <  NB501016:230:HFVJLAFXY:4:11607:8627:18072/151‑1 (MQ=60)
                                                                                                                        GGCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTC                           <  NB501016:230:HFVJLAFXY:3:21403:7310:10931/150‑1 (MQ=60)
                                                                                                                         GCGTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTG                         <  NB501016:230:HFVJLAFXY:1:21309:25142:13164/151‑1 (MQ=60)
                                                                                                                           GTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGCCGGTCCGGAAGCGCCGCACGCCATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGC                        >  NB501016:230:HFVJLAFXY:3:21511:24052:8339/1‑150 (MQ=60)
                                                                                                                           GTCAGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGC                        >  NB501016:230:HFVJLAFXY:1:21303:7825:4692/1‑150 (MQ=60)
                                                                                                                              AGCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGG                    <  NB501016:230:HFVJLAFXY:3:11411:14498:3803/151‑1 (MQ=60)
                                                                                                                               GCTCAACGTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGC                   >  NB501016:230:HFVJLAFXY:1:11204:20994:12049/1‑151 (MQ=60)
                                                                                                                                      GTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCTCAAGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGAAGGGCGTTA               >  NB501016:230:HFVJLAFXY:2:11110:1895:15242/1‑148 (MQ=60)
                                                                                                                                      GTGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGA            <  NB501016:230:HFVJLAFXY:2:11308:22065:11758/151‑1 (MQ=60)
                                                                                                                                       TGTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGG                                       <  NB501016:230:HFVJLAFXY:2:11208:5991:13733/123‑1 (MQ=60)
                                                                                                                                        GTTCGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATT          >  NB501016:230:HFVJLAFXY:1:21110:16149:3768/1‑151 (MQ=60)
                                                                                                                                           CGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTAac                                         >  NB501016:230:HFVJLAFXY:1:11108:22167:5461/1‑115 (MQ=60)
                                                                                                                                           CGAAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGAGTCTTCGGGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAAC       <  NB501016:230:HFVJLAFXY:2:21110:21473:1328/151‑1 (MQ=60)
                                                                                                                                             AAGGTAAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCG                                                              >  NB501016:230:HFVJLAFXY:3:21406:14662:6829/1‑94 (MQ=60)
                                                                                                                                                GTAAACTTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAA  <  NB501016:230:HFVJLAFXY:3:21404:11634:18880/151‑1 (MQ=60)
                                                                                                                                                  |                                                                                                                                                    
CGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAA  >  NC_000913/4188807‑4189101

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