Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 155,826 G→A intergenic (‑240/+59) unknown ← / ← unknown Type I restriction‑modification system, DNA‑methyltransferase subunit M (EC 2.1.1.72)/Transcriptional regulator, MerR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1155,8260GA100.0% 126.5 / NA 47intergenic (‑240/+59)unknown/unknownType I restriction‑modification system, DNA‑methyltransferase subunit M (EC 2.1.1.72)/Transcriptional regulator, MerR family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/24);  total (23/24)

TTTGTAGTTGGCGGCGTCGAGTTGCTGGCGGAGACGGTCGGCGGAAGACCAAAGTTTTTCGTCGAGGTCGTTTAGGAAGGCTTGTTGTGTGGCTGAATTTTCAACATAATCTTGGTGGCTTCTGGCTTCTGGCTTCTGGCTTCTGGCTTCTGCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCGTGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCGGGTTGTCTTGTAAT  >  1/155593‑156068
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
gattaagccaagcacaatgtgtttgtagttggcggcgtcgagttgctggcggagacggtcggcggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAAtt                                                                                                                                                                                                                                   >  1:693288‑M2/153‑251 (MQ=255)
      gccaagcacaatgtgtttgtagttggcggcgtcgagttgctggcggagacggtcggcggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATAttt                                                                                                                                                                                                                             >  1:593252‑M2/147‑251 (MQ=255)
       ccaagcacaatgtgtttgtagttggcggcgtcgagttgctggcggagacggtcggcggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATAttt                                                                                                                                                                                                                             <  1:831175‑M2/105‑1 (MQ=255)
            cacaatgtgtttgtagttggcggcgtcgagttgctggcggagacggtcggcggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACtt                                                                                                                                                                                                                       >  1:372546‑M2/141‑251 (MQ=255)
              caatgtgtttgtagttggcggcgtcgagttgctggcggagacggtcggcggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGCCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTaa                                                                                                                                                                                                                     >  1:393285‑M2/139‑251 (MQ=255)
                        gtagttggcggcgtcgagttgctggcggagacggtcggcggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATg                                                                                                                                                                                                           >  1:770658‑M2/129‑251 (MQ=255)
                                             ctggcggagacggtcggcggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATgagag                                                                                                                                                                                                                                              <  1:80567‑M2/88‑1 (MQ=255)
                                                gcggagacggtcgacggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTtcactc                                                                                                                                                                                    <  1:271097‑M2/146‑1 (MQ=255)
                                                       cggtcggcggaagaccaaagtttttcgtcgaggtcgtttaggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGc                                                                                                                                                                            <  1:365051‑M2/154‑1 (MQ=255)
                                                                                        tcgtttaggaaggcttgttgtgtggctgaattgtcaacataatcttggtggcttctggcttcggCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCGACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTAc                                                                                                                                           <  1:387328‑M2/187‑1 (MQ=255)
                                                                                              aggaaggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCAc                                                                                                                                                                                                  <  1:389694‑M2/132‑1 (MQ=255)
                                                                                                   ggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGc                                                                                                                                                                                                               <  1:309300‑M2/119‑1 (MQ=255)
                                                                                                   ggcttgttgtgtggctgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTAcc                                                                                                                                          <  1:337573‑M2/188‑1 (MQ=255)
                                                                                                                  tgaattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGc                                                                                                                 <  1:141309‑M2/213‑1 (MQ=255)
                                                                                                                    aattttcaacataatcttggtggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCa                                                                                                                <  1:605221‑M2/214‑1 (MQ=255)
                                                                                                                                         ggcttctggcttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTAcc                                                                                                                                          >  1:219605‑M2/16‑203 (MQ=255)
                                                                                                                                                  cttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTAt                                                                                 <  1:452043‑M2/245‑1 (MQ=255)
                                                                                                                                                   ttctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTAtt                                                                                <  1:393135‑M2/246‑1 (MQ=255)
                                                                                                                                                    tctgCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTAttt                                                                               <  1:132309‑M2/247‑1 (MQ=255)
                                                                                                                                                             aaTTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTa                                                                      >  1:493640/1‑251 (MQ=255)
                                                                                                                                                                     cTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTAtt                                                                                <  1:315106/233‑1 (MQ=255)
                                                                                                                                                                         gTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTAttt                                                                               <  1:363010/230‑1 (MQ=255)
                                                                                                                                                                                gAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAAtg                                                           >  1:503296/1‑243 (MQ=255)
                                                                                                                                                                                       taaaAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGt                                            >  1:259283/1‑251 (MQ=255)
                                                                                                                                                                                         aaaaaaTATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGt                                                                            <  1:758634/217‑1 (MQ=255)
                                                                                                                                                                                           aaaaTATACCTATTTTATAGTAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTAcc                                                                                                                                          >  1:269093/1‑153 (MQ=255)
                                                                                                                                                                                            aaaTATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGCGTTGCGCGGTGAGtttttt                                       >  1:159062/1‑251 (MQ=255)
                                                                                                                                                                                                    cTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAgg                                                                   >  1:303609/1‑215 (MQ=255)
                                                                                                                                                                                                      aTTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGc                                                                                                                                                             >  1:488063/1‑123 (MQ=255)
                                                                                                                                                                                                      aTTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGt                                                                                            >  1:464863/1‑188 (MQ=255)
                                                                                                                                                                                                      aTTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTAt                                                                                 >  1:614756/1‑199 (MQ=255)
                                                                                                                                                                                                      aTTTTATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCgg                             >  1:702887/1‑251 (MQ=255)
                                                                                                                                                                                                       ttttATAGGAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCgg                             >  1:788178/1‑250 (MQ=255)
                                                                                                                                                                                                               gAATTGTGCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGc                                                                                       <  1:566515/184‑1 (MQ=255)
                                                                                                                                                                                                                    gtgCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCgg                                                 <  1:353906/217‑1 (MQ=255)
                                                                                                                                                                                                                    gtgCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCAc                                    <  1:566808/230‑1 (MQ=255)
                                                                                                                                                                                                                    gtgCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCgg               >  1:642345/1‑251 (MQ=255)
                                                                                                                                                                                                                    gtgCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCgg               >  1:271480/1‑251 (MQ=255)
                                                                                                                                                                                                                    gtgCATGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCgg               >  1:830057/1‑251 (MQ=255)
                                                                                                                                                                                                                         tGTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCGGGTTGt          >  1:318061/1‑251 (MQ=255)
                                                                                                                                                                                                                          gTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTgttgcg                     <  1:727702/239‑1 (MQ=255)
                                                                                                                                                                                                                          gTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCGGGTTGt          >  1:312530/1‑250 (MQ=255)
                                                                                                                                                                                                                          gTTTAAAAGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCGGGTTGt          >  1:391561/1‑250 (MQ=255)
                                                                                                                                                                                                                                aaGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCGGGTTGTCTTGTaa   <  1:318394/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCGGGTTGTCTTGTaa   <  1:876366/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaGATAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGCTGaaaccggttttttg                                 <  1:195443/221‑14 (MQ=255)
                                                                                                                                                                                                                                   aTAAGCATGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCGGGTTGTCTTGTAAt  <  1:867525/249‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
TTTGTAGTTGGCGGCGTCGAGTTGCTGGCGGAGACGGTCGGCGGAAGACCAAAGTTTTTCGTCGAGGTCGTTTAGGAAGGCTTGTTGTGTGGCTGAATTTTCAACATAATCTTGGTGGCTTCTGGCTTCTGGCTTCTGGCTTCTGGCTTCTGCATAAAATTTGTCCTTTGTAAAGTGAGTAAATAAAAAAATATACCTATTTTATAGGAATTGTGCATGTTTAAAAGATAAGCGTGAGAGATTTTTTAATTATATTTCTACTTAAAATTGCAATGAAAAATCACCGCACTTTTCACTCACTCAAGCCCTCCAAAATCGAGCATTCCGGGCAGTCGTTACCTTGGCAGTTGTCGTGCCAACGTTGCAGTTCGGTGACCATTTTTTGTAACGCCTGTATTTGGTCGTTTAAGGTTTCAATGTGTTGCGCGGTGAGTTTTTTCACTTCGCGGCTGTTGCGTTGCGGGTTGTCTTGTAAT  >  1/155593‑156068

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