Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,775 C→T N360N (AAC→AAT)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,7750CT98.5% 196.6 / NA 66N360N (AAC→AAT) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/34);  minor base G (0/1);  total (31/35)
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.88e-01

CGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTC  >  NC_000913/4637539‑4638018
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
cgcgAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑CAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAAcc                                                                                                                                                                                                                                       >  2:456790/1‑251 (MQ=255)
     aGGTCCGCCACCGGAAGTGGTGGCTCGTTTTGCCGACAACATTCTGACGCAGAATGCGCGTATGCAGGCATTGGTAGAAGCGTTACTACGCCAGGCAAGACTGGAGAATCGTCATGAAGGCGGTCTGACTGTTGTTGATGTGGTGGCATTAGTTCGCCGCGTGTGCG‑AAGCGCGCGCCGTGCGTTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGt                                                                                                                                                                                                                                  <  1:309990/251‑1 (MQ=255)
      ggTCCGCCACCGGAAGTGGTGGCTCGGTTTAACGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGCTAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTTCATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGtt                                                                                                                                                                                                                                 <  2:375677/251‑1 (MQ=255)
           gcaaccgGAAGTGGTTGCTCGTTTTAACGAAAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTTTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTgg                                                                                                                                                                                                                            <  2:445172/247‑1 (MQ=255)
                cgGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCTTTACttgggccg                                                                                                                                                                                                                       >  2:248915/1‑244 (MQ=255)
                  gAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCCGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAACAAACATCACTCTACATGTTATGCCCCCTGCGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAggg                                                                                                                                                                                                                     >  2:130333/1‑250 (MQ=255)
                   aaGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCGGGCATTGGTAGAAACGTTACTCCTCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTGCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGGCAGCG‑AAGCGCGCTCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAAGGTTGCTGCTGAACCGTCGTTATTGGAGCAGgcg                                                                                                                                                                                                                    <  2:557048/251‑1 (MQ=255)
                        gtggCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTAATGTTGCTGCTGCACCGGCGTTACTGGAGCAGGCGCTggg                                                                                                                                                                                                               >  1:388765/1‑251 (MQ=255)
                         tggCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCAACGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACAGGCGTTACTGGAGCCGGCGCTgggg                                                                                                                                                                                                              >  2:466440/1‑250 (MQ=255)
                                ttttatcgaaaactttCTGACGCAATGTGCGCGAATTCAGGCATTGGTAGAAACGTTACTTCGCAAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCCTCCGCG‑GAGCGCGCACCGTGCCGGTTGCAGAAAAAAACATCACCCTACCTGTTCAGCCCACTGCGGTTACTTTTACTGCCGATCCGGCGTTACTGGATAAGGCGCTgggggatttct                                                                                                                                                                                                       >  2:361622/15‑243 (MQ=255)
                                ttttACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGCAACGTTACTACGCCAGGCAAGACTGGAGACTCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCACCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTTCATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGTGAATTTAc                                                                                                                                                                                                       >  2:140084/1‑251 (MQ=255)
                                  ttACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGCGCAGGCGCTGGGGAATTTACTg                                                                                                                                                                                                     >  1:238888/1‑251 (MQ=255)
                                  ttACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTAAGGAAGTCGTAATGACTGTTGTTGATGAGGCGGCATTATTTCGCCGCGTAAGCG‑GAGAGCGCACCGTGCAGTTGGAAGAAACAAACACAACTCTACATGTTTTGCCCCCTGCGGTTAAAGTTGCTGCTGAACAGGCGTTACTGGAGCAGGGGCTggggagattaacg                                                                                                                                                                                                     >  2:226168/1‑241 (MQ=255)
                                     ccgACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGAt                                                                                                                                                                                                  <  1:456790/249‑1 (MQ=255)
                                         caacaTTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACg                                                                                                                                                                                              >  1:602001/1‑251 (MQ=255)
                                           acaTTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGcc                                                                                                                                                                                            >  1:65720/1‑251 (MQ=255)
                                                  gACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTGGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTTAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCg                                                                                                                                                                                        <  1:421242/248‑1 (MQ=255)
                                                   aCGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTTGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGAttt                                                                                                                                                                                    <  1:302929/251‑1 (MQ=255)
                                                   aCGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTTGTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGGATTTACTGGTTAACGCCTTCGGttt                                                                                                                                                                                    >  2:231555/1‑251 (MQ=255)
                                                    cGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCTTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCAGTATTTCGCCGCGTCAGCG‑AAGCGCGGACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACGGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGTAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGAtttt                                                                                                                                                                                   <  2:275530/251‑1 (MQ=255)
                                                        aaaTGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTAccc                                                                                                                                                                               <  1:231555/251‑1 (MQ=255)
                                                                  aTGCAGGCATTGGTAGAAACGTTTCTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCgg                                                                                                                                                                     <  1:521880/251‑1 (MQ=255)
                                                                          aTTGGTAGAAACGTTACTACGCCAGGCAAGACGGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGGTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCGGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCAtaa                                                                                                                                                             <  2:85916/251‑3 (MQ=255)
                                                                                aGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTa                                                                                                                                                        <  2:447361/250‑1 (MQ=255)
                                                                                  aaaCGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGAG‑CAGCGCGCACCGGGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGGTACTGGCGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAgc                                                                                                                                                     >  1:280968/1‑251 (MQ=255)
                                                                                   aaCGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACGGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACCGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAgcg                                                                                                                                                    <  1:627929/251‑1 (MQ=255)
                                                                                     cGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGcc                                                                                                                                                  <  1:181235/251‑1 (MQ=255)
                                                                                              gCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGAt                                                                                                                                         <  1:46855/251‑1 (MQ=255)
                                                                                                            gcgaATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGGATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGGGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTcgc                                                                                                                           >  2:17545/3‑249 (MQ=255)
                                                                                                                aTCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCTCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTcgccctt                                                                                                                       <  1:396591/251‑7 (MQ=255)
                                                                                                                aTCGTCAGAAAGTCGATCTGACTGTTGTTGATGTGGCGGCCTTATTTCGCCGCGTCAGCG‑AAGCGGGAACCGTGCAGGTGACCGAAAAAAACATCACTCTACATTTTATGCCCACTGCGGTTAATGTTGCTGCTGAACCGGAGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGCCGAATAACGCTCAGAGCCGTAGTGGCTCAGGCACACGTccccctt                                                                                                                       >  2:379221/1‑245 (MQ=255)
                                                                                                                         aaGTCGTTCTGACTGTTGTTGATGTGGCGGCCTTATTTCGCCGCGTCAGCGCAAG‑GCGCACCGTGCAGTTGGCAGCACAAAACATACCTCTACAAGTTATGCCCACTGAGGTTAATGTTGCTGCTTCACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACCCTAAGCCCCGGAGTGGGTCAGGAACAGGTGGCCCTTAAGGTGCTg                                                                                                              >  2:81394/1‑251 (MQ=255)
                                                                                                                             cGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTTCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGAt                                                                                                           <  1:313093/250‑1 (MQ=255)
                                                                                                                                          gTTGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACGCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTa                                                                                              <  2:194143/250‑1 (MQ=255)
                                                                                                                                            tGATGTGGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTAtt                                                                                            <  2:248124/250‑1 (MQ=255)
                                                                                                                                                gtgGCGGCATTCTTTAGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGGGGATCAGGAACACGTCGCCCTTAAGGGGCTGGATACCGGTAGTGGTATTCCTg                                                                                        >  1:483271/1‑250 (MQ=255)
                                                                                                                                                gtgGCGGCATTATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGGGGTTAATGTTGCTGCTGAACCGGCGATACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGAACAGGAACACGTCGCCCTTAAGGAGCTGGATACCGGTAGGGGTATTCCTgg                                                                                       >  1:623454/1‑250 (MQ=255)
                                                                                                                                                gtgGCGGCAATATTTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGGAGTAAAAAACATCACTCTCCATGTTATGCCCCCTGAGGTTACTGTTGCTGCTGCACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGAATAACGCTAAGCGCCGAAGGGGGTCAGGAACACGTCGCCCATAAGGTGCTGGATACCGGTAGTGGGATTcatg                                                                                        >  1:363722/1‑247 (MQ=255)
                                                                                                                                                          tattTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCTCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCt                                                                              <  1:259783/250‑1 (MQ=255)
                                                                                                                                                          tattTCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCt                                                                              <  1:536973/250‑1 (MQ=255)
                                                                                                                                                              tcgccgcGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCa                                                                          <  1:248915/249‑1 (MQ=255)
                                                                                                                                                               cgccgcGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCAc                                                                         <  1:408735/250‑1 (MQ=255)
                                                                                                                                                                 ccgcGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAACACCTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGCTTTTACCCCCGAGAGCGGTCGCATAACGCTACGCGCCGAAGTGGCTCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACga                                                                       >  2:57617/1‑249 (MQ=255)
                                                                                                                                                                 ccgcGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTa                                                                      <  1:348878/251‑1 (MQ=255)
                                                                                                                                                                 ccgcGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGGAGTGGCATTCCTGACTCCGCGCTGTCACGt                                                                       >  2:244836/1‑250 (MQ=255)
                                                                                                                                                                  cgcgTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTAt                                                                     <  1:180501/251‑1 (MQ=255)
                                                                                                                                                                       cAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTAttttt                                                                 <  1:221011/250‑1 (MQ=255)
                                                                                                                                                                          cG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGACGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCGGACTACGCGCTGTCACGTATTTTTGAAc                                                             <  2:222434/251‑1 (MQ=255)
                                                                                                                                                                          cG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCAGCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAAc                                                             <  1:590521/251‑1 (MQ=255)
                                                                                                                                                                                cgcgcACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCttt                                                        >  1:265195/1‑251 (MQ=255)
                                                                                                                                                                                   gcACCGTGCAGTTGGCAGAAAACCACCTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGACGGGGATCAGGACCACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGACAGCTTTTAc                                                     >  2:397785/1‑251 (MQ=255)
                                                                                                                                                                                   gcACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTAc                                                     >  1:391613/1‑251 (MQ=255)
                                                                                                                                                                                     aCCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCCCTGAGGTTAATGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGGAATTTCCTGGATAACGCCATCGATTTTAACCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACAAGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTActc                                                   >  1:289193/1‑251 (MQ=255)
                                                                                                                                                                                      ccGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTActct                                                  <  2:388765/251‑1 (MQ=255)
                                                                                                                                                                                            aGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGCACGCTTTTACTCTTTGccg                                            <  1:473832/251‑2 (MQ=255)
                                                                                                                                                                                                  cAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCa                                      <  2:375002/251‑1 (MQ=255)
                                                                                                                                                                                                          aacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGc                               >  1:329851/4‑250 (MQ=255)
                                                                                                                                                                                                           acaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCa                              >  1:620401/3‑250 (MQ=255)
                                                                                                                                                                                                            caTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaa                             >  1:218688/2‑250 (MQ=255)
                                                                                                                                                                                                             aTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCaaaa                           <  1:130333/251‑1 (MQ=255)
                                                                                                                                                                                                                  tctacgtGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCTGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTACGCGCCGAAGTGGATAAGGCACACGTCGCCCTTAAGGTGCTGGATACCGGTAGGGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCACAAAagc                       >  2:593237/7‑250 (MQ=255)
                                                                                                                                                                                                                        tGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTc                >  2:176506/1‑251 (MQ=255)
                                                                                                                                                                                                                         gTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCt               <  1:54816/251‑1 (MQ=255)
                                                                                                                                                                                                                                  actgAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGc      <  1:503000/248‑1 (MQ=255)
                                                                                                                                                                                                                                  actgAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCCGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATACCGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGTCGCTGGAGACCGGTCGTGGTATTCCTGCATACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGAAGCGGTCTGGGGTTGGc      <  1:572714/248‑1 (MQ=255)
                                                                                                                                                                                                                                      aGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTc  >  1:223231/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
CGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCG‑AAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTC  >  NC_000913/4637539‑4638018

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