Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,775,854 G→A G610G (GGC→GGT)  unknown ← Type III restriction‑modification enzyme helicase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,775,8540GA98.7% 204.9 / NA 74G610G (GGC→GGT) unknownType III restriction‑modification enzyme helicase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/38);  minor base T (0/1);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAATAACTGTGCAAATGCGCCACCATCTGTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAGCCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAGCGCTTCGCTTGGACGGGCGGAATATTTGGC  >  1/1775607‑1776093
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
aCAAATAACTGTGCAAATGCGCCACCATCTGTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCg                                                                                                                                                                                                                                              >  1:236393/1‑251 (MQ=255)
  aaaTAACTGTGCAAATGCGCCACCATCTGTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCg                                                                                                                                                                                                                                              >  1:87461/1‑249 (MQ=255)
     tAACTGTGCAAATGCGCCACCATCTGTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCgctg                                                                                                                                                                                                                                           >  1:571264/1‑249 (MQ=255)
             cAAATGCGCCACCATCTGTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTcgc                                                                                                                                                                                                                                 <  1:571847/251‑1 (MQ=255)
                         caTCTGTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACaa                                                                                                                                                                                                                     >  1:501926/1‑251 (MQ=255)
                            cTGTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAAGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACgg                                                                                                                                                                                                                  <  1:333885/251‑1 (MQ=255)
                              gTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTa                                                                                                                                                                                                                >  1:716111/1‑251 (MQ=255)
                                    ccAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATgcgc                                                                                                                                                                                                          <  1:533044/251‑1 (MQ=255)
                                    ccAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATgcgc                                                                                                                                                                                                          <  1:678673/251‑1 (MQ=255)
                                              aaGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGa                                                                                                                                                                                                <  1:875054/251‑1 (MQ=255)
                                              aaGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGa                                                                                                                                                                                                <  1:654342/251‑1 (MQ=255)
                                                                cAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGa                                                                                                                                                                              >  1:897225/1‑251 (MQ=255)
                                                                  gATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGa                                                                                                                                                                                                >  1:183138/1‑231 (MQ=255)
                                                                  gATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGa                                                                                                                                                                                                >  1:368411/1‑231 (MQ=255)
                                                                  gATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGAtt                                                                                                                                                                            >  1:297976/1‑251 (MQ=255)
                                                                     aaTCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGAttt                                                                                                                                                                           >  1:174760/1‑249 (MQ=255)
                                                                      aTCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGa                                                                                                                                                                                                <  1:99177/227‑1 (MQ=255)
                                                                              cTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTaa                                                                                                                                                                                                               <  1:343636/204‑1 (MQ=255)
                                                                                             gcgcAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAat                                                                                                                                                                                                     <  1:192963/199‑1 (MQ=255)
                                                                                                   aGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATc                                                                                                                                            >  1:372245/1‑250 (MQ=255)
                                                                                                    gACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCa                                                                                                                                          >  1:699821/1‑251 (MQ=255)
                                                                                                     aCGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCa                                                                                                                                          >  1:89476/1‑250 (MQ=255)
                                                                                                      cGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCAcc                                                                                                                                        <  1:397373/251‑1 (MQ=255)
                                                                                                      cGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTCCTGTCGTTTTGGCGACAATTTGTGCCATATCCAcc                                                                                                                                        <  1:544927/251‑1 (MQ=255)
                                                                                                        atatAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGcc                                                                                                                                                 >  1:511726/1‑240 (MQ=255)
                                                                                                        atatAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCt                                                                                                                                       >  1:309916/1‑250 (MQ=255)
                                                                                                        atatAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCt                                                                                                                                       >  1:384177/1‑250 (MQ=255)
                                                                                                          ataGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCAt                                                                                                                                    >  1:863714/1‑251 (MQ=255)
                                                                                                          ataGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACATCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGcc                                                                                                                                                 >  1:362271/1‑238 (MQ=255)
                                                                                                             gTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGt                                                                                                                                                                    <  1:515741/216‑1 (MQ=255)
                                                                                                              tCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGAt                                                                                                                                 <  1:310067/250‑1 (MQ=255)
                                                                                                              tCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATc                                                                                                                                <  1:315997/251‑1 (MQ=255)
                                                                                                               cTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGa                                                                                                                                  <  1:471548/248‑1 (MQ=255)
                                                                                                                ttCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCg                                                                                                                              <  1:444140/251‑1 (MQ=255)
                                                                                                                                     aCCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTcgccgc                                                                                                                                                                                                                              >  1:808040/1‑134 (MQ=255)
                                                                                                                                               gaaGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATGTCGCGTAGGCCGGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAATCCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATGTCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAgg                                                                                                <  1:270476/248‑1 (MQ=255)
                                                                                                                                               gCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCcagcag                                                                                                                       >  1:876465/1‑227 (MQ=255)
                                                                                                                                               gCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCAt                                                                                                           <  1:257554/239‑1 (MQ=255)
                                                                                                                                                cAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGgt                                                                                               <  1:219345/250‑1 (MQ=255)
                                                                                                                                                cAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGgt                                                                                               <  1:30447/250‑1 (MQ=255)
                                                                                                                                                       gctgctTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATcc                                                                                                                                           <  1:821135/199‑1 (MQ=255)
                                                                                                                                                        ctgctTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGaa                                                                                       <  1:735001/250‑1 (MQ=255)
                                                                                                                                                             ttGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCt                                                                                 <  1:330204/251‑1 (MQ=255)
                                                                                                                                                               gCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATcc                                                                                                                                           <  1:155038/191‑1 (MQ=255)
                                                                                                                                                               gCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTtg                                                                               >  1:498639/1‑251 (MQ=255)
                                                                                                                                                                  cGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGGCCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTtgatg                                                                            <  1:416270/251‑1 (MQ=255)
                                                                                                                                                                    ttGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGa                                                                                                                                                                                                <  1:436344/133‑1 (MQ=255)
                                                                                                                                                                    ttGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGgt                                                                                               >  1:748845/1‑230 (MQ=255)
                                                                                                                                                                    ttGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTgatgat                                                                          <  1:407030/251‑1 (MQ=255)
                                                                                                                                                                      ggTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGAtt                                                                                                                                                                            >  1:492896/1‑151 (MQ=255)
                                                                                                                                                                      ggTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGAttt                                                                        >  1:192956/1‑251 (MQ=255)
                                                                                                                                                                      ggTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGAttt                                                                        >  1:51827/1‑251 (MQ=255)
                                                                                                                                                                      ggTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGAttt                                                                        >  1:398525/1‑251 (MQ=255)
                                                                                                                                                                       gTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGAtttt                                                                       >  1:61126/1‑251 (MQ=255)
                                                                                                                                                                            cAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTgg                                                                   >  1:13713/1‑250 (MQ=255)
                                                                                                                                                                              gATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGa                                                                                        >  1:697675/1‑227 (MQ=255)
                                                                                                                                                                              gATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCgg                                                                >  1:830925/1‑251 (MQ=255)
                                                                                                                                                                                     aaTGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCAt                                                         <  1:809520/251‑1 (MQ=255)
                                                                                                                                                                                             tttCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCCGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAAt                                                 >  1:91361/1‑251 (MQ=255)
                                                                                                                                                                                               tCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATaa                                               <  1:658466/251‑1 (MQ=255)
                                                                                                                                                                                                  cgTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGc                                            <  1:409371/251‑1 (MQ=255)
                                                                                                                                                                                                        cccGGTTGCAGGTGCGGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCGCCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTgaga                                      <  1:3065/251‑1 (MQ=255)
                                                                                                                                                                                                        cccGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCcagc                                                                                                                         <  1:188041/168‑1 (MQ=255)
                                                                                                                                                                                                                  ggTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAgcg                              >  1:238017/1‑249 (MQ=255)
                                                                                                                                                                                                                   gTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCcagc                                                                                                                         <  1:620783/157‑1 (MQ=255)
                                                                                                                                                                                                                             ccgACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAGCGCTTCGCTTGGAc                  <  1:770761/250‑1 (MQ=255)
                                                                                                                                                                                                                             ccgACATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAGCGCTTCGCTTGGAc                  <  1:35139/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTcc                                                                                  >  1:115654/1‑182 (MQ=255)
                                                                                                                                                                                                                                 cATCCAGTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATaa                                               >  1:156049/1‑217 (MQ=255)
                                                                                                                                                                                                                                       gTTTAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAGCGCTTCGCTTGGACGGGCGGAatat       >  1:342943/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tttAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAGCGCTTCGCTTGGACGGGCGGAatat       <  1:119666/250‑1 (MQ=255)
                                                                                                                                                                                                                                        tttAAACGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAGCGCTTCGCTTGGACGGGCGGAatat       <  1:520353/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCGGGTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAGCGCTTCGCTTGGACGGGCGGAATATTTggc  <  1:761621/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ggTGAAAACCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATcc                                                                                                                                           <  1:349575/110‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
ACAAATAACTGTGCAAATGCGCCACCATCTGTCCTGCCAAATCGTAAAGCAAATCGGCTTGGGTCAGATAATCAATATCTTCAAAATCAATCAGCGCAAAGACGATATAGTCTTCGGGGCGTTTTTCCTTCAAACCCAATTCGGCAGAAAGGCTGCTTTGCTCGTTGGTATGCAGATTGTGAATGGTGATTTCGCGTACGCCCGGTTGCAGGTGCAGACCGCCGACATCCAGTTTAAACGGGTGAAAGCCGCTGCTGACTTCGCCGCTTGGCACAACGGTAATGCGCGGAATATCGATGGTTTGCGCCGCCATGATTTCTGTCGTTTTGGCGACAATTTGTGCCATATCCACCTCATGATCCGCCAGCAGTTCGCCTTGCATCGGTTGCAAGGTCGTCTGAACTTCCTTGATGATTTTTTGGCGGATTTCATCGCTCAATAACGCTTGAGAAGTCAGCGCTTCGCTTGGACGGGCGGAATATTTGGC  >  1/1775607‑1776093

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