Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,647,434 A→T 5.1% pseudogene (203/279 nt) ydfX ← pseudogene, Qin prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,647,4340AT 5.1% 265.3 / 3.3 79pseudogene (203/279 nt)ydfXpseudogene, Qin prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (38/37);  new base T (2/2);  total (40/39)
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

GTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATTCCTG  >  NC_000913/1647288‑1647581
                                                                                                                                                  |                                                                                                                                                   
gTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCaaa                                                                                                                                                    >  1:405132/1‑148 (MQ=255)
   ttttGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAatat                                                                                                                                                 <  3:269781/148‑1 (MQ=255)
    tttGAGCGTAAATCAATACGGGCATTATCCTCTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATc                                                                                                                                                <  8:83881/148‑1 (MQ=255)
      tGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCtt                                                                                                                                              <  7:249658/148‑1 (MQ=255)
        aGCGTAAATCAGTACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCAc                                                                                                                                           <  5:191132/149‑1 (MQ=255)
        aGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCAc                                                                                                                                           <  5:362229/149‑1 (MQ=255)
             aaaTCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGAt                                                                                                                                      >  4:19124/1‑149 (MQ=255)
                tCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTaa                                                                                                                                   >  2:118312/1‑149 (MQ=255)
                  aaTACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAAc                                                                                                                                  <  5:14498/148‑1 (MQ=255)
                   aTACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGc                                                                                                                                >  8:120072/1‑149 (MQ=255)
                    tACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGc                                                                                                                                >  4:215924/1‑148 (MQ=255)
                        ggCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAgc                                                                                                                           <  5:112778/149‑1 (MQ=255)
                        ggCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAgc                                                                                                                           >  5:262032/1‑149 (MQ=255)
                           aTTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgct                                                                                                                        <  8:324704/149‑1 (MQ=255)
                               tCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTTAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcgc                                                                                                                     <  7:369497/148‑1 (MQ=255)
                               tCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGcc                                                                                                                               >  5:179815/1‑138 (MQ=255)
                               tCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGcc                                                                                                                               <  6:179815/138‑1 (MQ=255)
                                ccACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAgc                                                                                                                           >  8:58512/1‑141 (MQ=255)
                                ccACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAgc                                                                                                                           <  7:58512/141‑1 (MQ=255)
                                ccACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcgct                                                                                                                    >  1:347623/1‑148 (MQ=255)
                                 cACTCGTTCATGGACAGTCTTTTCATCATCCAGTGGTTGCATCCATAATTGCAGGCACTAGCGCCGTTAACTTTTCACGATTATCCCTGGTGTCGATAGCCTTCCAGCGTTCTAATATCTTCTCTCGATTAAAGCCAAGCTATCGCTGa                                                                                                                  <  8:91275/149‑1 (MQ=255)
                                 cACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGa                                                                                                                  <  5:301835/149‑1 (MQ=255)
                                      gTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAAt                                                                                                             <  6:262032/149‑1 (MQ=255)
                                          aTGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCATATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATcgcg                                                                                                         >  2:386953/1‑149 (MQ=255)
                                          aTGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCATATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATcgcg                                                                                                         >  4:128173/1‑149 (MQ=255)
                                             gACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctc                                                                                                                       <  2:292331/132‑1 (MQ=255)
                                             gACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctc                                                                                                                       >  1:292331/1‑132 (MQ=255)
                                                 ctctTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCATGCGCTCGCTGATCAATCGCGCCACCtt                                                                                                  <  1:118312/149‑1 (MQ=255)
                                                 ctctTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCtt                                                                                                  >  1:225911/1‑149 (MQ=255)
                                                  tctTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCtt                                                                                                  >  5:54440/1‑148 (MQ=255)
                                                       tcatcaTCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCatat                                                                                             <  2:405132/148‑1 (MQ=255)
                                                       tcatcaTCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCatat                                                                                             <  3:140713/148‑1 (MQ=255)
                                                       tcatcaTCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCatat                                                                                             <  2:225911/148‑1 (MQ=255)
                                                          gcatcCAGTGGAAGCGTGCATAATTGCAGGCTCTAGCGCCGGGAACTTTTCCCGCTTATCCCTGCTGGCGATACCCTTCCAGGGAGCAAAGATCTTCTCTCGATTCACCACACGCGCTCGCTCATCAATCGCGCCACCTTCAAATGTGa                                                                                         >  4:60997/2‑149 (MQ=255)
                                                               cAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGc                                                                                    <  2:347623/149‑1 (MQ=255)
                                                                         cTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCat                                                                                               <  7:163324/127‑1 (MQ=255)
                                                                         cTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCat                                                                                               >  8:163324/1‑127 (MQ=255)
                                                                          tCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCACTCGCGCCACCTTCATATGTGACACGCTGAACATCGAt                                                                         >  2:207039/1‑149 (MQ=255)
                                                                          tCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGAt                                                                         >  8:331307/1‑149 (MQ=255)
                                                                           ccATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATg                                                                        >  7:390722/1‑149 (MQ=255)
                                                                            cATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGAt                                                                         <  3:215924/147‑1 (MQ=255)
                                                                                aTTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCgg                                                                   <  1:77479/149‑1 (MQ=255)
                                                                                 ttGCAGGCACTAGCGCCGTGAACTTTTCCCGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCggg                                                                  >  8:336948/1‑149 (MQ=255)
                                                                                   gCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGgc                                                                 <  3:19124/148‑1 (MQ=255)
                                                                                    cAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGgcgc                                                               >  3:329622/1‑149 (MQ=255)
                                                                                       gCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGctct                                                            >  8:137572/1‑149 (MQ=255)
                                                                                         aCTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCttt                                                          >  4:222170/1‑149 (MQ=255)
                                                                                             gcgcCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCaaa                                                      <  6:54440/149‑1 (MQ=255)
                                                                                               gcCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGc                                                    >  7:343295/1‑149 (MQ=255)
                                                                                                  gTGAACTTTTCACGCTTATCACTGTTGTCGATAGCCTTCCAGCTTTCAAATATCTTCACTCGATTAACGACAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACAGAGATTTTCGGTCGCTCTTTCAAATcc                                                   >  6:395545/1‑147 (MQ=255)
                                                                                                   tGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTACAGCGTTCAAATATCTTCACTAGATTAACGCCAAGCGCTCGCTGATCAATCTCGCCACCTTCATATGTGACACGCTGAACATCTATGTTCGGGCGCTCTTTCAAAGcc                                                   >  8:44605/1‑146 (MQ=255)
                                                                                                        ttttCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGc                                           >  5:133475/1‑149 (MQ=255)
                                                                                                        ttttCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGc                                           >  2:215712/1‑149 (MQ=255)
                                                                                                        ttttCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGc                                           >  3:327258/1‑149 (MQ=255)
                                                                                                        ttttCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGc                                           >  6:275012/1‑149 (MQ=255)
                                                                                                           tCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGAGACGCTGAAGATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCtt                                         <  7:44605/148‑1 (MQ=255)
                                                                                                             aCGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCa                                       <  3:222170/148‑1 (MQ=255)
                                                                                                              cGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCATATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGt                                     <  3:128173/149‑1 (MQ=255)
                                                                                                              cGCTTATCCCTGGTGTCGATAGCCATCCAGCGTTCATATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGt                                     <  1:386953/149‑1 (MQ=255)
                                                                                                                 ttATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGAt                                  >  5:7084/1‑149 (MQ=255)
                                                                                                                   aTCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGAt                                  <  8:390722/147‑1 (MQ=255)
                                                                                                                    tCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTaa                               >  5:398291/1‑149 (MQ=255)
                                                                                                                    tCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTaa                               >  1:97273/1‑149 (MQ=255)
                                                                                                                          gtgtCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcgt                         <  1:389197/149‑1 (MQ=255)
                                                                                                                           tgtCGATCGCGTTCCAGCGTTCAAATATCTTCTCTCGATTAAGGCCAAGCGCTCGCTGATCAATGGCGCCACGTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAATCCCAGAATGCTTCAGTGATTAATAtcgtc                        <  3:60997/149‑1 (MQ=255)
                                                                                                                           tgtCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcgtc                        <  7:137572/149‑1 (MQ=255)
                                                                                                                               gATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACGCGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCt                    >  2:80946/1‑149 (MQ=255)
                                                                                                                               gATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCt                    >  7:24480/1‑149 (MQ=255)
                                                                                                                                 tAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGATCATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGc                  >  4:360751/1‑149 (MQ=255)
                                                                                                                                   gCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTc                <  7:331307/149‑1 (MQ=255)
                                                                                                                                   gCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCGCGCTGATCAATCGCGCCACCTGCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTc                <  1:207039/149‑1 (MQ=255)
                                                                                                                                     cTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGc                      >  4:373893/1‑141 (MQ=255)
                                                                                                                                     cTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGc                      <  3:373893/141‑1 (MQ=255)
                                                                                                                                         cAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGacac                                                                                      >  2:320407/1‑73 (MQ=255)
                                                                                                                                         cAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGacac                                                                                      <  1:320407/73‑1 (MQ=255)
                                                                                                                                            cGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGc                      >  3:249211/1‑134 (MQ=255)
                                                                                                                                            cGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGc                      <  4:249211/134‑1 (MQ=255)
                                                                                                                                             gTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCAtt      >  8:111399/1‑149 (MQ=255)
                                                                                                                                             gTTAAAATATCTTCACTAGATTAACGCCAAGCGCTCGCTGATCAATAGCGCCAACTTAATATGTTACACGCTGAACATCGATGTTCGGGCTCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCAtt      <  5:395545/149‑1 (MQ=255)
                                                                                                                                              ttCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATTc     >  5:400314/1‑149 (MQ=255)
                                                                                                                                                 aaaTATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCttt                                                          <  3:217407/93‑1 (MQ=255)
                                                                                                                                                 aaaTATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCttt                                                          >  4:217407/1‑93 (MQ=255)
                                                                                                                                                 aaaTATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATTCCTg  >  2:388279/1‑149 (MQ=255)
                                                                                                                                                 aaaTATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGATATTTCAAAGCCCAGAATGGTACAGTGATTAATATCGTCGCCTGCTCCTGTTTCATTCCTg  >  6:274651/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATTCCTG  >  NC_000913/1647288‑1647581

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