Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,300 A→T 6.0% intergenic (‑219/+20) waaH ← / ← tdh LPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,790,3000AT 6.0% 2121.2 / 83.6 590intergenic (‑219/+20)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (310/243);  new base T (21/15);  total (331/259)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-03

Only 100 of 598 total aligned reads displayed.

TAACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAAT  >  NC_000913/3790040‑3790571
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                
taACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                                                                                                                                                                                                    <  4:179187/275‑1 (MQ=255)
                                                                  aTCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGGTATTATTTTTACAGGTGAGGGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                                                                                                                                                                                                  >  2:126652/1‑275 (MQ=255)
                                                                                    tttttttttAAAAAATTATTAAGCCTGTGNNNTATTAATTTTACAGGGTGAAGAGTAACTTCAGGGTAAGGAGATAATGAGAACAACTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATagag                                                                                                                                                                                <  3:230070/268‑1 (MQ=255)
                                                                                                              gtgAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGAGTCAAAGCCTTTACGGAAATCATCGATAGAGAAACGATTGGTAATGATCGGCGAAAg                                                                                                                                                      >  4:225554/1‑275 (MQ=255)
                                                                                                                 agaTATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGaa                                                                                                                                                                                                                                                                                 >  7:33338/1‑149 (MQ=255)
                                                                                                                  gaTATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAAt                                                                                                                                                                                                                                                                                <  7:599213/149‑1 (MQ=255)
                                                                                                                                        agagTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                                                                                                                                                                                           >  7:470683/1‑149 (MQ=255)
                                                                                                                                         gagTCACCTCAGGGTAGAGAGATAGTGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACATTGACAGGACTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                                                                                                                                                                                           >  5:216730/1‑148 (MQ=255)
                                                                                                                                            tCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                                                                                                                                                                                                        >  6:140631/1‑148 (MQ=255)
                                                                                                                                            tCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                                                                                                                                                                                                        >  5:139633/1‑148 (MQ=255)
                                                                                                                                              aCTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGcc                                                                                                                                                                                                                                                                         >  7:289943/1‑129 (MQ=255)
                                                                                                                                                ttCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGcc                                                                                                                                                                                                                                                                         <  6:315237/127‑1 (MQ=255)
                                                                                                                                                 tCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCt                                                                                                                                                                                                                                                   >  6:442628/1‑148 (MQ=255)
                                                                                                                                                    gggTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                                                                                                                                                                                               >  8:52984/1‑149 (MQ=255)
                                                                                                                                                        agagagATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                                                                                                                                                                                           <  7:243366/149‑1 (MQ=255)
                                                                                                                                                             gaTAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACttt                                                                                                                                                                                                                                       <  5:366794/148‑1 (MQ=255)
                                                                                                                                                               tAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                                                                                                                                                                                                    <  8:405401/149‑1 (MQ=255)
                                                                                                                                                                aaTGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                                                                                                                                                                                                   >  8:733257/1‑149 (MQ=255)
                                                                                                                                                                aaTGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                                                                                                                                                                                                   <  8:662565/149‑1 (MQ=255)
                                                                                                                                                                 aTGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACa                                                                                                    >  1:18612/1‑275 (MQ=255)
                                                                                                                                                                          gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                                                                                                                                                                                                 <  2:503007/157‑1 (MQ=255)
                                                                                                                                                                                taAGGGATGACCACAAATGCCAAAACTTAATTAAACTGCCAGGATTTTACATATTTAAAAGATTAACCCCAGATAAAATGCCTGCAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCATAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTaaa                                                                                    <  4:319453/275‑1 (MQ=255)
                                                                                                                                                                                 aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                                                                                                                                                                                                 >  5:532553/1‑102 (MQ=255)
                                                                                                                                                                                    agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                     <  2:411631/175‑1 (MQ=255)
                                                                                                                                                                                     gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                                                                                                                                                                                               >  5:29563/1‑148 (MQ=255)
                                                                                                                                                                                     gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                                                                                                                                                                                               >  6:383907/1‑148 (MQ=255)
                                                                                                                                                                                        cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTcc                                                                                                                                                                                                                                     >  4:547793/1‑123 (MQ=255)
                                                                                                                                                                                          gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                                                                                                                                                                         <  8:270360/149‑1 (MQ=255)
                                                                                                                                                                                            acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGa                                                                                                                                                                                                                                 <  2:547795/123‑1 (MQ=255)
                                                                                                                                                                                              aaTTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGccc                                                                                                                                                                                                                                                                        >  2:356982/1‑82 (MQ=255)
                                                                                                                                                                                                     aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                                                                                                                                                                                                        <  6:610123/139‑1 (MQ=255)
                                                                                                                                                                                                     aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGGGGCCAGACTGAATCAGCACAGCCATCTTGGACCAGGTTTCAAACATCACACGCCGGTAAAAACCTTTACTGAAAAAGCCtt                                                               >  4:392085/1‑275 (MQ=255)
                                                                                                                                                                                                         ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                                                                                                                                                                                          >  8:100814/1‑149 (MQ=255)
                                                                                                                                                                                                         ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGAGCTGGAATACCAGCCCTTGTTCGTGTTAATTCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                                                                                                                                                                                           >  6:584919/1‑148 (MQ=255)
                                                                                                                                                                                                            tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat                                                                                                                                                                                        <  5:588922/148‑1 (MQ=255)
                                                                                                                                                                                                              ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                     >  7:72224/1‑149 (MQ=255)
                                                                                                                                                                                                              ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                     >  7:855716/1‑149 (MQ=255)
                                                                                                                                                                                                              ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACg                                                                                                                                                                           <  1:226471/159‑1 (MQ=255)
                                                                                                                                                                                                                   cTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga                                                                                                                                                                                 >  5:793878/1‑148 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                                                                                                                                                                                               >  4:436793/1‑85 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                                                                                                                                                                                               >  4:612200/1‑85 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCGTAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAAGATGGGTAATGATCGGCGAAAGATCGAGGCAAGACTGAATCAGCGCCGCCATCTTGTACCAGGTATCGAACATCCCACGACCGTAAATACCTTTAATGAACAAACCTTTAAAGACCACCTTgg                                                >  2:324943/1‑275 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca                                                                                                                                                                                         >  3:181645/1‑139 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATagag                                                                                                                                                                                >  6:891888/1‑148 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACGAGGTTTCAAACATCTCACGACCGTCAATACCTTTAATGAAAAAGCCTTTAAAGATCACTTTgg                                                >  2:193397/1‑275 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTgt                                                >  2:63366/1‑274 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTgg                                                >  2:663305/1‑275 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGATCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCg                                                                                                                                                  >  1:204220/1‑178 (MQ=255)
                                                                                                                                                                                                                    tGACAGGATTTTACATATTTAACAGATTATACTCAGATAAAAGGGCTGGAATACCAGCCCtt                                                                                                                                                                                                                                                                      >  1:48595/1‑62 (MQ=255)
                                                                                                                                                                                                                      aCAGGATTTTACATATTTAACAGATTATCCTCAGATAGAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAACGGCAAAGTCCTTCTTGAAATCATCGATAGAGAAACGATTGGTAATGATCTGCCAACGATCGCGGCCAGTCTTAATCAGCGCCGCCTCCTGGTACTAGGCTTAACACAACCCACTACCGGAACTACACTTAATCAAACAGCCTTTACAGAACACtttcgcc                                              >  4:355955/1‑271 (MQ=255)
                                                                                                                                                                                                                      aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaa                                                                                                                                                                              >  6:825074/1‑148 (MQ=255)
                                                                                                                                                                                                                      aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGCCTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACAtctc                                                                                                >  1:134092/1‑226 (MQ=255)
                                                                                                                                                                                                                      aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTAAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACAGTAAATACCTTTAAAGCACAACCCTTTAAAGATAAATTTGGta                                              >  1:107211/1‑274 (MQ=255)
                                                                                                                                                                                                                       cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAc                                                                                                                                                                            <  7:77502/149‑1 (MQ=255)
                                                                                                                                                                                                                        aGGATTTTACATATTTAACAGTTAAGCCTCAGAAAAAAGNGCTGGAATACCTGTCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAATCCCTTCTCGAATTCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGc                                                                                                                                              <  3:204219/178‑1 (MQ=255)
                                                                                                                                                                                                                         ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGa                                                                                                                                                                          <  8:488264/149‑1 (MQ=255)
                                                                                                                                                                                                                         ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCCTCTGGAAATCATCGATAGAGAAACg                                                                                                                                                                           >  5:674977/1‑148 (MQ=255)
                                                                                                                                                                                                                                 cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAGGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAAt                                                                                                                                                                  >  7:368505/1‑149 (MQ=255)
                                                                                                                                                                                                                                   tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                                                                                                                                                                         <  7:167656/108‑1 (MQ=255)
                                                                                                                                                                                                                                   tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGa                                                                                                                                                                >  7:400640/1‑149 (MQ=255)
                                                                                                                                                                                                                                    atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACt                                                                                                                                                                                                                                         >  5:431833/1‑75 (MQ=255)
                                                                                                                                                                                                                                         aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGc                                                                                                                                                           >  5:825460/1‑148 (MQ=255)
                                                                                                                                                                                                                                         aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCg                                                                                                                                                          >  8:378070/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCg                                                                                                                                                          >  8:167340/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCCACGGCCCGTAAATACCTTTAATGACCAAGCTTTTAAAGATAACTTTGGTCAAGTCGATGGAATTATCAg                           >  4:484214/1‑275 (MQ=255)
                                                                                                                                                                                                                                         aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAACCATCGATAGAGAAACGATGGGTAATGATCGGCGaa                                                                                                                                                        <  4:782784/151‑1 (MQ=255)
                                                                                                                                                                                                                                          cAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGc                                                                                                                                                                                                                                                    <  2:793904/58‑1 (MQ=255)
                                                                                                                                                                                                                                            gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaa                                                                                                                                                        <  5:881947/148‑1 (MQ=255)
                                                                                                                                                                                                                                             aTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaaa                                                                                                                                                       >  6:523758/1‑148 (MQ=255)
                                                                                                                                                                                                                                             aTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAg                                                                                                                                                      >  8:600146/1‑149 (MQ=255)
                                                                                                                                                                                                                                              ttATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAgc                                                                                                                                >  2:758103/1‑170 (MQ=255)
                                                                                                                                                                                                                                               tATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCg                                                                                                                                                  >  2:463637/1‑151 (MQ=255)
                                                                                                                                                                                                                                                aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGAt                                                                                                                                                    <  6:674977/148‑1 (MQ=255)
                                                                                                                                                                                                                                                aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGAt                                                                                                                                                    >  6:793844/1‑148 (MQ=255)
                                                                                                                                                                                                                                                  ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCg                                                                                                                                                  <  6:250363/148‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGa                                                                                                                                                 <  7:52984/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                     <  2:259030/112‑1 (MQ=255)
                                                                                                                                                                                                                                                   cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAg                                                                                                                                                <  7:118670/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                                                                                                                                                                                                     >  8:784824/1‑79 (MQ=255)
                                                                                                                                                                                                                                                    tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGTAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAgg                                                                                                                                               >  8:32848/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGa                                                                                                                                                                          <  7:670673/122‑1 (MQ=255)
                                                                                                                                                                                                                                                    tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAg                                                                                                                                                <  5:94391/148‑1 (MQ=255)
                                                                                                                                                                                                                                                     cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCgg                                                                                                                                                            <  5:164329/135‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGcc                                                                                                                                             >  8:225512/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGACGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCCCACGACAGCAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGAAATCACAGACGGCCGAATAc               >  4:298365/1‑274 (MQ=255)
                                                                                                                                                                                                                                                       gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaaa                                                                                                                                                       >  8:440966/1‑138 (MQ=255)
                                                                                                                                                                                                                                                       gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGcc                                                                                                                                             <  6:501530/148‑1 (MQ=255)
                                                                                                                                                                                                                                                        aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                                                                                                                                                                               >  3:48872/1‑113 (MQ=255)
                                                                                                                                                                                                                                                         tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAg           <  4:112402/275‑1 (MQ=255)
                                                                                                                                                                                                                                                           aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTaaa                                                                                    <  2:302072/200‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACt                                                                                                                                        >  6:196801/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                                                                                                                                                                                             <  4:731125/95‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaGGGATGGCAAACCAGCCCTTGTTCGGGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAg           <  4:663303/272‑1 (MQ=255)
                                                                                                                                                                                                                                                               ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACa                                                                     <  1:62216/211‑1 (MQ=255)
                                                                                                                                                                                                                                                                gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATc                                                                                                                                   >  7:80233/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCa                                            <  4:503648/235‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                                                                                                                                                                                               >  2:601397/1‑55 (MQ=255)
                                                                                                                                                                                                                                                                  tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGGTCGGCGAAAGATCGAGGCCAGACTGAATCAg                                                                                                                                 >  8:861983/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATgg                                                                                                                                                                       <  7:433561/110‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGAt                                                                                                                                                    >  5:499126/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                
TAACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAAT  >  NC_000913/3790040‑3790571

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