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 F2 I1 R1 15 39.5 1304580 98.1% 1279792 140.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,722,270 T→A I773N (ATC→AAC)  pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,722,2700TA98.2% 184.1 / NA 54I773N (ATC→AAC) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (26/27);  minor base C (1/0);  total (27/27)
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

CGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGT  >  NC_000913/2722125‑2722413
                                                                                                                                                 |                                                                                                                                               
cGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGa                                                                                                                                              <  7:113661/149‑1 (MQ=255)
cGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGa                                                                                                                                              <  3:45405/149‑1 (MQ=255)
 gctgctTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGAc                                                                                                                                             <  8:122183/149‑1 (MQ=255)
    gctTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTAc                                                                                                                                          >  8:54113/1‑149 (MQ=255)
     ctTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTAc                                                                                                                                          <  5:117109/148‑1 (MQ=255)
                     ttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGt                                                                                                                          <  1:11872/148‑1 (MQ=255)
                     ttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGAGCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGt                                                                                                                          <  4:34858/148‑1 (MQ=255)
                         cAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTg                                                                                                                       >  2:17085/1‑147 (MQ=255)
                           gaagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGc                                                                                                                   <  7:73693/149‑1 (MQ=255)
                            aagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTAATACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCg                                                                                                                  <  1:63466/149‑1 (MQ=255)
                             agaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCgg                                                                                                                 <  3:63282/149‑1 (MQ=255)
                                   gCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTacgac                                                                                                           >  6:167363/1‑149 (MQ=255)
                                      acaACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTa                                                                                                        <  2:91763/149‑1 (MQ=255)
                                      acaACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTa                                                                                                        >  1:90796/1‑149 (MQ=255)
                                           ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAg                                                                                                                    >  4:4873/1‑132 (MQ=255)
                                           ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAg                                                                                                                    <  3:4873/132‑1 (MQ=255)
                                           ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGAt                                                                                                   >  7:92919/1‑149 (MQ=255)
                                           ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGAt                                                                                                   >  6:131014/1‑149 (MQ=255)
                                                cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTAt                                                                                                       >  1:76438/1‑140 (MQ=255)
                                                cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTAt                                                                                                       <  2:76438/140‑1 (MQ=255)
                                                cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAc                                                                                              >  6:13072/1‑149 (MQ=255)
                                                cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAc                                                                                              <  5:131014/149‑1 (MQ=255)
                                                cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAc                                                                                              >  7:116658/1‑149 (MQ=255)
                                                cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAc                                                                                              <  7:26514/149‑1 (MQ=255)
                                                  aTTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACg                                                                                             >  3:61050/1‑148 (MQ=255)
                                                         tCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAgaga                                                                                      <  4:120975/148‑1 (MQ=255)
                                                             gAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTc                                                                                 >  3:110768/1‑149 (MQ=255)
                                                                  aCCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGt                                                                            >  6:112321/1‑149 (MQ=255)
                                                                  aCCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTCGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGt                                                                            >  7:5414/1‑149 (MQ=255)
                                                                  aCCAAAGAAGATATTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGCGAAATCGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGt                                                                            >  7:122858/1‑149 (MQ=255)
                                                                    cAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa                                                                          >  7:53721/1‑149 (MQ=255)
                                                                     aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAc                                                                         <  4:110768/149‑1 (MQ=255)
                                                                      aagaagATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACCCAGACCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                        >  6:42809/1‑149 (MQ=255)
                                                                          agaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTg                                                                    <  5:13072/149‑1 (MQ=255)
                                                                             tCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCg                                                                             <  5:102056/137‑1 (MQ=255)
                                                                             tCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCg                                                                             >  6:102056/1‑137 (MQ=255)
                                                                                  attaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTcc                                                            >  8:123386/1‑149 (MQ=255)
                                                                                       cGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGACAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTAGTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATcc                                                       >  8:124948/1‑149 (MQ=255)
                                                                                            cTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGa                                                    <  8:116658/147‑1 (MQ=255)
                                                                                            cTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGAt                                                   >  2:147588/1‑148 (MQ=255)
                                                                                                     gaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATg                                         >  4:17316/1‑149 (MQ=255)
                                                                                                              aTTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                <  8:92919/149‑1 (MQ=255)
                                                                                                              aTTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                <  1:147588/149‑1 (MQ=255)
                                                                                                                  aCCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTg                             <  5:112321/148‑1 (MQ=255)
                                                                                                                      aTGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAt                                   <  3:114273/138‑1 (MQ=255)
                                                                                                                      aTGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTg                        >  6:100885/1‑149 (MQ=255)
                                                                                                                      aTGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATAGATGCCGAAt                                   >  4:114273/1‑138 (MQ=255)
                                                                                                                                        tGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa      <  5:100885/149‑1 (MQ=255)
                                                                                                                                        tGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa      <  4:61050/149‑1 (MQ=255)
                                                                                                                                        tGACACAGAACGACTACGATCGGGAAATGGCGTTCGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa      <  8:5414/149‑1 (MQ=255)
                                                                                                                                          acacAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAgg    >  6:26505/1‑149 (MQ=255)
                                                                                                                                           cacaGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAggg   <  7:54113/149‑1 (MQ=255)
                                                                                                                                            acaGAACGACTACGATCGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGt  >  2:27125/1‑148 (MQ=255)
                                                                                                                                              aGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGt  >  2:4973/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGT  >  NC_000913/2722125‑2722413

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
ACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGT  >  NC_000913/2722124‑2722413
                                                                                                                                                  |                                                                                                                                               
ACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGAC                                                                                                                                             <  NB501016:228:HFK2FAFXY:2:11211:9236:7598/151‑1 (MQ=60)
ACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGAC                                                                                                                                             <  NB501016:228:HFK2FAFXY:4:21507:13862:14432/151‑1 (MQ=60)
     GCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACG                                                                                                                                         <  NB501016:228:HFK2FAFXY:3:21502:8666:17216/150‑1 (MQ=60)
                     TTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTT                                                                                                                         <  NB501016:228:HFK2FAFXY:1:11105:19409:8467/150‑1 (MQ=60)
                           AGAAGACGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCG                                                                                                                  <  NB501016:228:HFK2FAFXY:4:11612:4202:15273/151‑1 (MQ=60)
                            GAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTAATACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGG                                                                                                                 <  NB501016:228:HFK2FAFXY:1:11301:10337:1911/151‑1 (MQ=60)
                             AAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGT                                                                                                                <  NB501016:228:HFK2FAFXY:2:11307:1717:7563/151‑1 (MQ=60)
                                      CACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTAT                                                                                                       >  NB501016:228:HFK2FAFXY:1:11311:21474:14944/1‑151 (MQ=60)
                                           CTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATC                                                                                                  >  NB501016:228:HFK2FAFXY:4:21409:11950:13561/1‑151 (MQ=60)
                                           CTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGC                                                                                                                   <  NB501016:228:HFK2FAFXY:2:11103:13476:8600/134‑1 (MQ=60)
                                                GCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACG                                                                                             <  NB501016:228:HFK2FAFXY:3:21509:24753:3460/151‑1 (MQ=60)
                                                GCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACG                                                                                             >  NB501016:228:HFK2FAFXY:4:21509:19207:6934/1‑151 (MQ=60)
                                                GCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACG                                                                                             <  NB501016:228:HFK2FAFXY:4:11412:13355:15921/151‑1 (MQ=60)
                                                GCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATT                                                                                                      >  NB501016:228:HFK2FAFXY:1:11306:9805:3250/1‑142 (MQ=60)
                                                  AATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGG                                                                                            >  NB501016:228:HFK2FAFXY:2:11306:12295:8281/1‑150 (MQ=60)
                                                             CGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCG                                                                                >  NB501016:228:HFK2FAFXY:2:21206:24849:8211/1‑151 (MQ=60)
                                                                  CACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTCGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTC                                                                           >  NB501016:228:HFK2FAFXY:4:11403:16406:12181/1‑151 (MQ=60)
                                                                  CACCAAAGAAGATATTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGCGAAATCGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTC                                                                           >  NB501016:228:HFK2FAFXY:4:21512:9218:11367/1‑151 (MQ=60)
                                                                    CCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAC                                                                         >  NB501016:228:HFK2FAFXY:4:11602:20273:11110/1‑151 (MQ=60)
                                                                          AAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGC                                                                   <  NB501016:228:HFK2FAFXY:3:11406:3810:16165/151‑1 (MQ=60)
                                                                             ATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGT                                                                            <  NB501016:228:HFK2FAFXY:3:21408:15997:4224/139‑1 (MQ=60)
                                                                                                              AATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGC                               <  NB501016:228:HFK2FAFXY:1:21209:11844:16479/151‑1 (MQ=60)
                                                                                                                  TACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGT                            <  NB501016:228:HFK2FAFXY:3:21412:11015:14889/150‑1 (MQ=60)
                                                                                                                      CATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATT                                  <  NB501016:228:HFK2FAFXY:2:21208:17816:2589/140‑1 (MQ=60)
                                                                                                                                        ATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAG     <  NB501016:228:HFK2FAFXY:3:21407:6117:12997/151‑1 (MQ=60)
                                                                                                                                           ACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGC  <  NB501016:228:HFK2FAFXY:4:11602:17755:14649/151‑1 (MQ=60)
                                                                                                                                                  |                                                                                                                                               
ACGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGT  >  NC_000913/2722124‑2722413

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: