Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 11.2% intergenic (+31/‑186) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2610GT11.2% 285.4 / 23.2 90intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (37/43);  new base T (4/6);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-03

CGCTCTAACGACCTGACTGTAAACTGCCTGACCGGTAAGAAACTGACCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACT  >  NC_000913/4059991‑4060490
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
cGCTCTAACGACCTGACTGTAAACTGCCTGACCGGTAAGAAACTGACCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTCTCCGTATTCGTAAACGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAACGACGATTAATTTCGTTTTTTAATCTTACAAAACCTGCCAGCGag                                                                                                                                                                                                                                    >  1:111955/1‑274 (MQ=255)
                         gCCTGACCGGTAAGAAACTGACCAAACTGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTAGCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                                                                                                                                                           <  2:284418/274‑1 (MQ=255)
                              acgggTAAAAAACTGACCAAAAGGGCTTGTTCCTGTAAGGTAGAAGCCGTTGTTTTGGTTCCGCCTATCCGCATGAAGCTGGAACAAGCTCTGGAGTTCATCGATGAGTACGAACTTGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcgc                                                                                                                                                                                                     <  4:46307/272‑1 (MQ=255)
                                             aCCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTTACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCAGCTCCAACCGGGCACCGGAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcacccacatgcacgta                                                                                                                                                                                      >  4:245677/1‑262 (MQ=255)
                                               cAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                                                                                                    <  2:321250/275‑1 (MQ=255)
                                                 aCATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACGCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGca                                                                                                                                                                                                    <  2:192596/257‑1 (MQ=255)
                                                      ggTGCTTCCGTTACTGCCGATCCCTTGCTTCTGGTCACGCCCATCCGCCTGACTCTGGAACAAGCTCTGAAGTTCAACGATGATGACGAACGGTTAGAATTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAACACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                                                                                                             <  4:420502/274‑1 (MQ=255)
                                                       gTGCTTCCGGTACTGCCGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                                                                                                                                                            <  2:19903/275‑1 (MQ=255)
                                                       gTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTCTCTCCCTta                                                                                                                                                                             >  2:1744/1‑274 (MQ=255)
                                                                gTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                                                                                                                                                     >  2:158890/1‑273 (MQ=255)
                                                                  acgtaCGAAGCCGTCGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGATCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCAAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                                                                                                                                                 <  4:158888/271‑1 (MQ=255)
                                                                   cTGACAAAGCCGTGTTTCTGGTACCGCCGATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                                                                                                                                                 <  2:56691/274‑1 (MQ=255)
                                                                                       gTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGACAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATGAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCTATGGCAGGTTTTTTTTCGCCTGTCTTGAGCCAACACCAATGCGAg                                                                                                                                                                                        >  4:179420/1‑231 (MQ=255)
                                                                                        ttCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                                                                                                                                           >  2:168926/1‑275 (MQ=255)
                                                                                          ccGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTACACGTCACCTGACGGAAAACGAACGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGTTTACGTCAATAGAGAGCTTTGTc                                                                                                                                           >  2:332978/1‑273 (MQ=255)
                                                                                           cGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTCGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                                                                                                    <  2:179422/231‑1 (MQ=255)
                                                                                            gCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTg                                                                                                                                                                                                               >  4:159300/1‑203 (MQ=255)
                                                                                                aTCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                                                                                                                                                           <  2:159302/203‑1 (MQ=255)
                                                                                                 tCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                                                                                                                                           >  4:45970/1‑250 (MQ=255)
                                                                                                     cATGACTCTGTAACAAGGTCTGGAGTGCATGGTGGCGGCCGAACCGGTTGAATGTACTCCGACCCCGATCCATACTAGTAACCGTAACCTGACTGAAAACTTTAGTCGCCGCGCAAACCGCGCACCGAAAGACGATTAATTTGGTGTTTTAGTCATAAAAAACCTGGCAGCGATGGAAGGTTTTTTTTTGCCGGTCGTGAGCGCACACCAATGCGACTATCTCCCTTATAAGTCGGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                                                                                                              <  4:394912/275‑1 (MQ=255)
                                                                                                     cATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                                                                                                       <  2:45970/250‑1 (MQ=255)
                                                                                                               gAACAATGTCTTGTGTTCATGGATGCTACNAAACTGGTAAAAAAGACTCCGCCCCNTGACCCTATTCCTAAACGGCACCTGACGGAAAACGGTACTCGCCGCGCCAACCGCGCACCGAAAGACGAGTAATTTCTTTNTTTAGTCGTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACg                                                                                                                    <  3:90722/275‑1 (MQ=255)
                                                                                                                    aGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                                                                                                                                                     >  4:302186/1‑173 (MQ=255)
                                                                                                                     gCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTg                                                                                                                                                                                                               >  4:266386/1‑178 (MQ=255)
                                                                                                                        ctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                                                                                                                                                                 <  2:302190/173‑1 (MQ=255)
                                                                                                                         tGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                                                                                                                                                           <  2:266390/178‑1 (MQ=255)
                                                                                                                                 aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                                                                                                 <  5:177219/149‑1 (MQ=255)
                                                                                                                                  tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCa                                                                                                                                                                                                                                <  8:216340/149‑1 (MQ=255)
                                                                                                                                  tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                                                                                                                                                >  7:226458/1‑149 (MQ=255)
                                                                                                                                    gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                                                                                                              >  6:394842/1‑149 (MQ=255)
                                                                                                                                      tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtt                                                                                                                                                                                                                           <  6:31032/149‑1 (MQ=255)
                                                                                                                                       gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttt                                                                                                                                                                                                                          <  6:203761/149‑1 (MQ=255)
                                                                                                                                        atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                                                                                                                                                     >  2:423827/1‑138 (MQ=255)
                                                                                                                                            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                                                                                                 <  4:423822/138‑1 (MQ=255)
                                                                                                                                              aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                                                                                                                                                   >  7:53438/1‑149 (MQ=255)
                                                                                                                                                   gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                                                                                                                                                                                                              >  6:169862/1‑149 (MQ=255)
                                                                                                                                                   gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGa                                                                                >  2:158644/1‑275 (MQ=255)
                                                                                                                                                     aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                                                                                                                                                   <  7:205112/142‑1 (MQ=255)
                                                                                                                                                     aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                                                                                                                                                   >  8:205112/1‑142 (MQ=255)
                                                                                                                                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatga                                                                                                                                                                                                                   >  5:327396/1‑136 (MQ=255)
                                                                                                                                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatga                                                                                                                                                                                                                   <  6:327396/140‑5 (MQ=255)
                                                                                                                                                       aaGTGACTCCGACCCTTAACCGAATTCATAACCGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGAGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGc                                                                             <  2:253794/274‑1 (MQ=255)
                                                                                                                                                            aCTCCGAACTCTACCCCATTTCGAAAACGAGACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGACCCGAAAGACGATTAATTTCGTTTTTTAGTAATAAAAAACCGGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATACGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGATTTCTTTTTTAGCGTAAGCGTAAGATGCTTCATc                                                                       <  4:168924/275‑1 (MQ=255)
                                                                                                                                                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt                                                                                               >  2:431245/1‑250 (MQ=255)
                                                                                                                                                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATc                                                                       <  2:245681/274‑1 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                                                                                                                                                    <  6:81876/132‑1 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                                                                                                   >  8:119424/1‑149 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAATCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                                                                                                                                                    >  5:81876/1‑132 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTCTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                                                                                                                                                          >  6:340909/1‑142 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCACGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                                                                                                                                                          <  5:340909/142‑1 (MQ=255)
                                                                                                                                                                  aCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                                                                                                        >  2:90254/1‑156 (MQ=255)
                                                                                                                                                                      ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                                                                                                    <  4:90253/156‑1 (MQ=255)
                                                                                                                                                                           cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAGACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                                                                                                         <  8:38410/146‑1 (MQ=255)
                                                                                                                                                                           cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                                                                                                         >  7:38410/1‑146 (MQ=255)
                                                                                                                                                                            gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                                                                                                                                                     >  6:78927/1‑149 (MQ=255)
                                                                                                                                                                              aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAGTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTGTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTCTAAGTCTGTGGTTACGTCACTCGAGAGATTTGTCACAATTATCTGCAAAGTAATACGCCGTTAATTTCTTTCTTTTTTTGTGTAAGCGTAAGATCACTCATCTGGTTTAAACCAAACGGa                                                     >  4:317016/1‑275 (MQ=255)
                                                                                                                                                                               ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGCAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                                                                                                        >  2:434470/1‑207 (MQ=255)
                                                                                                                                                                                 cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                                                                                                  >  4:286856/1‑195 (MQ=255)
                                                                                                                                                                                    aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACa                                              <  2:388980/275‑1 (MQ=255)
                                                                                                                                                                                     aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                                                                                                                                                            <  7:119424/149‑1 (MQ=255)
                                                                                                                                                                                     aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                                                                                                              <  2:286860/195‑1 (MQ=255)
                                                                                                                                                                                         tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                                                                                                        <  8:3445/149‑1 (MQ=255)
                                                                                                                                                                                            ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                                                                                                                                                      <  6:254750/149‑1 (MQ=255)
                                                                                                                                                                                            ccTGACGGAAAACGATCGTCCCCGCGCCAACCGCGCACCGAAAGACGCTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGGGACCGCACACCAATGCGAGTATCCCCCTCATAAGTCCGTGGTCACGCCCATAGAGCGCCTTGCCCCAATTCCTCGCAACGCCATACGCCGCTACTTCCTTTCTTTTTTGGCGTAACCGTAAGATCCTCCATCGGCCTTAAACCAAAAGCATTaacacatgcca                                        >  4:124663/1‑265 (MQ=255)
                                                                                                                                                                                               gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                                                                                                  <  5:169862/149‑1 (MQ=255)
                                                                                                                                                                                                aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                                                                                                                                                 <  6:367123/149‑1 (MQ=255)
                                                                                                                                                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttt                                                                                                                                                                                                                          >  5:160549/1‑86 (MQ=255)
                                                                                                                                                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttt                                                                                                                                                                                                                          <  6:160549/86‑1 (MQ=255)
                                                                                                                                                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                                                                                                  >  3:145475/1‑174 (MQ=255)
                                                                                                                                                                                                         gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGCCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCAc                          <  2:317020/275‑1 (MQ=255)
                                                                                                                                                                                                          aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                                                                                                              <  1:145475/174‑1 (MQ=255)
                                                                                                                                                                                                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                                                                                                     <  8:53438/149‑1 (MQ=255)
                                                                                                                                                                                                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                                                                                                     <  8:226458/149‑1 (MQ=255)
                                                                                                                                                                                                                  cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCCATCCCACGTAGaa                   >  2:508533/1‑273 (MQ=255)
                                                                                                                                                                                                                        aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac                                                                                                                                         >  8:150653/1‑149 (MQ=255)
                                                                                                                                                                                                                                     aaGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaaa                                                                                                                            <  5:78927/149‑1 (MQ=255)
                                                                                                                                                                                                                                        aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTc                                                                                                                         <  7:32418/149‑1 (MQ=255)
                                                                                                                                                                                                                                          gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTAAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCAt                                                                                                                       <  7:368424/149‑1 (MQ=255)
                                                                                                                                                                                                                                            ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAAGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACt  >  2:76206/1‑264 (MQ=255)
                                                                                                                                                                                                                                                tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCg                                                                                                                 <  8:183130/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                                                                                                       <  5:396201/106‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                                                                                                       >  6:396201/1‑106 (MQ=255)
                                                                                                                                                                                                                                                        tttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGc                                                                                                         <  7:150653/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttttt                                                                                                >  5:208931/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAACTCATACGCCGTTAATTGCTTTCttttt                                                                                                <  6:208931/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                                                                                                       >  5:122843/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                                                                                                       <  6:122843/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc                                                                                            >  5:158623/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       ccTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgt                                                                                          >  7:429396/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         tGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa                                                                                        >  6:433576/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
CGCTCTAACGACCTGACTGTAAACTGCCTGACCGGTAAGAAACTGACCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACT  >  NC_000913/4059991‑4060490

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