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
A4 F81 I0 R1 263 75.1 3021180 96.0% 2900332 137.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,722,492 C→A 75.4% T847K (ACG→AAG)  pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,722,4920CA75.4% 162.4 / 63.3 97T847K (ACG→AAG) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (11/13);  new base A (35/38);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAC  >  NC_000913/2722349‑2722635
                                                                                                                                               |                                                                                                                                               
gATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGc                                                                                                                                            <  6:65672/149‑1 (MQ=255)
gATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGc                                                                                                                                            <  1:87915/149‑1 (MQ=255)
     ccGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCaaa                                                                                                                                        <  3:319046/148‑1 (MQ=255)
           cTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGt                                                                                                                                 <  4:71272/149‑1 (MQ=255)
                 aCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAAt                                                                                                                                    <  5:173443/140‑1 (MQ=255)
                 aCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAAt                                                                                                                                    >  6:173443/1‑140 (MQ=255)
                 aCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAtg                                                                                                                           <  8:121984/149‑1 (MQ=255)
                 aCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAt                                                                                                                            <  5:188598/148‑1 (MQ=255)
                        tGCCGAATTTTCTGTTCTGGTTCGCTGGGATCTCTAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGgcgc                                                                                                                    <  3:160549/149‑1 (MQ=255)
                           cGAATTTGCTGTACTGGTTCGCTCGGATCTCAAATGGTTAGTCTTAGGTCGACGCTTAATGGAAAAGTTGATTCCCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTa                                                                                                                  <  3:242357/148‑1 (MQ=255)
                           cGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAg                                                                                                                 <  1:207195/149‑1 (MQ=255)
                                    tGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGc                                                                                                        >  8:210579/1‑149 (MQ=255)
                                    tGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGc                                                                                                        >  3:179600/1‑149 (MQ=255)
                                      tACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTc                                                                                                      >  3:77673/1‑149 (MQ=255)
                                      tACTGGTTCGCTCGGATCTCAAAGGGTTAGGCATAGGTCGACGCTTAATGGAAAAGTTGATTACATATACGCGAGATCACGGACTACAACGTCTGAATGGTAATACGATGCCAAACAAACGTGCCATGGTGGCGCTACCCCGCAAGCTc                                                                                                      >  6:357698/1‑149 (MQ=255)
                                        cTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCgg                                                                                                    <  1:179734/149‑1 (MQ=255)
                                        cTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCgg                                                                                                    <  8:137290/149‑1 (MQ=255)
                                           gTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACTATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGtt                                                                                                 <  4:51637/149‑1 (MQ=255)
                                           gTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTgg                                                                                                                               <  7:370396/119‑1 (MQ=255)
                                           gTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTgg                                                                                                                               >  8:370396/1‑119 (MQ=255)
                                               gCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAAc                                                                                             >  2:188138/1‑149 (MQ=255)
                                                cTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACg                                                                                            >  3:161472/1‑149 (MQ=255)
                                                 tCGGATCTCAAAGGGTTAGGCTTAGGTCGCCGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGt                                                                                           >  7:270664/1‑149 (MQ=255)
                                                 tCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGGATTAAGATGCCAACCAATCGTGGAATGTTGGCGCTAGACCGCAAGCTCGGGGTTAACGt                                                                                           >  6:234897/1‑149 (MQ=255)
                                                   ggATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGtt                                                                                          >  5:7452/1‑148 (MQ=255)
                                                   ggATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGtt                                                                                          <  6:8169/148‑1 (MQ=255)
                                                   ggATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGtt                                                                                          >  5:246066/1‑148 (MQ=255)
                                                    gATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATc                                                                                                                                   >  5:223429/1‑106 (MQ=255)
                                                    gATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATc                                                                                                                                   <  6:223429/106‑1 (MQ=255)
                                                        tcAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATc                                                                                    >  6:251276/1‑149 (MQ=255)
                                                          aaaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGCTTACCTATACGCGAGATCACGGACTACAACGTCTGAAGGGTATTAAGATGCCAAACAATCGCGGCATGGTGGCGCTAGCCCGCAAGCTAGGGTTTACCGTTGATATc                                                                                    >  2:318750/1‑147 (MQ=255)
                                                          aaaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTGCGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATc                                                                                    >  1:15474/1‑147 (MQ=255)
                                                          aaaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGt                                                                                                                                 >  3:251610/1‑102 (MQ=255)
                                                          aaaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGt                                                                                                                                 <  4:251610/102‑1 (MQ=255)
                                                          aaaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGACAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATc                                                                                    >  4:5797/1‑147 (MQ=255)
                                                          aaaGGGTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAtggtg                                                                                                                        <  1:220055/110‑1 (MQ=255)
                                                          aaaGGGTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAtggtg                                                                                                                        >  2:220055/1‑110 (MQ=255)
                                                                        aGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAggg                                                                      <  4:77673/147‑1 (MQ=255)
                                                                        aGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAggg                                                                      <  5:251276/147‑1 (MQ=255)
                                                                        aGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACAGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAggg                                                                      <  6:7452/147‑1 (MQ=255)
                                                                             gACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAtggtg                                                                                                                        <  8:365225/92‑1 (MQ=255)
                                                                             gACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAtggtg                                                                                                                        >  7:365225/1‑92 (MQ=255)
                                                                               cGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTgg                                                             <  1:188138/149‑1 (MQ=255)
                                                                                gCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGTTCGTTggg                                                            >  1:76136/1‑149 (MQ=255)
                                                                                 cTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTGCGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGc                                                           <  2:15474/149‑1 (MQ=255)
                                                                                  ttAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCaa                                                                                                          <  5:188897/101‑1 (MQ=255)
                                                                                  ttAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCaa                                                                                                          >  6:188897/1‑101 (MQ=255)
                                                                                  ttAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCt                                                          >  3:69704/1‑149 (MQ=255)
                                                                                  ttAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCt                                                          >  5:260174/1‑149 (MQ=255)
                                                                                     attgAAAAGTTGAGTACCTATACGCGAGATCACTGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTAc                                                       <  6:120707/146‑1 (MQ=255)
                                                                                       ggAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGc                                                     >  2:384483/1‑149 (MQ=255)
                                                                                         aaaaGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCt                                                    >  5:351847/1‑148 (MQ=255)
                                                                                             gTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCggg                                                                                                   >  7:46057/1‑97 (MQ=255)
                                                                                             gTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCggg                                                                                                   <  8:46057/97‑1 (MQ=255)
                                                                                                 aTTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGTGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTg                                            <  4:84831/148‑1 (MQ=255)
                                                                                                     ccTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTAc                                                       <  4:203102/133‑1 (MQ=255)
                                                                                                     ccTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTAc                                                       >  3:203102/1‑133 (MQ=255)
                                                                                                      cTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATTCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCAAGCTCGAagagt                                                                       >  6:292685/1‑115 (MQ=255)
                                                                                                      cTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCTAAGAgg                                                                       <  5:292685/116‑1 (MQ=255)
                                                                                                       tataCGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAgcg                                     >  4:193218/1‑149 (MQ=255)
                                                                                                           cgcgAGATCACGGACTACAACGTCTGATTGGTATTAAGAGGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGa                                  <  4:200807/148‑1 (MQ=255)
                                                                                                           cgcgAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGa                                  <  4:161472/148‑1 (MQ=255)
                                                                                                                  tCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCAt                          <  2:189170/149‑1 (MQ=255)
                                                                                                                    cagttcTACAACGTCTTAATTGAATAAATAAGCCAAACAATCGTGTCATGGTTGCGCTAGCCAGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTTCCCAGCGCGAGTAATCAt                          <  4:91564/142‑1 (MQ=255)
                                                                                                                    aCGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCAt                          >  8:222224/1‑147 (MQ=255)
                                                                                                                    aCGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGa                        <  4:179600/149‑1 (MQ=255)
                                                                                                                      ggACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGc                                                           <  5:320933/112‑1 (MQ=255)
                                                                                                                      ggACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGc                                                           >  6:320933/1‑112 (MQ=255)
                                                                                                                           aCAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAgg                  >  6:59787/1‑148 (MQ=255)
                                                                                                                            cAACGTCTTAATGGTATTAAGATGACAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTAGGGTTTAACGTTGATATAAAGCTCGAAGAGGGGATCGTTGTGCTTACGCTAAATATTGCCCAGCGCGAGGAATCATGAgagggga                 >  3:91564/1‑142 (MQ=255)
                                                                                                                             aaCGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGtt                                                                                          >  5:80270/1‑74 (MQ=255)
                                                                                                                             aaCGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGtt                                                                                          <  6:80270/74‑1 (MQ=255)
                                                                                                                             aaCGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTGGCCCGCAAGCTCGGGTTTAACGTTGATATcc                                                                                   <  7:225514/81‑1 (MQ=255)
                                                                                                                             aaCGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTGGCCCGCAAGCTCGGGTTTAACGTTGATATcc                                                                                   >  8:225514/1‑81 (MQ=255)
                                                                                                                             aaCGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTAc               >  8:46039/1‑149 (MQ=255)
                                                                                                                              aCGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACt              <  3:5797/149‑1 (MQ=255)
                                                                                                                               cGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTg             <  6:260174/149‑1 (MQ=255)
                                                                                                                               cGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTg             <  3:193218/149‑1 (MQ=255)
                                                                                                                                gTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTg             >  7:191892/1‑148 (MQ=255)
                                                                                                                                 tCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTgg            >  7:238404/1‑148 (MQ=255)
                                                                                                                                   tGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGaaa         >  4:261681/1‑149 (MQ=255)
                                                                                                                                     aaTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAgg                  >  4:87053/1‑138 (MQ=255)
                                                                                                                                     aaTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAgg                  <  3:87053/138‑1 (MQ=255)
                                                                                                                                     aaTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGANATg       >  5:36455/1‑149 (MQ=255)
                                                                                                                                     aaTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATg       <  4:69704/149‑1 (MQ=255)
                                                                                                                                     aaTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATg       >  8:209137/1‑149 (MQ=255)
                                                                                                                                      aTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGt      <  6:351847/149‑1 (MQ=255)
                                                                                                                                       tGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGTAATCATGAGTAAGGTACTGGAAATGtt     <  5:357698/149‑1 (MQ=255)
                                                                                                                                       tGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGATCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGtt     <  5:234897/149‑1 (MQ=255)
                                                                                                                                       tGGTATTAAGATGCCAAACAATCGNTGCNTGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGtt     <  6:36455/149‑1 (MQ=255)
                                                                                                                                         gTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAACTCGGGTTTAACGTTGATATCCAGCTCGAAGAgg                                                                       >  6:337657/1‑81 (MQ=255)
                                                                                                                                         gTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAACTCGGGTTTAACGTTGATATCCAGCTCGAAGAgg                                                                       <  5:337657/81‑1 (MQ=255)
                                                                                                                                          tATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAc  <  8:238404/149‑1 (MQ=255)
                                                                                                                                           aTTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGa                                  >  4:34668/1‑116 (MQ=255)
                                                                                                                                           aTTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGa                                  <  3:34668/116‑1 (MQ=255)
                                                                                                                                            ttAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGa                                                                    <  5:292196/81‑1 (MQ=255)
                                                                                                                                            ttAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGa                                                                    >  6:292196/1‑81 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAC  >  NC_000913/2722349‑2722635

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ 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
CGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTG  >  NC_000913/2722348‑2722633
                                                                                                                                                |                                                                                                                                             
CGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCC                                                                                                                                         <  NB501016:230:HFVJLAFXY:1:11204:14123:18739/151‑1 (MQ=60)
     TCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAAC                                                                                                                                     <  NB501016:230:HFVJLAFXY:2:21306:17060:2281/150‑1 (MQ=60)
                 AACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATC                                                                                                                                 <  NB501016:230:HFVJLAFXY:3:11609:5746:15950/142‑1 (MQ=60)
                 AACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATG                                                                                                                         <  NB501016:230:HFVJLAFXY:3:11612:26774:11834/150‑1 (MQ=60)
                        ATGCCGAATTTTCTGTTCTGGTTCGCTGGGATCTCTAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCT                                                                                                                 <  NB501016:230:HFVJLAFXY:2:11309:25974:16389/151‑1 (MQ=60)
                           CCGAATTTGCTGTCCTGGTGCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAcgattgtttggcatcttaataccatt                                                                                                              <  NB501016:230:HFVJLAFXY:4:21406:25598:19551/151‑27 (MQ=60)
                           CCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAATGGTTAGTCTTAGGTCGACGCTTAATGGAAAAGTTGATTCCCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAG                                                                                                               <  NB501016:230:HFVJLAFXY:2:21202:9632:11393/150‑1 (MQ=60)
                           CCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGC                                                                                                              <  NB501016:230:HFVJLAFXY:1:21103:7419:8043/151‑1 (MQ=60)
                                    CTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCT                                                                                                     >  NB501016:230:HFVJLAFXY:2:21101:23694:9508/1‑151 (MQ=60)
                                      GTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCG                                                                                                   >  NB501016:230:HFVJLAFXY:2:11205:18682:3494/1‑151 (MQ=60)
                                        ACTGGTGCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGG                                                                                                 <  NB501016:230:HFVJLAFXY:1:11310:6358:10882/151‑1 (MQ=60)
                                           GGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGC                                                                                                                            <  NB501016:230:HFVJLAFXY:4:21609:14014:6819/121‑1 (MQ=60)
                                                GCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGT                                                                                         >  NB501016:230:HFVJLAFXY:2:11309:21416:19746/1‑151 (MQ=60)
                                                 CTCGGATCTCAAAGGGTTAGGCTTAGGTCGCCGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTT                                                                                        >  NB501016:230:HFVJLAFXY:4:21502:9964:6449/1‑151 (MQ=60)
                                                 CTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGA                                                                                                                                             <  NB501016:230:HFVJLAFXY:4:21507:10425:5250/98‑1 (MQ=60)
                                                   CGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTG                                                                                       >  NB501016:230:HFVJLAFXY:3:11402:23863:8845/1‑150 (MQ=60)
                                                   CGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTG                                                                                       >  NB501016:230:HFVJLAFXY:3:21411:25548:10561/1‑150 (MQ=60)
                                                    GGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCG                                                                                                                                >  NB501016:230:HFVJLAFXY:3:21407:10873:5358/1‑108 (MQ=60)
                                                          CAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTGCGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCC                                                                                 >  NB501016:230:HFVJLAFXY:1:11103:11458:12713/1‑149 (MQ=60)
                                                          CAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTG                                                                                                                              >  NB501016:230:HFVJLAFXY:2:21204:3370:10418/1‑104 (MQ=60)
                                                          CAAAGGGTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGG                                                                                                                     <  NB501016:230:HFVJLAFXY:1:21105:23257:11898/112‑1 (MQ=60)
                                                                        TAGGTCGACGCTTAATGGAAAAGTTTATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGG                                                                   <  NB501016:230:HFVJLAFXY:3:21412:26219:10049/149‑1 (MQ=60)
                                                                             CGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATAACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGG                                                                                                                     >  NB501016:230:HFVJLAFXY:4:21608:19455:6410/1‑94 (MQ=60)
                                                                               ACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGG                                                          <  NB501016:230:HFVJLAFXY:1:11312:15231:2537/151‑1 (MQ=60)
                                                                                CGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGTTCGTTGGGG                                                         >  NB501016:230:HFVJLAFXY:1:11202:2788:13710/1‑151 (MQ=60)
                                                                                  CTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTT                                                       >  NB501016:230:HFVJLAFXY:2:11203:13647:9725/1‑151 (MQ=60)
                                                                                  CTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTT                                                       >  NB501016:230:HFVJLAFXY:3:21503:15293:3387/1‑151 (MQ=60)
                                                                                  CTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAG                                                                                                       <  NB501016:230:HFVJLAFXY:3:11612:26642:12921/103‑1 (MQ=60)
                                                                                         GAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTA                                                 >  NB501016:230:HFVJLAFXY:3:21608:21789:7964/1‑150 (MQ=60)
                                                                                             AGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGT                                                                                                >  NB501016:230:HFVJLAFXY:4:11409:8599:2976/1‑99 (MQ=60)
                                                                                                     ACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACG                                                    >  NB501016:230:HFVJLAFXY:2:21106:9480:6758/1‑135 (MQ=60)
                                                                                                      CCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCTAAGAGGG                                                                    <  NB501016:230:HFVJLAFXY:3:21509:26571:7403/118‑1 (MQ=60)
                                                                                                           ACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGTCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGG                              <  NB501016:230:HFVJLAFXY:3:21410:25425:13203/151‑1 (MQ=60)
                                                                                                              CGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCC                                                                                                                                         >  NB501016:230:HFVJLAFXY:1:11302:11497:16812/1‑41 (MQ=60)
                                                                                                                      CGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCT                                                        <  NB501016:230:HFVJLAFXY:3:21602:24143:14207/114‑1 (MQ=60)
                                                                                                                            ACAACGTCTTAATGGTATTAAGATGACAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTAGGGTTTAACGTTGATATAAAGCTCGAAGAGGGGATCGTTGTGCTTACGCTAAATATTGCCCAGCGCGAGGAATCATGAGaggggaa              >  NB501016:230:HFVJLAFXY:2:11208:1432:2086/1‑143 (MQ=60)
                                                                                                                             CAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTG                                                                                       >  NB501016:230:HFVJLAFXY:3:11504:11526:5357/1‑76 (MQ=60)
                                                                                                                             CAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTGGCCCGCAAGCTCGGGTTTAACGTTGATATCCA                                                                                <  NB501016:230:HFVJLAFXY:4:21405:7497:18831/83‑1 (MQ=60)
                                                                                                                              AACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTG           <  NB501016:230:HFVJLAFXY:2:11102:1674:4651/151‑1 (MQ=60)
                                                                                                                               ACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGG          <  NB501016:230:HFVJLAFXY:2:21104:19163:5806/151‑1 (MQ=60)
                                                                                                                                CGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGG          >  NB501016:230:HFVJLAFXY:4:11611:24319:16759/1‑150 (MQ=60)
                                                                                                                                 GTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGA         >  NB501016:230:HFVJLAFXY:4:21408:12140:4358/1‑150 (MQ=60)
                                                                                                                                     GAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAATCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGT               <  NB501016:230:HFVJLAFXY:2:11207:20455:2989/140‑1 (MQ=60)
                                                                                                                                     GAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGANATGT    >  NB501016:230:HFVJLAFXY:3:11407:18875:20386/1‑151 (MQ=60)
                                                                                                                                       ATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGTAATCATGAGTAAGGTACTGGAAATGTTG  <  NB501016:230:HFVJLAFXY:3:21609:13919:8843/151‑1 (MQ=60)
                                                                                                                                       ATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGATCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTG  <  NB501016:230:HFVJLAFXY:3:21409:6562:8994/151‑1 (MQ=60)
                                                                                                                                         GGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAACTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGG                                                                    <  NB501016:230:HFVJLAFXY:3:21605:4557:15144/83‑1 (MQ=60)
                                                                                                                                           TATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAG                               <  NB501016:230:HFVJLAFXY:2:11108:26563:4257/118‑1 (MQ=60)
                                                                                                                                            ATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGAT                                                                 <  NB501016:230:HFVJLAFXY:3:21509:26118:5711/83‑1 (MQ=60)
                                                                                                                                                |                                                                                                                                             
CGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTG  >  NC_000913/2722348‑2722633

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