Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,362 G→A intergenic (+496/+260) rtcR → / ← glpG DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,3620GA90.5% 249.2 / 10.8 97intergenic (+496/+260)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base A (57/29);  total (62/35)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

CGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATT  >  NC_000913/3560116‑3560606
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
cGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGAt                                                                                                                                                                                                                                                     >  1:692082/1‑251 (MQ=255)
  gaccTGGGAGATAAAACCGAGCGCCGCTTCGCCAAGGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTGTTCCCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGGGCGTGCGCTACCGCTGCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATgc                                                                                                                                                                                                                                                   <  2:315725/248‑1 (MQ=255)
   aaCTGGGAGATAAAATCGAGCTTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATgcg                                                                                                                                                                                                                                                  <  1:544062/251‑1 (MQ=255)
   aaCTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATgcg                                                                                                                                                                                                                                                  >  2:124871/1‑251 (MQ=255)
   aaCTGGGAGATAAAATCGAGCGGCGCTTCGCCAATGATCCCGCCATCGCAGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAACTCTTCTTTTACCATCACCGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATgcg                                                                                                                                                                                                                                                  >  1:728165/1‑251 (MQ=255)
      ggggAAATAAAATCGATGGTCGTTTCCCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGctc                                                                                                                                                                                                                                               <  2:755684/250‑1 (MQ=255)
       gggAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCCCTTTGCGGACGATGCGctct                                                                                                                                                                                                                                              >  2:107289/1‑251 (MQ=255)
         gagaTAAAATCGAGCGTCTCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCCCCGGCGAGAATGATGCGAAAATCTTCTTTTACCACCAACGT‑GTTAGCACCGTTGAGATTGTGGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTTCCGCTTCCGGCTTGGTGCCGATATCTATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCttt                                                                                                                                                                                                                                            <  1:374399/251‑1 (MQ=255)
            aaaaaaTCGATCGTCGATCCGCCAATGAACCCGCCATCGCGGCTGCGAAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTGACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCTCTTCCGGCGTGGTGCCGATATCGATAAACAGCTTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTttc                                                                                                                                                                                                                                         <  1:371820/249‑1 (MQ=255)
             aaaaaTCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTttct                                                                                                                                                                                                                                        <  1:271471/250‑1 (MQ=255)
                    gAGCGTCGCTTCGCCAGTGATACCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTGTTACCATCCACGT‑GTTAGCGACGTTGAGATTGTTGGTCACATTGCGCAAATTGCGGAGATTGAGCAGTGCGTGCACTACCGCTTCTGGCGTGGTCCCGATCTCGATAAACAGCGTCGTGCCATTGGGGGTTTGC‑TCAGCCACTTTTAGGACGATGCGCTCTTTttctactggcg                                                                                                                                                                                                                                 >  2:489603/1‑244 (MQ=255)
                    gAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGCTAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCt                                                                                                                                                                                                                                 >  2:676933/1‑251 (MQ=255)
                        gTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGt                                                                                                                                                                                                                             >  2:668513/1‑251 (MQ=255)
                          cGCTTCGCCAATGATCCCGCCAGCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGGGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTgg                                                                                                                                                                                                                           <  1:6135/251‑1 (MQ=255)
                           gCTTCGCCAATGATCCAGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGc                                                                                                                                                                                                                          >  2:665209/1‑251 (MQ=255)
                              tCGCCAATGATCCCGCCATAGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGTTTTGC‑TCCGCCACTTTGCGGCCGATGCGCTCTTTTTCTTCGGTATTGGTGGCCtt                                                                                                                                                                                                                       >  2:720300/1‑251 (MQ=255)
                                cccaaggaccccGCCATCGCGGTGGCGCAATTCGCCACCGGCGAAAAAGAGTCAAAAACCTTCTTTGACCATCAACGT‑GTTAGCAACGTTGGGATTGTGGGTGACAAGGCGCAAATTGCTGTGATTGAGCAGTGCGGGCGCTACCGCTTCCGGCTTGTGGCCGATATCGATAAACAGCGTCGAGCCATTGGGGGTTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGc                                                                                                                                                                                                                     <  2:444002/242‑1 (MQ=255)
                                  cAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGAAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGTCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCg                                                                                                                                                                                                                    >  1:452559/1‑250 (MQ=255)
                                   aaTGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGAt                                                                                                                                                                                                                  >  1:769828/1‑251 (MQ=255)
                                    aTGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGGGCAGTGCGTGCGCTACCGCTTCCGGCGTGGGGCCGATATCGATAAACAGCGTCGAGCCTTTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                                 <  1:468356/251‑1 (MQ=255)
                                    aTGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                                 <  1:317428/251‑1 (MQ=255)
                                      gATCCCGCCATCGCGGCTGCGTAATTCGCCAACGGCGAGAATGAGGCGAAAATCTTCTTTTACCAGCAACGG‑GTTAGCAACGTTGAGATTGTGGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTTCCGCTGCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGt                                                                                                                                                                                                               <  1:518808/251‑1 (MQ=255)
                                         cccGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGcc                                                                                                                                                                                                            >  1:614316/1‑251 (MQ=255)
                                              cATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACggg                                                                                                                                                                                                       >  2:766677/1‑250 (MQ=255)
                                               aTCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCg                                                                                                                                                                                                      >  2:606967/1‑251 (MQ=255)
                                                 cgcgGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGGTAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCgtg                                                                                                                                                                                                    >  2:152486/1‑251 (MQ=255)
                                                 cgcgGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCgtg                                                                                                                                                                                                    <  1:235010/251‑1 (MQ=255)
                                                    ggCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGtt                                                                                                                                                                                                  >  1:743698/1‑250 (MQ=255)
                                                       ggctgaacgtgCCACCGGCGAGAAAGATGCGAAAAGCTTCGTTTACCATCAACGT‑TTTAGCAACGTTGAGATTGGTGGGGACAATGCGCAAATGGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTTCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAcc                                                                                                                                                                                              <  1:550558/241‑1 (MQ=255)
                                                        gcgacactCGCCACTGGCGAGACTGAGTCGAAATTCTTCTTTTACCAT‑AACGTTGTTAGCAACGTGGAGATTGGTGGTGACAATGCGCAAATGGCTGTGATTGAGAAGTGCGGGCTCTACCGCTTCCGGCGTGGGGCCGATATCTATAAACAGTGGCGCGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAc                                                                                                                                                                                               <  2:43128/242‑1 (MQ=255)
                                                         cGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTGCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCAGCCACTTTTCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCgt                                                                                                                                                                                            >  2:62932/1‑250 (MQ=255)
                                                            aatgcGCCAACGGCGAGAATGAGGCGAAAATCTTTTTTTACCATCAACGT‑GTTAGAAACGTTGAGATTGTTGGGGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTCCCGGCGTGTGGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACt                                                                                                                                                                                         <  1:759691/247‑1 (MQ=255)
                                                            aaTTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACt                                                                                                                                                                                         >  2:384455/1‑251 (MQ=255)
                                                             aTTCGCCACCGGCGAGAAGGATGCGAAAATCTTCTTGTACCATCAACGG‑GTTAGCAACGTTGAGATGGGTGGTGACAATGCGCAAATTGCTGGGATTGAGCAGTGCGTGCGCTACCGCTGCGGGCGTGGTGCCGATATCGAAAAACAGCGTCGAGCCATTGGGGATTGGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTg                                                                                                                                                                                        <  2:192726/251‑1 (MQ=255)
                                                               tCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAATGT‑GTTAGCAACGTTGTGATTGTTGGTTACAATGCGCAAATTTCTGTGATTGAGCAGTTCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTgg                                                                                                                                                                                       <  1:430035/250‑1 (MQ=255)
                                                                    cccggcaaaaaTGATGAGAAAATCTTTATTTCCCCTCAACGT‑TTTATCAACGTGGAGTTTGTTGGTGACAATGCGCAATTTGCGTTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCTTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCa                                                                                                                                                                                 <  2:563301/242‑1 (MQ=255)
                                                                    cccGGCGAGAATGATGCGAAAATCTTCTTTTACCACCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGTTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCa                                                                                                                                                                                 <  1:668150/250‑1 (MQ=255)
                                                                    aCCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGCTTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTTGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGc                                                                                                                                                                                  >  2:667701/1‑250 (MQ=255)
                                                                            gaATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGcgccgc                                                                                                                                                                          >  1:562615/1‑250 (MQ=255)
                                                                               tgatgCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTATCCGGACTGGCAGGCCGCGCCGCACCg                                                                                                                                                                      >  2:426333/1‑251 (MQ=255)
                                                                                 atgCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGcc                                                                                                                                                                    >  1:361066/1‑251 (MQ=255)
                                                                                     gAAAATCTTCTTTTACCATCAACGA‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCatg                                                                                                                                                                 >  1:587872/1‑250 (MQ=255)
                                                                                        aaTCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATgg                                                                                                                                                             >  2:24717/1‑251 (MQ=255)
                                                                                        aaTCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATgg                                                                                                                                                             >  1:301283/1‑251 (MQ=255)
                                                                                            ttcttTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cgc                                                                                                                                                         >  2:567660/1‑248 (MQ=255)
                                                                                                      gtCAAGGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCCAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGACATTGGGGATTTTC‑TCCGCCACTTTGCGGACGATGCGATCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGAGCCACGGCGTGATCACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAgggt                                                                                                                                             >  2:178552/2‑249 (MQ=255)
                                                                                                      aTCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGtt                                                                                                                                             >  1:443679/1‑251 (MQ=255)
                                                                                                        cAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGtttt                                                                                                                                           >  1:435467/1‑251 (MQ=255)
                                                                                                        cAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGtttt                                                                                                                                           >  1:767843/1‑251 (MQ=255)
                                                                                                         aaCGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGAACCGCCATGATGG‑CGCA‑GGATCAGGTTTTg                                                                                                                                          >  2:636339/1‑251 (MQ=255)
                                                                                                         aaCGC‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGAAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATCTCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCAGGGTGGCCTTGCGCTCGTGCCACGCCGTGTTAACCGACCAGGACGGCAGCGCAGCACCGCCCTGATGG‑CGCA‑GGCTCATGTTTTg                                                                                                                                          >  2:479096/1‑251 (MQ=255)
                                                                                                           acgcgtTAGCAAAGTTGAAATTGTTTGTGACAATGCGCAAATTGCTGAGTTGGGGCAGTGCGTGCGCTACCGCTTCATGCGCGGTGCCGATATCGATAAACAGCGTGGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACTGTGCGATATTTTTCTTCGGTCTGGGAGGCCTTGCGATCGTGCCACACCGTGTTAACCGAACTGGACGGCAACGCCGCACCCCCTCGACGG‑C‑CACGGAACAGGtttttc                                                                                                                                         >  2:538614/5‑249 (MQ=255)
                                                                                                                     aacggggagatttgtgGTGCCAATGCGCAAATTGCTGTGATGGAGCAGGTCGGGCGCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGGGATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                                 <  2:299378/157‑1 (MQ=255)
                                                                                                                     aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                                 >  1:84533/1‑171 (MQ=255)
                                                                                                                     aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                                 >  1:299378/1‑171 (MQ=255)
                                                                                                                     aaCGTTGAGATGTTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATAGCGATAAACAGCGGCGAGCCTGTGGGGATGTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                                 <  2:84533/171‑1 (MQ=255)
                                                                                                                        gTTGAGATTGTTGGTGACAATGCGCAACTTGCTGTGATTGCGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGAGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTTCTCCGCCAGCTCAt                                                                                                                            >  1:322054/1‑251 (MQ=255)
                                                                                                                                  ttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGGTCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGCTGG‑CGCC‑GGATCAGGTTTTGCTCCGCCAGCACATTGAGGTcgcg                                                                                                                  >  2:108771/1‑251 (MQ=255)
                                                                                                                                    gttGACAATGCGCAAATTGTTGTGTTTGAGCAGTGCGGGCGCTACCGTTCCCGGCGGGGTGCAGATATCGATAAACGGCGTCGAGCCATGGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                                <  1:567660/249‑1 (MQ=255)
                                                                                                                                    ggTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                                >  1:60698/1‑251 (MQ=255)
                                                                                                                                         cAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                            >  1:544035/1‑250 (MQ=255)
                                                                                                                                          aaTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCCGGTTTTTCTCCGCCAGCTCATTGAGGTCgcggcgcat                                                                                                            >  2:70614/1‑246 (MQ=255)
                                                                                                                                              cgcAAAATGCTGTGATTGAGCAATGCGTGCTCTACCGCTTCCGGAGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGAGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCATTGCGATCGTGCAACGGCGTGTTAACCGAACAGGAAGGCAGCACAGCACCGCCATGATGGAG‑CA‑GACTCAGGTTTTGATACGCAATCTCATTAGGGTCGCGGCAACTAGTCtt                                                                                                      >  2:565344/1‑250 (MQ=255)
                                                                                                                                                  aaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATTGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCGCCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                  >  1:767604/1‑251 (MQ=255)
                                                                                                                                                  aaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                  >  1:7596/1‑251 (MQ=255)
                                                                                                                                                     tGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGACCTGGAAGGCAGCGCCGCACCGCAATGACGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                               >  2:88596/1‑251 (MQ=255)
                                                                                                                                                               gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                     >  1:550512/1‑251 (MQ=255)
                                                                                                                                                                      ggtggcgcTCCCGCTCCCGGGGTGGTGCCGATACCGAAAACCGGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                              <  1:107289/247‑1 (MQ=255)
                                                                                                                                                                       cgtgcgCTCCCGCTCCCGGCGTGGTGCCGATATAGATAAACAGCGTCGGGCCTTTGGGGTTTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                             <  1:665209/251‑1 (MQ=255)
                                                                                                                                                                           cgcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                                         <  1:24717/251‑1 (MQ=255)
                                                                                                                                                                                    ttCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTACCAGCTCTTCGGTAc                                                                >  2:736917/1‑251 (MQ=255)
                                                                                                                                                                                    ttCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCTGCTCTTCGGtcc                                                                >  2:278744/1‑249 (MQ=255)
                                                                                                                                                                                     cccGGGGTGGTCCCGATACCGATAAACCGCGTCGAGCCTTTGGGGTTTTGC‑CCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                               <  1:384455/250‑1 (MQ=255)
                                                                                                                                                                                      ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAACTGCTCTACCAGCTCTTCGGTACTg                                                              >  2:538320/1‑251 (MQ=255)
                                                                                                                                                                                       cGGCGTGGGGCCGATATCGATAAACAGGGTCGGGCCTTGGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                             <  1:241166/251‑1 (MQ=255)
                                                                                                                                                                                          cGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGGGAAATGCTCTACCAGCTCTTCGGTACTGCCAt                                                          >  2:99109/1‑251 (MQ=255)
                                                                                                                                                                                                gCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                    >  1:49515/1‑251 (MQ=255)
                                                                                                                                                                                                 ccGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                                   >  1:387636/1‑251 (MQ=255)
                                                                                                                                                                                                   gATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                                 >  1:129959/1‑251 (MQ=255)
                                                                                                                                                                                                      atCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGAACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGATCCTTGAGGTCGCGGCGAATAGTCTGAGGGATGACGGAGAAATGCTCTACCAACTCTTCGGTACTGACATAACCctgcattt                                              >  2:609641/1‑247 (MQ=255)
                                                                                                                                                                                                       tCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                             >  1:617617/1‑251 (MQ=255)
                                                                                                                                                                                                        cataaaaCAGCGTGGGGCCATGGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                             <  1:667701/246‑1 (MQ=255)
                                                                                                                                                                                                         gATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                                                                                         >  1:3120/1‑173 (MQ=255)
                                                                                                                                                                                                          aTAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                          >  2:708799/1‑251 (MQ=255)
                                                                                                                                                                                                          aTAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                          >  1:434707/1‑251 (MQ=255)
                                                                                                                                                                                                           aaaccagcgtggggccttgggggATTTGC‑TCCGCCACTTGGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                         <  1:441242/232‑1 (MQ=255)
                                                                                                                                                                                                                  gtgtgggcgcatGGGGGTTTGCCTCCGCCACTTTGGGGCCGAGGCGCTCTTTTTCTGCGGTCTGGGTGTCCTGGCGACCGGGCCACGGGGTGTTAACCGAACTGGAAGGCAGGGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                                   <  2:191411/240‑1 (MQ=255)
                                                                                                                                                                                                                      gAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGTGGTCGCGGCGACTAGACTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACATATCCATCCTGCTTAACCAGCTCgtcaat                              >  2:332437/1‑246 (MQ=255)
                                                                                                                                                                                                                       ggCCATTGGGGATTTGC‑CCCGCCACTTGGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                             <  1:676933/250‑1 (MQ=255)
                                                                                                                                                                                                                        cccATTGGGGATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                            <  1:124871/250‑1 (MQ=255)
                                                                                                                                                                                                                              gggggttttcccccGCCTCTTTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGTTCGTGCCCCGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTtg                       <  1:88596/239‑1 (MQ=255)
                                                                                                                                                                                                                                 gATTTGC‑TCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGATGTTTAACCAGTTCGATAATACCGTTGTGAc                   >  1:518516/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGGAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGCCGTTGtgtttt         >  1:50780/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cccACTTTGCGGACGATGCGCTCTTTTTCTTCTGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg         <  1:720300/250‑1 (MQ=255)
                                                                                                                                                                                                                                               aCTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTgtttgttt      >  1:216372/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ttGCGGGCTATGCGCTCTTTTTTGCCTGGCTGGGGGGCCTTGCGTCCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCCCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCa    <  2:587872/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   tGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAtt  >  1:764366/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
CGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGT‑GTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGC‑TCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CA‑GGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATT  >  NC_000913/3560116‑3560606

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