Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,522,187 C→T A189T (GCG→ACG)  hofM ← DNA utilization protein HofM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,522,1870CT94.1% 178.7 / ‑5.9 68A189T (GCG→ACG) hofMDNA utilization protein HofM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/39);  minor base G (3/1);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATAT  >  NC_000913/3521946‑3522420
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
gggCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCa                                                                                                                                                                                                                                   <  1:178660/251‑1 (MQ=255)
         attaatgcgCGGTTTCACCAGTACGCCTTCCCAAGCGCCAGCCCCGGCGCGATGGCAAAGTCGCCGCCTGGGGTCGGTAGCGGCGGCTGAGGAACAGTAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATGTCGCGACGCTTGCGGGATCA‑ACGGAAATCGCGGCCGCCAGCGCCTACGCGCTAGTCATCCCCC‑CCGCCAGTGTGCGCCCCCAGGTAGAGCTCGTCGTCCTCAGCCACTGTTC‑Gtt                                                                                                                                                                                                                         <  2:412741/243‑1 (MQ=255)
               tGGCGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCAc                                                                                                                                                                                                                    <  1:151830/250‑1 (MQ=255)
                  cGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGCCCCCCCCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACgcca                                                                                                                                                                                                                <  1:147558/251‑1 (MQ=255)
                     gTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCCCCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTAGAGCGCGTCTTCCACAGCCACTGTTC‑GTTATCACgccaggc                                                                                                                                                                                                             <  2:254990/251‑1 (MQ=255)
                     gTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACCGCCTCCCAGGGATCAAATCCACCTTCGCCACCTATCGCGCCGCCTTCGGGATCA‑ACGGATCACGCTGCCGCCCGCTCCTTCGCGCTAGTCCTCCCCC‑CCGCCAGTTTGCGCCCCCAGCTCTAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCCCgccaggc                                                                                                                                                                                                             >  2:188942/1‑251 (MQ=255)
                       tcatcaGTCCGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACgccaggcc                                                                                                                                                                                                            <  1:120399/250‑1 (MQ=255)
                                cGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGATGAAGAACAGAAACGGCCTCACAGGGATCAAATCCGCCTTCCCCACAAATCACGACGCTTTCGGGATCA‑ACGGATAACGATGCCGCCCGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAAATTGCGCCCCCAGCTATAGAGCGTCGTCCACAACCACCGTTCCGTCA‑CACCCCAGGCCAGaactttct                                                                                                                                                                                                  >  1:19014/1‑244 (MQ=255)
                                 gCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCCGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGGTTGCGCCCCCAGCTATAGCGCGTCGGCCACAGCCACTGTTC‑GTTATCCCTCCAGGCCAGACATTGCtg                                                                                                                                                                                                 >  1:101214/1‑251 (MQ=255)
                                  ccTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAGAAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCC‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCtg                                                                                                                                                                                                 >  1:197129/1‑250 (MQ=255)
                                         aaGCGCAAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGCCGCACAGAAACGGCCTCCCAGGGATCCAATCCCCCTTCGCCACATATCGCGACGCTTTCGGGATCC‑ACGGACAACGCTGCCGCCAGCTCCTTCGCGCTACTCATCCCCC‑CCCCCCGTTTGCGCCCCCAGCTACAGCCCCTCGTTCACAGCCACCGGTC‑GTTATCACccccggcccgaaaatgctggtgcggcg                                                                                                                                                                                         >  2:193974/1‑232 (MQ=255)
                                             gcCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGAGACGCTTTCGGGATCC‑ACGGATAACGCTGCCGCCCGCTCCTTCGTGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATCGCGCGTCGTCCCCAGCCACCGTTC‑GTTATCAAGCCAGGCCAGACCTTGCTGATGAGAAGGTa                                                                                                                                                                                      >  1:53678/1‑250 (MQ=255)
                                                     cAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATACCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCa                                                                                                                                                                             <  1:277165/251‑1 (MQ=255)
                                                        cgcgATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGCCGAACAGAAACGGCCTCCCAGGGAGCAAATCAGCCTTCGCCACATATCGCGCCGCTTTCGGGGTCA‑ACGGATAAAGCTGCCGCCCGCTCCTTCGCGCCAGTCATCACCA‑CCCCCCGTTTGCGCCCCCACCTATAGCGCGTCGTCCCCAACCACCTTTC‑GTTATCACTCCAGGCCCGACATTGCTGGTGAGAAGGGAAAAAAGGCAggg                                                                                                                                                                          >  2:76188/1‑250 (MQ=255)
                                                          cgATGGCAAATTCTCCACCGGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAAt                                                                                                                                                                        <  1:251593/251‑1 (MQ=255)
                                                           gATGGCAAAGTCTCCACCGGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCACCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATc                                                                                                                                                                       <  1:188942/251‑1 (MQ=255)
                                                              ggCAAAGTCACCCCCTGGCGTCGGTAGCGGCGGCTGACGAACATAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCAGCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGc                                                                                                                                                                     <  1:363510/250‑1 (MQ=255)
                                                                cAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAAAAGAAACGGCCTCCCCGGGATCAAAACCGCCTTCGCCACATATCGCGACGCTTTCGGGGTCA‑CCGTATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑ACGCCAGTTTGAGCCCCCAGCTATAGCCCGTCGGCCACAGCCCCCGTTC‑GTTTTCACCCCCGGCCAGCCATTGCTGATGAGAAGGTaaaaaaagaaagaacagctgt                                                                                                                                                                  >  2:315547/1‑236 (MQ=255)
                                                                cAAAGTCCCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCGCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGt                                                                                                                                                                  <  1:182810/251‑1 (MQ=255)
                                                                   aGTCTCCACCGGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGAACAAATCCGCCTCCGCCACATACCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAAt                                                                                                                                                               <  2:46766/251‑1 (MQ=255)
                                                                      ctcCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GGTATCACGCCAGGCCAGACATTGCTGATGTGAAGGTAAAAAAGGCAGGAATCGCTGTAATGcc                                                                                                                                                            >  1:4663/1‑250 (MQ=255)
                                                                               gcggctGTAGCGGCGGCTGACGAACAAAACCGCCCTCCCGGGGATCAAACCCGCCGTCGCCACATATCGCGCCGCTTTCGGGACCA‑ACGGAGAACGCTGCCGCCAGCTCCTTCGCGCTAGTCACCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATcc                                                                                                                                                   <  2:293292/251‑1 (MQ=255)
                                                                               gcggcgGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATACCGCGACGCTTGCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCGCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATcc                                                                                                                                                   <  1:23128/251‑1 (MQ=255)
                                                                                cggcggTAGCGGCGGCTGACGAACAGAAACGGCCCCCCCGGGATCAAATCCGCCTTCGCCCCATATCTCGACGCTTTCGGGATCA‑AAGGAAAAAGCTGCCGCCAGCTCCCTCGCGCTAGGCATCCCCA‑CCGCCAGTTTTCGCCCCCCACTATAGAGCGTCGGCCACAGCCACTGTTC‑GTTTTCCCGCCAGGACAGACATTGCTGATGGGAAGGGACAAAAAGCAGGAAACGCCGTAATGCACTGGCAtacc                                                                                                                                                  >  2:314474/1‑248 (MQ=255)
                                                                                  gcggTGGCGGCGGCTGCCGAACAGAAACGGCCTCCAAGGGATCAAACCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAAGGCTGCCGCCAGCTCCTTCGCGCTAGTCACCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCggg                                                                                                                                                <  2:380764/251‑1 (MQ=255)
                                                                                    ggTAGCGGCGGCTGCCGAACAGAAACGGCCCCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATACCGCGGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGTCAGTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGt                                                                                                                                              <  2:212455/251‑1 (MQ=255)
                                                                                    ggTAGCGGCGGCTGAAGAACAGAAAAAGCCTCCCCGGGAACAAATCCGCCTTCGCCACATATCGAGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAACTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAACCACTGTTC‑GTTATCACGCCAGGCCAGCCATTGCTGATGAGAAGGTAAAAAAAGCAGGAATAGCTGTAATGCACTGGCATCCggggg                                                                                                                                              >  2:342618/1‑250 (MQ=255)
                                                                                       aGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAAGCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGa                                                                                                                                            <  1:294482/250‑1 (MQ=255)
                                                                                         cggcggCTGACGAACAGAAACCGCCTCCCAGGGCTCACATCAGCCTTCGCCACATCTCGACACGCTTTCGGGATCC‑ACGGCTAACCCTGCCGCCAGAACCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCCGACATTGCTGATGAGCAGGTAACAAAGGCAGGAATCGCTGTAATGCCCCGGCATCCGGGGTGATcg                                                                                                                                         >  2:411067/1‑251 (MQ=255)
                                                                                                   cGAACAGAAACGGCCCCCCAGGGAACAAATCCGCCTTCGCCCCAAATCGCGACGCTTTCGGGATCA‑ACGGATATCGCTGCCGCCAGCTCCTTCGCGCATGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCGCTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATg                                                                                                                               <  2:31052/251‑1 (MQ=255)
                                                                                                      aCAGAAACGGCCTCCCAGGGATCAAACCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAac                                                                                                                            <  1:149177/251‑1 (MQ=255)
                                                                                                         gAAACGTCCTCCCAGGGATCAAATCCGCCGTCGCCACTAACCGCGACGCTTTCGGTATCA‑ACGGAGAACGCTTCCGCCAGCTCCTTCGCGCTAGTCACCCCCA‑CCGCCAGGTTGCGCCCCCAGCTATAGCGCGTCGGCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATTCACTGGCATCCGGGGTGATCGCGCTCACATGAACACg                                                                                                                         <  2:340245/251‑1 (MQ=255)
                                                                                                            aCGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCa                                                                                                                       >  1:260486/1‑250 (MQ=255)
                                                                                                              ggCCTCCCAGGGATCAAATCCGCCTTCGCCACATACCGCGACGCTTTCGGGATCA‑ACGAATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCACCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAAcc                                                                                                                    <  1:205878/251‑1 (MQ=255)
                                                                                                                            caaacccactttcctcacaTACCGCGACGTTTTCGG‑AACACACGGTTAACTCTGCCGCCAGGTCCTTCGCTCTAAACACCCCCT‑CCGCCAGTTTGCGCCCCCAGATATAGCTCGTCGTCCACAGCCACTGTTC‑GGTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCTATGGCATCCGGGGTGATCGCGCTCACATGGACACGCAACCTTTCTGCCagagga                                                                                                      <  2:12970/236‑3 (MQ=255)
                                                                                                                             aaaccctcctgcGCCACATAGCGCGACGCTTTCGTGTTCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTaa                                                                                                     <  1:266329/240‑1 (MQ=255)
                                                                                                                                      ttCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTg                                                                                            >  2:359572/1‑251 (MQ=255)
                                                                                                                                       gtgccctctttcgtgactcTTTCGGGCTCA‑CCGGATCACGCTGCAGCCCTCTCCTTCGCGCTAGGCA‑CCCCAGCCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCGCTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCGGTAATGCACTGGCGTCCGGGGTGATCGCGCTCACAGGAACACGCAACCTTGCTTCCAGAGTAAGCAGCGtt                                                                                             <  2:27700/231‑1 (MQ=255)
                                                                                                                                         gccccaTATCGGGACGCTTGCGGGATCA‑ACGGATAACGCGGCCGCCAGCTCCTTCGCGCTAGTCACCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCGTTCGGCCAGAGTAAGCGGCGGTGCCa                                                                                         <  2:179357/247‑1 (MQ=255)
                                                                                                                                                   gctaCGCTTTCGGGATCA‑ACGGATACCTCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTctttgct                                                                               <  1:359249/248‑1 (MQ=255)
                                                                                                                                                   gcgACGCTTTCGGGAGCA‑ACGGATAACGCGGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTGCTGCCAGAGTAAGCAGCGTTGCCAGCTctttgct                                                                               <  2:81174/251‑1 (MQ=255)
                                                                                                                                                           tgcGGGACCA‑ACGGAAAACGCTGCCGCCAGCTCCTCCGCGCTAGTCACCCCCA‑CCGCCAGGTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGTGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGGTGCCAGCTCTTTGCTTTgcgc                                                                         <  2:197648/247‑1 (MQ=255)
                                                                                                                                                             ggggtccA‑ACGGATAACGCTGCCGCCACCTCCTTCGCTCTAGTCACCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTTTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTTCTCGGCa                                                                     <  2:247491/245‑1 (MQ=255)
                                                                                                                                                             cGGGATCC‑ACGGACACCGCTGCCGCCAGCGCCTTCGCGCTAGTAACCCCCT‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCa                                                                     <  2:39973/251‑1 (MQ=255)
                                                                                                                                                               ggATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGt                                                                   >  1:111451/1‑251 (MQ=255)
                                                                                                                                                                 aTCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTc                                                                  >  1:176353/1‑250 (MQ=255)
                                                                                                                                                                   cA‑ACGTATAACGCTGCCGCCACCTCCTTCGCGCTAGTTATCCCCG‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCTGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCGCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACg                                                               <  1:379825/251‑1 (MQ=255)
                                                                                                                                                                   cA‑ACGGAAAAAGCTGCAGCCAACTCCTTCGCGCTCGTCATACCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTGATCACGCCAGGCCAGACATTGCTGATGAGACGGTACAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGACCACGCAACATTTCTGCCATCGTAAGCAGCGTTGCCTGCACTTCGCTTTGCGCGGCTGTCACg                                                               >  2:156856/1‑251 (MQ=255)
                                                                                                                                                                    a‑aCGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCACCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGt                                                              <  1:51145/251‑1 (MQ=255)
                                                                                                                                                                        ggATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTtat                                                           >  1:187630/1‑251 (MQ=255)
                                                                                                                                                                           tAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAgg                                                        >  1:325341/1‑251 (MQ=255)
                                                                                                                                                                            aaCGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGc                                                       >  1:277460/1‑251 (MQ=255)
                                                                                                                                                                                      ccAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGAGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtgagg                                             >  1:295284/1‑250 (MQ=255)
                                                                                                                                                                                            cctccgggcTAGTCACCCCCT‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTTTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCt                                       <  2:364629/244‑1 (MQ=255)
                                                                                                                                                                                            ccTTCGCGCTAGACATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCAGCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCt                                       <  2:197129/251‑1 (MQ=255)
                                                                                                                                                                                               tCGCGCTAGACATCCCCA‑CCGCCAGTTGGCGCCCCCAGATATATCGCGTCGTCCACAGCCACTGTTC‑GGTATCACGCCAGGCCAGACAGTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTGTGCTTTGCGCGGCAGTCACGTTATAGGCGGTACTGAGTGAGTCTTc                                     <  2:57827/250‑1 (MQ=255)
                                                                                                                                                                                                cgcgCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTACGCAGCGTTGCCGGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtgcggggagcttcgc                                   >  2:128933/1‑238 (MQ=255)
                                                                                                                                                                                                  cgcTAGTCATACCCA‑CAGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGCCATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCCCATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGAGAGTCTTTGCTg                                 >  2:399062/1‑251 (MQ=255)
                                                                                                                                                                                                   gcTAGTCATCCCCA‑CAGCCAGTTTGCGACCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCCGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCggggcatggtgcgtctttcctgt                                >  2:346918/1‑231 (MQ=255)
                                                                                                                                                                                                      aGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGCGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACtagaacgggagagccgcgggg                             >  2:63303/1‑231 (MQ=255)
                                                                                                                                                                                                      aGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACt                                                 <  1:63303/231‑1 (MQ=255)
                                                                                                                                                                                                         cATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCg                           <  1:141589/250‑1 (MQ=255)
                                                                                                                                                                                                             cccA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGTCCACCGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGCAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGCGGGAGTCTTTGCTGGAATcacaacg                      >  2:398837/1‑245 (MQ=255)
                                                                                                                                                                                                                       aGTTTGCGCCCCCAGCTATAGCGCGTCGTCCCCAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGATGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGACATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGAGCGGCAGTCACGTTATAGGCGGGACAGGGGTTGTCTTCGCTGTAATccgagctcaaggagtc             >  1:355519/1‑236 (MQ=255)
                                                                                                                                                                                                                        gtgtgctcccccAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGTCATTGCTGATGAGAAGGTAAAAAAGGCAGGAACCGCTTTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAATCGCAGGTAGTcc            <  2:399808/244‑1 (MQ=255)
                                                                                                                                                                                                                               cccccAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAAGCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCa     <  1:163611/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccccAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GGTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCat    <  1:13659/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGCTATAGCGCGTCGTCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGGATGGGTGAGGCTTCGCTGGAAACAAAGCGCACGGGGTCCGGGACCatat  >  2:236247/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
GGGCAAAAAATTAATTGGCGGGTTCATCAGTACGCCTTCCCAAGCGCCAGCCCCAGCGCGATGGCAAAGTCTCCACCTGGCGGCGGTAGCGGCGGCTGACGAACAGAAACGGCCTCCCAGGGATCAAATCCGCCTTCGCCACATATCGCGACGCTTTCGGGATCA‑ACGGATAACGCTGCCGCCAGCTCCTTCGCGCTAGTCATCCCCA‑CCGCCAGTTTGCGCCCCCAGCTATAGCGCGTCGCCCACAGCCACTGTTC‑GTTATCACGCCAGGCCAGACATTGCTGATGAGAAGGTAAAAAAGGCAGGAATCGCTGTAATGCACTGGCATCCGGGGTGATCGCGCTCACATGAACACGCAACCTTTCTGCCAGAGTAAGCAGCGTTGCCAGCTCTTTGCTTTGCGCGGCAGTCACGTTATAGGCGGGACTGAGTGAGTCTTCGCTGTAATCGAAGCGCAGGGAGTCCGGATCCATAT  >  NC_000913/3521946‑3522420

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