Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,455 +G intergenic (+589/+167) 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,4551.G98.7% 252.7 / NA 78intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (32/45);  minor base T (0/1);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATG‑CGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGC  >  NC_000913/3560209‑3560696
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
tcttTTACCATCAACGTGTTAGCAACGTGGAGATTGTTGGTGACAATC‑CGCAAATTGCTGTAATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCTTCGAGCCATTGGGGATTTC‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGc                                                                                                                                                                                                                                                    <  1:583929/247‑1 (MQ=255)
    taaCCATCAACGTGTTAGCAACGTTGAGATTGTGGGTGACAATG‑CGCAAATTTCTGTGATTGAGCAGGGCGGGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGGGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAt                                                                                                                                                                                                                                             <  1:503544/248‑1 (MQ=255)
         aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATG‑CGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGtt                                                                                                                                                                                                                                       >  2:23337/1‑251 (MQ=255)
            aaCGTGTTAGCAACGTGGAGATGGTGGGTGACAAGG‑CGCAAATTGCTGTGATGGAGCAGGGCGGGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTGG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                                                                                                                                    <  1:425450/251‑1 (MQ=255)
                  ttAGCAACTTTGAGTTTGTGGGTGACAATC‑GGCAAATTGCTGTGATTGAGCGGTGCGTGCGCTACCGCTCCCGGCGTGGGGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                                                                                                                                               <  1:528106/250‑1 (MQ=255)
                    aGCAACGTTGAGATTGTTGGTGACAATG‑CGCAAATTGCTGTGCTTCAGCAGCGCGTGCGCTACCACTTCCGGCGTGGCGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATTCGCTCTTTTTCTTCGGTCCGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGTCGCAGGATCAGGTTTTGCTCCGccc                                                                                                                                                                                                                            >  2:45002/1‑250 (MQ=255)
                         cgtggagaTTGTTGGTGA‑AATGGCGCAAATGGCTGGGATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGGGCCGATACCGATAAACAGCGTCAAGCCTTTGGGGATTTG‑CTCCGCCACT‑TGGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                                                                                                                                       <  1:311453/247‑1 (MQ=255)
                          gTTGAGATTTTGGGTGCAAATG‑CGCAATTGGCTGTGATTGAGCAGGGCGTGCGCTACCGTTTCCGGCGTGGTGCCGATATCGATAAACGGCGTCGAGCCATTGGGGTTTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGGGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAt                                                                                                                                                                                                                      <  1:322122/251‑1 (MQ=255)
                                 ttgttgGTGACAATG‑CGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTc                                                                                                                                                                                                               >  2:265360/1‑251 (MQ=255)
                                   gtgggTGACAAGG‑CGAAAAGGGCTGGGATTGAGCAGGGGGTGCGCTCCCGTTTCCGGCGTGGGGCGGATATGGATAAACAGCGTCGAGCCATGGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                                                                                                                              <  1:436801/247‑1 (MQ=255)
                                              tG‑CGCCAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                                                                                                                                  >  2:566796/1‑251 (MQ=255)
                                              tG‑CGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                                                                                                                                  >  2:512740/1‑251 (MQ=255)
                                                              gATTGCGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                     >  2:219466/1‑249 (MQ=255)
                                                               attgagaagttggtccgttACCGCTTCCGGCTGGGTCCCGATATCAAAAAACGGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                                                                                                                                  <  1:232037/233‑1 (MQ=255)
                                                               aTTGAGCAGTGAGTGCCATACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCCCG‑GCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCAAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                     >  2:617099/1‑248 (MQ=255)
                                                                  gAGCAGGGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAACCAGCGTCGACCCATTGGGGTTTTG‑CTCCGCCACT‑TTGCGGACGAGGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCAGG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                                                                                                                                <  1:580931/250‑1 (MQ=255)
                                                                   aGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                                                                                                                              >  1:492385/1‑251 (MQ=255)
                                                                     cAGTGCGTGCGTTCCGGCTTCGGGCGGGGTGCCGATATCGATAAACAGCGCCGAGCTATGGGGGATTGG‑CCCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTGCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                                                                                                                            <  1:226830/251‑1 (MQ=255)
                                                                       gtgcgtgcGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGCTCACGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCAGGCTGACGGAGAAATGctcc                                                                                                                                                                          >  2:120161/1‑250 (MQ=255)
                                                                          cgtgcgCTAACGCTTCTGGCGAGGAGCCCATATCGATCAACAGACTCGAGCCATTGGGGATTTG‑CTCCGGCACT‑TTGCGGACGATGCGCTCTTTTTCGTCGGTCTGGGTGG‑CCATGCGATCGTGCCACG‑GCGTGTTCACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCCCCTTGAGGACGCGGCGAATACTACGCGGCCTGACGGAGAAATGCTCTaac                                                                                                                                                                       >  2:626834/1‑249 (MQ=255)
                                                                            ggcgcTACCGCTGCGGGCGTGGTGCGGATATCGATAAACAGCGTCGAGCCATTGGGGATGTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGGGCCACG‑GGTTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                                                                                                                                     <  1:566796/250‑1 (MQ=255)
                                                                                  aCCGTTTCCGGCGTGGTGCCGATTGCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGGGG‑TCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCATCTCTTc                                                                                                                                                               <  1:489106/251‑1 (MQ=255)
                                                                                         cgggggtggggCCGATATCGATAAACAGCGCCGAGCCGTTGGGGATTTG‑TTCCGCCACT‑TTGCGGACGATGCTCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGAGCGTGCCACG‑GCTTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                                                                                                                        <  1:381795/246‑1 (MQ=255)
                                                                                         ccgttggggtgCCGATATCGATAACCATCGTCGAGCCATTGGGGATTTG‑CGCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCTCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                                                                                                                         <  1:300429/243‑1 (MQ=255)
                                                                                                ggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTAATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                                                                                                                                 >  2:260094/1‑251 (MQ=255)
                                                                                                 gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                                >  1:375418/1‑251 (MQ=255)
                                                                                                   ccctatatCTATAAACAGGGTCGAGCC‑TGTGGGATTTGGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTCCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                              <  1:300095/247‑1 (MQ=255)
                                                                                                      aatagtgATAAACAGCTTCGCGCCATTGGGGGTTTG‑CTCCGCC‑CTGTGGTGGAGGATGCGCTTGTTTTCTTCGTTC‑GGGTGGGCCTTGCGTTCGGGCCACG‑GCTTTTTACCCGAACCGGAAGTCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                                                                                                                           <  2:518891/245‑1 (MQ=255)
                                                                                                        tatCGCTAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                                                                                                                         >  2:103203/1‑251 (MQ=255)
                                                                                                          tCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                                                                                                                            >  2:483432/1‑214 (MQ=255)
                                                                                                          tCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                                                                                                                            <  1:483432/214‑1 (MQ=255)
                                                                                                                aaCAGCGTCGGGCCATTGGGGATTGG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                                                                                                                 <  1:160698/251‑1 (MQ=255)
                                                                                                                  cagggtctggcttgtGGGGTTTTG‑CTCCGCCTCT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                               <  2:490499/237‑1 (MQ=255)
                                                                                                                     cGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGACGATGCAATCTTTTTCTTCGGTCTGGGAGG‑CCTTGCGACAATGCCACG‑GCGTATTACCCGAAATGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGACCAAGTTTTGCTCCGCCAGCTCATTGAGGTCACGGCGAATAGTCTGCGGGCTGCCGGAGACATGCTCTACCAGCTCTTCGGAACTGACATAACCCTGCTGTTTAACCCGTTCg                                                                                                                             >  1:287606/1‑250 (MQ=255)
                                                                                                                          gtgccatgggtgattgctccccccACT‑TGGCGGACGAGGCGCTCTTTTTCTCCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGGGTGAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata                                                                                                                          <  2:30693/227‑1 (MQ=255)
                                                                                                                           gCCATGGGTGATTTG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCGTCGGTCTGGGTGG‑CCTGGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                      <  1:173770/251‑1 (MQ=255)
                                                                                                                                      tttG‑CTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGTGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt                                                                                                           >  1:364069/1‑251 (MQ=255)
                                                                                                                                            cccGCCACT‑TTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCCCG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTgtgt                                                                                                      <  1:217475/250‑1 (MQ=255)
                                                                                                                                             ccGCCACT‑TTGCGGACGATGCGCTTTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt                                                                                                     <  1:410639/251‑1 (MQ=255)
                                                                                                                                              cGCCACT‑TTGCGGACGTGGCCCTCTTTTTCTTCGGTCGGGGTGG‑CCTTGCGCTCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                                    <  2:492385/251‑1 (MQ=255)
                                                                                                                                                 cctctttggggaggATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGGGACCGTGCCACG‑GCGGGTGAACCGAACTGGGAGGCAGCGCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgt                                                                                                  <  2:517364/238‑1 (MQ=255)
                                                                                                                                                    t‑ttGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGAGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                               >  2:103107/1‑250 (MQ=255)
                                                                                                                                                    attggCGGACGATGCGCTCTTGTCCTCCGGCCGGGGTGG‑CCTTCCGATGGGGCCACG‑GCTTGTTAACCAACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                               <  1:249070/247‑1 (MQ=255)
                                                                                                                                                            acGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTCCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                        >  1:224496/2‑250 (MQ=255)
                                                                                                                                                             cGATGCGCTCTTTTTCTGCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTTTTACCCGAACTGGAAGTCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTTTTCGGTACTGACATAACCCTGCTGTTTATCCAGTTCGATAATACCGTTGTGACGTTGCGTTTGTTTCATTTATaaa                                                                                      <  1:373238/251‑1 (MQ=255)
                                                                                                                                                                 gcgcTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATccc                                                                                  >  2:454092/1‑251 (MQ=255)
                                                                                                                                                                  cgcTCTTTTTCTGCGGGCGGGGGGG‑CCGTGCTTTCGTGCCACG‑GCGTGTTACCCGAACGGGAAGGCGGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                                 <  2:504604/251‑1 (MQ=255)
                                                                                                                                                                    gtctTTTTCTTGGGTTTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTGAACCAAACTGGAAGGCAGCGCCGCGCCGCGATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATATCCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                                 <  2:510349/248‑1 (MQ=255)
                                                                                                                                                                      ctTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTgg                                                                               <  1:98627/249‑1 (MQ=255)
                                                                                                                                                                       tttttCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGGGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGaa                                                                             <  1:255080/250‑1 (MQ=255)
                                                                                                                                                                        ttttCTGCGGCCTGGTTGG‑CTTGGCGACCGTGCC‑CGTGCGTGTTAACCGAACTGGAAGGCAGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCTGCGAATAGACTGCGGGCTGACGGAGAAATGCTCGTCCAGCTCTTCGGTACTGACATATGCTTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGGTGTGGTTGTTTCTGTTAGAAATCCCTGGAat                                                                            <  2:255197/250‑1 (MQ=255)
                                                                                                                                                                          ttcttcGGTCTGGGTGG‑CCTTGCGTTCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATGTATAAATCCCTGTAatta                                                                          <  1:495890/250‑1 (MQ=255)
                                                                                                                                                                          ttctcaggTCTGGGTGT‑CCTTGCGATCGGGCCACG‑TTGTGTTAACCGAATTGGAGGGCAGCCCCGCACCGCCATGATGGCGCATGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTACGTTGTGTTTGTTTCATGTATTAATCCCTGGAatta                                                                          <  2:584120/244‑1 (MQ=255)
                                                                                                                                                                              gcGGTCTGGTTGG‑CCTTGCTAGCGGGCCACG‑GCTTTTTAACCAATCTGGAAGTCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATATCCCTTCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttt                                                                       <  1:71311/248‑1 (MQ=255)
                                                                                                                                                                                ggTCTGGGTGG‑CGTTGCGTCCGGGCCACG‑GCGGGTTACCGGATCGGGAAGGCAGCGCAGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                   <  2:224681/251‑1 (MQ=255)
                                                                                                                                                                                     gggTGG‑CCTTGCGATCGGGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttcg                                                              <  1:31482/251‑1 (MQ=255)
                                                                                                                                                                                      ggTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                                                                                                                         <  1:231822/159‑1 (MQ=255)
                                                                                                                                                                                      ggTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                                                                                                                         >  2:231822/1‑159 (MQ=255)
                                                                                                                                                                                       gTGG‑CCTTGCGAGCGGGCCACG‑GGGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgcg                                                            <  1:238625/251‑1 (MQ=255)
                                                                                                                                                                                          ggacTTGCGATCGGGCCACG‑GCGTGTTAACAGAACTGGACGGCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAACTCCCTGGAATTATTTTCGTTTTCgcgc                                                           >  2:548978/4‑250 (MQ=255)
                                                                                                                                                                                            cccttgccgtcctgccccgtcGTGTTAACCTATCTGGAAGGCCGCTCCACACTGCCATGATGCCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGTGCAGACGGAGAAATGCTCTTCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAGCCCTGCAATTATTTTCGTTTTCTCGCAt                                                         <  2:109510/231‑1 (MQ=255)
                                                                                                                                                                                                gCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCACTGGAATTATTTTCGTTTTCGCGCATTGAgg                                                    >  2:564606/1‑250 (MQ=255)
                                                                                                                                                                                                  gATCGTGCCACG‑GCGTGTTAACCGGACTGGCAGGCAGAGCCGCACCGCCATGATGGCGCCGGATAAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTTTGAAGTTTTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCg                                                   >  2:37616/1‑250 (MQ=255)
                                                                                                                                                                                                           ccG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                                         <  2:470765/250‑1 (MQ=255)
                                                                                                                                                                                                             g‑gCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAg                                       <  2:224496/251‑1 (MQ=255)
                                                                                                                                                                                                                    ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCGTTGAGCGAATCAAAAAAAGCCATCg                                 >  1:357680/1‑251 (MQ=255)
                                                                                                                                                                                                                            aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAATACCGTTGTGAAGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAccc                         >  1:272015/1‑251 (MQ=255)
                                                                                                                                                                                                                              tGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCAc                       >  2:422139/1‑251 (MQ=255)
                                                                                                                                                                                                                              tGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgcg                                                            >  1:79853/1‑214 (MQ=255)
                                                                                                                                                                                                                               gtaaggtggcgctgcgccGCCATGATGGCGCAGGATCAGGGTTTGCTCCGTCAGCTCAGTGATGTCGCGGCTAATAGTCTGCGGGCTGACGGATAAATGCTCTCCCAGCTCTTCGGTTCTGTCATAACCCTGCTGATTAACCATTGCGATAATACCGTTGTGAAGTGGTGGTGGTTTCATTAATAATTCCCCGGTATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAATCCATCGCTAAAGCCACg                      <  2:617219/235‑1 (MQ=255)
                                                                                                                                                                                                                               ggAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg                      >  2:179910/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaGGCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGc                    <  1:6660/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gcagcgcagcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTaa                 <  2:122428/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg            >  1:521393/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gccgcACCGCCATCATGGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGACTTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg            >  2:9738/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ccccGCCATGATGGCTCAGGATCTGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTTTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAt        <  2:396127/249‑1 (MQ=255)
                                                                                                                                                                                                                                             ccccGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTTCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAGTTATAAATCCCTGGAATGATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTATCCCGGCGAt        <  2:609912/249‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTACCCCGGCGAtg       >  1:105339/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   cATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGc  >  1:276012/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
TCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATG‑CGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACT‑TTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGG‑CCTTGCGATCGTGCCACG‑GCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGC  >  NC_000913/3560209‑3560696

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