Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 947,942 T→C I108V (ATT→GTT)  unknown ← 3‑deoxy‑manno‑octulosonate cytidylyltransferase (EC 2.7.7.38)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1947,9420TC98.4% 174.0 / NA 64I108V (ATT→GTT) unknown3‑deoxy‑manno‑octulosonate cytidylyltransferase (EC 2.7.7.38)
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (31/32);  total (31/33)
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

TGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAATGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAA  >  1/947699‑948172
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
tGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGa                                                                                                                                                                                                                                 >  1:564111/1‑251 (MQ=255)
        tgACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCgg                                                                                                                                                                                                                         <  1:201377/251‑1 (MQ=255)
         gACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCg                                                                                                                                                                                                                          >  1:825075/1‑249 (MQ=255)
              aaaTATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTtcgtc                                                                                                                                                                                                                   <  1:542001/251‑1 (MQ=255)
                atatTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTtcgtcg                                                                                                                                                                                                                  >  1:314633/1‑250 (MQ=255)
                atatTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGc                                                                                                                                                                                                                 <  1:201993/251‑1 (MQ=255)
                 tatTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGcc                                                                                                                                                                                                                <  1:458575/251‑1 (MQ=255)
                    tCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTg                                                                                                                                                                                                             <  1:107544/251‑1 (MQ=255)
                       gTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAAt                                                                                                                                                                                                          >  1:62773/1‑251 (MQ=255)
                         aTTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGt                                                                                                                                                                                                        <  1:61940/251‑1 (MQ=255)
                             aGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGa                                                                                                                                                                                                     >  1:72584/1‑250 (MQ=255)
                              gCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGt                                                                                                                                                                                                        >  1:84992/1‑246 (MQ=255)
                                          ccGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCAtc                                                                                                                                                                                        <  1:797493/250‑1 (MQ=255)
                                                 gTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAAttt                                                                                                                                                                                            >  1:174263/1‑239 (MQ=255)
                                                              aaTTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGttttt                                                                                                                                                                    >  1:222238/1‑250 (MQ=255)
                                                              aaTTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGttttt                                                                                                                                                                    >  1:633870/1‑250 (MQ=255)
                                                               aTTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTcc                                                                                                                                                                  <  1:388352/251‑1 (MQ=255)
                                                                 tGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTccac                                                                                                                                                                <  1:312701/251‑1 (MQ=255)
                                                                    tCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTccac                                                                                                                                                                <  1:122648/248‑1 (MQ=255)
                                                                                 gAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAAATTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGt                                                                                                                                                >  1:358297/1‑251 (MQ=255)
                                                                                         gcgcGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAg                                                                                                                                                                         >  1:502934/1‑218 (MQ=255)
                                                                                                  aaaTAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATtgtg                                                                                                                               >  1:365693/1‑251 (MQ=255)
                                                                                                   aaTAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTgg                                                                                                                              >  1:317196/1‑251 (MQ=255)
                                                                                                    aTAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGt                                                                                                                             >  1:466936/1‑251 (MQ=255)
                                                                                                    aTAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGt                                                                                                                             >  1:133523/1‑251 (MQ=255)
                                                                                                         gCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAAGTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGaa                                                                                                                                    <  1:130103/239‑1 (MQ=255)
                                                                                                                 ccGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCa                                                                                                                 <  1:235030/250‑1 (MQ=255)
                                                                                                                     ttttGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGcac                                                                                                            <  1:721049/251‑1 (MQ=255)
                                                                                                                                  cACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGTGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGcac                                                                                                            <  1:879398/238‑1 (MQ=255)
                                                                                                                                  cACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGccc                                                                                                  <  1:200924/248‑1 (MQ=255)
                                                                                                                                   aCTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACCTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGaa                                                                                              >  1:106069/1‑251 (MQ=255)
                                                                                                                                    cTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGaaa                                                                                             >  1:490006/1‑251 (MQ=255)
                                                                                                                                      tttACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAAtt                                                                                           >  1:290586/1‑251 (MQ=255)
                                                                                                                                                             ttCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAAtt                                                                    <  1:538875/251‑1 (MQ=255)
                                                                                                                                                              tCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTTCCACCaaat                                                                    >  1:255042/1‑248 (MQ=255)
                                                                                                                                                              tCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGGCACCCAtt                                                                    >  1:55726/1‑250 (MQ=255)
                                                                                                                                                                  tCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAATGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTa                                                               <  1:445441/251‑1 (MQ=255)
                                                                                                                                                                  tCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTa                                                               <  1:350340/251‑1 (MQ=255)
                                                                                                                                                                   cgcgTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCacac                                                                                                          <  1:468809/207‑1 (MQ=255)
                                                                                                                                                                   cgcgTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTAt                                                              <  1:895635/251‑1 (MQ=255)
                                                                                                                                                                           caATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGGTTCCACCAATTCGTTTTCGGTGGCg                                                      >  1:351418/1‑251 (MQ=255)
                                                                                                                                                                              tttttACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCCTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAGTTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGccccccactccttttcgggggcggta                                                   >  1:443080/1‑227 (MQ=255)
                                                                                                                                                                                    ccGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGttttc                                                             >  1:242670/1‑232 (MQ=255)
                                                                                                                                                                                    ccGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGAcg                                             <  1:139232/251‑1 (MQ=255)
                                                                                                                                                                                    ccGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGAcg                                             <  1:149433/251‑1 (MQ=255)
                                                                                                                                                                                                 cATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCgg                                <  1:199111/251‑1 (MQ=255)
                                                                                                                                                                                                   tGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTc                                                                   <  1:827147/214‑1 (MQ=255)
                                                                                                                                                                                                   tGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGGCGCGGTCGGCACCGGAt                              <  1:45459/251‑1 (MQ=255)
                                                                                                                                                                                                   tGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGAt                              <  1:603107/251‑1 (MQ=255)
                                                                                                                                                                                                   tGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATACACACTTCCGCCCCGAAATTTCTTGCGGCGGTTTCCACCACTTCGTTTTCGGGGGCGGTAATGACGCCGTCGGcccccggt                              >  1:448830/1‑244 (MQ=255)
                                                                                                                                                                                                   tGTTCACAGGGTATTTTGCAAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGAt                              <  1:805891/251‑1 (MQ=255)
                                                                                                                                                                                                    gTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTTTCGGTGGCGGTAATGACGCGGGCGGCACCGGAtt                             >  1:821645/1‑251 (MQ=255)
                                                                                                                                                                                                    gTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAAATCCTTTTTGGTGGGGATAATGAAGCGGGCGGCACCGGAtt                             >  1:868895/1‑251 (MQ=255)
                                                                                                                                                                                                           gggTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGTTATCGGGGGCGATAATGACGCGGGCGGCACCGGGTTGttttgg                      >  1:343817/1‑248 (MQ=255)
                                                                                                                                                                                                            ggTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGc                      <  1:823464/250‑1 (MQ=255)
                                                                                                                                                                                                             gTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCCCCAATTCGTTATCGGGGGGGGTAATTACGCGGGCGGCACCGGGTTGTTTgggtt                    >  1:806276/1‑248 (MQ=255)
                                                                                                                                                                                                             gTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGTTATCGGGGGGGGTAATGACGCGGGCGGCACCGGGTTGttttggtt                    >  1:781112/1‑246 (MQ=255)
                                                                                                                                                                                                                         gTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCCCCAATTCGTTATCGGGGGGGGTAATGACGCGGTCGGCACCCGGTTGTTTTGCTTGTTTTAAAACgg        >  1:759217/1‑250 (MQ=255)
                                                                                                                                                                                                                         gTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCCATTCGTTATCGGTGGCGATAATGACGCGGGCGGCACCGGATTGTTTGGCTTGTTCaaaaaac         >  1:90981/1‑248 (MQ=255)
                                                                                                                                                                                                                         gTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCaaaaa           <  1:603961/248‑1 (MQ=255)
                                                                                                                                                                                                                         gTTTTCCGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAAc          <  1:415414/249‑1 (MQ=255)
                                                                                                                                                                                                                               cGCCACTTGTTTCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGaa  <  1:501495/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aCTTGTTGCACGATAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGtt     <  1:868176/244‑1 (MQ=255)
                                                                                                                                                                                                                                               aTAACGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGa                                               <  1:278985/190‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
TGCCGATGTGACGTAAATATTCGGTATTCAGCCGTAACTTTGCCGTATCGTCCGTTTGTGCGAATTGATCGCGATCCCAGGGAATCACGGCGCGGGAGAAATAAAGCACATAGCCGTTTTTGTCGGTCAGCACTTTTACCACGTTGGGATTAAACAATTCTTTCGCGTCATCAATTTTTACCGCAAGGCTTGCCATGTTCACAGGGTATTTTGCCAGGTTTTCCGCCACTTGTTTCACGATAATGGGCGGAATGAGCGGTTCGTCGCCTTGAATGTTGACGATAATTTCATCGTCGGGCAGCACCAGTTTTTCCACCACTTCCGCCAGGCGTTCCGTGCCGGAATTGTGGTTTTTGGCAGTCATGCACACTTCCGCCCCGAAATTTCTTGCGGCGGTTGCCACCAATTCGTTATCGGTGGCGATAATGACGCGGTCGGCACCGGATTGTTTGGCTTGTTCAAAAACGTGTTGAA  >  1/947699‑948172

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