Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,880,033 A→G 23.3% S290G (AGT→GGT)  selA → L‑seryl‑tRNA(Sec) selenium transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,880,0330AG23.3% 302.5 / 33.1 179S290G (AGT→GGT) selAL‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand):  ref base A (21/109);  new base G (37/6);  total (60/119)
Fisher's exact test for biased strand distribution p-value = 1.03e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

CGTCAGGCGGTGAATGAAAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCAAATCCAGGC  >  NC_007779/3879799‑3880268
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
    aGGCGGTGAATGAAAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCGGGAGTTGATTGCGGCGGGGGTCCGGCTGGTGCGTTtctc                                                                                                                                                                                                                                          >  1:340750/1‑238 (MQ=255)
           gaatgaaAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGAc                                                                                                                                                                                                                                  <  1:1238443/238‑1 (MQ=255)
                aaaaTACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGGCTGGTGAGGTTCTCCGGCGACACgttg                                                                                                                                                                                                                           >  1:1276905/1‑239 (MQ=255)
                      ccGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCGGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGAATGCGGCGGGGGGTCGTCTGGGTGAGGTTACCGGCGACAACGTGCTAggg                                                                                                                                                                                                                     >  2:1232460/1‑238 (MQ=255)
                        gCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGTCAGTCTGGGGAGGTTTTCCGGCGGCAAGTTGTTAGGCgg                                                                                                                                                                                                                   >  1:174206/1‑239 (MQ=255)
                         cACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGGCAGGCTGGTGAGTTTTTCCGGCGAGAAATTGGTAGGTggg                                                                                                                                                                                                                  >  1:224339/1‑239 (MQ=255)
                              ttGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGTGTGGCTGGGGGGTTTTCCCCGGGACGAGTTGTTTGGGGGGgcgcg                                                                                                                                                                                                             >  2:407604/1‑238 (MQ=255)
                                  tgaAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGGTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCgcaggc                                                                                                                                                                                                          <  1:43478/238‑1 (MQ=255)
                                  tgaAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCgcaggc                                                                                                                                                                                                          <  2:788735/238‑1 (MQ=255)
                                  tgaAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCgcaggc                                                                                                                                                                                                          <  2:1464263/238‑1 (MQ=255)
                                  tgaAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCgcaggc                                                                                                                                                                                                          <  2:1731977/238‑1 (MQ=255)
                                  ggaAAGTACATACCAGTAACTACAGCATGCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCgcaggc                                                                                                                                                                                                          <  2:874793/237‑1 (MQ=255)
                                      aGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGGCAGTCTGGTGAGGTTCTCCCGCGGCAAGATGTTTTGGGGGGCGcaggcaggga                                                                                                                                                                                                     >  1:547978/1‑237 (MQ=255)
                                        tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGGCCGGCTGGGGGGTTTCTCCCGGGCCAAGTTGGTAGGGGGGCCGCCGGCAGGGatt                                                                                                                                                                                                   >  1:1461198/1‑239 (MQ=255)
                                        tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGGCCAGCTGGGGGGGTTTTCCGGCCGACAAGTTGTAAGAGGggcggcggccggagtt                                                                                                                                                                                                    >  1:1367377/1‑229 (MQ=255)
                                        tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGGGGGGTTTTCCGGCGGCAAGGTTTTTGGGGGGCCGccggcagggaat                                                                                                                                                                                                   >  1:1807705/1‑235 (MQ=255)
                                        tACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGGCAGTCTGGGGGGGTTTTCCGGCGGCAAGTTGTTAGGCGGGGCGCCGGGAGGGatt                                                                                                                                                                                                   >  1:1156450/1‑239 (MQ=255)
                                             gccAGTAACTACAGCATTCAGGGGTTCACCAAGGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                              <  2:14316/238‑1 (MQ=255)
                                             gccAGTAACTACAGCATGCAGTGGTTCACCAAAGCGATAGATGGAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTGAGGCGGGTGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGTTCGCAGGCAGGAATTAttgt                                                                                                                                                                                              <  2:1114599/238‑1 (MQ=255)
                                             aCCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                              <  1:937002/239‑1 (MQ=255)
                                             aCCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                              <  2:1385868/239‑1 (MQ=255)
                                               cAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGTGGGGCTGGTGGGTGTCTTCCGGCGCCAGGTGTTTGGGGGGGCGCCGGCGGGGATTAttgttt                                                                                                                                                                                            >  2:1133613/1‑238 (MQ=255)
                                               cAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGGCCGTCTGGGGGGGTTTTCCGGCGGCCAGATGTTGGGGGGGGCGCCGGCAGGGAGTAtttttg                                                                                                                                                                                            >  1:256261/1‑239 (MQ=255)
                                                  tAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGGCCGGCTTGGGGGGTTTTTCGGCGGCCAAGTGGTAGGGGGGGCGCCGGGGGGAATTATTTTTTGTa                                                                                                                                                                                        >  1:1289653/1‑239 (MQ=255)
                                                  tAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGGCCGTCTGGGGGGTTTTTCCGGCGGCAAGTTGTTAGGCGGGCCGCCGGGCGGGATTATTTTTTGtt                                                                                                                                                                                        >  1:595041/1‑238 (MQ=255)
                                                   aaCTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGAc                                                                                                                                                                                                                                  <  1:1628822/198‑1 (MQ=255)
                                                   aaCTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTa                                                                                                                                                                                        <  1:422158/238‑1 (MQ=255)
                                                        cAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGGCCGTCTGGGGAGGTTTTCCCGCGGCAAGTTGTTTGGCGGGGCGCCGGGAGGGATTATTTTTTGTaaaaaaa                                                                                                                                                                                  >  1:1660197/1‑238 (MQ=255)
                                                         aGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGGCAGGCTGGGGGGTTTTTCCCGCGGCCAGGTGTTAGGCGGGGCGCCGGGAGGGAGTATTTTTGGGgaaaaaaa                                                                                                                                                                                 >  1:868317/1‑237 (MQ=255)
                                                           cATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGGCAGTCTGGGGGGTTTTTCCCGCGGCAAGTTGTTTGGGGGGGCGCCGGCAGGGATTATTTTTTGTGAAAAAaaga                                                                                                                                                                               >  1:1125411/1‑239 (MQ=255)
                                                            aTTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAag                                                                                                                                                                                  <  1:639844/235‑1 (MQ=255)
                                                             ttCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                              <  1:1393725/223‑1 (MQ=255)
                                                              tCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGGCAGGCTGGGGGGGTTCTCCGGCGGCCAGTTGTTAGGCGGGCCGCCGGCCGGGAATATTGTTTGTAAAAAAGAgatga                                                                                                                                                                            >  1:984916/1‑239 (MQ=255)
                                                               cAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAgatgat                                                                                                                                                                           <  1:105754/239‑1 (MQ=255)
                                                                aGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGGGGCAGGCTGGGGGGGTTTTTCCGCGGCAAAGTTGTAGGCGGGGCGCCGGCCGGGAATATTGTTGGTGAAAAAAAAGTGATc                                                                                                                                                                          >  1:1643534/1‑239 (MQ=255)
                                                                  gggTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGGGGTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGGATTATTGTTGGTAAAAAAGAGATGAGcg                                                                                                                                                                         >  1:878576/1‑238 (MQ=255)
                                                                    tttacccAAGGGGATGGATGAGGCGAGACGGGTGGCGCGCGGAAAAGAGCTGGATGTTCCCGTGGGGACTGATTGAGGCAGTGGCTGGCTGGTGGTTCTTAGCCAGTAGGTTTTGCCGAAAGGGCAAATGCCGCAGGAGTTGATTGCGGCGGGCGGCGGTCTGGTGATTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGATTTATTGTTGGTAAAAAAGAATTGATcgcc                                                                                                                                                                       <  2:938786/233‑1 (MQ=255)
                                                                     ttCACCAAAGCTATGGATGAGGCGGAACTGGGGGCGCTCGGCAAAGAGTTGGTTGTTCCCGTGGTGACTGATTTGGGAGGTGGCTCGCTGGTCGTTCTTGGCCAGTCGGTTTTGCGGAAGGAGCCAATGCCGCAGGAGTTGATTGCGGGGGGCGTCAGTCTGGTGATTTTCTCCGGCGAAAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATcgccc                                                                                                                                                                      <  2:984916/238‑1 (MQ=255)
                                                                       cACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGGGGGGTTCTCCCGCGGCAAGTTGTTTGGCGGGCCGCCGGCAGGGATTATTTTTTGGAAAAAAAAAATGATagccccc                                                                                                                                                                    >  1:1675134/1‑236 (MQ=255)
                                                                       cACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGGCTGGGGGGGTTTTCCCGCGGCAAGTTGTTTGGGGGGGCGCCGGGAGGGATTATTTTTTGGGAAAAAAAGAAGATcgccccc                                                                                                                                                                    >  1:1078078/1‑236 (MQ=255)
                                                                                gATAGATGAAGGGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                          <  1:1359226/239‑1 (MQ=255)
                                                                                gATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTCCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                          <  1:191619/239‑1 (MQ=255)
                                                                                 atggaggagggggACCTGGGGGCGCTGGGCAAAGGGTTGGTTTTCCCCGTGGTGACTGTTTGAGGCTGTGGCTCGCTGGTCGATCTTAGCCAGTCGGGTTTGCCAAAAGGGCCAATGCCGCAGGGGTGGATTGCGGCGGGCGTCAGTCGGGTGAGTTTCTCGGGCGACAATTTGTTAGGCGGGCCGCAGGCAGGAATTATGGTTGGAAAAAAAGAGATGATCCCCCGCCTGCAAAGcc                                                                                                                                                          <  2:788164/227‑1 (MQ=255)
                                                                                 aTAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                          <  1:834916/238‑1 (MQ=255)
                                                                                 aTAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTGGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                          <  1:1253798/238‑1 (MQ=255)
                                                                                 aTAGATGAAGCGGAACTGGTGGCGCTCGGAAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGcc                                                                                                                                                          <  1:710487/238‑1 (MQ=255)
                                                                                  tAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCAc                                                                                                                                                        <  1:705071/239‑1 (MQ=255)
                                                                                     aTGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGAc                                                                                                                                                                                                                                  <  1:120527/164‑1 (MQ=255)
                                                                                              aaCTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCg                                                                                                                                             <  1:19535/238‑1 (MQ=255)
                                                                                                       cgcTCGGAAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTAcg                                                                                                                                    <  1:25784/238‑1 (MQ=255)
                                                                                                        gcTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGGCCGTCTGGGGGGTTTCTCCGGCGGCAAGTTGTTAGGCGGGGCGCCGGCAGGGATTATTGTTTGGGAAAAAAAGATGATCGCCCGCCCGCCAAGACACCCGCCGCAGCGGGCATTAcgc                                                                                                                                   >  1:576537/1‑238 (MQ=255)
                                                                                                           gggCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTGGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCgg                                                                                                                               <  1:1535315/238‑1 (MQ=255)
                                                                                                            ggaaaaGGGCGGGATGTTCCGGTGGTAGCTGATTTGGGAGGTGGCGGGCGGGTGGATTTTGGCAAGTGCGGTTTGCCGAAAGAGCCAATGCCGCAGGGTTTGATGGGGGGGGGCGTCAGTCTGGTGATTTTCCCGGGCGAAAAGTTGTTAGGGGGGCCGCAGGAAGGAATTATTGTTGGAAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCgg                                                                                                                               <  2:1687505/235‑1 (MQ=255)
                                                                                                            agacgggggtgggATGTTCCGTGGGTGACGGTTTTGGGCGGGGGCTCGCGGGTGGATTTTAGCCGGGGCGGTTTGCCAAAAGAGCCAAGGCCGAGGGAGTGGATTGGGGGGGGGGTAGGTCTGGTGATTTTCTCGGGCGAAAAGTTGTTAGGCGGGCCGCGGGAAGGAATTATTGTGGGAAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCgg                                                                                                                               <  2:1041206/227‑1 (MQ=255)
                                                                                                             gCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGa                                                                                                                              <  1:1773984/238‑1 (MQ=255)
                                                                                                              cAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCg                                                                                                                                                                                                               <  1:849706/157‑1 (MQ=255)
                                                                                                                          gtttCCGGTGGGGACTTGTTTGGGCAGGGGCTGGCGGGCGGTTCTGAGCCAGTACGGTTTGCGAAAAGAGCCAATGCCGCGGGAGTTGATTGCGGGGGGCGTCGGTCTGGTGGTTTTCTCCGGCGAAAAGTTGTTAGGCGGGCCGCAGGCAGGAATTTTTGTTGGTAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTcg                                                                                                                <  2:1660197/237‑1 (MQ=255)
                                                                                                                             tCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGaa                                                                                                                     >  1:1041206/1‑230 (MQ=255)
                                                                                                                                 gTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCt                                                                                                         <  1:978725/239‑1 (MQ=255)
                                                                                                                                   aGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTAGGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTgg                                                                                                       <  1:1446469/239‑1 (MQ=255)
                                                                                                                                           aTTTAGGCGGTGGCTCGCTGGTCGATCTTGGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGCCCGCAGGCAGGAATTATTGTTGTTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTAcg                                                                                                                                    <  1:788148/202‑1 (MQ=255)
                                                                                                                                           aTTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTAttgt                                                                                                                                                                                              <  1:352018/145‑1 (MQ=255)
                                                                                                                                              tAGGCGGTGTCTGGCTGGTCGATTTTAGCCAGTCCGGTTTGCCGAAAGACCCAAGGCCGCGGGAGTTGATGGGGGGGGGCGTCAGTCTGGTGAGTTTCTCCGGCGAAAAGTTGTTGGGCGGCCCGCAGGCAGGAATTATTGTTGGAAAAAAAGAGAGGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTg                                                                                           <  2:525802/239‑1 (MQ=255)
                                                                                                                                                   aGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTc                                                                                       <  1:407604/238‑1 (MQ=255)
                                                                                                                                                     tGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGGGGGTTTTTCCGGCGGCAAGTTGTTAGGCGGGGCGCAGGGAGGGATTTTTTTTGGTAAAAAAAAGATGATCGCCCGCCCGCCAAACCACCCGCCGaa                                                                                                                                                >  1:1550425/1‑180 (MQ=255)
                                                                                                                                                     gggCTCGCTGGTCGATCTTGGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCgg                                                                                                                               <  1:1153362/196‑1 (MQ=255)
                                                                                                                                                      ggCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGGGAGGTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGGATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCCGCAAAGCCACCCGCTGAAGCGTGCATTAcgc                                                                                                                                   >  1:939409/1‑192 (MQ=255)
                                                                                                                                                      ggCGCGCGGGGCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCt                                                                                      <  1:1231565/236‑1 (MQ=255)
                                                                                                                                                         gctgtGGTCATTCTTAGCGGGTACGGTTGGCGGAAAGAGCCAATGCCGCGGGGGTTGTTGGCGGGGGGCGTCAGTCTGGTGAGTTTCCCCGGCGAAAGGTTGTTAGGGGGCCCGCAGGAAGGAATTTTTGTTGGTAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCtttattt                                                                                <  2:1807705/235‑1 (MQ=255)
                                                                                                                                                          cgcgGGTCGATCTTAGCCAGTCCGGTTTGCCGAAAGAGCCAATGCCGCAGGGGTGGATGGCGGCGGGCGTCAGTCTGGTGAGTTTCCCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCtttattta                                                                               <  1:307561/239‑1 (MQ=255)
                                                                                                                                                          cGCTGGTCGATCTTAGCCAGTACGGTTGGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCtttattta                                                                               <  1:1018136/239‑1 (MQ=255)
                                                                                                                                                               gTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCg                                                                                                                                                                                                                                        <  1:1397310/85‑1 (MQ=255)
                                                                                                                                                                    tCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTCCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTc                                                                                       <  1:727696/221‑1 (MQ=255)
                                                                                                                                                                      ttAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGct                                                                   <  1:1251244/239‑1 (MQ=255)
                                                                                                                                                                        ggCCAGGGCGGTTGGCCGAAAGGGCCATGGCCGGAGGAGTTTATTGCGGCGGGCGGCGGTCTGGTGGGTTTCTCGGGCGACAAGTTGTTAGGCGGGCCGAAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGctc                                                                  <  2:1289653/237‑1 (MQ=255)
                                                                                                                                                                        aGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGctct                                                                 <  1:948118/239‑1 (MQ=255)
                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGGGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCCGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGAACCCCGCGGCGCTGGGAGCCACGTTGCGTCttt                                                                                    >  1:225211/1‑214 (MQ=255)
                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGGTAAAATGACCCTCGCGGCGCTGGGAGCCACGGTGCGGCTTTTTTTACACCCTGAAGCTCTGAGGGa                                                           >  1:770544/1‑239 (MQ=255)
                                                                                                                                                                              tACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGGGGGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCCGGCAGGGATTATTGTTGGTAAAAAAGAGATGAACGCCCGCCTGCAAAACCACCCGCTGAAGCGTGCATTAc                                                                                                                                     >  1:1104664/1‑166 (MQ=255)
                                                                                                                                                                               aCGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGAAAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGaa                                                          <  1:380084/239‑1 (MQ=255)
                                                                                                                                                                                      gCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGGGGGGTTCTCCCGCGACAAGGTGGTAGGCGGGGCGCAGGCAGGGAATATTTTTGGGTAAAAAAAGAAGGTCGCCCGCCCCCCAACCCACCCCCCGAAACGTGGATATAGCGCGGGGAAAAATAACCCCGCCGCGCCGGGAGCACCCGTTGCGCCTTTTTTATACCCCTAAGCCCTTAGCGAAAAAATa                                                    >  2:728184/1‑238 (MQ=255)
                                                                                                                                                                                      gCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGGGCGGCTGGTGGGGTTTTCCCGCGACCAGTTGTTAGGCGGGGCGCGGGCCGGGCGTATTTTTGGTGAAAAAAAGGTGATCGGCCGCCCGCCAAGGCACCCCCCGAAGGGGGGATTACGCGCGGGGAAAATAAACCCCGCGGCGCCGGGAGCCCCGGTGCGTCTTTATTTATACCCCGAAGCTCCGACTGAAAAATTaa                                                   >  2:362473/1‑238 (MQ=255)
                                                                                                                                                                                          aaaGGGCAAGTGCCGAGGGAGTGTGTTGCGGCGGGGGTAGTTTTGGTAGTTTTCTCCGGCGAAAGTTGGTTGGGGGGGCCGCAGGCAGGATTTTTTTTTGGAAAAAAAGAGAGGACCCCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGa                                                <  2:269501/238‑1 (MQ=255)
                                                                                                                                                                                           aaGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGCCTGGTGAGTTTCTCCGGCGACAAGTTGTTGGGGGGGCCGCAGGCAGGAATTATTGTGGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGAcc                                              <  1:1117041/239‑1 (MQ=255)
                                                                                                                                                                                              agCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTc                                                                                       <  1:233899/195‑1 (MQ=255)
                                                                                                                                                                                              agCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAG‑‑‑CCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACAcc                                                                           <  1:170846/204‑1 (MQ=255)
                                                                                                                                                                                               gcaatcgCCGCGGGGTTTGATGGGGGGGGGCGTCGGTCTGGGGGTTTTTGCGGGCAAAAGGTGGTTGGGGGGGCCGCAGGAAGGAATTATTGTTGGAAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTg                                           <  2:770544/232‑1 (MQ=255)
                                                                                                                                                                                                ccAGTGCCGCAGGAGTTGATTGCGGGGGGCGGCGGTCTGGTGAGTTTCTCCGGCGAAAATTTGTTAGGCGGGCCGCAGGCAGGAATTATGGTTGGTAAAAAAGAGATGGTCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTgcg                                         <  1:357194/239‑1 (MQ=255)
                                                                                                                                                                                                ccAAGGCCGCAGGAGTTGATGGGGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGAAAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTGGGAAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTgcg                                         <  1:1827176/239‑1 (MQ=255)
                                                                                                                                                                                                             aGTTGATTGCGGCGGGCTTCAGTCTGGTGATTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcg                             <  1:57672/238‑1 (MQ=255)
                                                                                                                                                                                                             aGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGAAAAAAAGAGATGACCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcg                             <  1:1374420/238‑1 (MQ=255)
                                                                                                                                                                                                                   ttGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGAAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcgcagcgc                       <  1:130877/238‑1 (MQ=255)
                                                                                                                                                                                                                   ttGCGGCGGGCGCCAGATTGGTGAGTTTCTCCGGAGACAAGTTGTTAGGCGCGCCGGAGGCAGAATTTATGGTTGAAAAAAAAGAGATGCCCACCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCcgcagcgca                      <  1:448743/239‑1 (MQ=255)
                                                                                                                                                                                                                               tcattttGGTGATTTTCGCGGGCGAAATGTTGTGAGGGGGGCCGCGGGAGGGAATTTTGGTTGCAAAAAAAGGCAGGATCCCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCa           <  2:256261/232‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCAAATc       <  1:1774390/237‑1 (MQ=255)
                                                                                                                                                                                                                                    cTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGACTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAACGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCAAATcc      <  1:1399511/238‑1 (MQ=255)
                                                                                                                                                                                                                                    aggtgtgtgtTTTGCCGGCGACAGGTTAGTAGGCGGGCGGTGGGCAGAAATTTTTGTGGAGAAAAAAAGAATGATCTCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCAAATc       <  2:1430953/234‑1 (MQ=255)
                                                                                                                                                                                                                                       ggTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCAAATCCAg    >  2:196662/1‑238 (MQ=255)
                                                                                                                                                                                                                                         tGATTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCAAATCCAGGc  <  1:792599/238‑1 (MQ=255)
                                                                                                                                                                                                                                         tGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCAAATCCAGGc  >  1:1519441/1‑238 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                
CGTCAGGCGGTGAATGAAAATACCGCACTGTTGATGAAAGTACATACCAGTAACTACAGCATTCAGGGGTTCACCAAAGCGATAGATGAAGCGGAACTGGTGGCGCTCGGCAAAGAGCTGGATGTTCCCGTAGTGACTGATTTAGGCAGTGGCTCGCTGGTCGATCTTAGCCAGTACGGTTTGCCGAAAGAGCCAATGCCGCAGGAGTTGATTGCGGCGGGCGTCAGTCTGGTGAGTTTCTCCGGCGACAAGTTGTTAGGCGGGCCGCAGGCAGGAATTATTGTTGGTAAAAAAGAGATGATCGCCCGCCTGCAAAGCCACCCGCTGAAGCGTGCATTACGCGCGGATAAAATGACCCTCGCGGCGCTGGAAGCCACGTTGCGTCTTTATTTACACCCTGAAGCTCTGAGTGAAAAATTACCGACCCTGCGCCTGCTTACCCGCAGCGCAGAGGTCATTCAAATCCAGGC  >  NC_007779/3879799‑3880268

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