Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 953,254 G→A P45S (CCA→TCA)  unknown ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1953,2540GA100.0% 142.1 / NA 48P45S (CCA→TCA) unknownhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/27);  total (21/27)

AGAGATTTGTTTTTCGGTTAAGGAAATAATGACTTCATTTGGCAACATTGACGATTGCAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGGGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAACACATTTTTATTAGGGATAAACGTTATAATAATTTTTAGGAGTTCATTATGCAAAAGCACAAAATTATCTTAGGT  >  1/953024‑953498
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
agagaTTTGTTTTTCGGTTAAGGAAATAATGACTTCATTTGGCAACATTGACGATTGTAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACt                                                                                                                                                                                                                                   >  1:417093/1‑250 (MQ=255)
     tttGTTTTTCGGTTAAGGAAATAATGACTTCATTTGGCAACATTGACGATTGTAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCAttt                                                                                                                                                                                                                             <  1:459793/251‑1 (MQ=255)
      ttGTTTTTCGGTTAAGGAAATAATGACTTCATTTGGCAACATTGACGATTGTAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCAtttt                                                                                                                                                                                                                            >  1:757393/1‑251 (MQ=255)
          ttttCGGTTAAGGAAATAATGACTTCATTTGGCAACATTGACGATTGTAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGAtt                                                                                                                                                                                                                        <  1:187544/251‑1 (MQ=255)
                 ttAAGGAAATAATGACTTCATTTGGCAACATTGACGATTGTAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAacac                                                                                                                                                                                                                  <  1:560458/250‑1 (MQ=255)
                                     tttGGCAACATTGACGATTGTAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGgat                                                                                                                                                                                             <  1:559307/251‑1 (MQ=255)
                                          cAACATTGACGATTGTAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGt                                                                                                                                                                                        <  1:130349/251‑1 (MQ=255)
                                          cAACATTGACGATTGTAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGt                                                                                                                                                                                        <  1:246653/251‑1 (MQ=255)
                                                              aTGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGgat                                                                                                                                                                                             >  1:787796/1‑226 (MQ=255)
                                                                                cAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGa                                                                                                                                                                                                                          >  1:594442/1‑179 (MQ=255)
                                                                                  aaGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGgatg                                                                                                                                                                                            >  1:355709/1‑207 (MQ=255)
                                                                                  aaGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAAttt                                                                                                                                                 >  1:601100/1‑250 (MQ=255)
                                                                                    gAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAacac                                                                                                                                                                                                                  >  1:11107/1‑183 (MQ=255)
                                                                                    gAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGgat                                                                                                                                                                                             <  1:237883/204‑1 (MQ=255)
                                                                                    gAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTaaaaa                                                                                                                                                                                   <  1:318885/214‑1 (MQ=255)
                                                                                                cGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATttgattg                                                                                                                                  <  1:679679/251‑1 (MQ=255)
                                                                                                   cAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTgaga                                                                                                                               >  1:843570/1‑251 (MQ=255)
                                                                                                      tataATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCGTTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTtc                                                                                                                                                         >  1:23563/1‑222 (MQ=255)
                                                                                                         aaTGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACat                                                                                                                         <  1:712148/251‑1 (MQ=255)
                                                                                                              cacTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGAATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGc                                                                                                                    <  1:740333/251‑1 (MQ=255)
                                                                                                                        ttCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTAc                                                                                                               >  1:199239/1‑246 (MQ=255)
                                                                                                                           aTTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAGAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAg                                                                                                       <  1:119281/251‑1 (MQ=255)
                                                                                                                           aTTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAg                                                                                                       <  1:315253/251‑1 (MQ=255)
                                                                                                                                      aaTGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTc                                                                                            <  1:775107/251‑1 (MQ=255)
                                                                                                                                                aTCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGt                                                                                                                                                                                        <  1:211828/149‑1 (MQ=255)
                                                                                                                                                aTCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCa                                                                                                                   >  1:166531/1‑218 (MQ=255)
                                                                                                                                                  cTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTtaat                                                                                <  1:865550/251‑1 (MQ=255)
                                                                                                                                                          aTTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATcac                                                                                                                           <  1:12426/200‑1 (MQ=255)
                                                                                                                                                                          ttGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACata                                                                                                                        <  1:578356/187‑1 (MQ=255)
                                                                                                                                                                                   cATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACat                                                                                                                         >  1:653820/1‑177 (MQ=255)
                                                                                                                                                                                        aCTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAacac                                                                          <  1:68001/219‑1 (MQ=255)
                                                                                                                                                                                            aaTTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGt                                                                                                                                                                                        <  1:574882/105‑1 (MQ=255)
                                                                                                                                                                                             aTTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAACACATTTTTATTAGGGATAAACGTTATAATAAtt                                           <  1:583927/245‑1 (MQ=255)
                                                                                                                                                                                                    gtgAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACat                                                                                                                         >  1:526431/1‑160 (MQ=255)
                                                                                                                                                                                                      gAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGt                                                                                                                                                                                        >  1:764440/1‑95 (MQ=255)
                                                                                                                                                                                                      gAAATAAAGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACatat                                                                                                                       >  1:778676/1‑160 (MQ=255)
                                                                                                                                                                                                           aaaGTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAat                                                                                                                                                                <  1:277465/114‑1 (MQ=255)
                                                                                                                                                                                                              gTGCTTAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCAt                                                                                                                  >  1:783400/1‑157 (MQ=255)
                                                                                                                                                                                                                   tAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTaaaaa                                                                                                                                                                                   >  1:682552/1‑87 (MQ=255)
                                                                                                                                                                                                                   tAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACatat                                                                                                                       >  1:117729/1‑147 (MQ=255)
                                                                                                                                                                                                                   tAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACata                                                                                                                        <  1:638547/146‑1 (MQ=255)
                                                                                                                                                                                                                   tAATTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAACACATTTTTATTAGGGATAAACGTTAtaataa                                             <  1:120556/221‑1 (MQ=255)
                                                                                                                                                                                                                     aTTATACTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGt                                                                                                                                                                                        >  1:66171/1‑80 (MQ=255)
                                                                                                                                                                                                                        ataCTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTctct                                                                                                                                                                    <  1:213653/97‑1 (MQ=255)
                                                                                                                                                                                                                        ataCTGTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAACACATTTTTATTAGGGATAAACGTTATAATAATTTTTAGGAGTTCATTATGCAAAAGCACAAAATTa          >  1:387161/1‑251 (MQ=255)
                                                                                                                                                                                                                             gTCCTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAACACATTTTTATTAGGGATAAACGTTATAATAATTTTTAGGAGTTCATTATGc                         <  1:321457/231‑1 (MQ=255)
                                                                                                                                                                                                                               ccTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAACACATTTttat                                                                  >  1:220976/1‑188 (MQ=255)
                                                                                                                                                                                                                                cTTTTGAGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAACACATTTTTATTAGGGATAAACGTTATAATAATTTTTAGGAGTTCATTATGCAAAAGCACAAAATTATCTTAGGt  <  1:342229/251‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
AGAGATTTGTTTTTCGGTTAAGGAAATAATGACTTCATTTGGCAACATTGACGATTGCAACAATGCTATGGCGATGTTGGCAAAGAAATAGAGTGCCGACAATATAATGACACTTTGTTTTTCATTTTGTATGAAATGGAAGCTATCTTTAATGATTTGGGTAAAGTTAATTGGCATTTCATCTACTCAATTTGATGTGAAATAAAGTGCTTAATTATACTGTCCTTTTGGGATAAAAACATAAAAAACTCCATTTTGATTTAACACAATGTTTTCCTCACGCAGGATGATGTAAAAATGATGTAGATTCTCTAAATTATTTTCTCTAATTTCCAAAAATTTGATTGAGATCACATATTGCATTACTCATAAAGTGGTATATTTCACCGCGCTTAATAAACACATTTTTATTAGGGATAAACGTTATAATAATTTTTAGGAGTTCATTATGCAAAAGCACAAAATTATCTTAGGT  >  1/953024‑953498

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