Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 838,183 G→T intergenic (‑258/+7) ybiJ ← / ← ybiI DUF1471 family putative periplasmic protein/DksA‑type zinc finger protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913838,1830GT93.3% 186.0 / 4.0 60intergenic (‑258/+7)ybiJ/ybiIDUF1471 family putative periplasmic protein/DksA‑type zinc finger protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (28/28);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GCCATAGCAGCAACAACAGTATTGATAGTTTTCATAATTAATCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCC  >  NC_000913/837891‑838468
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                             
gCCATAGCAGCAACAACAGTATTGATAGTTTTCATAATTAATCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGtt                                                                                                                                                                                                                                                                                       >  4:276601/1‑301 (MQ=255)
  cATAGCAGCAACAACAGTATTGATAGTTTTCATAATTAATCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTa                                                                                                                                                                                                                                                                                      >  3:339032/1‑300 (MQ=255)
               acaGTATTGATAGTTTTCATAATTAATCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTgg                                                                                                                                                                                                                                                                            >  2:174799/1‑297 (MQ=255)
                    aTTGATAGTTTTCATAATTAATCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTc                                                                                                                                                                                                                                                                   <  4:296669/301‑1 (MQ=255)
                    aTTGATAGTTTTCATAATTAATCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTc                                                                                                                                                                                                                                                                   <  3:33594/301‑1 (MQ=255)
                           gTTTTCATAATTAATCTCTCGCAGGTGATCTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCTTGTTCTTTTACTATAGACAAAATCGAAGAGACTTGACCGCCTTGTTCAAATAAATTGATTGATATCTAAACTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACTGCGTGGCAGTCATAGCAGCGTTAAAGGTCCGGAATGACACGCCAACTACAAATTATACTCGCTATGCTTTTCGCCATTTGTCGTTACCGTACTTGTCTGTCGTTCGACCctt                                                                                                                                                                                                                                                            >  3:58421/1‑300 (MQ=255)
                                 aTAATTAATCGCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTTGTCGATGTGATGAGTATCACTTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGAt                                                                                                                                                                                                                                                      <  2:171312/301‑1 (MQ=255)
                                   aattaatCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATta                                                                                                                                                                                                                                                    <  1:186083/301‑1 (MQ=255)
                                   aattaatCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGAGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATta                                                                                                                                                                                                                                                    <  4:451749/301‑1 (MQ=255)
                                          ctctcGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATcc                                                                                                                                                                                                                                               >  4:347141/1‑299 (MQ=255)
                                          ctctcGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATcc                                                                                                                                                                                                                                               <  3:347135/299‑1 (MQ=255)
                                                       tttttGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGttt                                                                                                                                                                                                                                >  3:181062/1‑301 (MQ=255)
                                                        ttttGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGtttt                                                                                                                                                                                                                               >  4:105879/1‑301 (MQ=255)
                                                           tGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTg                                                                                                                                                                                                                             <  3:276597/300‑1 (MQ=255)
                                                               taAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAg                                                                                                                                                                                                                                                                                         <  3:201235/236‑1 (MQ=255)
                                                               taAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAg                                                                                                                                                                                                                                                                                         >  4:201238/1‑236 (MQ=255)
                                                                aaGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGCTTATATCCTGTATAAGCTGGTTTTTGTAAAtt                                                                                                                                                                                                                       >  1:64733/1‑300 (MQ=255)
                                                                 aGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAACTACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATACGCTGGTTTTTGTAAATCt                                                                                                                                                                                                                      >  3:189082/1‑301 (MQ=255)
                                                                 aGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCt                                                                                                                                                                                                                      <  4:181063/301‑1 (MQ=255)
                                                                  gTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTg                                                                                                                                                                                                                                                                                            >  3:217125/1‑230 (MQ=255)
                                                                  gTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCtt                                                                                                                                                                                                                     <  3:212394/301‑1 (MQ=255)
                                                                      cGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGc                                                                                                                                                                                                                                                         >  3:9898/1‑261 (MQ=255)
                                                                      cGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGc                                                                                                                                                                                                                                                         <  4:9898/261‑1 (MQ=255)
                                                                                  gATGTGATGAGTATCCCGTTTTTGTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATg                                                                                                                                                                                                     <  3:310179/301‑1 (MQ=255)
                                                                                  gATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAAtt                                                                                                                                                                                                     >  4:455282/1‑300 (MQ=255)
                                                                                  gATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTATTGGATTATATCCTTTCTAAGCTGGtagatcggaagagcgtctggagggaaaggg                                                                                                                                                                                                     >  3:273013/1‑272 (MQ=255)
                                                                                  gATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGACTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGt                                                                                                                                                                                                                                  <  4:273017/272‑1 (MQ=255)
                                                                                              aTCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATGAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTcc                                                                                                                                                                                                                <  4:220576/278‑1 (MQ=255)
                                                                                              aTCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCTCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTccagatcggaagagcgtcgattcgg                                                                                                                                                                                         >  3:220573/1‑278 (MQ=255)
                                                                                                 aCGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCAc                                                                                                                                                                                      <  3:105878/301‑1 (MQ=255)
                                                                                                             gagagaTAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGc                                                                                                                                                                          <  4:339038/301‑1 (MQ=255)
                                                                                                               gagaTAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACCGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCtt                                                                                                                                                                        >  4:274981/1‑301 (MQ=255)
                                                                                                                                                  aataaatTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCGGCGGGATGGGGGCACCGAActct                                                                                                                                     >  3:452753/1‑301 (MQ=255)
                                                                                                                                                      aatTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACCGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTcac                                                                                                                                 <  3:274977/301‑1 (MQ=255)
                                                                                                                                                                                 aaCCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCg                                                                                                                                                         >  3:393566/1‑250 (MQ=255)
                                                                                                                                                                                 aaCCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCg                                                                                                                                                         <  4:393574/250‑1 (MQ=255)
                                                                                                                                                                                                         aaCTTGTTTTATACCGCGTGGCAGTGACAGCAGCGTTAAAGGCACGGAATGACACGCAGAATACAAATTTTACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCg                                                                              <  3:455272/301‑1 (MQ=255)
                                                                                                                                                                                                            ttGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCAt                                                                           >  4:203279/1‑301 (MQ=255)
                                                                                                                                                                                                                 ttATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTc                                                                       <  4:189083/300‑1 (MQ=255)
                                                                                                                                                                                                                    taCCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATATTCAAtt                                                                   >  3:126925/1‑301 (MQ=255)
                                                                                                                                                                                                                     aCCGCGTGGCAGTCACAGCAGCGTGAAAGGAACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTg                                                                  <  2:128087/301‑1 (MQ=255)
                                                                                                                                                                                                                      ccGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTTTGCGATTATATCCTGtat                                                                                                                                                                                                                                          >  4:433563/1‑132 (MQ=255)
                                                                                                                                                                                                                      ccGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGtat                                                                                                                                                                                                                                          <  3:433553/132‑1 (MQ=255)
                                                                                                                                                                                                                                gTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGGACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTCCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGTCGCCGGGAATTTCACCCCGCGCCCGGGCAATCGCATCTTCAATTGTACTGTTgtt                                                        >  2:69375/1‑298 (MQ=255)
                                                                                                                                                                                                                                     aggagcGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTc                                                  <  4:452763/298‑1 (MQ=255)
                                                                                                                                                                                                                                          cGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCCTCAATTGTACTGTTGATCTGTTCGTTgg                                             >  3:6546/1‑300 (MQ=255)
                                                                                                                                                                                                                                                 ggCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATgg                                                                                                                                                   <  4:272144/192‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATgg                                                                                                                                                   >  3:272140/1‑192 (MQ=255)
                                                                                                                                                                                                                                                 ggCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCAcc                                                                                               <  3:304972/244‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCAcc                                                                                               >  4:304977/1‑244 (MQ=255)
                                                                                                                                                                                                                                                   cACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCAt                                    >  3:261490/1‑301 (MQ=255)
                                                                                                                                                                                                                                                     cGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCg                                  <  1:175029/301‑1 (MQ=255)
                                                                                                                                                                                                                                                             cacGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGccc                           >  4:414355/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 aaaagTACAAATTATACGCGCTATGCATCTGGCCATTTTTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAAACTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGGGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGTGTCATCGTTAGCCCCACCg                      <  2:64654/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aCAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGcc                <  3:414345/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGcc                            <  4:278121/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGcc                            >  3:278117/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                          gcgtgtgctcttccgatcgctcGCCATTTGTAGTTAACGTAACTGGCTGTCTTGCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTGTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGTGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCata             <  4:353928/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     cTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTcccc  >  3:132640/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             ctcGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACATCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGTGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGGCATCGTTAGCCCAACCGGCTGCCata             >  3:353921/1‑282 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                             
GCCATAGCAGCAACAACAGTATTGATAGTTTTCATAATTAATCTCTCGCAGGTGATTTTTGTATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCC  >  NC_000913/837891‑838468

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: