Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,803,523 A→C intergenic (+21623/+31861) unknown → / ← unknown tRNA‑Asp‑GTC/tRNA‑Lys‑CTT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,803,5230AC95.7% 131.7 / ‑5.7 47intergenic (+21623/+31861)unknown/unknowntRNA‑Asp‑GTC/tRNA‑Lys‑CTT
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (24/21);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATCATCGTGAGTTTTATCAAATTCACGAATCGCATCTTGAAGATAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCACATGAGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGT  >  1/1803279‑1803763
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
atcatcGTGAGTTTTATCAAATTCACGAATCGCATCTTGAAGATAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAatat                                                                                                                                                                                                                                             >  1:89030/1‑251 (MQ=255)
          gTTTTATCAAATTCACGAATCGCATCTTGAAGATAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCAc                                                                                                                                                                                                                                    >  1:750095/1‑250 (MQ=255)
                 cAAATTCACGAATCGCATCTTGAAGATAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATAc                                                                                                                                                                                                                            <  1:833772/251‑1 (MQ=255)
                           aaTCGCATCTTGAAGATAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAAc                                                                                                                                                                                                                 <  1:416054/251‑1 (MQ=255)
                            aTCGCATCTTGAAGATAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACa                                                                                                                                                                                                                >  1:535618/1‑251 (MQ=255)
                                      gAAGATAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGt                                                                                                                                                                                                      >  1:93244/1‑251 (MQ=255)
                                          aTAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAg                                                                                                                                                                                                  <  1:163003/251‑1 (MQ=255)
                                              gTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTAt                                                                                                                                                                                               <  1:293807/250‑1 (MQ=255)
                                                ttGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCta                                                                                                                                                                                            <  1:54215/251‑1 (MQ=255)
                                                   ccATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCtata                                                                                                                                                                                          <  1:265736/250‑1 (MQ=255)
                                                    cATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATc                                                                                                                                                                                        <  1:251399/251‑1 (MQ=255)
                                                    cATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATc                                                                                                                                                                                        <  1:657651/251‑1 (MQ=255)
                                                    cATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATc                                                                                                                                                                                        <  1:601993/251‑1 (MQ=255)
                                                      tttATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCAt                                                                                                                                                                                      <  1:497826/251‑1 (MQ=255)
                                                                                        acaacaTTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATc                                                                                                                                                                                        >  1:351754/1‑215 (MQ=255)
                                                                                                   aaagagccaggCTTGAGACTTTAATAGGATCCCCCCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGTTAAAACAAATAACAACAAAAATTAAAATACTGCTATCCCCCACCACAACAATCATATGCAAATATAATTTTCACCAATAAACATTTT‑AAAACAAACCTTCCGTTGAGTATCTATATCATAAAAAGAAGTTTTCCTGAAACCGCCATTTTTTGAAAGTGTGagat                                                                                                                                         <  1:288336/241‑1 (MQ=255)
                                                                                                   aGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATAc                                                                                                                                                                                                                            >  1:573048/1‑169 (MQ=255)
                                                                                                    gCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAAt                                                                                                                                                                                                                         >  1:589437/1‑171 (MQ=255)
                                                                                                     aaaTACCAGTTTGAGCCGTTAATAGGTTTCCCCCAAAAAGTAACGAAAGTAGTGTTTTTTTCATATTGCTCCTCCTTTTTATCGTTTAAACAAATAACAACAAAAATTAAAATACTCCTATCCCCCCCCAAAACAATCATATGCGAATAAAAATTTCCCCAATATAAAATTT‑AAAACAAATCGTCCGTTGATTATCTATATAATAAAAAAAAGTTTTCCTGCAACCGCCATTTTTTGAAAGTGTGagat                                                                                                                                         <  1:868719/248‑1 (MQ=255)
                                                                                                             gCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTAtt                                                                                                                               >  1:222559/1‑251 (MQ=255)
                                                                                                                        tAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTc                                                                                                                                                                                                         <  1:91433/166‑1 (MQ=255)
                                                                                                                             ggATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATc                                                                                                                                                                                        <  1:19916/178‑1 (MQ=255)
                                                                                                                                  cccccacAAAGAAACGCAAGTATTTTCTTTTTAATAATGCTCCTCCCTTTAATCGATTAAACAAAAAAAAACAAAAATTAAAATTCTGCTATACCCCGCCACAACAATCATATGCGAAAATAAATTTCACCAATATA‑AATTTGGAAAAAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAAATTTCCTACAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAact                                                                                                          <  1:9863/248‑1 (MQ=255)
                                                                                                                                  ccaccaAAAAGAAACGAAAGTAGTTTTTTTTTCATAATGCTCCCCCTTTTTATCGGTTAAACGAAAAAAAAAAAAAATAAAAATACTGCTATACCCCCCCACAACAATCATATGAGAATATAAATTCCACAAAAAAACAATTT‑GAAACAAATCGTCCGTTGAGTATTTATATAATAAAAAAAAGATTTCCGGAAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAact                                                                                                          <  1:540316/251‑1 (MQ=255)
                                                                                                                                          aaGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCaaa                                                                                                  >  1:522467/1‑251 (MQ=255)
                                                                                                                                          aaGTAACCAAAGTGGTGTCTTTTTCATATTGCTACTCCCTTTTATCGTTTAAACGAAAAACAAAAAAAATAAAAATTCTGCTATCCCCCCCCACAACAATCATAGGCGAATATAAATTTCACAAATATACAATTA‑AAAACAAATGGCCCGGTGATTATCAATATCATAAAAAAAAGTTTCCCGGCAACCGCCATTTTTTGAAAGTGTGagat                                                                                                                                         <  1:187939/212‑1 (MQ=255)
                                                                                                                                            gTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTAttt                                                                                                                              >  1:755433/1‑221 (MQ=255)
                                                                                                                                              aaCGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCta                                                                                                                                                                                            <  1:407293/157‑1 (MQ=255)
                                                                                                                                                    agtagtGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTaa                                                                                        >  1:144002/1‑251 (MQ=255)
                                                                                                                                                            cTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGcc                                                                                                                                                             <  1:582233/174‑1 (MQ=255)
                                                                                                                                                            cTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAAtatttta                                                                                >  1:78656/1‑251 (MQ=255)
                                                                                                                                                                              cccTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTAt                                                                                                                                                                                               >  1:126792/1‑122 (MQ=255)
                                                                                                                                                                              cccTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCtata                                                                                                                                                                                          <  1:451216/127‑1 (MQ=255)
                                                                                                                                                                                cTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCGATATACAATTT‑GAAACAAATCGTcc                                                                                                                                                                                                        >  1:622417/1‑111 (MQ=255)
                                                                                                                                                                                                          aaaaaTTAAAATACTGCTATCCCCCACCACAACAATAATATGCAAAAATAAATTCCACAAATATAAAATTT‑AAAAAAAATCGTCCGTTGAGTATCTATATCATAAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTc                                  <  1:373157/251‑1 (MQ=255)
                                                                                                                                                                                                             aaTTGAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTAttt                                                                                                                              >  1:674773/1‑156 (MQ=255)
                                                                                                                                                                                                                 gAAATACTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATc                                                                                                                                                                                        >  1:53562/1‑94 (MQ=255)
                                                                                                                                                                                                                  aaaTACTTCTATCCCCCACCACAAAAATCATAGGCAAATATAAATTTCACCAATATACAATTT‑GAAAAAAATCGTCCTTTGAGTATCTATATCATAAAAAGAAGTTTTCCTGAAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGaa                          <  1:812807/251‑1 (MQ=255)
                                                                                                                                                                                                                       cTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGagatag                                                                                                                                       >  1:232114/1‑137 (MQ=255)
                                                                                                                                                                                                                       cTGCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATc                                                                                                                                    >  1:719460/1‑140 (MQ=255)
                                                                                                                                                                                                                         gCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAAt                                                                                                                                                                                                                         >  1:619305/1‑54 (MQ=255)
                                                                                                                                                                                                                         gCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTa                   >  1:846799/1‑251 (MQ=255)
                                                                                                                                                                                                                         gCTATACCCCACCACAACAATCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTa                   >  1:867489/1‑251 (MQ=255)
                                                                                                                                                                                                                                       caacaaTCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCtata                                                                                                                                                                                          >  1:270431/1‑70 (MQ=255)
                                                                                                                                                                                                                                       caacaaTCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTa     >  1:798291/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aaaTCATATGAGAATATAAATTTCACAAATATACAATTT‑GAAACAAATCGTCCGTGGAGTATCTATATCATAAAAAGAAGTTTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGt  <  1:792644/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaTCATATGCGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCta                                                                                                                                                                                            >  1:321544/1‑64 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
ATCATCGTGAGTTTTATCAAATTCACGAATCGCATCTTGAAGATAGGTTTGCCATTTATCGGAGAAGGCAACCATCGGAGCGCCAATTACAACATTTTTAGCAATAACAGCTTGAGACGTTAATAGGATTCCACCACAAAGTAACGCAAGTAGTGTCTTTTTCATAATGCTACTCCCTTTTATCGGTTAAACGAATAACAACAAAAATTGAAATACTGCTATACCCCACCACAACAATCACATGAGAATATAAATTTCACCAATATACAATTT‑GAAACAAATCGTCCGTTGAGTATCTATATCATTAAAAGAAGATTTCCTGCAACCGCCATTTTTTGAAAGTGTGAGATAGATCTTATTTTTTGGAAGGAAAATTTAACTACTGCAAAAATAAACTAAATATTTTATTTTTAGCAAAAAAATAGGAGTCCAAAGACTCCTATGTTCATGCTCAATACGAATTCGTTACGCAACCTCTTCTACGT  >  1/1803279‑1803763

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