Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,443 C→G intergenic (+577/+179) 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,4430CG91.1% 153.7 / 11.2 57intergenic (+577/+179)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base G (23/28);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01

AATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCG  >  NC_000913/3560193‑3560685
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
aaTGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCAc                                                                                                                                                                                                                                                       >  2:544127/1‑251 (MQ=255)
   gatgCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCATCGCCGCAgcgc                                                                                                                                                                                                                                                    >  1:73728/1‑251 (MQ=255)
        gAAAATCTTCTTTTACCATCAACGTGTGAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGCTAAACAGCGTCGAGCCCTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTCACCGAACTGGTAGGCGGCGCCGCGGCGCCatgg                                                                                                                                                                                                                                               >  2:102565/1‑250 (MQ=255)
                   tttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGcgcag                                                                                                                                                                                                                                       <  1:236161/248‑5 (MQ=255)
                     tacaatcagcgggtTAGCAACGTGTGGAGTGTTGGTGACAAGGCGAAAATTGCTGTGAGTGAGCTGGGGCTGCGCTACCGCTTCCGGCGTGGGGCCGATAACGATAAACGGCGTCGGGCTATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                   <  2:125171/238‑1 (MQ=255)
                      aCCATCAACGTGTTAGCAACTTTGAGTTTTTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCa                                                                                                                                                                                                                                  <  1:419508/250‑1 (MQ=255)
                       ccaacaccGTGTTAGCAACGTTGAGATTGTTGGTGACAATCCGCAATTTGCTGTGATTGAGCAGTGCGTCCGCTACCGCTTCCGGCGCGCTGCCGATATCGATACACAGCGTCGAGCCATTGGGGATTGGCTCCGCCACGTTGCTGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGAGCCACGGCGTGTTAACCGAACCGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAgg                                                                                                                                                                                                                                <  1:469455/244‑1 (MQ=255)
                         aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTACGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGtt                                                                                                                                                                                                                              >  2:355374/1‑251 (MQ=255)
                          tCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCAGCGCCATGATGGCGCAGGATCAGGttt                                                                                                                                                                                                                             >  2:156461/1‑251 (MQ=255)
                              cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCt                                                                                                                                                                                                                         <  1:243105/251‑1 (MQ=255)
                               gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                        >  1:93170/1‑251 (MQ=255)
                                    aGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGTCGCAGGATCAGGTTTTGCTCCGccc                                                                                                                                                                                                                   >  1:179540/1‑250 (MQ=255)
                                              gagaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCCTGCTGGCGCAGGATCCGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                           >  2:347638/1‑250 (MQ=255)
                                                 aTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAATGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTgcggg                                                                                                                                                                                                       >  2:233906/1‑247 (MQ=255)
                                                         tGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGAGGGCGATGCGCTCTTTTTCTTCGGTATGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCAGCGCCATGATGGCGCAGGAACAGGTTTTGATCCGCCAGCCCATTGAGGTagcgccg                                                                                                                                                                                               >  2:373265/1‑247 (MQ=255)
                                                         tGACAATGCGCAAATTGCTGTGATTGAGCAGGGCGTGCGATACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCCTTGGGGATTTGCTCCGCCACTTCGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTACCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCCGCGCCATGATGGCGCAGGATCATGTTTTGCTCCGCCAGCTCATGAGGGTCgcggcg                                                                                                                                                                                               >  2:263970/1‑250 (MQ=255)
                                                                             ggATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCAATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                                                                                                                          <  1:200101/250‑1 (MQ=255)
                                                                              gATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                           >  1:205292/1‑249 (MQ=255)
                                                                                                         ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGAGCCGCAGCGCCATGATGGCGCAGGCTCAGGTTTTGCTCCGCCCGCTCATTGGGGTCGGGGCGAATAGTCTGAGGGGTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                                                                                                              >  2:320456/1‑251 (MQ=255)
                                                                                                             ggtggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                          <  1:324521/250‑1 (MQ=255)
                                                                                                                 gtgCCGATCTCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                      <  1:431887/251‑1 (MQ=255)
                                                                                                                     cGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGTTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgc                                                                                                                                  >  1:363557/1‑251 (MQ=255)
                                                                                                                           cGATAAACAGCGTGGAGCCATTGGGGATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTa                                                                                                                            <  1:210637/251‑1 (MQ=255)
                                                                                                                               aaaCAGCGTCGCGCCATTGGGGATTTGCTCCGCCACTTTGAGGGCGATGCGCTCTTTTTCTTCGGTCTGGGCGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGGAGGCGGCGCCGCAGCGCCATGCTGGCGCAGGATCAGGTTTTACTCCGCCCGCTCATTGAGGTCGCGGCGACTACTCTGCGGGCTGGCGGTGACATGCTCTAACAGCTATTCGGATCTGCTCTAACCCTGCTGTTTAAcc                                                                                                                         >  2:529087/1‑250 (MQ=255)
                                                                                                                                  ccgtggtcgtcatgtGGGGTTTTTCTCCGCCACCTGGCGGGTGATGAGCTCTTTTTTTGCTGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGATGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGCGGTACTGACATAACCCTTCTGTTTAACCAGtg                                                                                                                     <  2:17354/237‑2 (MQ=255)
                                                                                                                                    ggctgggttcagtGGGGAGTTTGCGCGGCCACTTTGGGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                    <  2:420356/239‑1 (MQ=255)
                                                                                                                                         gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                              <  1:56180/251‑1 (MQ=255)
                                                                                                                                            cctttGGGGTTTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTGCGTGCTGGGTGGCCTTGCGTTCGTGCCACGGCGTGTTACCCGAACTGGAAGGCTGCGCCGCACCGCCATGATGGCGCAGGATCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                          <  1:233512/248‑1 (MQ=255)
                                                                                                                                            ccATTGGGGATTTGCTCCTCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaaatccccaatgg                                                                                                                                                           <  2:61781/203‑12 (MQ=255)
                                                                                                                                            ccATTGGGGATTTGCTCCTCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaaatccccaatgg                                                                                                                                                           >  1:61781/1‑192 (MQ=255)
                                                                                                                                            ccATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                          <  1:363787/251‑1 (MQ=255)
                                                                                                                                                         tGCTACGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt                                                                                               <  1:355374/250‑1 (MQ=255)
                                                                                                                                                          gccccgccgCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg                                                                                              <  1:389998/247‑1 (MQ=255)
                                                                                                                                                           cagcggcattgtgtgggCGGTGCTCTTTTTGTCTTCGGTCTGGTGGGCTCTGCGCTCGGCCCAGGGCGTGTTAACTGACTTTGAAGGCAGCCCCGCACCGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                            <  2:395182/237‑1 (MQ=255)
                                                                                                                                                              cgccccgttGCGGGCGATGCGCTCTTTTCCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                         <  2:214480/244‑1 (MQ=255)
                                                                                                                                                               gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGTTAATAACGTTGTGACGTTGtgtttg                                                                                        >  2:130685/1‑251 (MQ=255)
                                                                                                                                                                        cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAAGTCGCGGCGAATAGTCTGAGGGCTGACGCAGCAATGCTCTACCAGCTCTTCGGAACTGACATAAACCTGCTGTTTAACCAGTTAGATAATACACGTGTGACGTtggtatgatttaatttt                                                                               >  2:425073/1‑236 (MQ=255)
                                                                                                                                                                        cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAACCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                >  1:117499/1‑250 (MQ=255)
                                                                                                                                                                           gCGATGCGCTCGTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaa                                                                            <  1:181373/251‑1 (MQ=255)
                                                                                                                                                                            cgttGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTTCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                           <  2:105605/248‑1 (MQ=255)
                                                                                                                                                                            cggggcggttctgttcttcGGGCTGGGTGGCCTTGCGAACGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCATGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTTTTTATCCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                           <  2:53739/238‑1 (MQ=255)
                                                                                                                                                                            cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaaa                                                                           <  1:481487/251‑1 (MQ=255)
                                                                                                                                                                             gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCTGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTTCCAGCTCTTCGGTACTGCCATAACCCTGCTGTTTAACCAGGTCGATAATTCCGTTGTGACGTTGTGTTTGTTTCATTTATAAAt                                                                          >  2:43189/1‑251 (MQ=255)
                                                                                                                                                                                gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                        <  1:544127/250‑1 (MQ=255)
                                                                                                                                                                                             tcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttt                                                            <  1:320456/249‑1 (MQ=255)
                                                                                                                                                                                             tcGGGCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCATCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAtttt                                                           <  2:73728/250‑1 (MQ=255)
                                                                                                                                                                                                  cTGGGTGGCCTTTCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGATTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACAGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGATCGATAATCCCGTTGTGACGTTGTGTTTGTTTCTTTTATAAATCCCTGGAATTAttttcgtttt                                                     >  1:80560/1‑251 (MQ=255)
                                                                                                                                                                                                   tGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGTATTAttttcgttttc                                                    >  1:527396/1‑251 (MQ=255)
                                                                                                                                                                                                    gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCCCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                     <  1:74919/249‑1 (MQ=255)
                                                                                                                                                                                                        ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCAACGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACAGTTGTGACGTTGTGTTTGTTTCATTTTTAAATCCCTGGAATTATTTTCGTTATCgcgc                                                >  2:230480/1‑250 (MQ=255)
                                                                                                                                                                                                           cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTg                                            >  1:522930/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCGTGTTAACCGAACTGTAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTCACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGc                           >  2:508011/1‑251 (MQ=255)
                                                                                                                                                                                                                               gtgtTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGATGTTTAACCAGTTCGCTAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAAAAAAAGcc                          >  1:248644/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gtTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATc                       >  1:306297/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gtTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTGCGCGCATTGAGCGAATCAACAAAAGCCAt                        <  1:156461/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aaCTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAcc               <  1:83857/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ggAAGGCCGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCAc            <  1:43189/250‑1 (MQ=255)
                                                                                                                                                                                                                                              gAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACgg          <  1:276439/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      cgccgcAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg  <  1:288187/251‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
AATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCG  >  NC_000913/3560193‑3560685

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