Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,332,310 C→T G321G (GGG→GGA)  basS ← sensor histidine kinase BasS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,332,3100CT95.3% 137.1 / ‑4.2 43G321G (GGG→GGA) basSsensor histidine kinase BasS
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (40/1);  total (40/3)
Fisher's exact test for biased strand distribution p-value = 3.32e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAACCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGGCGCTCTCT  >  NC_000913/4332064‑4332556
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
cACTTAATCGCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCGTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAAtccc                                                                                                                                                                                                                                                     <  2:28518/250‑1 (MQ=255)
 aCTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGAGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑GTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAg                                                                                                                                                                                                                                                   <  2:260730/251‑1 (MQ=255)
     aaTCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGAAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGAcc                                                                                                                                                                                                                                                >  1:536135/1‑250 (MQ=255)
             aCGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGcc                                                                                                                                                                                                                                       >  1:440009/1‑251 (MQ=255)
                cgcATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAAAAGAAATTATGAAAGTTTAACTACAGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATCTCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCAt                                                                                                                                                                                                                                     >  1:49451/1‑250 (MQ=255)
                              cagtttAACGGAGGAGAGTGCAATGAAAAAACGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGTCCC‑GCGTGCCGGCAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCACCTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGCCTGTCCATAc                                                                                                                                                                                                                      >  2:544438/5‑251 (MQ=255)
                                     aCGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACaa                                                                                                                                                                                                                >  1:159708/1‑250 (MQ=255)
                                       ggaggaGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTa                                                                                                                                                                                                                   >  2:451741/1‑245 (MQ=255)
                                           gagaGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCttt                                                                                                                                                                                                         >  2:181863/1‑251 (MQ=255)
                                               gTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTc                                                                                                                                                                                                     >  1:305706/1‑251 (MQ=255)
                                                ggCAATGAAAACCCGTGTTTATAAAATTTTACATTCCTTTTCCAGCGTTCTGGTGCGCAGCTGTTTTTTTTATAT‑TGGTTTTCCCACGTAGTGATTTTTCTTCA‑CCCTTCCCAGGCCCCGCGTGCCGGAACTCTCTCGCCTG‑TTTTGCAA‑AAAAACCGCCCCT‑TGAGGCATCTGTGAGGGGGGGTTAACAA‑TATTTAACCCCAGCCCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCa                                                                                                                                                                                                    <  1:542128/250‑1 (MQ=255)
                                                      gAAAAACCGGGTTTATGACAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGCAAAACTGCCC‑G‑TGATGCAACTGTGTACTGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTccc                                                                                                                                                                                              >  2:113947/1‑251 (MQ=255)
                                                        aaaaCCGTGTTTATGAACGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCg                                                                                                                                                                                            >  1:205511/1‑251 (MQ=255)
                                                          aaCCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGCAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAAAGGCTGTCCATACGTACAAACGCTTTACTCAACTCCACGc                                                                                                                                                                                           >  1:272655/1‑250 (MQ=255)
                                                                       aaaGTTTAACTACCGTGTTCCGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCtcatc                                                                                                                                                                              >  1:184161/1‑250 (MQ=255)
                                                                       aaaGTTTAACTACCGTGTTCAGCGTGCGGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCACGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGCTGCAACTGTGTAATGCGGCTGACAC‑TACTTAATCCCAGACCACTACCGCCATAACGGCTGTCCATCCGTACAAACGCTTTACTCAACTCCCCGCATTTACTCtcatc                                                                                                                                                                              >  2:565543/1‑250 (MQ=255)
                                                                        aaGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCACTCCCGCCATAACGGCTGTCCATACGTACAAACCCTTTACTCACCTCCCCGCATTTACTCtcatca                                                                                                                                                                             >  2:95646/1‑250 (MQ=255)
                                                                                    cGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCtt                                                                                                                                                                >  2:322077/1‑251 (MQ=255)
                                                                                     gtgtTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTc                                                                                                                                                               >  1:584006/1‑251 (MQ=255)
                                                                                           aGCGTGCTGGTGGTCAGCCGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTc                                                                                                                                                         >  2:375573/1‑251 (MQ=255)
                                                                                                gCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGACAAACTGCCC‑G‑TGCTGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACTTACAAACGCTTTACTCAACTCACCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTg                                                                                                                                                    >  1:15700/1‑251 (MQ=255)
                                                                                                 cTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGc                                                                                                                                                   >  1:174388/1‑251 (MQ=255)
                                                                                                        tCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCACGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTTTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGAcc                                                                                                                                            >  1:331325/1‑251 (MQ=255)
                                                                                                                       atatTTGGTTT‑CCCACGTACTTATTTTGCTGCAGCCGGCCCCCGGCCG‑GGGTGCCGGAGGTCTCT‑GCCGGGTTTGGCAG‑AAAAACCCCCG‑GGTGATCCACCTGTGTAATGCGCCT‑ACAAATCCGTACCCCCACACCACTCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCAGCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCAtcttct                                                                                                                             <  1:375573/251‑1 (MQ=255)
                                                                                                                          tCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGCCCGCTCCGTCATCTTCTTGc                                                                                                                          >  1:113145/1‑251 (MQ=255)
                                                                                                                                        gTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCACCTGcaat                                                                                                                                                >  1:187133/1‑213 (MQ=255)
                                                                                                                                                       cAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTtgct                                                                                                      >  1:141159/1‑243 (MQ=255)
                                                                                                                                                        aGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCgg                                                                                             >  1:492340/1‑251 (MQ=255)
                                                                                                                                                        aGCCGTACCCAGGCAC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTAATCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTATTCATAATGTTTCTGCCTTGccg                                                                                             >  2:442739/1‑249 (MQ=255)
                                                                                                                                                               cccAGGCCC‑GCGTGCCGGCAGTCTCTTGCCGG‑TTTTGCAGGAAAACCTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACGAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCCTCTTCAACTGCAATGACCGCTcc                                                                                                                                       >  1:10941/1‑202 (MQ=255)
                                                                                                                                                               cccAGGCCC‑GCGTGCCGGCAGTATCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACCA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAACCGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCATTGCGGGCTGTa                                                                                       >  2:539271/1‑250 (MQ=255)
                                                                                                                                                                             gCCGGCAGTCTCTTGCCGG‑CTTTGCAGGACACACTGCCC‑G‑TGAAGCAACTGTGTAATGCGGCTGACAC‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGtttt                                                                         >  2:79387/1‑251 (MQ=255)
                                                                                                                                                                              ccGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATAAGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCCTCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGTTGGGCGTTTTc                                                                        >  2:475449/1‑251 (MQ=255)
                                                                                                                                                                                ggAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACGGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGCTGGGAGTTTTCt                                                                       >  1:480554/1‑250 (MQ=255)
                                                                                                                                                                                                      cAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATcc                                                 >  1:118737/1‑251 (MQ=255)
                                                                                                                                                                                                        ggAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCgc                                               >  1:206815/1‑251 (MQ=255)
                                                                                                                                                                                                           aaaaCTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCgcagc                                            >  2:264611/1‑251 (MQ=255)
                                                                                                                                                                                                             aaCTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCgcagc                                            >  1:576153/1‑249 (MQ=255)
                                                                                                                                                                                                                tGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGt                                       >  1:362567/1‑251 (MQ=255)
                                                                                                                                                                                                                   cc‑G‑TGATGCCACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACAAATCCCGCCATAACGGCTGACCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTAATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCCGCAGGGTTGc                                    >  2:25122/1‑251 (MQ=255)
                                                                                                                                                                                                                           tGCAACTGTGTAATGCGGCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCAcc                              >  2:8072/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ggCTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTcg              >  1:110603/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cTGACAA‑TACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTcggc            >  2:380907/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        ttAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGGCGctctct  >  1:391027/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
CACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT‑GCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCC‑GCGTGCCGGAAGTCTCTTGCCGG‑TTTTGCAGGAAAAACTGCCC‑G‑TGATGCAACTGTGTAATGCGGCTGACAA‑TACTTAACCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGGCGCTCTCT  >  NC_000913/4332064‑4332556

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