Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,223,638 T→C intergenic (‑10/‑190) iclR ← / → metH DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,6380TC97.8% 143.1 / ‑5.8 46intergenic (‑10/‑190)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (12/33);  total (13/33)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTA  >  NC_000913/4223393‑4223873
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           
aGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCCCGGGTCGCGGGGAATGG‑TAGGACCATGAC‑AGGC‑TCcttttt                                                                                                                                                                                                                                        >  1:277472/1‑251 (MQ=255)
      ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTAt                                                                                                                                                                                                                                  <  1:277929/251‑1 (MQ=255)
         accggACGCCCGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCTGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGt                                                                                                                                                                                                                               <  1:24774/247‑1 (MQ=255)
             gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTTATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCTTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCGTTGTCTGTATCGTTGaa                                                                                                                                                                                                                           <  2:513689/251‑1 (MQ=255)
              aCGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATGGGGTAACCCTGCTTGTTTCGCCAGTTCCTTGAGTGCCCCGCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGGCGCTGTTGCGGTGGCAGCGGCTGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCGGTATCGTGGaaa                                                                                                                                                                                                                          <  1:264237/251‑1 (MQ=255)
                 cACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCGCGTTTTGCGGG‑GATGGGGGCGGCCATTACCAG‑C‑TCCTTTTTTTGTATCGTGGTAAttc                                                                                                                                                                                                                       >  1:236780/1‑249 (MQ=255)
                    gAAACCCTGCTGTTGCATCGTTGTTAGTAGGTGGTGGGTCGTGGAATTTGGTAACCCGGCTTGTTGCGCTAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCGCGTTTCGCGGG‑AATGGGTGCTACCATGTC‑AGTC‑TCCTTTTTCTGTATCGTGGAAAtcattt                                                                                                                                                                                                                    <  1:376602/251‑1 (MQ=255)
                    gAAACCCTGCTGTTGCAGCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAAtcattt                                                                                                                                                                                                                    <  2:128687/251‑1 (MQ=255)
                     aaaCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTTCCATTGGATTCGGCAATCCACTCCAGTAAATTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTGTCTGCCGCGTTTCGCGGG‑AAGGGGTTCGACCAGGAC‑AGTC‑TCCGTTTTCTGTATCGTGGAAAtcatttt                                                                                                                                                                                                                   <  2:30879/251‑1 (MQ=255)
                         ccTGCTGTTGCATCGTTGTTAGCAGACGATGGGTCGTGGAATTGTGGAACCACGCTGGTTGCGCCAGTCCCGTAAGTGCCGCACTGCCATTGGATTCGGCAATCCACTCCAGGAATGTCAGGCCAGGCGTAAAAAGCTGAACCGTTCCAGTCGCTGGAGCGTTGGCAACGGCGGGTTTTCTGCCGCTTTTCGCGGG‑AATGTGTGCGACCATGGC‑AGTC‑TCCTGTTTCTGTATCGTGGAGAagattttcatt                                                                                                                                                                                                               <  2:17431/251‑1 (MQ=255)
                                    aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTCGCGGG‑GATGGGGGGGACCATGAC‑AGGCGTCTTTTTT‑TGTATCGGGGAGATCatttttattttttttttt                                                                                                                                                                                                      >  1:389936/1‑243 (MQ=255)
                                       gTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTCATTGTTAGCTa                                                                                                                                                                                                 <  1:506365/251‑1 (MQ=255)
                                               cAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGCTTTCTGCCGCGTTTAGCGGG‑GATGGGTGCGGCCCTGAC‑AGTC‑TCCTTTTTCTGTTTCTTGGGAATCATTTTCATTTTTATTGTTAGCTTATGGAAt                                                                                                                                                                                          >  1:375594/1‑250 (MQ=255)
                                                                aaTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGTTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCg                                                                                                                                                                        <  2:375894/251‑1 (MQ=255)
                                                                     ggTAACCCGGCTTGTTGCTCCAGTTCCTTGAGTGCCACACTGCCATTGGATTCGGCAATCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTTCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAGAGTTACTGAACGGATCCGATGAg                                                                                                                                                                   <  2:279919/251‑1 (MQ=255)
                                                                          cccGGCTTGTTGCGCCATTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAAt                                                                                                                                                              <  1:149462/251‑1 (MQ=255)
                                                                                 tgttgCGCCAGTGCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCTCTCCAGTAATTTCAGGCCATGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGTTGGTACCATGAC‑AGTG‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGa                                                                                                                                                         <  1:253210/249‑1 (MQ=255)
                                                                                      cgcCAGTCCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGGCGCGGTTTGCGGG‑GATGGGGGCGGCCATGAC‑AGGC‑TCCTTTTTTTTTTTCGTGGGAATCATTTTTCTTTTTTTTTTTATCTAATGCAATAAGTACTTCACTGGTCCCATTAGTTAATGTTGataaaaat                                                                                                                                                  >  2:153696/1‑244 (MQ=255)
                                                                                                   aGTGCCACACTGCCATTTGATTCGGCAATCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTTTATCGTGGAAATCACTTTCATTTTTATTGTTATCTAATGCAATAGTTTCTGAACTGATCCGATGTGTTAATGTTGAACAAATCTCATGTTGCgtgg                                                                                                                                     <  1:431152/251‑1 (MQ=255)
                                                                                                      gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg                                                                                                                                   <  2:381070/250‑1 (MQ=255)
                                                                                                      gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGtggtgg                                                                                                                                  <  1:380512/251‑1 (MQ=255)
                                                                                                              gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCttt                                                                                                                           <  2:113988/250‑1 (MQ=255)
                                                                                                                    ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcac                                                                                                                     <  1:202534/250‑1 (MQ=255)
                                                                                                                         cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGATCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGc                                                                                                               <  2:33888/251‑1 (MQ=255)
                                                                                                                          ggCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaact                                                                                                                                                                              <  1:420030/187‑1 (MQ=255)
                                                                                                                           gCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGt                                                                                                             <  1:239436/251‑1 (MQ=255)
                                                                                                                                         tAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGGGCGACCATGAC‑AGTC‑TCCTTTTTCTGTTTCGTGGAAAACATTTTCATTTTTTTTGTTTGCTAATACAATAGTTACTGGACTGATCCGCTGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGGCGCTTTTAACACCGATTCGTTTATGCCCGTATgg                                                                                               >  1:136232/1‑251 (MQ=255)
                                                                                                                                                aGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGG‑GATGGGGGGGGCCCTGAC‑CGTC‑TCCTTTTTTTGTATCGTGGGAATCCTTTTTCTTTTTTTTGTTTGCTTATTCAATAGTTTCTGAACTGGTTCGGTGAGTTTATGTTTAACAAATCTCATTTTTCGTGGGGGTCGCTTTTTCCACCAATGCGTTTTTGCCCGTATGGgttgttt                                                                                        >  2:166412/1‑250 (MQ=255)
                                                                                                                                                            aaaGACTGAACCGTTCCAGGCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGTTGCTACTATGTC‑AGTC‑TCCTTTTTCTGTATCGTTTAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATTCGTTTATGCCAGTATGGTTTGTTTAAGTTttatta                                                                             <  2:375594/250‑1 (MQ=255)
                                                                                                                                                               gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATc                                                                         <  1:166412/251‑1 (MQ=255)
                                                                                                                                                                 cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCTCGGG‑GATGGGTGCGACCCTGAC‑AGTC‑TCCTTTTTCTGTATCGTTGGAATCATTTTCATTTTTTTTTTTTACTTATGCAATAGTTACTGAACTTATCCCGTGAGTTAATTTTTAACAAATCTCATTTTGCGTGGTGGTCGCTTTTTCCACAGATGCGTTTAttcc                                                                                                      >  1:364440/1‑217 (MQ=255)
                                                                                                                                                                  tGAACCTGTCCAGTCGCTGGTGCGGTGGCCACGGCGGGTTTTTTTCCGCGGTTTGAGGG‑GATGGGTGCGGCCATTGC‑AGTC‑TCCCTTTTTTGTATTGGGGGAATCATTTTTATTTTTTTTGTTAGCTAATGCCAAAGTTTATGCAATGGTCCCCTTGGTTAATGTTGGACAAATCCCCTGTTGCGTGGGGGTCCGTTTTTCCACAAATGCGTTTATGCCAGTATGGTTTTTTGAAtttttttttaaattgg                                                                      >  2:364038/1‑240 (MQ=255)
                                                                                                                                                                      ccTGTCCAGTCGTTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCACGGG‑AATGGTTGCGACCATGAC‑AGTC‑TCCGTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGAAATAGTTACTGAACTGAGCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTAGCTTTTACCACAGATGCGTTTAGGCCAGTATGGTTTGGTGAATTTTTATTAGACCTGGGTTg                                                                  <  2:185095/251‑1 (MQ=255)
                                                                                                                                                                         gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                                                <  1:84851/250‑1 (MQ=255)
                                                                                                                                                                          tcaaGTCGCTGTTGCTGTGGCACAGGCTGGGTTTCTGCCGCGTGTCGCGGG‑AATGGGGGCGACCATTAC‑AGTC‑TCCTTTTCCTGTATCGTGTAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTTGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTGTATTAAATCTGTGTTGAGc                                                               <  1:377519/247‑1 (MQ=255)
                                                                                                                                                                           ccAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCgt                                                             <  1:23659/251‑1 (MQ=255)
                                                                                                                                                                                 gCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTTGCGGG‑AATGGGGTAGCCCCTGAC‑AGTC‑TCCTTTTTCTGTATCTGGGGAATCATTTTCATTTTTTTTTTTAGCTAATGTAATAGTTAATGAACTGATCCGATGAGTTAATGTTGAAAAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGGGCGGGTCggg                                                       >  1:15736/1‑251 (MQ=255)
                                                                                                                                                                                             gCAACGGCGGGGTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGTCGATGAT‑AGTC‑TTCTTTTTCTGTATCGTGGAAAGTATTTTCATTTTTATTGGTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGGTAATGTTGAACAAAGCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGag                                           <  2:44352/251‑1 (MQ=255)
                                                                                                                                                                                                        ttttCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTTATCCGATGGGTTAATGTTGAACAAATCGCATGGTATTTGTTGGTCGGTTTTTCCACATATTCGTTTTAGCGAGTATGGCGTGTTGAATATGTATATAATCTGGGTTGAGCGTGTTGTGAGTAAGGGTGCGCAGCAAAGTg                                 <  2:255597/250‑1 (MQ=255)
                                                                                                                                                                                                        ttttATGCCGCGTTTCGCGGG‑CATGGGTGCGACCCTGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATCGTTACTGAACTGATCCGATGAGTTAATGTTGAGCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGATCGAGTAGGGCGCAAGTGTGAGCAGCAAAGTgg                                >  2:269084/1‑251 (MQ=255)
                                                                                                                                                                                                           tCTGCCGCGTTTAGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGAAGCAAAGTGGaac                             >  2:361671/1‑251 (MQ=255)
                                                                                                                                                                                                                cgcgTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGACTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtg                        >  1:274520/1‑251 (MQ=255)
                                                                                                                                                                                                                         gggg‑AATGGGGTCGACCAGGAC‑AGGC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGGTTGTTGAATTTTTAGTAAATCTGGTTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAg               <  2:64054/250‑1 (MQ=255)
                                                                                                                                                                                                                           gg‑AATGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt             <  1:7574/251‑1 (MQ=255)
                                                                                                                                                                                                                               aTGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaaa          <  1:198927/251‑1 (MQ=255)
                                                                                                                                                                                                                                tGGGTGCGACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt         >  1:511492/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gACCATGAC‑AGTC‑TCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATTGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTa  <  2:236780/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           
AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACCATGAC‑AGTC‑TCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTA  >  NC_000913/4223393‑4223873

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