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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,3620GA96.0% 141.8 / NA 51intergenic (+496/+260)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (23/25);  minor base C (0/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCA  >  NC_000913/3560119‑3560604
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aaCTGGGAGATAAAATCGAGCGTCGCTGCGCCAATGACCCCGCAATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTCCCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGGGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGGCGAGCCATTGGGGATTTGTTCCGCCACTTTGCGGACGATgcg                                                                                                                                                                                                                                               <  2:280166/251‑1 (MQ=255)
           aaaaTCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTttctt                                                                                                                                                                                                                                    >  1:229814/1‑251 (MQ=255)
                cGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCACCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTc                                                                                                                                                                                                                               >  2:245431/1‑251 (MQ=255)
                 gAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCGGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCt                                                                                                                                                                                                                              <  1:484757/251‑1 (MQ=255)
                    cGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGAAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCCAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTggg                                                                                                                                                                                                                           >  1:149382/1‑251 (MQ=255)
                    cGTCGATTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCCACGTGTTAGCACCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATACACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTggg                                                                                                                                                                                                                           >  2:203499/1‑251 (MQ=255)
                            cGCCAAAGATCCCGCCCTCGCTGCGGCGTAATTCGCCACCTGCGAGAATGATGCGAAAATCGTCTTTTCCCAACAACGTGTTAGCAACGTTGAGATTGTTGGTGACAAGGCGCAATTTGCTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGTGGCCGATAGCGATAAACGGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCtt                                                                                                                                                                                                                    <  2:314128/250‑1 (MQ=255)
                             gCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCGCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGc                                                                                                                                                                                                                  <  1:313618/251‑1 (MQ=255)
                             gCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCACCAACGTGTAAGCAACGTTGAGATTGTGGGTGACAAGGCGCAAATTGCTGTGATTGAGCAGTGCGGGCGCTACCGCTTCCGGCGTGGTGCCAATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGc                                                                                                                                                                                                                  <  2:136063/251‑1 (MQ=255)
                                 aTGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                              <  1:421010/251‑1 (MQ=255)
                                   gATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGt                                                                                                                                                                                                            <  1:468837/251‑1 (MQ=255)
                                         gccaccctggCTGCGTAATGCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCACCAACTTGTTAGCAATGTGTAGATTGTTGGTGACAAGGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCTAGATCGATAAACGGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTTCCACg                                                                                                                                                                                                      <  2:443718/243‑1 (MQ=255)
                                           cATCCCGGCTGCGTAATTCGCCACCGGCGAGAATGAGGCAAAAATCTTCTTTTACCACAAACGTGTTAGCAACGTTGAGATGGTTGGTGACAATGCGCAAATTGCTGTTATTGAGCAGGGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGc                                                                                                                                                                                                    <  2:293107/251‑1 (MQ=255)
                                            accgcgGCTGCGTAATTCGCCACCGGCGGGAATGATGCGAAAATCTTTTTTTACCATCAAGGGGGTAGCAACTTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGGTCTCGCTTCCGGCGGGGGGCCGATATTGATAAACGGTGTCGAGTCATTGGGGATTTGCTCCGCCACTTTGCGTACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCg                                                                                                                                                                                                   <  1:457102/249‑1 (MQ=255)
                                                       gTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCACCGTGTTAGCAACGTGGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGATCAGTGCGTGCGCTACCGTTTCCGGCGTGGTGCCGATATCGATAAACAGCGCCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGaa                                                                                                                                                                                        <  2:136375/251‑1 (MQ=255)
                                                       gTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGa                                                                                                                                                                                         >  1:267677/1‑250 (MQ=255)
                                                         aaTTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACt                                                                                                                                                                                      >  1:92261/1‑251 (MQ=255)
                                                             cGCCACCGGGGAGAATGAGGCGAAAATCTTCTTTTACCAGAATCGTGTGAGCAACGTTGAGATTGTTGGTGCCAATGCGCAAATTGCTGTGATTGAGCAGGGCGTGCGATACCGCTTCCGGCGTGGTGCCGATATCGAAAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGaa                                                                                                                                                                                  <  2:137965/251‑1 (MQ=255)
                                                                  ccGGGGAGAAGGATGCGAAAACCTTCTTTTCCCATAAACGTGTTAGAATGTTTGAGATTGTTGGTGACAATGCGCTAATTGCTGTGATTGAGCAGTGCGGGCGCTACCGCTTCCGGCGGGGGGCCGATATCGATAAACAGCGCCGAGCCATTGGGGATTGGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAg                                                                                                                                                                             <  2:229814/251‑1 (MQ=255)
                                                                      cGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGcgcc                                                                                                                                                                         >  1:152417/1‑251 (MQ=255)
                                                                                    aaaTTTTCTTTTACCAACAACGGTTTAGAAACGTTGGGTTTGTGGGTGAAAAGGCGAAAATTGCGGTGATTGAGCAGGGGGGGCGCTCCCGCTCCCGGGGGGGGGCCGATATCGAAAAACAGGGGCGAGCATTGGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCatgatg                                                                                                                                                           <  2:38888/251‑1 (MQ=255)
                                                                                            ttttACAATCAACGAGTTAGCAACGTTGAGATTGTTGGCGACAATGCGCAAATTGCTGTGATTGAGCAGTGAGTGCGCTACCGCTTCCGGCGAGGTGCCGATATCGAAAAACAACGTCGAGCCATTGGGGATGGGCTCGGCCACGTTGCGGACGATGCGCTCTTTTTATTCGGCATGTGTGGCCTTGCGATCGTACCAAGGCGTGTTAACCGAACTGGGAGGCAGCGCCGCACCCACATGATGGcgcaggg                                                                                                                                                   >  2:258242/1‑245 (MQ=255)
                                                                                                  cATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAATTTGCTGTGATTGAGCAGTGAGTGCGCTACCGCTTACGGCGTGGTGCCGATATCGAGAACAAGCGTCGAGCCATTGGGGCTTTGCTCCGCCACTTTGCGGACGGTGCGCTCTGTTTCTTCGGTCTGGGTGACCTTGAGATCGTGCCACGGCGAGTTAACAGGACTGGAAGGCAGACCCGCACCGCCaaaaaggagcaggatccagg                                                                                                                                             >  2:278535/1‑232 (MQ=255)
                                                                                                         gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                                      >  1:438611/1‑251 (MQ=255)
                                                                                                           gtTAGCAACTTTGATATTGTTGGTGACAATGCGCAAATTGCTGTGATTGGGCAGTGCGTGCGCTACCGCTTCGGGCGGGTGGCCGATATCGATAAACGGCTTCGAGCAATTGGTGTTTTGCTCCGCCACTTTGCGGACGATGCGCGCTTTTTCTTCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCTAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                                                    <  1:199071/251‑1 (MQ=255)
                                                                                                                 aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                              >  1:34261/1‑251 (MQ=255)
                                                                                                                   cGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                             >  1:403270/1‑250 (MQ=255)
                                                                                                                         agaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGCTATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTTCTCCGCCCGCTCATTg                                                                                                                         >  2:85535/1‑249 (MQ=255)
                                                                                                                                tgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                >  2:445859/1‑251 (MQ=255)
                                                                                                                                   tGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                              >  1:212766/1‑250 (MQ=255)
                                                                                                                                      cAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                           >  2:85024/1‑250 (MQ=255)
                                                                                                                                        aTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATCTCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                                        >  1:205599/1‑251 (MQ=255)
                                                                                                                                                        gATTGAGCAGTGCGTGAGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGGGGTAGCGGCGACTAGTCTGCGGGCTGACgg                                                                                          >  2:102939/1‑249 (MQ=255)
                                                                                                                                                                   gcgtgcgCTCCCGCTCCCGGCTGGGTCCCGATATCAATAACCAGCGTCGGCCCATTGGGGGTTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                            <  1:409653/251‑1 (MQ=255)
                                                                                                                                                                         gcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGATATGCTCTACAAGctc                                                                      >  2:318568/1‑251 (MQ=255)
                                                                                                                                                                                   ccgggcgggggacagATATGGAAAAACAGGGTCGACCATTGGGGGTTTGGTCCGGCCACTTTGCGGACGATGCGCTCTTTTCCTTCGGTGTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTACCCGAACGGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                               <  2:236679/235‑1 (MQ=255)
                                                                                                                                                                                   ccGGCGGGGTGCATATATCAAAAAACAGCGTCGGGCCATTGGGGTTTTGCCCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                             <  1:479921/250‑1 (MQ=255)
                                                                                                                                                                                   ccGGCGGGGGGCCAATATCGATAAACGGCGCCGGGCAATGGGGGTTTTGCTCCGCCACTTGGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTGGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                            <  1:245431/251‑1 (MQ=255)
                                                                                                                                                                                            tgCCGATAACGCTAAACAGCGTCAAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGAAGCGCTCTTTTTCTTCGGTCTGGGTGGCCTAGCGATCGTGCCACGGCGCGTTAACCGAACTGGAAGGCATCGCGGCACCACCAAGATGGCGCAGGATCAGGATTTGCTCCGCCAGCTCATTGAGGGCACGACGAATCATCTGCGGGCTGACGGAGAAATGCTCTAACAAGTCTTCGGTAcggccatcaccc                                                   >  2:436504/1‑240 (MQ=255)
                                                                                                                                                                                             ggctgattacgaaaagcggCGTCGAGCCATTGGGGTTTTCCCCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                   <  1:445859/233‑1 (MQ=255)
                                                                                                                                                                                               ccGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAACGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGACTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAATGACATAACCctt                                                 >  1:285036/1‑250 (MQ=255)
                                                                                                                                                                                                  atatCGAAAAACAGCTTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                              <  1:27026/251‑1 (MQ=255)
                                                                                                                                                                                                       gataaccggctgcgggCCATGGGGGTTTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                                         <  1:411853/236‑1 (MQ=255)
                                                                                                                                                                                                       gATAAACCGAGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGAGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCGGACGGGGAACATGTATACCAGCTCTTCGGTACTGAAATAACCCTGCTGTTTa                                          >  2:56022/1‑250 (MQ=255)
                                                                                                                                                                                                       gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTTCTGTTTa                                          >  2:122443/1‑250 (MQ=255)
                                                                                                                                                                                                                   cGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGACGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGGTGGCGAAGGACCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGACTATTCTCCGGGCTGCCGTGGAAATGCACTACCAACTCTTCGTTACTGACCTAACCCTGCTGTTTAACCAGTTCGata                              >  2:247380/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gccccGCCACTTTGCGGACGAGGCTCTCTTTTTCTGCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGAACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg             <  1:122443/247‑1 (MQ=255)
                                                                                                                                                                                                                                         ccgtcACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATGGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg           <  1:385351/245‑1 (MQ=255)
                                                                                                                                                                                                                                           ccctcTTTGCGGACGATGCGCTTTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTACGTTGtgtttg       <  1:425714/247‑1 (MQ=255)
                                                                                                                                                                                                                                            ccccTTTGCGGCCGAGGCGCTCTTTCGCTTCGGTCTGGGTGTCCTGGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGGCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt        <  1:203499/246‑1 (MQ=255)
                                                                                                                                                                                                                                                 tggCGGCCGATGCGCTCTTTTTCTTGGGGTGGGGTGGCCTGGCGAGCGTGCCCCGGGGTGTTACCCGAACTGGAAGGCTGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTTGTGTTTGTTTCa  <  2:205599/248‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
AACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCA  >  NC_000913/3560119‑3560604

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