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

TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCT  >  NC_000913/3560206‑3560697
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
tcttctTTTACCATCATCGTGTTAGCAACGTTGAGATTGTTGGTGACACTGCGCAAATTGCGGTGATTAAGCAGTGCGTGCGATACCGCTTCCGGCGTGATGCCGATATCGATAAACAGCGTCGAGCCATTGGGTTTTTGCACCGCCACTTTGCGGGCGATGCGCTCTTTTTCGTCGGGCTGGGTGGCATTGCGGTCGTGCCACGGCGTGTTAAACGAACTGGA‑GGGCAGCGCCGGACCGCCATGA‑‑TGGcg                                                                                                                                                                                                                                                    >  2:53043/1‑250 (MQ=255)
 cttcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGcgc                                                                                                                                                                                                                                                   >  2:128971/1‑249 (MQ=255)
 cttcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGcgc                                                                                                                                                                                                                                                   >  1:258261/1‑249 (MQ=255)
      tttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGAt                                                                                                                                                                                                                                              >  1:477292/1‑251 (MQ=255)
          ccATCAACGTGTTAGCAAGTTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCTTGGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAgg                                                                                                                                                                                                                                          <  1:416169/251‑1 (MQ=255)
            atccacgagtTAGCAACGTTGAGATTGTTGGTGACAATGCACAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATAGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCACCACCGCCCTGA‑‑TGGCGCAGGATCAgggt                                                                                                                                                                                                                                        >  2:540939/9‑249 (MQ=255)
               aaCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTg                                                                                                                                                                                                                                     >  2:48774/1‑251 (MQ=255)
                               tGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGC‑‑TGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTg                                                                                                                                                                                                                     >  2:530741/1‑251 (MQ=255)
                                                       aaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGTTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCTGACAGCTCCTTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                                                                                                             >  2:127505/1‑251 (MQ=255)
                                                              ggtttgaagaagggggtgcctgcccGCTCCCGGCGTGTTCCCGATATCGTTAACCTCCTCCGACCCATTGGGGATTTTCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTGTGGGTGGCTTTGCGAGCGTGCCACGGCTTGTTAACGGAACTGGA‑AGGAAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTTCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                      <  1:327228/228‑1 (MQ=255)
                                                                 aTTGAGCAGTGAGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGAGTCGAGTCATTGGGGATTTGCTCCGACACTTTGCGGGCGCTGCGCTCTTTTTCTACGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACAGAACTGGA‑AGGCAGCGACGAACCGCCATGACTTCGC‑GAG‑AACAGGTTTTGCTGCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGAGCTGACgg                                                                                                                                                                                      >  1:407337/1‑248 (MQ=255)
                                                                    gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGCTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGACACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCCGGTTTTGCTCCGCCAGCTCAATGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                                                                                                                >  2:454222/1‑251 (MQ=255)
                                                                          tgcgtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACGTTGCGGGCGATGCGCTCTTTTTCTTCGGTATGGGTGGCCTTTCGATCGTGCCACGGCGTGTTAACCGCACTGGA‑CGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCCGGTTTTGCTACGCCAGCTCAAGGCGGCCCGGGCCAATAGCCTGCGGGCGGACGGAGaaaggcctt                                                                                                                                                                           >  2:380666/1‑244 (MQ=255)
                                                                             gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGAGTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCAAATAGTCTGCGGGATGACGGAGAACTGCTCTAcc                                                                                                                                                                        >  1:24113/1‑250 (MQ=255)
                                                                               gcgctagcgttcccGGCGTGGTGCCGATATCGATAAACGGCGTCGAGCCATTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                                                                                                                                     <  1:56567/239‑1 (MQ=255)
                                                                                    accgtttccggttgggtgCTGATACCGAAAAACGGCGTCGCGCCTTGGGGGTTTTGCCCCGCCACTTGGCGGGCGAGGCGCTCTTTTTCTTCGGTCTGGTGGGCCTTGCGAGCGTGCCCCGGCGTGTTAACCGAACGGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGCTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                 <  1:319785/236‑1 (MQ=255)
                                                                                    aCCGCTTCCGGCGTGGCGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                                >  1:132579/1‑251 (MQ=255)
                                                                                             gggggggggCCGATATCAATAAACGGCGTGGGGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                                                                                                                       <  1:409212/243‑1 (MQ=255)
                                                                                              gggtggggcagatacctaTAAACAGCGTCGGGCAATTGGGTATTTGCTCCGCCTCTGTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                                                                                                                                      <  1:438205/234‑1 (MQ=255)
                                                                                                  gttgcggatttcgaaaaacgccttggAGCCTTTGGGTATTTGCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGTTCGGGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                                                                                                                                  <  1:192065/226‑1 (MQ=255)
                                                                                                   gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGTGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                                 >  1:166760/1‑251 (MQ=255)
                                                                                                     gcctatatGGATAAACAGGGTCGAGCCATGGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                               <  1:112037/247‑1 (MQ=255)
                                                                                                                       cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCGGCGGGCCGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGTCATAACCCTGCTGTTTAACCAGTTCg                                                                                                                              >  2:181491/1‑250 (MQ=255)
                                                                                                                        gccgtcccATGGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                                                                                                                            <  1:47163/245‑1 (MQ=255)
                                                                                                                         gggacccttgggggTTTTGCTCCGCCACTTTGCGGGCTATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGGGTGTTAACCGAACGGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata                                                                                                                           <  1:194122/241‑1 (MQ=255)
                                                                                                                          cGAGCCATTGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGA‑AGGCAGCGCCGAACCGCCATGG‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATACTCTGCGGGCTGCCGGCGAAATGCTCTAACAGCTCTTCGGTACTNACATAACCCTGCTGTTTAACCAGTTCGaaa                                                                                                                           >  2:289074/1‑248 (MQ=255)
                                                                                                                                 ggggggATTTCCCCCGCCACTTTGCGGGCGAGGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCGGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                                     <  1:289074/247‑1 (MQ=255)
                                                                                                                                         ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt                                                                                                           >  1:429149/1‑251 (MQ=255)
                                                                                                                                          tGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg                                                                                                          >  1:95038/1‑251 (MQ=255)
                                                                                                                                           gCTCCGCCACTTTGCAGGCGATGCGCTCTCTTTCTTCGGTCTGGGTGGCCTTGCGATCGAGCCACGGCGTGTTAACCGAACTGAA‑AGGCAGCGCCGCACCGCCATGA‑‑CGGCGCAGGATCAGGGTATGCTCCGCCAACTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCACCTCTTCGGTACTGCCATAACCCAGCTGTTTAAACACTTCGATAATACAGCTGTGACggtgg                                                                                                         >  2:486242/1‑247 (MQ=255)
                                                                                                                                            ccccgcctcttggggggCGATGCGCTCTTTTTCTGCGGTCGGGGTGGCTTTGCGATCGTGCCACGGCGTTTTAACCGACCTGGA‑GGGCAGCGCCGCTCCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                        <  1:261769/237‑1 (MQ=255)
                                                                                                                                            ccccGCCTCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                        <  1:531869/249‑1 (MQ=255)
                                                                                                                                              gcGCCACTTTGCGGGCTATGCGCTTTTTCTCTTGGGTCGGTGTGGCTTTGCGAGCGTGCCAGGGCGTTGTAACCTAACTGGA‑AGCCAGCGCCGCCCCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAAAGCCGTTGTTACGTTgtgt                                                                                                       <  2:485010/249‑1 (MQ=255)
                                                                                                                                               ggccccTTTGCGGCGGATTCGTTCTTTTTGTNCGGTCGGGGGGGCCTGGCGAGCGTCCCTCGGCGTGTTAACCGACCTGGA‑AGGCGGCCCCGCACCGCCATGA‑‑TGGCGCAGTATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGATTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                                     <  2:277432/246‑1 (MQ=255)
                                                                                                                                                gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                                    >  1:477589/1‑251 (MQ=255)
                                                                                                                                                   aCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGC‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTgtttgttt                                                                                                 >  2:131433/1‑251 (MQ=255)
                                                                                                                                                      tggCGGGCGATGCGCTTTTTTTCTGGGGGCNGGGTGGCCTTGCTATCGTGCCACGGCGTGTTAACCGACTTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCa                                                                                               <  2:291061/248‑1 (MQ=255)
                                                                                                                                                        gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                            <  1:529990/251‑1 (MQ=255)
                                                                                                                                                          ggtcgatgaccattatttCTCCGGTCGTGGTGGCCTTGCGACCGTACCAGGGCGATTTACCGGAACTGAA‑AGGCAGCACCCGACCGCTATGA‑‑TGGCGCAGGATCATGTTTTGCTCCTCCAGCTGATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATATCCCTGCTGTTTAACCAGTTCTATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                          <  1:334310/236‑1 (MQ=255)
                                                                                                                                                            gggATGCGCTCTTTTTCTGGGGGCTGTTGGGCCGTGCGGTCGGGCCACGGCGTTTTAACCGAACTGGA‑AGGCCGCGCCGCTCCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGTCATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaa                                                                                        <  2:299361/249‑1 (MQ=255)
                                                                                                                                                                      tttttcCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTgg                                                                                <  1:348174/248‑1 (MQ=255)
                                                                                                                                                                        cctacccgccgctctccctacctGCGATCTTGCCACGGCGTGTTAACCGAACGGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatt                                                                            <  2:9524/229‑1 (MQ=255)
                                                                                                                                                                          ttcttgtttttGGGTGGCCTTGCGACCGGGCCACGGGGTGTTAACCGAACGGGA‑GGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatta                                                                           <  1:524216/240‑1 (MQ=255)
                                                                                                                                                                          ttcttgggcttggggggacTTGCGATCGTGCCACGGCGTGTTAACCGAAC‑GGAGAGGCAGCGCGGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatta                                                                           <  1:302061/232‑1 (MQ=255)
                                                                                                                                                                                ggTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATTA‑‑TGGCGCAGGAACAGGTTTTGCTCCGCCCGCTCATTGAGGACGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCCGCTCTTCGGTCCTGACCTAACCCTGCTGTTTAACAAGTTCGATAATCCCGTTGTGACGTAGCGTTTGGTTCATTTATAAATCCCTGGAATTAttttcgg                                                                    >  2:170466/1‑250 (MQ=255)
                                                                                                                                                                                   ggggtgggCCTGGCGGTCGTCCCACGGGGGGTTAACCGAACTGGA‑AGGCGGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                                 <  2:468850/245‑1 (MQ=255)
                                                                                                                                                                                   cTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCAAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                                 >  1:399285/1‑251 (MQ=255)
                                                                                                                                                                                    gggttgcgcttgggtgcggcccccggcggGTTAACCGACTTGGT‑AGGCGGCGCCGACCCGCCATGA‑‑TGGCGCAGGAACAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACTGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttc                                                                <  2:418119/228‑1 (MQ=255)
                                                                                                                                                                                      ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑GGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgc                                                              <  1:1205/251‑1 (MQ=255)
                                                                                                                                                                                               gCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatt                                                                            >  1:9524/1‑228 (MQ=255)
                                                                                                                                                                                                  aTCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGa                                                   >  1:52438/1‑250 (MQ=255)
                                                                                                                                                                                                    cGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAt                                                 <  1:134395/250‑1 (MQ=255)
                                                                                                                                                                                                    cGTGCCACGGCGTGTTAACAGAACTGGA‑AGGCAGCGCCGAACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGATCCGCCAGCACATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACACTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTGTATTTCATTTATAAATCCCTGGAATTATTATCGTGTTCACGCATTGAGCGAATc                                                >  1:420470/1‑251 (MQ=255)
                                                                                                                                                                                                     gTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                                >  1:417292/1‑250 (MQ=255)
                                                                                                                                                                                                       gCCACGGCGTGTTAACCGAACGGGA‑AGGCAGCGCCGCCCCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCTTTGTGTCGTTGTGGTTGTTGCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaac                                             <  2:220597/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaCCGAACTGGA‑AGGCGGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCt                                <  1:521084/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGAACTGAA‑AGGCAGCGTGGCCCCGCCATGA‑‑TGGCGCAGGATCTGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTaaa                             <  1:564353/251‑1 (MQ=255)
                                                                                                                                                                                                                        gAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAc                            >  1:412036/1‑251 (MQ=255)
                                                                                                                                                                                                                           cgtgggagccaggtcccccccGCCATGA‑‑TGGCGCAGGTTCAGGTTGTGCGTCTCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTTCTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTTGGGTTTGTTTCATTTATAAATCCCTGGAATTATTTGCGTTTTCGCTCACTGAGCGAATCAACAAAAGCCATCGCTAAAccc                          <  2:457237/232‑1 (MQ=255)
                                                                                                                                                                                                                               a‑aGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGTCGTTGTGTTTGTTTCAtt                                                                                             <  2:490856/179‑1 (MQ=255)
                                                                                                                                                                                                                               a‑aGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAtt                                                                                             >  1:490856/1‑179 (MQ=255)
                                                                                                                                                                                                                                 aGGCAGCGCCGCACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCt                    <  1:185778/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCAGCGCCGAACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa                   <  1:447384/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cgcggcACCGCCATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg              <  1:131433/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccATGA‑‑TGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgtg    <  1:128971/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   atga‑‑tgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCt  >  1:397351/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGA‑AGGCAGCGCCGCACCGCCATGA‑‑TGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCT  >  NC_000913/3560206‑3560697

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