Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,733,959 T→A 5.6% S80S (TCT→TCA xylH → D‑xylose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,733,9590TA 5.6% 286.0 / 5.3 89S80S (TCT→TCAxylHD‑xylose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (44/40);  new base A (3/2);  total (47/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 = 7.98e-01

TGATGATTGCAGCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCGGTGCCTGGAACGGA  >  NC_000913/3733811‑3734103
                                                                                                                                                    |                                                                                                                                                
tgatgaTTGCAGCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTct                                                                                                                                                  >  5:314648/1‑149 (MQ=255)
 gatgatTGCAGCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCACGTAACGTATCAAACCAGTTACTCCAGCCCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGGTCGTCATAATTTctg                                                                                                                                                 >  3:259798/1‑149 (MQ=255)
  atgatTGCAGCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCTACCTGTTAGGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGTATGGTGTTCGTCATAATTTctgc                                                                                                                                                <  6:197703/149‑1 (MQ=255)
   tgatTGCAGCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTctgct                                                                                                                                               <  8:165149/149‑1 (MQ=255)
      ttGCAGCTATCATCGCAATCATGCTGTTCTTTACCTGGCCCACCGATGGTGCCTACTTAAGCGCCCGTATCGTCTCCAACCTGGTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGGTCGTCATAATTTCTGCTGa                                                                                                                                             >  1:224441/1‑148 (MQ=255)
         cAGCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATTGTGTTCGTCATAATTTCTGCTGAAATc                                                                                                                                         <  3:115683/149‑1 (MQ=255)
         cAGCTATCATCGCAATCATGCTCTTCTTTACATGGACCACCGATGGTGCCTAATTAAGCGCCCGTAACGTCTAAAACCTGTTACGCCAGACCGCGATTAACGGCATCCTAGCGGTAGGAATGGTGTTCGTCATAATTTATGCTGaaaaa                                                                                                                                         >  2:152544/1‑147 (MQ=255)
           gCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGa                                                                                                                                       <  1:232517/149‑1 (MQ=255)
           gCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCGCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGa                                                                                                                                       <  3:17654/149‑1 (MQ=255)
            cTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGAc                                                                                                                                      >  6:154480/1‑149 (MQ=255)
              atcatcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGTATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGAc                                                                                                                                      <  4:133390/147‑1 (MQ=255)
              atcatcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGTATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGAc                                                                                                                                      >  3:133390/1‑147 (MQ=255)
               tcatcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTAACGGCATCGTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCtt                                                                                                                                   >  3:155269/1‑149 (MQ=255)
                catcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCttt                                                                                                                                  >  2:4258/1‑149 (MQ=255)
                 atcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCttt                                                                                                                                  <  7:226063/148‑1 (MQ=255)
                 atcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCttt                                                                                                                                  <  5:195170/148‑1 (MQ=255)
                  tcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGa                                                                                                                                       >  4:195991/1‑142 (MQ=255)
                  tcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGa                                                                                                                                       <  3:195991/142‑1 (MQ=255)
                  tcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCAGCAGAAAtt                                                                                                                                         <  4:288560/140‑2 (MQ=255)
                  tcGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCAGCAGAAAtt                                                                                                                                         >  3:288560/1‑139 (MQ=255)
                   cGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCg                                                                                                                               >  7:240744/1‑149 (MQ=255)
                   cGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCg                                                                                                                               >  7:143660/1‑149 (MQ=255)
                      aaTCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCg                                                                                                                            >  3:167988/1‑149 (MQ=255)
                              tGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAat                                                                                                                     >  5:55204/1‑148 (MQ=255)
                               gTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAAt                                                                                                                                          <  3:75945/126‑1 (MQ=255)
                               gTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAAt                                                                                                                                          >  4:75945/1‑126 (MQ=255)
                                   tttACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGCCCTTTCCGTCGGCTCAATGATggg                                                                                                               >  4:171997/1‑149 (MQ=255)
                                   tttACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATggg                                                                                                               >  8:40936/1‑149 (MQ=255)
                                    ttACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATgggg                                                                                                              >  6:311654/1‑149 (MQ=255)
                                     tACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATggg                                                                                                               <  4:307744/147‑1 (MQ=255)
                                       ccTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTctgct                                                                                                                                               >  7:267967/1‑113 (MQ=255)
                                       ccTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTctgct                                                                                                                                               <  8:267967/113‑1 (MQ=255)
                                         tGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGt                                                                                                          >  7:231234/1‑148 (MQ=255)
                                         tGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACAGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGt                                                                                                          >  1:103472/1‑148 (MQ=255)
                                          ggACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCAGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGt                                                                                                          >  6:263960/1‑147 (MQ=255)
                                              caccGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTg                                                                                                    <  5:154480/149‑1 (MQ=255)
                                              caccGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTg                                                                                                    >  8:170098/1‑149 (MQ=255)
                                              caccGATGGTGCCTACTTAAGCGCCAGTAACGTATCCAACCTGTTACGCCAGACAGAGATTACCGGCATCCTCGAGGTAGGAATGGTTTTAGTCATAATTTCTGCTGAAATCGCCCTTTCCGTCGGCTCAATTATGGGGCTTTTAGGTg                                                                                                    >  1:248426/1‑149 (MQ=255)
                                               accGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTgg                                                                                                   >  5:190766/1‑149 (MQ=255)
                                                ccGCTGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGc                                                                                                  <  5:90970/149‑1 (MQ=255)
                                                ccGATGGTGCCTACTTAAGCGCCAGTCACGTCTCGATCCTGTTACGCAAGAGCGCGATTACCGGCATCCTCGCGGTAGGATTGGTGTTCGTCATAATTTATGCTTAAATCTACTTTTCCTTCTGCTCAATGTTGGGGCTGTTAGGTGGc                                                                                                  <  2:248426/149‑1 (MQ=255)
                                                   aTGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTc                                                                                               <  2:73241/149‑1 (MQ=255)
                                                       tGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGATTGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTAGGCTCAATGATGGGGCTTTTAGGTGGcttagctg                                                                                           >  2:100186/1‑142 (MQ=255)
                                                        gCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCgcggc                                                                                          >  4:67233/1‑149 (MQ=255)
                                                           tACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGAt                                                                                       >  4:41862/1‑149 (MQ=255)
                                                             cTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGAttt                                                                                     <  6:190766/149‑1 (MQ=255)
                                                             cTTAAGCGCCCGTAACGTCTCCAACCTGCTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGTGCTGTTAGGTGGCGTCGCGGCGAt                                                                                       >  7:306356/1‑147 (MQ=255)
                                                             cTTAAGAGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGAt                                                                                       <  6:208890/147‑1 (MQ=255)
                                                              ttAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTtg                                                                                    >  1:165295/1‑149 (MQ=255)
                                                                aaGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGTCTGTTAGGTGGCGTCGCGGCGATTtgt                                                                                   >  1:86482/1‑148 (MQ=255)
                                                                aaGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCAGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTtgt                                                                                   <  5:263960/148‑1 (MQ=255)
                                                                 aGCGCCCGTAACGTCTCCAACCTGCTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGa                                                                                 <  8:306356/149‑1 (MQ=255)
                                                                     cccGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTc                                                                             <  1:4258/149‑1 (MQ=255)
                                                                         tAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGt                                                                         <  8:231234/149‑1 (MQ=255)
                                                                           aCGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGt                                                                         >  6:146490/1‑147 (MQ=255)
                                                                           aCGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTa                                                                       <  6:163919/149‑1 (MQ=255)
                                                                             gTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATTGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAgg                                                                     <  5:192079/149‑1 (MQ=255)
                                                                              tctcCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGc                                                                    <  5:311654/149‑1 (MQ=255)
                                                                                  cAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGc                                                                >  3:130499/1‑149 (MQ=255)
                                                                                   aaCCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGcc                                                               <  3:67233/149‑1 (MQ=255)
                                                                                         ttACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGTCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGc                                                          <  2:86482/148‑1 (MQ=255)
                                                                                         ttACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGc                                                          <  7:19792/148‑1 (MQ=255)
                                                                                              ccAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGt                                                                              >  5:297039/1‑123 (MQ=255)
                                                                                              ccAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGt                                                                              <  6:297039/123‑1 (MQ=255)
                                                                                                 gACCGCGATTACCTGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACca                                                 <  4:155269/149‑1 (MQ=255)
                                                                                                   ccGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACcatc                                               <  4:130499/149‑1 (MQ=255)
                                                                                                    cgcgATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACcatca                                              <  3:201690/149‑1 (MQ=255)
                                                                                                    cgcgATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACcatca                                              >  7:24428/1‑149 (MQ=255)
                                                                                                    cgcgATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTAGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTCCcatca                                              >  5:201746/1‑149 (MQ=255)
                                                                                                     gcgATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACcatcat                                             >  5:39165/1‑149 (MQ=255)
                                                                                                      cgATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCAtt                                            <  7:170098/149‑1 (MQ=255)
                                                                                                        aTTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATtgt                                          <  6:55204/149‑1 (MQ=255)
                                                                                                               taaTCCTCGCGGTAGGAAAGGTCTACGTCATAATTTCTGCTGAAATCGACATATCCGACGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGCCGATTTGTGACGTCTGGATAGCCTGGCCTTTGACACTTACAATCAATGTGACGCTg                                   >  8:68758/3‑149 (MQ=255)
                                                                                                                   ccTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTc                               >  5:34361/1‑149 (MQ=255)
                                                                                                                    cTCGCGGAAGGAATGGTGTTCGTCATAATTTCTGCAGAAATCGACCTTTCCGTCGGCTCAATGAAGGGGCAGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTAAGCCCCTTACCATCATTGTGACGCTggat                                >  4:261592/1‑145 (MQ=255)
                                                                                                                     tCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGAttt                                                                                     <  5:237405/93‑1 (MQ=255)
                                                                                                                     tCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGACGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGAttt                                                                                     >  6:237405/1‑93 (MQ=255)
                                                                                                                       gcgGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTgg                            <  6:34361/148‑1 (MQ=255)
                                                                                                                       gcgGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTgg                            <  2:165295/148‑1 (MQ=255)
                                                                                                                              gAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCtt                    <  6:39165/149‑1 (MQ=255)
                                                                                                                                  ggTGTTCGTCATAATTTCTGCTGAAATTGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGa                                        >  3:48672/1‑125 (MQ=255)
                                                                                                                                  ggTGTTCGTCATAATTTCTGCTGAAATTGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGa                                        <  4:48672/125‑1 (MQ=255)
                                                                                                                                   gtgtTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCgg               <  5:117928/149‑1 (MQ=255)
                                                                                                                                   gtgtTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCg                <  3:41862/148‑1 (MQ=255)
                                                                                                                                    tgtTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCGGt              <  8:143660/149‑1 (MQ=255)
                                                                                                                                         gTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACCTCTGGTTAGGCAGGCCTTTGCCACTTACCATCAGTGTGACGCGGGTTCTGGGACTGCTTCTCGGTGCCTg         >  2:208868/1‑149 (MQ=255)
                                                                                                                                          tCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGc                                     >  2:187861/1‑120 (MQ=255)
                                                                                                                                          tCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGc                                     <  1:187861/120‑1 (MQ=255)
                                                                                                                                          tCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCGGTGCCTg         >  3:165406/1‑148 (MQ=255)
                                                                                                                                           gtcgaTTTCAGCAGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTTTGACGCTGGTTCTGGGACTGCTTCTCGGTGCCTgg        >  4:189846/5‑148 (MQ=255)
                                                                                                                                            aTAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCGGTGCCTGGaa      >  1:176838/1‑149 (MQ=255)
                                                                                                                                             tattttCTGCTGAAATCGACCTTTCAGTCGGATCAATGATTGGGCTGTTAGGTGGCGTCGCGGCGATTAGTGACGTCTGGTAAGGCTGGCCTTTGCCACTTACCCTCATTGTGACGCTCTTTCTGGGACTGCTTGTCGGTGACTGTAAc     >  6:335944/4‑149 (MQ=255)
                                                                                                                                             tAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCGGTGCCTGGAAc     <  6:283799/149‑1 (MQ=255)
                                                                                                                                              aaTTTCTGCTGGAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCGGTGCCTGGAACg    >  3:282664/1‑149 (MQ=255)
                                                                                                                                                tttCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGAGTTGTGACGTCTGGTTAGGCTGGCCTTTGCCCCTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCGGTGCCTGGAACGGa  >  8:296946/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGATGATTGCAGCTATCATCGCAATCATGCTGTTCTTTACCTGGACCACCGATGGTGCCTACTTAAGCGCCCGTAACGTCTCCAACCTGTTACGCCAGACCGCGATTACCGGCATCCTCGCGGTAGGAATGGTGTTCGTCATAATTTCTGCTGAAATCGACCTTTCCGTCGGCTCAATGATGGGGCTGTTAGGTGGCGTCGCGGCGATTTGTGACGTCTGGTTAGGCTGGCCTTTGCCACTTACCATCATTGTGACGCTGGTTCTGGGACTGCTTCTCGGTGCCTGGAACGGA  >  NC_000913/3733811‑3734103

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: