Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,631,200 A→G 5.0% pseudogene (1086/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,2000AG 5.0% 250.9 / 2.6 79pseudogene (1086/1284 nt)ydfJpseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base G (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGCACAGTTATGAATCACAAT  >  NC_000913/1631055‑1631343
                                                                                                                                                 |                                                                                                                                               
cAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATTGTGCAGAAAATACCCGCCAGGATAGGa                                                                                                                                              <  3:298317/149‑1 (MQ=255)
cAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGa                                                                                                                                              <  1:171612/149‑1 (MQ=255)
cAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGa                                                                                                                                              <  3:116997/149‑1 (MQ=255)
  gCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGAcc                                                                                                                                            <  8:236343/149‑1 (MQ=255)
  gCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGAcc                                                                                                                                            <  3:171076/149‑1 (MQ=255)
           aGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTTCAGAAAATACCCGCCAGTATAGGACCAAAACCgg                                                                                                                                    >  4:21856/1‑148 (MQ=255)
            gcgcACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCATGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGaa                                                                                                                                  <  5:10553/149‑1 (MQ=255)
                aCTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCaaaa                                                                                                                                        <  3:23509/139‑1 (MQ=255)
                aCTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCaaaa                                                                                                                                        >  4:23509/1‑139 (MQ=255)
                aCTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATACGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGa                                                                                                                              >  4:133190/1‑149 (MQ=255)
                 cTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGAt                                                                                                                             <  7:171715/149‑1 (MQ=255)
                  ttAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGAt                                                                                                                             >  7:91472/1‑148 (MQ=255)
                  ttAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGAt                                                                                                                             >  1:41999/1‑148 (MQ=255)
                    aaaTCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGAt                                                                                                                             >  1:162852/1‑146 (MQ=255)
                     aaTCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAg                                                                                                                          <  7:280031/148‑1 (MQ=255)
                         aCGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGaccac                                                                                                                     <  2:41999/149‑1 (MQ=255)
                           ggTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGaccacc                                                                                                                    >  6:382213/1‑148 (MQ=255)
                             tCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAAAATGCCATCATCATAATGGCGATCGCATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAAt                                                                                                                 >  5:98973/1‑149 (MQ=255)
                               tttCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGa                                                                                                                              <  2:401492/134‑1 (MQ=255)
                               tttCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGa                                                                                                                              >  1:401492/1‑134 (MQ=255)
                                ttCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATGGGACCAAAACCGGAAGCGATAAGACCACCAATTTc                                                                                                              >  7:95611/1‑149 (MQ=255)
                                   aCCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGAAACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCttt                                                                                                           <  5:272889/149‑1 (MQ=255)
                                   aCCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATACGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCttt                                                                                                           <  3:133190/149‑1 (MQ=255)
                                        tGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGaaa                                                                                                      >  7:268052/1‑149 (MQ=255)
                                            aTTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGTCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGc                                                                                                  <  6:250454/149‑1 (MQ=255)
                                             ttttCAGCGCAGAGATTAAACCAATCACTGAATATGCCAAGATCATAATGGCGATCGCATACCAGGATTCCGTCATGGTGCAGAAAATACACGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGc                                                                                                  >  6:358975/1‑148 (MQ=255)
                                                  aGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCTGATACCAGTATTCCGTCATGGTGCAGAAAATACCCTCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCAt                                                                                               <  5:151277/146‑1 (MQ=255)
                                                  aGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCTCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCAt                                                                                               >  6:151277/1‑146 (MQ=255)
                                                   gcgcAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCgg                                                                                           <  4:381921/149‑1 (MQ=255)
                                                    cgcAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCggg                                                                                          <  2:228681/149‑1 (MQ=255)
                                                           aTTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCg                                                                                   >  8:292795/1‑149 (MQ=255)
                                                               aaCCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGttt                                                                               >  1:360801/1‑149 (MQ=255)
                                                                 ccAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGttttt                                                                             <  8:3721/149‑1 (MQ=255)
                                                                  cAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTa                                                                            <  3:21856/149‑1 (MQ=255)
                                                                  cAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTa                                                                            <  3:309720/149‑1 (MQ=255)
                                                                  cAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTa                                                                            >  2:204727/1‑149 (MQ=255)
                                                                     tCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAg                                                                         <  2:134515/149‑1 (MQ=255)
                                                                       aCTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATGGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGcc                                                                       >  5:149728/1‑149 (MQ=255)
                                                                           aaTATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGaa                                                                    >  2:297673/1‑148 (MQ=255)
                                                                            atatGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACa                                                                  >  5:295149/1‑149 (MQ=255)
                                                                                gCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTc                                                              >  5:293119/1‑149 (MQ=255)
                                                                                 ccATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATGGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCt                                                             <  6:149728/149‑1 (MQ=255)
                                                                                  cATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTa                                                                            <  5:121335/133‑1 (MQ=255)
                                                                                  cATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTa                                                                            >  6:121335/1‑133 (MQ=255)
                                                                                  cATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTg                                                            >  8:112986/1‑149 (MQ=255)
                                                                                   aTGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGc                                                           >  6:310485/1‑149 (MQ=255)
                                                                                    tGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGc                                                           >  5:386926/1‑148 (MQ=255)
                                                                                            aTGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAAt                                                  <  1:204727/149‑1 (MQ=255)
                                                                                            aTGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAAt                                                  <  6:293119/149‑1 (MQ=255)
                                                                                               gCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGtt                                               <  2:360801/149‑1 (MQ=255)
                                                                                                     ggATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCa                                          <  1:297673/148‑1 (MQ=255)
                                                                                                     ggATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAACGTTTTCCa                                          <  8:268052/148‑1 (MQ=255)
                                                                                                       aTACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCagag                                       <  5:310485/149‑1 (MQ=255)
                                                                                                          ccAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCa                                     <  5:224084/148‑1 (MQ=255)
                                                                                                          ccAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCa                                     >  7:111459/1‑148 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCTCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGCGTAAAGCGGCTTTTTCAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGTGCaaa                                   <  6:98973/149‑1 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGCCCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCaaa                                   >  5:271617/1‑149 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAg                                                                                                   >  6:277952/1‑85 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAg                                                                                                   <  5:277952/85‑1 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCg                                                                                   <  1:261064/101‑1 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCg                                                                                   >  2:261064/1‑101 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTa                                                                            <  7:279892/108‑1 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTa                                                                            >  8:279892/1‑108 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTc                                            <  6:90292/140‑1 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTc                                            >  5:90292/1‑140 (MQ=255)
                                                                                                           cAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCaaa                                   >  3:305549/1‑149 (MQ=255)
                                                                                                                ttCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATc                              <  7:311289/149‑1 (MQ=255)
                                                                                                                  ccGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATcc                             <  7:112986/148‑1 (MQ=255)
                                                                                                                        tGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATTGGCATGCGGGTATAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCACAGCAcc                      >  7:204352/1‑149 (MQ=255)
                                                                                                                           tGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCCCCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGca                   >  6:142603/1‑149 (MQ=255)
                                                                                                                            gCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGca                   >  5:406500/1‑148 (MQ=255)
                                                                                                                             cAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGcaca                 <  2:162852/149‑1 (MQ=255)
                                                                                                                              aGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGCACAg                <  5:382213/149‑1 (MQ=255)
                                                                                                                                aaaaTACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGCACAGtt              <  6:386926/149‑1 (MQ=255)
                                                                                                                                  aaTACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGCACAGTTAt            >  5:74817/1‑149 (MQ=255)
                                                                                                                                   aTACCCGCCAGGATGGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGCACAGTTATg           <  8:95611/149‑1 (MQ=255)
                                                                                                                                    tACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGCACAGTTATGa          <  7:292795/149‑1 (MQ=255)
                                                                                                                                        cGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCCCCGCACAGTTATGAATca      >  5:232743/1‑149 (MQ=255)
                                                                                                                                              gATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCAAGCACCGCACAGTTATGAATCACAAt  >  5:14420/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGAGATTAAACCAATCACTGAATATGCCATGATCATAATGGCGATCGGATACCAGGATTCCGTCATGGTGCAGAAAATACCCGCCAGGATAGGACCAAAACCGGAAGCGATAAGACCACCAATTTCTTTAGAAATAGCCATCCGGGTAAAGCGGTTTTTACAGCCGAACATTTCTGCCATGGTAATGTTTTCCAGAGCAAATAATCCCAGCACCGCACAGTTATGAATCACAAT  >  NC_000913/1631055‑1631343

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