Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,620,278 C→A 42.6% R226S (CGT→AGT)  deoB → phosphopentomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,620,2780CA42.6% 47.1 / 85.1 75R226S (CGT→AGT) deoBphosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base A (18/14);  total (41/34)
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.53e-01

GATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAA  >  NC_000913/4620133‑4620424
                                                                                                                                                 |                                                                                                                                                  
gATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTa                                                                                                                                                 >  6:165444/1‑149 (MQ=255)
  ttGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTAcc                                                                                                                                               <  8:294758/149‑1 (MQ=255)
         ccATGAAGAAACTTTCTGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTAATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTaa                                                                                                                                          <  7:260896/147‑1 (MQ=255)
          cATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGCCCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTACCCg                                                                                                                                       >  4:7336/1‑149 (MQ=255)
                 aaaCTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAc                                                                                                                                <  7:34312/149‑1 (MQ=255)
                   aCTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCt                                                                                                                              <  6:2872/149‑1 (MQ=255)
                     tttCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTgg                                                                                                                            >  8:28858/1‑149 (MQ=255)
                           tCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTg                                                                                                                      <  8:122814/149‑1 (MQ=255)
                               gATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAgcc                                                                                                                  >  2:17859/1‑149 (MQ=255)
                               gATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAgcc                                                                                                                  <  5:1908/149‑1 (MQ=255)
                               gATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAgcc                                                                                                                  >  1:33716/1‑149 (MQ=255)
                               gATAAACTCTACGAACGGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATCTCGGGCGGGTTATCGCTGTTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCGGTTGAgcc                                                                                                                  <  8:273321/149‑1 (MQ=255)
                                   aaCTCTACGAACTGTGCGAAATCGCACGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTACAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCa                                                                                                              >  8:17307/1‑149 (MQ=255)
                                    aCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAg                                                                                                             >  1:153740/1‑149 (MQ=255)
                                        tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGtt                                                                                                                       >  7:78211/1‑135 (MQ=255)
                                        tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGtt                                                                                                                       <  8:78211/135‑1 (MQ=255)
                                        tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCacc                                                                                                         >  7:49915/1‑149 (MQ=255)
                                          cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCac                                                                                                          <  2:267389/146‑1 (MQ=255)
                                            aaCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGGTCGTCCGTTTATCGGCGACAAAGCCGGTATCTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCaccgacc                                                                                                     <  5:124348/149‑1 (MQ=255)
                                                 tgCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgct                                                                                                <  2:153740/149‑1 (MQ=255)
                                                   cGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctgc                                                                                              <  3:18930/149‑1 (MQ=255)
                                                   cGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctgc                                                                                              <  7:28858/149‑1 (MQ=255)
                                                   cGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctg                                                                                               >  3:83234/1‑148 (MQ=255)
                                                    gAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctg                                                                                               >  3:225608/1‑147 (MQ=255)
                                                    gAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtg                                                                                                  <  6:386040/144‑1 (MQ=255)
                                                    gAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtg                                                                                                  >  5:386040/1‑144 (MQ=255)
                                                      aaTCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGa                                                                                           <  4:259157/149‑1 (MQ=255)
                                                        tCGCCCGTGAAGAGCTGACCACCGGCGGCTACAATATCGGGCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGaaa                                                                                         <  5:127802/149‑1 (MQ=255)
                                                        tCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGaaa                                                                                         >  3:219207/1‑149 (MQ=255)
                                                          gCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATAGCTCGTCCGTTTATCTGCTACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACt                                                                                       >  6:122879/1‑149 (MQ=255)
                                                            ccGTGAAGAGCTGACCAACGGCGTCTACAATATCGGTCTTGTTATCGCTCGTCCGTTTATCGGAGACAAAGCCGGAAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTgg                                                                                     >  6:70540/1‑149 (MQ=255)
                                                              gTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGtt                                                                                   >  1:346250/1‑149 (MQ=255)
                                                                   gagCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTtgatg                                                                               >  6:237340/1‑148 (MQ=255)
                                                                         aCCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAaca                                                                        >  4:29991/1‑149 (MQ=255)
                                                                         aCCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAac                                                                         <  4:83234/148‑1 (MQ=255)
                                                                           cAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAaca                                                                        >  5:54164/1‑147 (MQ=255)
                                                                           cAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACg                                                                      >  5:194548/1‑149 (MQ=255)
                                                                            aaCGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTg                                                                                      >  3:102893/1‑132 (MQ=255)
                                                                            aaCGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTg                                                                                      <  4:102893/132‑1 (MQ=255)
                                                                                     tACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAggtggt                                                            >  1:99907/1‑149 (MQ=255)
                                                                                        aaTATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTc                                                         >  5:194153/1‑149 (MQ=255)
                                                                                            tCGGGCGGGGTAGCGCGCGTCCGTTTAGCGGCGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAACCACGGCCAGGTGGTTTCTGTc                                                     <  5:294075/149‑1 (MQ=255)
                                                                                               gTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGCCCAGGTGGTTTCTGTCgg                                                   >  6:273866/1‑148 (MQ=255)
                                                                                                   tgtTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTaaaa                                              <  3:29991/149‑1 (MQ=255)
                                                                                                      tATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTg                                           >  7:234064/1‑149 (MQ=255)
                                                                                                       aTCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGc                                          >  7:67407/1‑149 (MQ=255)
                                                                                                        tCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGTTAAAATTGCg                                         <  7:181817/149‑1 (MQ=255)
                                                                                                        tCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCg                                         <  4:219207/149‑1 (MQ=255)
                                                                                                        tCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGc                                          <  6:194548/148‑1 (MQ=255)
                                                                                                        tCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCg                                         >  7:111060/1‑149 (MQ=255)
                                                                                                            tCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACa                                     >  2:231700/1‑149 (MQ=255)
                                                                                                               tCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACGGTGCTGCAGAAACTGGTTGATTAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCt                                  <  8:49915/149‑1 (MQ=255)
                                                                                                               tCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAc                                      >  4:109443/1‑145 (MQ=255)
                                                                                                               tCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAc                                      <  3:109443/145‑1 (MQ=255)
                                                                                                                ccGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGACGCCAGCACCGACCGTGCAGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTa                                 >  7:284875/1‑149 (MQ=255)
                                                                                                                 cGTTTATCGTCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGCCCTTGCTGTTGAGCCGCCAGCACCGACCGTGCTGAAGAACCTGGTGGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTa                                 <  7:164927/148‑1 (MQ=255)
                                                                                                                       tCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAAc                          >  8:249388/1‑149 (MQ=255)
                                                                                                                             aCAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAACCACGGCCAggtgg                                                             >  5:58633/1‑108 (MQ=255)
                                                                                                                               aaaGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGt                                                      >  7:276505/1‑113 (MQ=255)
                                                                                                                               aaaGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGt                                                      <  8:276505/113‑1 (MQ=255)
                                                                                                                               aaaGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGAAGAAACTGGATGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTAt                  >  7:228472/1‑149 (MQ=255)
                                                                                                                               aaaGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAAtt                                            >  5:109284/1‑123 (MQ=255)
                                                                                                                               aaaGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAAtt                                            <  6:109284/123‑1 (MQ=255)
                                                                                                                                aaGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATc                 <  7:72270/149‑1 (MQ=255)
                                                                                                                                   ccGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCAc               >  1:184417/1‑148 (MQ=255)
                                                                                                                                      gTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCaa            >  8:205071/1‑148 (MQ=255)
                                                                                                                                         aCTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTa                                 >  6:212633/1‑124 (MQ=255)
                                                                                                                                         aCTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTa                                 <  5:212633/124‑1 (MQ=255)
                                                                                                                                            tCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      >  1:267945/1‑148 (MQ=255)
                                                                                                                                            tCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACTGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      >  6:83073/1‑148 (MQ=255)
                                                                                                                                             ccAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACAGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      <  6:54164/147‑1 (MQ=255)
                                                                                                                                             ccAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      <  2:81904/147‑1 (MQ=255)
                                                                                                                                             ccAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      <  2:33716/147‑1 (MQ=255)
                                                                                                                                               aGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGaaa  <  5:316060/149‑1 (MQ=255)
                                                                                                                                               aGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGa    >  5:363307/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAA  >  NC_000913/4620133‑4620424

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