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,3620GA97.9% 126.4 / NA 49intergenic (+496/+260)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (18/29);  minor base C (0/1);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGGCGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTG  >  NC_000913/3560118‑3560599
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
gAACTGGGAGATAAAATCGAGCGTCGCTTTGCCAAGGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TATCAACGTTGAGAGTGT‑TGGTGACAATGCGCAAAG‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGAGAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATgc                                                                                                                                                                                                                                            <  1:114889/251‑1 (MQ=255)
      ggagataaattctggcgcctcttggCCAATGATCCCGCCACCGCGTCGTCGTAATCCGCTCCCGGCGAGAATTATGCGAAAAGCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTTT‑TGTTGACAATGCGCAAAT‑GGCTGGGATGGAGCAGGGCGTGCGCTACCTATGCCGGGGTGGTGCCGATATTGATAAACAAAGT‑CGAGCCATTGGGTATTTGCTCCGCCCCT‑TTGCGGACTTTGCGctctc                                                                                                                                                                                                                                      <  2:36821/227‑2 (MQ=255)
      ggAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCtt                                                                                                                                                                                                                                      <  1:50416/251‑1 (MQ=255)
                          cTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGCGCCATTGGGGATTTGCTCCGCCACT‑TTGAGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGcc                                                                                                                                                                                                                 >  2:103689/1‑251 (MQ=255)
                           ttCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCt                                                                                                                                                                                                                >  2:58096/1‑251 (MQ=255)
                             cGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCTCCGGCGAGAATTATGCGAAAATCTTCTTTTACCATCAAGGTGT‑TGGCAACGTTGAGATTGT‑TGGTGACAAGCCCCAAAT‑TGCTGTGATTGAGCAGTGCTGGCGCTACCGTCTCCGGCGTGGTCGCGATATCGATATACTTCGG‑CGGCCCATTGGGGATTTGCTCCGTGACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCtt                                                                                                                                                                                                               <  1:252266/250‑1 (MQ=255)
                                  aTGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTTACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATGGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                         <  2:73220/251‑1 (MQ=255)
                                  aTGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                         <  1:47077/251‑1 (MQ=255)
                                      tCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGCTGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGc                                                                                                                                                                                                     >  2:92123/1‑251 (MQ=255)
                                           ccATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGAGGCGAAAATCTTCTTGTACCATCAACGTGT‑TAGCAACGTGGAGATTGG‑TGGTAACAAGGCGCAAAT‑TGCTGTGATTGAGCAGGGCGGGCGCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGT‑CGAGTCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCTCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACgg                                                                                                                                                                                                <  1:108388/251‑1 (MQ=255)
                                                  ggCTGCGTAATTCGCCCCCGGCGAGAATGATGCGAAAATCTTCTTTTCCCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTa                                                                                                                                                                                         <  2:84670/251‑1 (MQ=255)
                                                  ggCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGtt                                                                                                                                                                                          >  1:201900/1‑250 (MQ=255)
                                                      gCGTAATTCGCCACCGGCGAGAAGGAGGCGAAAATCTTCTTTTACCAGCAGCGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCg                                                                                                                                                                                     <  2:16811/251‑1 (MQ=255)
                                                          aaTTCGCCACCGGTGAGAATCAGGCGAAAATCTTCTTTTACCATCAACG‑GTGTAGCAACGTTGAGATTGT‑TGGGGACAATGCGCAAAT‑TGCTGTGATGTAGCAGTGCGGGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCTACT‑TTGCGGACGATGCGCGCTTTGTC‑GTCGGTCTGGGTGGCCTTGCGATCGGGCCACGTCGTGGTAACCGAACt                                                                                                                                                                                 <  2:27762/251‑1 (MQ=255)
                                                                    cGGCGAAAAGGATGCGAAAATCTTCTTTGACCATCAACGTGT‑TAGCACCGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAgc                                                                                                                                                                       <  2:206885/251‑1 (MQ=255)
                                                                        gagaATGATGCGGAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGCTTTGCTCCGCCATT‑TTGCGGACGAGGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGCCGTGTTAACCGAACTGGAAGGCAGcgccg                                                                                                                                                                   >  2:74868/1‑251 (MQ=255)
                                                                            atgatgCGAAAATCTTCTTTTCCCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAcc                                                                                                                                                               <  1:193235/251‑1 (MQ=255)
                                                                            atgatgCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATGGT‑TGGGGACAATGCGCAAAT‑TGCTGGGATTGAGCAGTGCGTGCGCGACCGCGTCCGGGGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAcc                                                                                                                                                               <  1:130401/251‑1 (MQ=255)
                                                                                tgCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTTGGTGGCCTTGCGGTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGcc                                                                                                                                                            >  2:63219/1‑250 (MQ=255)
                                                                                      aaTCTCCTTTTCACATCAACGCGT‑TAGGAACGTTGAGATTGT‑TGGTGACAATGCGCACAT‑TGCTGTGATTGCGCAGTGCGTGCGCTACCGCTTACGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGCTTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGGCTGGGTGCCCTCGCGATCGTGCCACGGCGTGTTCACCGAACTGGAAGGCAGCGCCGCACCGCCATGATgg                                                                                                                                                     >  2:269470/1‑251 (MQ=255)
                                                                                              tttACCACAAACGTGT‑TAGCAACTTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TCCGGGGATTGAGCAGGGCGTGCGCTACCGCTGCCGGCGTGGGGCCGATATGAAAAAACAGCGT‑CGGGCCATGGGGGATTGGCTCCGCCACT‑TTGCGGACGAGGCGCTCTTTTTC‑TGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCTCCGCCATGATGG‑cgcag                                                                                                                                               <  1:225629/248‑5 (MQ=255)
                                                                                                      cAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtttt                                                                                                                                    >  1:203692/1‑251 (MQ=255)
                                                                                                             gtatgCAACGTTGGGATTGTGTGG‑GATAATGCCCAAAT‑TCCTGTGATTGAGCAGGTCTTGCCCTCCCGCTTCCGTCGTGGTGCCGATATCGATACACACCGT‑GTAGCCATTGGGGGTTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCTCGGCGTGTTATCCGATGTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCg                                                                                                                              <  1:22563/247‑1 (MQ=255)
                                                                                                                 gCAACGTTGAGATTGT‑GGGTGAAAATGCG‑AAATGTGCTGGGTTTGAGCATTGCGGGCGCTACCGCTTCCGGGTTGGTGCCGATATCGGTAAACAGGGT‑GGGGCCAGGGGGGATTTGCTCCGCCTCG‑GTGCGGACGGTGCTCTCTTTTGT‑TTCTGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAGGGCATCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAg                                                                                                                          <  1:243136/251‑1 (MQ=255)
                                                                                                                         gagaTTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGCGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATTATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                   >  1:138858/1‑250 (MQ=255)
                                                                                                                            aTTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGCGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAGCAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGGAGGCCTGGCCGCACCGCCATGGTGGAG‑CAGGATCAGGTTTTGCTCCGCCAGCACATTGAggg                                                                                                               >  2:62077/1‑250 (MQ=255)
                                                                                                                                   ggTGACAATGCGCAAAT‑TGCTGTGATGGAGCTTTGCGTGCGCTCCCGCTTCCGGGGTGTTTCCGATTCCTATAAACAGCGT‑CGGGCCATGGGGTATGTGCCCCCCCACT‑TTGCGGACTATGCGCTCTTTTTC‑TTCGGTCGGGGTGGCCTTGCGATCTTGCCACGGCGTTTTAACCGAACTGGAGGGCAGCCCCGCACCGGCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                         <  1:237224/251‑1 (MQ=255)
                                                                                                                                    gTGACAATGCGAAAAT‑TGATGTGATTGAGCAGTGCGTTCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGAGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                        >  2:233972/1‑251 (MQ=255)
                                                                                                                                     tGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                        >  2:218271/1‑250 (MQ=255)
                                                                                                                                       aCAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                     >  1:106640/1‑251 (MQ=255)
                                                                                                                                                               gAGCAGTGGGGGCGCTCCCCCTTCCGGCGGGGTGCCGATATCGATAAACAGCGT‑CGGGCCATTGGGGATTGGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                              <  1:155554/251‑1 (MQ=255)
                                                                                                                                                                       cgtgcggtattgcttcgggtgttgtgCCGTTACCGGTAAACAGCGT‑CGATGGATTGGGAATTTGCCCCGCGACT‑TTGGGGACGGGGCGCTTTTTTTC‑TTCGTGCTGGGTGTTCTCGCGATCGTTCCCCGCCGTTTCAGCCGAACTGTAATGCAGCGCGGCACGGCCGTGATGG‑CGCAGGATTAGGTTTTGCTCCGCCGGCTCGTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                       <  1:29123/227‑1 (MQ=255)
                                                                                                                                                                                        ggCGGGGTGCCGATATCTATAACCAGCGT‑CGAGCCATTGGGGTTTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGGCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                     <  1:159534/251‑1 (MQ=255)
                                                                                                                                                                                        ccttttttcccGATATCGAAAACCAG‑GTGCGACCCATTGGGGATTGGCCCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                     <  2:146329/242‑1 (MQ=255)
                                                                                                                                                                                               tgCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                     >  1:146329/1‑244 (MQ=255)
                                                                                                                                                                                               tgCCGATATAGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                              >  2:94007/1‑251 (MQ=255)
                                                                                                                                                                                                      atcaaTAAACACCGT‑GGGGCCAGTGGGGATTTGCCCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                       <  2:201900/247‑1 (MQ=255)
                                                                                                                                                                                                            aaaCAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGCTGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCAACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCCGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCACTGCGGTACTGACATAACCCTGATGCTTAAcc                                  >  1:32141/1‑250 (MQ=255)
                                                                                                                                                                                                                 ggcttagatctttcggggATTTCCTCCGCCACT‑TGGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTCTCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                             <  2:204411/237‑1 (MQ=255)
                                                                                                                                                                                                                    ccggggcctgtGGGGATTTGTCCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                          <  2:240226/241‑1 (MQ=255)
                                                                                                                                                                                                                        ggcattgggggATTTGCTCCGCC‑CTGTTGCGGACGATGCGCTTTTTTTC‑TTGGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                      <  1:92123/244‑1 (MQ=255)
                                                                                                                                                                                                                         gCCATTGGGGCTTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                     >  1:75672/1‑251 (MQ=255)
                                                                                                                                                                                                                         gCCATTGGGGATTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                     <  1:58096/251‑1 (MQ=255)
                                                                                                                                                                                                                              gggggtgtggCTCCGCCATT‑GTGCGGACGGTGCGCTCTTTTTT‑TGCGGTCTGGGTGGCCTTGGGATCGGGCCACGGCGTGTTAACCGAACGGGAAGGCAGGGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt                 <  2:203692/241‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTTGCTCCGCCACT‑TTGCGGACGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGCTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg           >  1:114051/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cccGCCCCT‑TTGCGGACGATGCGCTTTTTTGC‑TTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt     <  1:63219/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cGCCACT‑TTGCGGACGTTGCGCACTTTTTC‑TTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGATGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt   <  1:226745/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ccctcT‑TTGCGACCGATGCGCT‑TGTTTCATTCGGTTTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg  <  2:172527/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 t‑ttGCGGGCGATGCGCTCTTTTTC‑TTCGGTCTGCGTGGCCTAGCGATCGTGCCACGGTGTGTTAATCCAACTGGGAGGCAGCGCCGCACCGACATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCAGCGCGCTGACGGAGACATGCTCTACCAGCCCTTCGGTACTGACATAACCCTGCTGTTTAAcccatt                              >  2:240200/1‑214 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
GAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGT‑TAGCAACGTTGAGATTGT‑TGGTGACAATGCGCAAAT‑TGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGT‑CGAGCCATTGGGGATTTGCTCCGCCACT‑TTGCGGGCGATGCGCTCTTTTTC‑TTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTG  >  NC_000913/3560118‑3560599

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: ‑