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.0% 179.3 / ‑4.4 50intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (24/25);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑C‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT  >  NC_000913/3560206‑3560692
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
tcttctTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcg                                                                                                                                                                                                                                                >  1:101837/1‑250 (MQ=255)
   tcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                                                                                                                                               >  1:266096/1‑247 (MQ=255)
   tcttTTACCACAAACGTGTGACCAAGGTTGAGTTTGTGGGTGAAAAGGCGAAAATTGCGGTGATTGAGAAGTGCGGGCGCGCCCTCTTCGGGGGTGGTGCCGATATCGATAAACAG‑CGGCGTGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCGTGCGATCATGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGGc                                                                                                                                                                                                                                                 <  2:147447/247‑1 (MQ=255)
    cttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgcagg                                                                                                                                                                                                                                            >  1:259002/1‑246 (MQ=255)
      tttACCATCAACGTGTAAGCAACGTTGAGTTGGTTGGGGACAAGGGGCAAATGCCTGTAATGGGGCGGTGCGTGCGCTACCGCTCCGGGCGTGGGGCGAATATGGAAAAACAG‑CGTCGGGCAATGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGa                                                                                                                                                                                                                                          <  2:140567/250‑1 (MQ=255)
        gcccATAAACGTTTTAGCAACGTTGAGTTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGGGCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAG‑CGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATc                                                                                                                                                                                                                                        <  1:415657/248‑1 (MQ=255)
            aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCCTGCTGGCGC‑AGGAACAGGtt                                                                                                                                                                                                                                   >  2:189507/1‑251 (MQ=255)
            aTCAAAGTGTTAGCAACGTTGAGCTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGACATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGtt                                                                                                                                                                                                                                   >  2:423884/1‑251 (MQ=255)
                               ggagaTTGTGGGTGACAATGCGCAAATTGTTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATGGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                <  1:434549/250‑1 (MQ=255)
                                               aatgcgaaattttctgggttggAGCAGTGCGTGCGCTCCCGCTTCCGGGGTGGTGCCGTTATCGATAAACGG‑CGTCGAGCCATTGGGTATTTGCCCCGCCTCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCTCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                                                                                                                                <  1:417508/230‑1 (MQ=255)
                                                aTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                                                                                                                               >  2:339888/1‑251 (MQ=255)
                                                    gaaaaTTGCTGGTATTGAGCAGGGCGGGCGCTCCCGCTCCGGGCGGGGTGCCGATATCGATAAACGG‑CGTGGAGCATTTGGGGTTTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                                                                                                                           <  1:241983/249‑1 (MQ=255)
                                                      aaaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                         >  1:270805/1‑251 (MQ=255)
                                                           gttgtgATGGAGCGGTGCGTGCGCTACCGCGTCGGGCGGGGTGCCGATATGGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                                                                                                                     <  1:248516/248‑1 (MQ=255)
                                                                gagtGAGCGGTGCGGGCGTTCCCGTTTCCGGGGGGGTGCTGATAGCGATAAACGG‑CGTCGAGCCATGGGGGATTTGCTCCTCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGTTGGCTTTGCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                                                                                                                                <  1:423884/247‑1 (MQ=255)
                                                                aattgagaagggcgggccttcccGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGCCTAGCAATGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCGCTTTTTCTTCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                                                                                                               <  1:343709/230‑1 (MQ=255)
                                                                 aTTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                 >  2:338160/1‑248 (MQ=255)
                                                                        aGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                                                                                                                       >  2:217821/1‑251 (MQ=255)
                                                                                  cTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctct                                                                                                                                                             >  1:211339/1‑251 (MQ=255)
                                                                                      cgcttcaagCGTGGTGCCGTTATCGCTAAACAG‑AGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTCCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGAAGCGCCGCACCGCCATGATGACGC‑AGGATCCGGTTTTGCTCCGCCAGCTAATTGCGGTCGCGGCGAATAACATGCGGGCTGACGGCGCAATGCTCTACCCGCTCTTctg                                                                                                                                                         >  2:425459/9‑249 (MQ=255)
                                                                                            gggCGTGGTGCCGATATCGATAACCAG‑CGCCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                                                                                                                   <  1:32343/250‑1 (MQ=255)
                                                                                             ggCGTGGTGCCGCTAACGATAAACAG‑CGTCGAGCCATTGGGGATCTGCTCCGCCCCTTTGAGGGAGATACGCTCTTATACTTCGGTCTGGGAGGCCTTGCGAACGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGACGCACCGCCATGATGGCGC‑AGGATCAGGCTTTGCTCCGACAGCTCATTGAGGTCGCGGCGAATAGTCTACGGGCTGACGGAGAAATGATCTACCAGCTCTTCGGTACTGAc                                                                                                                                                  >  2:287147/1‑251 (MQ=255)
                                                                                              gtgtggtgCGGATATCAATAA‑CAGGCGTCGGTCCTTTGGGGTTTTGCTCCGCCTCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGTGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGATGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                                                                                                                                 <  1:440713/249‑1 (MQ=255)
                                                                                                      ccgaatatcatgaaagagcttCGAGCCATTGGGGTTTTGCGCCGCTACTTTGGGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGAACTGTAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                          <  1:434553/231‑1 (MQ=255)
                                                                                                      ccGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAACGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACTGACATAACCctg                                                                                                                                         >  2:117537/1‑251 (MQ=255)
                                                                                                      ccGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGACTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACTGAGCAATGCTCTACCCGCTCTTCGGAACTGAAATAACCctc                                                                                                                                         >  2:416399/1‑250 (MQ=255)
                                                                                                        gATATCGAAAAACGG‑CGTCGAGCCTTTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCGCCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                                                                                                                       <  1:217821/251‑1 (MQ=255)
                                                                                                          tatCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                                                                                                                     >  2:473194/1‑251 (MQ=255)
                                                                                                              gtataacgggccttgggcaaTTGGGGTGTTGCTCCTCCGGTTTGCGGGCGAGGCGCTCTTTTCCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCTACCGGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTa                                                                                                                                  <  2:296547/232‑1 (MQ=255)
                                                                                                              gATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                                                                                                        >  1:301715/1‑196 (MQ=255)
                                                                                                                    cAG‑CGACGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGTGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                           >  2:151825/1‑251 (MQ=255)
                                                                                                                      accgtagcccgttcggggCTTTGCCCCGCCATTTTGCGGGGGTGGCCCGCTTTTTCTTCGTTCTGTTGGGCCGTTCGACCGTCCCACGGCGGGTTACCCGAGCTGTAAGGCCTCGCCTCTCCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                          <  2:210600/237‑1 (MQ=255)
                                                                                                                               ccATTGGGGTTTTGCTCCGCCATTTTGCGGGCGTTGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                                 <  1:339888/251‑1 (MQ=255)
                                                                                                                                   tGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                                                                                                                   >  1:114240/1‑181 (MQ=255)
                                                                                                                                                                  gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATccc                                                                              <  1:400904/251‑1 (MQ=255)
                                                                                                                                                                    gcTCTTTTTATTCGGTCTGGGTGGCCTTGCGATCGTGACAAGGCGCGTTAACCGAACTGGCAGGCAGCCCCGCACCGCCCTGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTATGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACCTAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTg                                                                            >  2:231668/1‑251 (MQ=255)
                                                                                                                                                                       ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTgg                                                                           <  1:267462/249‑1 (MQ=255)
                                                                                                                                                                         ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatt                                                                       >  2:169288/1‑251 (MQ=255)
                                                                                                                                                                            tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                                     <  1:100479/250‑1 (MQ=255)
                                                                                                                                                                            tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGAGTTCACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattatt                                                                    >  2:96921/1‑251 (MQ=255)
                                                                                                                                                                             cttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGCATTAttt                                                                   >  2:13166/1‑251 (MQ=255)
                                                                                                                                                                               ccacgacgcccttcgttctgATCGGGCCACGGGGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttt                                                                   <  2:202818/230‑1 (MQ=255)
                                                                                                                                                                                   tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                       >  1:419991/1‑225 (MQ=255)
                                                                                                                                                                                      gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttt                                                             <  1:338160/248‑1 (MQ=255)
                                                                                                                                                                                         gggCCTTGGGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCGGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                                       <  2:101837/250‑1 (MQ=255)
                                                                                                                                                                                                  gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttt                                                                   >  1:202818/1‑230 (MQ=255)
                                                                                                                                                                                                            gggggtgtTAACCGACCGGGAAGGCCGACCCGCCCCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCATCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaaa                                    <  1:473626/247‑1 (MQ=255)
                                                                                                                                                                                                                gtgtTAAACGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                 >  2:469054/1‑250 (MQ=255)
                                                                                                                                                                                                                      ccaaaatgggagtgcagcgcgtgcccGCCATGATGGCGC‑AGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAGTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTa                          <  2:92312/227‑1 (MQ=255)
                                                                                                                                                                                                                            cTGGAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAccc                     >  1:144212/1‑250 (MQ=255)
                                                                                                                                                                                                                              ggAAGGCAGCGCCGCACCGCCATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg                  <  1:96921/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cagCCATGATGGC‑GGAGGATCAGGGTGTGGTCCGGCAGCGTAGTGAGGTCGCGGCGAATAGGCTGCGGGATGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGGTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTGTCGCTCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgt  <  1:277200/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                           >  2:233172/1‑95 (MQ=255)
                                                                                                                                                                                                                                                 ccATGATGGCGC‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                           <  1:233172/95‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAG‑CGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑C‑AGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGT  >  NC_000913/3560206‑3560692

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