Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,088,199 C→A 60.3% intergenic (+341/‑83) metK → / → galP S‑adenosylmethionine synthetase/D‑galactose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,088,1990CA60.3% 38.3 / 81.6 58intergenic (+341/‑83)metK/galPS‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base C (15/8);  new base A (21/14);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCG  >  NC_000913/3087926‑3088483
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                            
cgcATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCt                                                                                                                                                                                                                                                                                >  2:510079/1‑288 (MQ=255)
                                     ttcttcCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTAACTTAATTc                                                                                                                                                                                                                                           >  2:446812/1‑288 (MQ=255)
                                         tcCTGATCACGTTTCACCGCAGATTATCATCACAGCTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATAGTTTCATATTATTATTCCGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGACGTACTCACCTAGCTTAATTCACAAt                                                                                                                                                                                                                                      <  1:253066/289‑1 (MQ=255)
                                             gATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATaaaa                                                                                                                                                                                                                                  >  1:578866/1‑289 (MQ=255)
                                             gATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTAcg                                                                                                                                                                                                                                                                 >  2:307077/1‑258 (MQ=255)
                                             gATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTAcg                                                                                                                                                                                                                                                                 <  1:307077/258‑1 (MQ=255)
                                                  cGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGAGTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATaa                                                                                                                                                                                                                             <  2:675189/289‑1 (MQ=255)
                                                                    catcaCAACTGAAACCGATTAAACCCACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGTGATTACACTGATGTGATTTGCTTCCCATCTTTTTACGTCGTACTTACCTATCTTAATTCACAATAAAAAATAACCATTTTGGTTGGcatca                                                                                                                                                                                                           >  2:378512/1‑289 (MQ=255)
                                                                                 aCCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTaa                                                                                                                                                                                              >  2:611675/1‑289 (MQ=255)
                                                                                              accaCAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAAtt                                                                                                                                                                                                                                            <  2:201101/230‑1 (MQ=255)
                                                                                              accaCAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAAtt                                                                                                                                                                                                                                            >  1:201101/1‑230 (MQ=255)
                                                                                                    acaGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAAc                                                                                                                                                                           >  1:59481/1‑289 (MQ=255)
                                                                                                        aCAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAg                                                                                                                                                                        >  1:5268/1‑288 (MQ=255)
                                                                                                             gATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATaa                                                                                                                                                                                                                             >  2:730667/1‑230 (MQ=255)
                                                                                                             gATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATaa                                                                                                                                                                                                                             <  1:730667/230‑1 (MQ=255)
                                                                                                                        ttttCATATTATTATTCGGTTTTCACAGCTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTCCTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGt                                                                                                                                                       >  1:257325/1‑289 (MQ=255)
                                                                                                                                            tttCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTgg                                                                                                                                   <  2:5268/289‑1 (MQ=255)
                                                                                                                                              tCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCg                                                                                                                                 >  1:189893/1‑289 (MQ=255)
                                                                                                                                                     tgttACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTa                                                                                                                           <  2:189893/288‑1 (MQ=255)
                                                                                                                                                                        aaaGTCTTAATTGCAGATAACATCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGTGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat                                                                                                          <  2:500080/286‑1 (MQ=255)
                                                                                                                                                                                     cAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat                                                                                                          <  2:783768/273‑1 (MQ=255)
                                                                                                                                                                                     cAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat                                                                                                          >  1:783768/1‑273 (MQ=255)
                                                                                                                                                                                      aGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGct                                                                                                                                       >  2:771267/1‑243 (MQ=255)
                                                                                                                                                                                      aGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGct                                                                                                                                       <  1:771267/243‑1 (MQ=255)
                                                                                                                                                                                          aaCAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGgcgc                                                                                        >  2:407233/1‑286 (MQ=255)
                                                                                                                                                                                               cGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCa                                                                                       <  1:257585/282‑1 (MQ=255)
                                                                                                                                                                                               cGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCa                                                                                       >  2:257585/1‑282 (MQ=255)
                                                                                                                                                                                                         atgatgATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGa                                                                      >  2:225450/1‑289 (MQ=255)
                                                                                                                                                                                                         atgatgATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGa                                                                      >  2:753016/1‑289 (MQ=255)
                                                                                                                                                                                                          tgatgaTATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGCGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGCTTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGTCGTTTATTGCAgtt                                                                     >  2:92695/1‑287 (MQ=255)
                                                                                                                                                                                                             tgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTTCTTCACATCTTTTTACGTCGTACGCACCTATCTTAATTCACAATAAAAAAGAACCAGATGGGAGGGCATCATGCCTGACGCTAAAAAACATGGTCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATTGGTGTAATTGCTGGCGCACTGCCGTGTATTGCAGATGat                                                                  <  1:353422/289‑2 (MQ=255)
                                                                                                                                                                                                             tgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATg                                                                    <  1:510079/287‑1 (MQ=255)
                                                                                                                                                                                                             tgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGaa                                                                  <  1:225450/289‑1 (MQ=255)
                                                                                                                                                                                                                tataACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGAtt                                                                                                                            >  1:283885/1‑228 (MQ=255)
                                                                                                                                                                                                                tataACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTGTTCGTCTGCTTCCTTGCCGCTCTGGCGGGAtt                                                                                                                            <  2:283885/228‑1 (MQ=255)
                                                                                                                                                                                                                          ttattTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACt                                                      >  2:445715/1‑288 (MQ=255)
                                                                                                                                                                                                                          ttattTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTAAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCCCTCCCGTTTCTTGCAGATGAATTCCAGATTACtt                                                     >  1:337319/1‑289 (MQ=255)
                                                                                                                                                                                                                                tCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGcaca                                               <  1:611675/289‑1 (MQ=255)
                                                                                                                                                                                                                                       cTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGa                                        >  2:753578/1‑289 (MQ=255)
                                                                                                                                                                                                                                                tAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGt                               <  2:59481/289‑1 (MQ=255)
                                                                                                                                                                                                                                                tAACTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGt                               <  1:753578/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  aCTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCa                                                                                       >  2:755140/1‑231 (MQ=255)
                                                                                                                                                                                                                                                  aCTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCa                                                                                       <  1:755140/231‑1 (MQ=255)
                                                                                                                                                                                                                                                  aCTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTActct                                                                                                                       >  1:265622/1‑198 (MQ=255)
                                                                                                                                                                                                                                                  aCTAAGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTActct                                                                                                                       <  2:265622/198‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGATAAATGTTAGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTTTTGCAGATGACTTCAACATTACTTCGCACACGCAAGAATGGGTCGTAAGCt                          >  1:502086/1‑288 (MQ=255)
                                                                                                                                                                                                                                                             tGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCatgatg                  >  2:255783/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                 aGTGTAAGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCg              >  2:57343/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  gtgtAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat                                                                                                          <  2:14093/196‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gtgtAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGatat                                                                                                          >  1:14093/1‑196 (MQ=255)
                                                                                                                                                                                                                                                                  gtgtAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGca                                                 >  1:507917/1‑253 (MQ=255)
                                                                                                                                                                                                                                                                  gtgtAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGca                                                 <  2:507917/253‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCCGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCg         >  2:805317/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                       aGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCCTGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGCCGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGGCTACTATTTGGCCTGGATATCGGTGTAATTGCTGTCGGCCTGCCGTTTATCGCAGATGAATTCCCGATTCTTTCGCACACGCAAGAATGTGTCGTATCCTCCATGATGTTTGGTGCgg        >  1:66954/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                       aGCGATTACAATGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCg         >  2:525735/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTTCGGCAGTcg  >  2:415249/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                             tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTcg  >  2:625588/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                             tACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTcg  >  1:274570/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                            
CGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCTGCTTCCTTGCCGCTCTGGCGGGATTACTCTTTGGCCTGGATATCGGTGTAATTGCTGGCGCACTGCCGTTTATTGCAGATGAATTCCAGATTACTTCGCACACGCAAGAATGGGTCGTAAGCTCCATGATGTTCGGTGCGGCAGTCG  >  NC_000913/3087926‑3088483

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