Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,631,023 C→T 6.6% pseudogene (1263/1284 nt) ydfJ ← pseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,631,0230CT 6.6% 439.9 / 14.4 136pseudogene (1263/1284 nt)ydfJpseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (65/62);  new base T (4/5);  total (69/67)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 139 total aligned reads displayed.

GCCTAAATGAACAATTCTTGGAGCCAGGTTATTAAGATCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAA  >  NC_000913/1630830‑1631277
                                                                                                                                                                                                 |                                                                                                                                                                                                                                                               
   tAAATGAACAATTCTTGGAGCCAGGTTATTAAGATCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAgc                                                                                                                                                                             >  2:499538/1‑275 (MQ=255)
                                              gggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACACTTCTTTAAATATAACTATGGCAAGCTATATCTTTTGTTATATGAATAAAAATCCCCTCCCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGCAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCCTTTTAAGCGCAGAGCTTAAACCAATCACTGAATATACACTGACCAAAATgg                                                                                                                                  >  4:609152/1‑275 (MQ=255)
                                              gggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcg                                                                                                                                                                                                                                                                <  5:771527/149‑1 (MQ=255)
                                               ggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcgc                                                                                                                                                                                                                                                               >  6:612171/1‑149 (MQ=255)
                                                   tgtCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGc                                                                                                                                                                                                                                                            <  7:922328/148‑1 (MQ=255)
                                                      cGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCtt                                                                                                                                                                                                                                                        <  6:159154/149‑1 (MQ=255)
                                                       gCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGTGCAGCTCtt                                                                                                                                                                                                                                                        <  7:759912/148‑1 (MQ=255)
                                                       gCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCtt                                                                                                                                                                                                                                                        <  7:30415/148‑1 (MQ=255)
                                                       gCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCtt                                                                                                                                                                                                                                                        <  8:342002/148‑1 (MQ=255)
                                                       gCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTa                                                                                                                                                                                                                                                       <  5:1055273/149‑1 (MQ=255)
                                                             gCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTacaaca                                                                                                                                                                                                                                                  >  8:80907/1‑148 (MQ=255)
                                                             gCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAacac                                                                                                                                                                                                                                                 >  5:239669/1‑149 (MQ=255)
                                                                gctAACAAATTATTTCCCATAACAATTCGTTAAATATAAATTTGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTTAGAGAAGGGATTAAGGGCTTACAGACTCCCGGAAGTTTCCGCAGCTCTTACAACAcgc                                                                                                                                                                                                                                               >  8:400175/3‑148 (MQ=255)
                                                                   aaCAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGtt                                                                                                                                                                                                                                           <  6:586745/149‑1 (MQ=255)
                                                                    aCAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTg                                                                                                                                                                                                                                          >  5:685989/1‑149 (MQ=255)
                                                                      aaaTTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGAt                                                                                                                                                                                                                                        >  5:550812/1‑149 (MQ=255)
                                                                      aaaTTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGAt                                                                                                                                                                                                                                        <  6:389547/149‑1 (MQ=255)
                                                                        attattTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATc                                                                                                                                                                                                                                       >  7:808119/1‑148 (MQ=255)
                                                                          tattTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCtt                                                                                                                                                                                                                                     <  7:162692/148‑1 (MQ=255)
                                                                                ccATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAg                                                                                                                                                                                                                               <  8:22468/148‑1 (MQ=255)
                                                                                ccATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGc                                                                                                                                                                                                                              >  5:546380/1‑149 (MQ=255)
                                                                                   tAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGt                                                                                                                                                                                                                            >  7:501864/1‑148 (MQ=255)
                                                                                   tAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGt                                                                                                                                                                                                                            >  7:695274/1‑148 (MQ=255)
                                                                                    aaCAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTc                                                                                                                                                                                                                           >  7:35780/1‑148 (MQ=255)
                                                                                    aaCAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCt                                                                                                                                                                                                                          >  5:890461/1‑149 (MQ=255)
                                                                                         ttCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTGCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGc                                                                                       <  2:720581/275‑1 (MQ=255)
                                                                                              taaatataaataTGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCa                                                                                                                                                                                                                 >  7:794499/1‑148 (MQ=255)
                                                                                                    taaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAAt                                                                                                                                                                                                          <  6:552289/149‑1 (MQ=255)
                                                                                                     aaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGGAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGTGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAAt                                                                                                                                                                                                          <  8:779509/148‑1 (MQ=255)
                                                                                                      aatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCa                                                                                                                                                                                                        >  6:549163/1‑149 (MQ=255)
                                                                                                        tatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAc                                                                                                                                                                                                       <  8:529693/148‑1 (MQ=255)
                                                                                                        tatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAc                                                                                                                                                                                                       <  8:351796/148‑1 (MQ=255)
                                                                                                            gCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACctc                                                                                                                                                                                          >  4:65727/1‑157 (MQ=255)
                                                                                                              aaGCTATATGTTTTGTTATAGGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCGTACAACACGCGGTTGATCATCAGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCAATTTCAGCGCAGAGATTAAACCAATAACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAAGACCCGCCAGGATAGGACCAAAACCGGa                                                                  <  2:609152/275‑1 (MQ=255)
                                                                                                               aGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCtt                                                                                                                                                                                                                                     >  7:784028/1‑111 (MQ=255)
                                                                                                               aGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCtt                                                                                                                                                                                                                                     <  8:784028/111‑1 (MQ=255)
                                                                                                               aGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATg                                                                                                                                                 >  4:749541/1‑195 (MQ=255)
                                                                                                                gCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGc                                                                                                                                                                                      <  2:65727/157‑1 (MQ=255)
                                                                                                                gCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTCGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCa                                                                                                                                                                                              >  6:723726/1‑149 (MQ=255)
                                                                                                                  tataTGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCAc                                                                                                                                                                                             >  8:290196/1‑148 (MQ=255)
                                                                                                                   atatGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACct                                                                                                                                                                                           >  6:710912/1‑149 (MQ=255)
                                                                                                                   atatGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCAt                                                                                                                                             <  2:749541/195‑1 (MQ=255)
                                                                                                                     atGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACctct                                                                                                                                                                                         <  6:890461/149‑1 (MQ=255)
                                                                                                                           tGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTa                                                                                                                                                                   >  2:192212/1‑165 (MQ=255)
                                                                                                                           tGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTTCAGAACATACCCGCCAGGAAAGGACCAAAACCGGAAACAAAAACacca                                                     >  4:196766/1‑275 (MQ=255)
                                                                                                                           tGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATa                                                                                            >  1:10817/1‑236 (MQ=255)
                                                                                                                              tataTGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGCCTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCAtttt                                                                                                                                                                                 <  8:1024235/148‑1 (MQ=255)
                                                                                                                              tataTGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATAc                                                                                                                       <  1:233800/206‑1 (MQ=255)
                                                                                                                              tataTGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATAc                                                                                                                       >  3:233800/1‑206 (MQ=255)
                                                                                                                               atatGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAAcc                                                                                                                                                               <  4:192212/165‑1 (MQ=255)
                                                                                                                               atatGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGc                                                                                       <  4:217598/237‑1 (MQ=255)
                                                                                                                                 atGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAg                                                                                                                                                                              >  8:705886/1‑148 (MQ=255)
                                                                                                                                  tGAATACAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGAtt                                                                                                                                                                    <  4:764085/157‑1 (MQ=255)
                                                                                                                                  tGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAAttt                                             >  2:328379/1‑275 (MQ=255)
                                                                                                                                   gAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAgcgc                                                                                                                                                                           >  6:223362/1‑149 (MQ=255)
                                                                                                                                    aaTAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAgcgc                                                                                                                                                                           >  8:432043/1‑148 (MQ=255)
                                                                                                                                     aTAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTcc                                                                                                                                                                                                                                   >  6:190228/1‑91 (MQ=255)
                                                                                                                                     aTAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTcc                                                                                                                                                                                                                                   <  5:190228/91‑1 (MQ=255)
                                                                                                                                     aTAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCag                                                                                                                                                                         >  5:83051/1‑149 (MQ=255)
                                                                                                                                      tAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGGTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCag                                                                                                                                                                         >  8:1045902/1‑148 (MQ=255)
                                                                                                                                        aaaaTCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCagag                                                                                                                                                                       <  7:432043/148‑1 (MQ=255)
                                                                                                                                        aaaaTCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAgaga                                                                                                                                                                      <  6:246757/149‑1 (MQ=255)
                                                                                                                                        aaaaTCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAgaga                                                                                                                                                                      >  5:582792/1‑149 (MQ=255)
                                                                                                                                         aaaTCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGAt                                                                                                                                                                     >  6:844861/1‑149 (MQ=255)
                                                                                                                                             ccccTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTaaa                                                                                                                                                                 <  5:131323/149‑1 (MQ=255)
                                                                                                                                             ccccTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTaa                                                                                                                                                                  >  8:157875/1‑148 (MQ=255)
                                                                                                                                             ccccTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTaa                                                                                                                                                                  <  7:376087/148‑1 (MQ=255)
                                                                                                                                             ccccTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTaa                                                                                                                                                                  <  8:501864/148‑1 (MQ=255)
                                                                                                                                             ccccTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTaa                                                                                                                                                                  <  7:688553/148‑1 (MQ=255)
                                                                                                                                              cccTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTaaa                                                                                                                                                                 <  8:808119/148‑1 (MQ=255)
                                                                                                                                              cccTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCa                                                                                     >  2:167499/1‑224 (MQ=255)
                                                                                                                                                  ctcCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCaa                                                                                                                                                             <  7:400175/148‑1 (MQ=255)
                                                                                                                                                    ccGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGAt                                                                                                                                                                                                                                        <  8:380496/71‑1 (MQ=255)
                                                                                                                                                       gTAAGAGAAGGGATTAAGGGTTTACAGACTTCTCGAAGGTTGCGCAGCTCTTACAACACGCGGTTGCTCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAAGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCACTTTCTTTAGAAATAGCCATCCgg                         >  2:648156/1‑274 (MQ=255)
                                                                                                                                                          agagAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATa                                                                                            >  3:233344/1‑205 (MQ=255)
                                                                                                                                                          agagAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATa                                                                                            <  1:233344/205‑1 (MQ=255)
                                                                                                                                                            agaAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAatat                                                                                                                                                  >  6:725272/1‑149 (MQ=255)
                                                                                                                                                             gaAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATg                                                                                                                                                 <  6:854799/149‑1 (MQ=255)
                                                                                                                                                              aaGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGTGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGc                                                                                                                                                <  5:407214/149‑1 (MQ=255)
                                                                                                                                                               aGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGTGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATg                                                                                                                                                 >  4:707913/1‑147 (MQ=255)
                                                                                                                                                               aGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCa                                                                                                                                                              >  2:605270/1‑134 (MQ=255)
                                                                                                                                                                gggATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCa                                                                                                                                              >  5:720612/1‑149 (MQ=255)
                                                                                                                                                                   aTTAAGGGTTTACAGACTTCTGGAAGGTTGTGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCAt                                                                                                                                             <  2:707913/147‑1 (MQ=255)
                                                                                                                                                                   aTTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCa                                                                                                                                                          <  4:605270/134‑1 (MQ=255)
                                                                                                                                                                   aTTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCAt                                                                                                         >  4:167548/1‑183 (MQ=255)
                                                                                                                                                                       aGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTg                                                                                                     <  2:167548/183‑1 (MQ=255)
                                                                                                                                                                         ggTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATa                                                                                                                                      <  7:554552/148‑1 (MQ=255)
                                                                                                                                                                         ggTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATa                                                                                                                                      <  7:878648/148‑1 (MQ=255)
                                                                                                                                                                           tttACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAAt                                                                                                                                    <  7:290196/148‑1 (MQ=255)
                                                                                                                                                                           tttACAGCTTTTGGAAGGTTGCGCAGCTTTTACAACACGCGTTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCGTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCg    <  4:499538/274‑1 (MQ=255)
                                                                                                                                                                            ttACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTaaa                                                                                                                                                                                                           >  4:778461/1‑76 (MQ=255)
                                                                                                                                                                            ttACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGa   >  4:129767/1‑275 (MQ=255)
                                                                                                                                                                            ttACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGa   >  2:341020/1‑275 (MQ=255)
                                                                                                                                                                             tACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTGTTTCACCTCTGGCATTTTCAGCGACGAGATTAAACCAATAACTGAATATGCCCTGATCATAATGGCAATCGGATACCAGGCTTCCGTCATGTTGAGGAACATACCCGCCAGGATAGGACCCAAACCGGAAGCGAGAAGCCACACAATTTCTTTAGAAATAACCATCCGGGAAAAGCGGTTTTTCCACCCTaa  >  4:379604/1‑275 (MQ=255)
                                                                                                                                                                             tACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGaa  >  2:84837/1‑275 (MQ=255)
                                                                                                                                                                                aGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAc                                                                                                                                                                                                       <  2:778461/76‑1 (MQ=255)
                                                                                                                                                                                aGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAatat                                                                                                                                                  <  5:540978/129‑1 (MQ=255)
                                                                                                                                                                                aGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAatat                                                                                                                                                  >  6:540978/1‑129 (MQ=255)
                                                                                                                                                                                         gAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCa                                                                                                                     <  5:35793/149‑1 (MQ=255)
                                                                                                                                                                                              ttGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGAt                                                                                                                 >  8:281422/1‑148 (MQ=255)
                                                                                                                                                                                                 |                                                                                                                                                                                                                                                               
GCCTAAATGAACAATTCTTGGAGCCAGGTTATTAAGATCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAA  >  NC_000913/1630830‑1631277

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