Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,450,568 A→T intergenic (+142/‑150) unknown → / → unknown hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,450,5680AT95.2% 109.1 / ‑0.4 42intergenic (+142/‑150)unknown/unknownhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (18/22);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

ACATTTAAAAAAGTAGCCGGTTACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGTAAATCCAAAAATGATACTAGGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑TTTTCCGCTTGC‑GGAAAATGCAGGGCGAAGGCTTGCAGCCCGGTAAACAATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCGCGCTTAAATTTAAGCGAGAA  >  1/1450457‑1450807
                                                                                                               |                                                                                                                                                                                                                                                 
atctaattttgttaagacggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTc                                                                                                                     <  1:735103/217‑1 (MQ=255)
atctaattttgttaagacggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttga                                                                                                       >  1:897988/22‑248 (MQ=255)
     attttgttaagacggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttgagcaag                                                                                                  >  1:47361/17‑243 (MQ=255)
               gacggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAAc                                                                                                                                                                                  <  1:598161/156‑1 (MQ=255)
               gacggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGAc                                                                                                                                                       <  1:356024/183‑1 (MQ=255)
               gacggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTa                                                                                                                   <  1:645917/219‑1 (MQ=255)
               gacggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttgagcaagctcaaaaatg                                                                                        >  1:631333/7‑233 (MQ=255)
                  gtatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGAt                                                                                                                                                                                                                                  >  1:469852/4‑111 (MQ=255)
                  ggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAATAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttgagcaagctcaaaaatgtag                                                                                     >  1:92644/4‑230 (MQ=255)
                  ggatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGAc                                                                                                                                                       >  1:126598/4‑186 (MQ=255)
                   gatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTg                                                                                                                                             >  1:710266/3‑195 (MQ=255)
                   gatACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑tttttgagcaagctcaaaaatgtagg                                                                                    >  1:60513/3‑229 (MQ=255)
                    atACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCAt                                                                                                                                                                                                        >  1:216121/2‑135 (MQ=255)
                     tACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTg                                                                                                                                             <  1:573317/193‑1 (MQ=255)
                     tACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTAc                                                                                                                                      >  1:740583/1‑200 (MQ=255)
                          gtgtTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTg                                                                                                                                                                                                                             >  1:501948/1‑108 (MQ=255)
                                    atcTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCaaaa                                                                                                                                                                    <  1:858884/155‑1 (MQ=255)
                                           acacTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAAc                                                             <  1:828510/251‑1 (MQ=255)
                                           acacTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTCTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAAc                                                             <  1:639394/251‑1 (MQ=255)
                                            cacTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAAcc                                                            <  1:119271/251‑2 (MQ=255)
                                             acTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAAcca                                                           <  1:283417/251‑3 (MQ=255)
                                              cTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGAc                                                                                                                                                       <  1:700016/158‑1 (MQ=255)
                                                  taAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTc                                                      <  1:21991/251‑1 (MQ=255)
                                                    aaaCGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAAc                                                             >  1:882242/1‑242 (MQ=255)
                                                      aCGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATg                                                  <  1:871823/251‑1 (MQ=255)
                                                       cgcgACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGt                                                 <  1:796521/251‑1 (MQ=255)
                                                         cgACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTa                                               <  1:624089/251‑1 (MQ=255)
                                                                gTCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCggg                                         >  1:315815/1‑250 (MQ=255)
                                                                 tCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGatat                                                                                                                                                 >  1:119379/1‑145 (MQ=255)
                                                                 tCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGt                                                 <  1:92661/241‑1 (MQ=255)
                                                                 tCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGc                                        <  1:251811/250‑1 (MQ=255)
                                                                 tCACAGATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCa                                       >  1:321878/1‑251 (MQ=255)
                                                                      gATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCa                                                                                                                                                                                                         <  1:120165/84‑1 (MQ=255)
                                                                      gATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCaaaaa                                                                                                                                                                   >  1:435792/1‑122 (MQ=255)
                                                                      gATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGa                                                                                                                                            <  1:472060/145‑1 (MQ=255)
                                                                      gATGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCatttt                                   >  1:791560/1‑250 (MQ=255)
                                                                       aTGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACt                                                                                                                         <  1:595573/163‑1 (MQ=255)
                                                                       aTGTAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCattttta                                 >  1:18069/1‑251 (MQ=255)
                                                                          tAAGTGTAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAGCCATGCTCGATGTTATTCGGGCatttttattt                              <  1:721810/251‑1 (MQ=255)
                                                                             gtgtAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTg                                                                                                                                             >  1:645560/1‑137 (MQ=255)
                                                                               gtAAAAAATGGCGATGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAg                                                                                                                                               <  1:602292/133‑1 (MQ=255)
                                                                                            aTGTGATTGTTCTTGGTGGTCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTgcgcgtacctaaatt            >  1:212832/1‑241 (MQ=255)
                                                                                               tgattgttcttggtggtCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTgcgcgtacc                  <  1:57031/231‑5 (MQ=255)
                                                                                               tgattgttcttggtggtCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCGTACCTAAATTTaa         >  1:137407/12‑251 (MQ=255)
                                                                                                    gttcttggtggtCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTgcgcgtacctaaat             <  1:275550/236‑10 (MQ=255)
                                                                                                    gttcttggtggtCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCGTACCTAAATTTAATCgag    >  1:717953/7‑251 (MQ=255)
                                                                                                    gttcttggtggtCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCGTACCTAAATTTAATCgag    >  1:322230/7‑251 (MQ=255)
                                                                                                      tcttggtggtCTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCGTACCTAAATTTAATCGAGaa  >  1:203027/5‑251 (MQ=255)
                                                                                                               |                                                                                                                                                                                                                                                 
ACATTTAAAAAAGTAGCCGGTTACAGGTGTTAATCAATCTCCGACACTTATAAAACGCGACATTGTCACAGATGTAAGTGTAAAAAATGGCGATGTGATAGTTTTGGGTGGACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGTAAATCCAAAAATGATACTAGGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑TTTTCCGCTTGC‑GGAAAATGCAGGGCGAAGGCTTGCAGCCCGGTAAACAATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCGCGCTTAAATTTAAGCGAGAA  >  1/1450457‑1450807

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