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

CGCTCTAACGACCTGACTGTAAACTGCCTGACCGGTAAGAAACTGACCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACT  >  NC_000913/4059991‑4060490
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
cGCTCTAACGACCTGACTGTAAACTGCCTGACCGGTAAGAAACTGACCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTCTCCGTATTCGTAAACGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAACGACGATTAATTTCGTTTTTTAATCTTACAAAACCTGCCAGCGag                                                                                                                                                                                                                                    >  1:111955/1‑274 (MQ=255)
                         gCCTGACCGGTAAGAAACTGACCAAACTGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTAGCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGt                                                                                                                                                                                                           <  2:284418/274‑1 (MQ=255)
                              acgggTAAAAAACTGACCAAAAGGGCTTGTTCCTGTAAGGTAGAAGCCGTTGTTTTGGTTCCGCCTATCCGCATGAAGCTGGAACAAGCTCTGGAGTTCATCGATGAGTACGAACTTGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgcgc                                                                                                                                                                                                     <  4:46307/272‑1 (MQ=255)
                                             aCCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTTACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCAGCTCCAACCGGGCACCGGAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGcacccacatgcacgta                                                                                                                                                                                      >  4:245677/1‑262 (MQ=255)
                                               cAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                                                                                                    <  2:321250/275‑1 (MQ=255)
                                                 aCATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACGCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGca                                                                                                                                                                                                    <  2:192596/257‑1 (MQ=255)
                                                      ggTGCTTCCGTTACTGCCGATCCCTTGCTTCTGGTCACGCCCATCCGCCTGACTCTGGAACAAGCTCTGAAGTTCAACGATGATGACGAACGGTTAGAATTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAACACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                                                                                                             <  4:420502/274‑1 (MQ=255)
                                                       gTGCTTCCGGTACTGCCGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                                                                                                                                                            <  2:19903/275‑1 (MQ=255)
                                                       gTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTCTCTCCCTta                                                                                                                                                                             >  2:1744/1‑274 (MQ=255)
                                                                gTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                                                                                                                                                     >  2:158890/1‑273 (MQ=255)
                                                                  acgtaCGAAGCCGTCGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGATCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCAAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                                                                                                                                                 <  4:158888/271‑1 (MQ=255)
                                                                   cTGACAAAGCCGTGTTTCTGGTACCGCCGATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                                                                                                                                                 <  2:56691/274‑1 (MQ=255)
                                                                                       gTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGACAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATGAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCTATGGCAGG‑TTTTTTTTCGCCTGTCTTGAGCCAACACCAATGCGAg                                                                                                                                                                                        >  4:179420/1‑231 (MQ=255)
                                                                                        ttCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                                                                                                                                           >  2:168926/1‑275 (MQ=255)
                                                                                          ccGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTACACGTCACCTGACGGAAAACGAACGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGTTTACGTCAATAGAGAGCTTTGTc                                                                                                                                           >  2:332978/1‑273 (MQ=255)
                                                                                           cGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTCGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                                                                                                    <  2:179422/231‑1 (MQ=255)
                                                                                            gCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTg                                                                                                                                                                                                               >  4:159300/1‑203 (MQ=255)
                                                                                                aTCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGt                                                                                                                                                                                                           <  2:159302/203‑1 (MQ=255)
                                                                                                 tCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                                                                                                                                           >  4:45970/1‑250 (MQ=255)
                                                                                                     cATGACTCTGTAACAAGGTCTGGAGTGCATGGTGGCGGCCGAACCGGTTGAATGTACTCCGACCCCGATCCATACTAGTAACCGTAACCTGACTGAAAACTTTAGTCGCCGCGCAAACCGCGCACCGAAAGACGATTAATTTGGTGTTTTAGTCATAAAAAACCTGGCAGCGATGGAAGG‑TTTTTTTTTGCCGGTCGTGAGCGCACACCAATGCGACTATCTCCCTTATAAGTCGGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                                                                                                              <  4:394912/275‑1 (MQ=255)
                                                                                                     cATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                                                                                                       <  2:45970/250‑1 (MQ=255)
                                                                                                               gAACAATGTCTTGTGTTCATGGATGCTACNAAACTGGTAAAAAAGACTCCGCCCCNTGACCCTATTCCTAAACGGCACCTGACGGAAAACGGTACTCGCCGCGCCAACCGCGCACCGAAAGACGAGTAATTTCTTTNTTTAGTCGTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACg                                                                                                                    <  3:90722/275‑1 (MQ=255)
                                                                                                                    aGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttttt                                                                                                                                                                                                                     >  4:302186/1‑173 (MQ=255)
                                                                                                                     gCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTg                                                                                                                                                                                                               >  4:266386/1‑178 (MQ=255)
                                                                                                                        ctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                                                                                                                                                                 <  2:302190/173‑1 (MQ=255)
                                                                                                                         tGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGt                                                                                                                                                                                                           <  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)
                                                                                                                                      tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tt                                                                                                                                                                                                                           <  6:31032/149‑1 (MQ=255)
                                                                                                                                       gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttt                                                                                                                                                                                                                          <  6:203761/149‑1 (MQ=255)
                                                                                                                                        atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                                                                                                                                                     >  2:423827/1‑138 (MQ=255)
                                                                                                                                            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                                                                                                 <  4:423822/138‑1 (MQ=255)
                                                                                                                                              aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                                                                                                                                                   >  7:53438/1‑149 (MQ=255)
                                                                                                                                                   gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGt                                                                                                                                                                                                              >  6:169862/1‑149 (MQ=255)
                                                                                                                                                   gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGa                                                                                >  2:158644/1‑275 (MQ=255)
                                                                                                                                                     aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                                                                                                                                                   <  7:205112/142‑1 (MQ=255)
                                                                                                                                                     aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                                                                                                                                                   >  8:205112/1‑142 (MQ=255)
                                                                                                                                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatga                                                                                                                                                                                                                   >  5:327396/1‑136 (MQ=255)
                                                                                                                                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatga                                                                                                                                                                                                                   <  6:327396/140‑5 (MQ=255)
                                                                                                                                                       aaGTGACTCCGACCCTTAACCGAATTCATAACCGGCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGAGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGc                                                                             <  2:253794/274‑1 (MQ=255)
                                                                                                                                                            aCTCCGAACTCTACCCCATTTCGAAAACGAGACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGACCCGAAAGACGATTAATTTCGTTTTTTAGTAATAAAAAACCGGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATACGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGATTTCTTTTTTAGCGTAAGCGTAAGATGCTTCATc                                                                       <  4:168924/275‑1 (MQ=255)
                                                                                                                                                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt                                                                                               >  2:431245/1‑250 (MQ=255)
                                                                                                                                                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATc                                                                       <  2:245681/274‑1 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                                                                                                    <  6:81876/132‑1 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                                                                                                   >  8:119424/1‑149 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAATCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                                                                                                    >  5:81876/1‑132 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTCTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTg                                                                                                                                                                                                          >  6:340909/1‑142 (MQ=255)
                                                                                                                                                              tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCACGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTg                                                                                                                                                                                                          <  5:340909/142‑1 (MQ=255)
                                                                                                                                                                  aCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                                                                                                        >  2:90254/1‑156 (MQ=255)
                                                                                                                                                                      ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                                                                                                    <  4:90253/156‑1 (MQ=255)
                                                                                                                                                                           cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAGACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                                                                                                         <  8:38410/146‑1 (MQ=255)
                                                                                                                                                                           cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                                                                                                         >  7:38410/1‑146 (MQ=255)
                                                                                                                                                                            gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                                                                                                                                                     >  6:78927/1‑149 (MQ=255)
                                                                                                                                                                              aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAGTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TGTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTCTAAGTCTGTGGTTACGTCACTCGAGAGATTTGTCACAATTATCTGCAAAGTAATACGCCGTTAATTTCTTTCTTTTTTTGTGTAAGCGTAAGATCACTCATCTGGTTTAAACCAAACGGa                                                     >  4:317016/1‑275 (MQ=255)
                                                                                                                                                                               ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGCAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                                                                                                        >  2:434470/1‑207 (MQ=255)
                                                                                                                                                                                 cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                                                                                                  >  4:286856/1‑195 (MQ=255)
                                                                                                                                                                                    aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGC‑AAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACa                                              <  2:388980/275‑1 (MQ=255)
                                                                                                                                                                                     aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                                                                                                                                                            <  7:119424/149‑1 (MQ=255)
                                                                                                                                                                                     aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                                                                                                              <  2:286860/195‑1 (MQ=255)
                                                                                                                                                                                         tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                                                                                                        <  8:3445/149‑1 (MQ=255)
                                                                                                                                                                                            ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                                                                                                                                                      <  6:254750/149‑1 (MQ=255)
                                                                                                                                                                                            ccTGACGGAAAACGATCGTCCCCGCGCCAACCGCGCACCGAAAGACGCTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGGGACCGCACACCAATGCGAGTATCCCCCTCATAAGTCCGTGGTCACGCCCATAGAGCGCCTTGCCCCAATTCCTCGCAACGCCATACGCCGCTACTTCCTTTCTTTTTTGGCGTAACCGTAAGATCCTCCATCGGCCTTAAACCAAAAGCATTaacacatgcca                                        >  4:124663/1‑265 (MQ=255)
                                                                                                                                                                                               gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                                                                                                                                                  <  5:169862/149‑1 (MQ=255)
                                                                                                                                                                                                aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                                                                                                                                                 <  6:367123/149‑1 (MQ=255)
                                                                                                                                                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttt                                                                                                                                                                                                                          >  5:160549/1‑86 (MQ=255)
                                                                                                                                                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttt                                                                                                                                                                                                                          <  6:160549/86‑1 (MQ=255)
                                                                                                                                                                                                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                                                                                                  >  3:145475/1‑174 (MQ=255)
                                                                                                                                                                                                         gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGCCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCAc                          <  2:317020/275‑1 (MQ=255)
                                                                                                                                                                                                          aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa                                                                                                                              <  1:145475/174‑1 (MQ=255)
                                                                                                                                                                                                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                                                                                                     <  8:53438/149‑1 (MQ=255)
                                                                                                                                                                                                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                                                                                                                                                     <  8:226458/149‑1 (MQ=255)
                                                                                                                                                                                                                  cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCCATCCCACGTAGaa                   >  2:508533/1‑273 (MQ=255)
                                                                                                                                                                                                                        aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac                                                                                                                                         >  8:150653/1‑149 (MQ=255)
                                                                                                                                                                                                                                     aaGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaaa                                                                                                                            <  5:78927/149‑1 (MQ=255)
                                                                                                                                                                                                                                        aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTc                                                                                                                         <  7:32418/149‑1 (MQ=255)
                                                                                                                                                                                                                                          gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTAAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCAt                                                                                                                       <  7:368424/149‑1 (MQ=255)
                                                                                                                                                                                                                                            ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAAGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACt  >  2:76206/1‑264 (MQ=255)
                                                                                                                                                                                                                                                tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCg                                                                                                                 <  8:183130/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                                                                                                       <  5:396201/106‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                                                                                                       >  6:396201/1‑106 (MQ=255)
                                                                                                                                                                                                                                                        tttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGc                                                                                                         <  7:150653/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttttt                                                                                                >  5:208931/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAACTCATACGCCGTTAATTGCTTTCttttt                                                                                                <  6:208931/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                                                                                                       >  5:122843/1‑109 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                                                                                                       <  6:122843/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc                                                                                            >  5:158623/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       ccTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgt                                                                                          >  7:429396/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         tGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa                                                                                        >  6:433576/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
CGCTCTAACGACCTGACTGTAAACTGCCTGACCGGTAAGAAACTGACCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACT  >  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: ‑