Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 889,247 C→T T117T (ACC→ACT)  unknown → Transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1889,2470CT98.7% 208.5 / ‑4.8 75T117T (ACC→ACT) unknownTranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (33/41);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

CGCGCAGTTTGACGGTCTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGCGAAAGAATTAACGCAACTGACGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACCTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTACAATTGGAAAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCGGTGCCTTTTCGCTTAGATTTTTTTGATGATGA  >  1/889004‑889474
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   
cgcgCAGTTTGATGGTTTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGc                                                                                                                                                                                                                              <  1:96925/251‑1 (MQ=255)
cgcgCAGTTTGATGGTTTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGc                                                                                                                                                                                                                              <  1:625003/251‑1 (MQ=255)
cgcgCAGTTTGATGGTTTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGc                                                                                                                                                                                                                              <  1:237532/251‑1 (MQ=255)
cgcgCAGTTTGATGGTTTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGc                                                                                                                                                                                                                              <  1:141273/251‑1 (MQ=255)
cgcgCAGTTTGATGGTTTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGc                                                                                                                                                                                                                              <  1:139884/251‑1 (MQ=255)
    cAGTTTGATGGTTTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACg                                                                                                                                                                                                                          >  1:601373/1‑251 (MQ=255)
    cAGTTTGATGGTTTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACg                                                                                                                                                                                                                          >  1:254216/1‑251 (MQ=255)
        ttGATGGTTTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGtt                                                                                                                                                                                                                      <  1:164593/251‑1 (MQ=255)
             ggtttAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGcc                                                                                                                                                                                                                  >  1:701136/5‑250 (MQ=255)
             ggtttAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGcc                                                                                                                                                                                                                  >  1:42624/5‑250 (MQ=255)
              gtttAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGccc                                                                                                                                                                                                                 >  1:463466/4‑250 (MQ=255)
              gtttAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCccg                                                                                                                                                                                                                >  1:661224/4‑251 (MQ=255)
              gtttAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCccg                                                                                                                                                                                                                >  1:470996/4‑251 (MQ=255)
              gtttAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTAAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCccg                                                                                                                                                                                                                >  1:450179/4‑251 (MQ=255)
                    ccGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGt                                                                                                                                                                                                            >  1:740596/1‑249 (MQ=255)
                         gtggtgACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAATGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTAtt                                                                                                                                                                                                      <  1:436793/250‑1 (MQ=255)
                            gtgACCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTgc                                                                                                                                                                                                   <  1:285476/250‑1 (MQ=255)
                               aCCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTgcagca                                                                                                                                                                                               >  1:596105/1‑251 (MQ=255)
                               aCCAATGACACCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTgcagc                                                                                                                                                                                                <  1:633776/250‑1 (MQ=255)
                                      acacCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGt                                                                                                                                                                                         <  1:581635/250‑1 (MQ=255)
                                      acacCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGt                                                                                                                                                                                         <  1:198065/250‑1 (MQ=255)
                                      acacCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGt                                                                                                                                                                                         <  1:151299/250‑1 (MQ=255)
                                      acacCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTa                                                                                                                                                                                        <  1:220330/251‑1 (MQ=255)
                                       cacCCGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAt                                                                                                                                                                                       >  1:355818/1‑251 (MQ=255)
                                         cccGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttg                                                                                                                                                                                     >  1:717776/1‑251 (MQ=255)
                                          ccGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttgt                                                                                                                                                                                    <  1:832728/251‑1 (MQ=255)
                                          ccGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttgt                                                                                                                                                                                    <  1:111340/251‑1 (MQ=255)
                                          ccGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttgt                                                                                                                                                                                    <  1:850166/251‑1 (MQ=255)
                                          ccGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttgt                                                                                                                                                                                    <  1:348723/251‑1 (MQ=255)
                                          ccGAAGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttgt                                                                                                                                                                                    <  1:17860/251‑1 (MQ=255)
                                             aaGTGCGGTGCGTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGGCTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttgttg                                                                                                                                                                                  <  1:541702/250‑1 (MQ=255)
                                               gtgcggtgcgTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGAtt                                                                                                                                                                               >  1:378112/1‑251 (MQ=255)
                                               gtgcggtgcgTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGAtt                                                                                                                                                                               >  1:85694/1‑251 (MQ=255)
                                                  cggtgcgTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATtaat                                                                                                                                                                            <  1:473061/251‑1 (MQ=255)
                                                   ggtgcgTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTgc                                                                                                                                                                                                   <  1:282417/227‑1 (MQ=255)
                                                     tgcgTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATtaataa                                                                                                                                                                          <  1:650386/250‑1 (MQ=255)
                                                       cgTTTGGAAAAAGAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAgg                                                                                                                                                                       <  1:207399/251‑1 (MQ=255)
                                                                   gAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGt                                                                                                                                                           >  1:725181/1‑251 (MQ=255)
                                                                   gAATTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGt                                                                                                                                                           >  1:836198/1‑251 (MQ=255)
                                                                    aaTTAAAGCAACTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGt                                                                                                                                                           <  1:631229/250‑1 (MQ=255)
                                                                             aaCTGGCGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAAtct                                                                                                                                                 <  1:465051/251‑1 (MQ=255)
                                                                                       ttCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttgt                                                                                                                                                                                    <  1:830520/206‑1 (MQ=255)
                                                                                              ttattGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGa                                                                                                                                >  1:277842/1‑251 (MQ=255)
                                                                                                                         gATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTgcag                                                                                                                                                                                                 >  1:629631/1‑159 (MQ=255)
                                                                                                                             gggAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCg                                                                                                                                                                 >  1:18087/1‑187 (MQ=255)
                                                                                                                              ggAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTa                                                                                                 <  1:665448/250‑1 (MQ=255)
                                                                                                                                  aCCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGa                                                                                            <  1:234515/251‑1 (MQ=255)
                                                                                                                                  aCCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGa                                                                                            >  1:34864/1‑251 (MQ=255)
                                                                                                                                  aCCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTG‑CCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGaa                                                                                           >  1:452732/1‑251 (MQ=255)
                                                                                                                                   cccTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTGTTGGaa                                                                                           <  1:424178/251‑1 (MQ=255)
                                                                                                                                    ccTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAAc                                                                                          <  1:871290/251‑1 (MQ=255)
                                                                                                                                       gCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTAttgt                                                                                                                                                                                    <  1:496464/158‑1 (MQ=255)
                                                                                                                                       gCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATg                                                                                       >  1:532448/1‑251 (MQ=255)
                                                                                                                                                  gTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGaa                                                                                                                                                    <  1:265387/179‑1 (MQ=255)
                                                                                                                                                  gTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATgcg                                                                            >  1:44288/1‑251 (MQ=255)
                                                                                                                                                             catcaAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATgcggtg                                                                         <  1:827429/243‑1 (MQ=255)
                                                                                                                                                                    atatTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATgcggtgcgcggtg                                                                  <  1:218237/243‑1 (MQ=255)
                                                                                                                                                                     tattatTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGc                                                                                                                          <  1:550759/186‑1 (MQ=255)
                                                                                                                                                                     tattatTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAAc                                                                                          >  1:680547/1‑218 (MQ=255)
                                                                                                                                                                     tattatTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACCTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTgg                                                         >  1:869456/1‑251 (MQ=255)
                                                                                                                                                                             ccGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGtt                                                  >  1:25659/1‑250 (MQ=255)
                                                                                                                                                                                  cgcTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGAt                                            >  1:865375/1‑251 (MQ=255)
                                                                                                                                                                                  cgcTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGAt                                            >  1:824628/1‑251 (MQ=255)
                                                                                                                                                                                    cTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATgg                                          >  1:732360/1‑251 (MQ=255)
                                                                                                                                                                                      tGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTCCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGa                                        <  1:217955/251‑1 (MQ=255)
                                                                                                                                                                                             tttGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCgg                                 <  1:577446/251‑1 (MQ=255)
                                                                                                                                                                                                gTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGcc                                                                                                                                          <  1:307390/143‑1 (MQ=255)
                                                                                                                                                                                                     agttagtttaGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCa                                                                                                     >  1:315661/10‑175 (MQ=255)
                                                                                                                                                                                                       ttGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCGGTGCCTTTTCg                       >  1:285857/1‑251 (MQ=255)
                                                                                                                                                                                                        tGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGGAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCGGTGCCTTTTCGc                      <  1:204782/251‑1 (MQ=255)
                                                                                                                                                                                                        tGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCGGTGCCTTTTCGc                      <  1:765932/251‑1 (MQ=255)
                                                                                                                                                                                                              aGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTAt                                                                                                >  1:322180/1‑171 (MQ=255)
                                                                                                                                                                                                                    aaCAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCGGTGCCTTTTCGCTTAGAttttttt          >  1:559771/1‑251 (MQ=255)
                                                                                                                                                                                                                      cAAGGCATTTTCATTTTACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCGGTGCCTTTTCGCTTAGATTTTTTtga        <  1:125823/251‑1 (MQ=255)
                                                                                                                                                                                                                            attttcattttACCTGTCGCCACTTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTGCAATTGGAGAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCGGTGCCTTTTCGCTTAGATTTTTTtgatgatga  >  1:93556/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   
CGCGCAGTTTGACGGTCTAACCGTGGTGGTGACCAATGACACCCGAAGTGCGGTGCGTTTGGCGAAAGAATTAACGCAACTGACGCCTTCCGAATTATTGCAAAACGTGACTTTTTTCCCCGATTGGGAAACCCTGCCTTATGATAGTTTTTCGCCTCATCAAGATATTATTTCCGCGCGCTTGAGCGCTTTGTTTCATTTGCAACAGCGCAAACAAGGCATTTTCATTTTACCTGTCGCCACCTTAATGCAACGGGTTTGCCCGCCGTCTTATTTGCAGCATAACGTATTGTTGATTAATAAAGGCGATCGTTTTGTCATCGAATCTCTTCGCCTACAATTGGAAAAAGCGGGTTATCGTGCGGTGGAGCAGGTATTGGAACATGGCGAATATGCGGTGCGCGGTGCGTTGTTGGATTTGTTCCCGATGGGAAGTGCGGTGCCTTTTCGCTTAGATTTTTTTGATGATGA  >  1/889004‑889474

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑