Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,072,982 C→T 5.8% E344K (GAA→AAA)  rutA ← pyrimidine oxygenase, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,072,9820CT 5.8% 217.4 / 2.3 69E344K (GAA→AAA) rutApyrimidine oxygenase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

CGGTCGAGCGGTTAAGGTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTC  >  NC_000913/1072840‑1073127
                                                                                                                                              |                                                                                                                                                 
cGGTCGAGCGGTTAAGGTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCAcc                                                                                                                                             <  2:286803/149‑1 (MQ=255)
 ggTCGAGCGGTTAAGGTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCa                                                                                                                                            <  8:33043/149‑1 (MQ=255)
 ggTCGAGCGGTTAAGGTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCa                                                                                                                                            <  7:192220/149‑1 (MQ=255)
          gTTAAGGTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCtt                                                                                                                                   <  1:376213/149‑1 (MQ=255)
           ttAAGGTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTg                                                                                                                                  >  8:237824/1‑149 (MQ=255)
              aGGTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGa                                                                                                                               <  2:262817/149‑1 (MQ=255)
               ggTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAACTCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGAc                                                                                                                              >  6:39039/1‑149 (MQ=255)
                gtcgtcATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGAGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATAGTCTAAGGTTAACAGCACGCCTTAGGCACCAGGCACGCTTGCGACt                                                                                                                             <  4:9827/149‑1 (MQ=255)
                   gtcATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTc                                                                                                                           <  5:219504/148‑1 (MQ=255)
                    tcatcaTGCCACCTCCTGCCTCAGCGCAGGTAGATGGGCGCGGCATTGCATCACTCGTGGAAGGCGCGCGCCGAAGGGGTCGACTCCCGACAGAAACTCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCAt                                                                                                                         <  3:186647/149‑1 (MQ=255)
                       tcatGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCt                                                                                                                       >  3:311282/1‑148 (MQ=255)
                         atGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTaa                                                                                                                     <  4:117242/148‑1 (MQ=255)
                         atGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGAATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTaa                                                                                                                     >  2:249930/1‑148 (MQ=255)
                          tGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACa                                                                                                                   <  5:181625/149‑1 (MQ=255)
                          tGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGCCATTGCAAAAGTGGTTGAATCCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCCGCACCAGGCACGCTAGCGACTTCATCTAACa                                                                                                                   >  8:270930/1‑149 (MQ=255)
                             cACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATgc                                                                                                                <  6:66446/149‑1 (MQ=255)
                                 tcctGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACa                                                                                                                   >  4:51437/1‑142 (MQ=255)
                                 tcctGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACa                                                                                                                   <  3:51437/142‑1 (MQ=255)
                                 tcctGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATgcgcgc                                                                                                            >  3:315278/1‑149 (MQ=255)
                                   ctGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATgcgcgc                                                                                                            <  2:375046/147‑1 (MQ=255)
                                   ctGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATgcgcgc                                                                                                            <  8:376215/147‑1 (MQ=255)
                                      aGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTTGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGacac                                                                                                       >  4:143361/1‑149 (MQ=255)
                                      aGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTTGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGacac                                                                                                       >  4:148573/1‑149 (MQ=255)
                                      aGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGaca                                                                                                        <  8:113457/148‑1 (MQ=255)
                                       gTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACt                                                                                                      >  3:132706/1‑149 (MQ=255)
                                        tCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAAGCACGCATGCCACTTCATCTAACATGCGCGCGACACTg                                                                                                     >  2:178602/1‑149 (MQ=255)
                                           gcgcAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCg                                                                                                  >  8:91838/1‑149 (MQ=255)
                                             gcAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGt                                                                                                 >  3:26388/1‑148 (MQ=255)
                                                   aGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAAcc                                                                                          >  3:37110/1‑149 (MQ=255)
                                                        ggCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACt                                                                                      <  2:269699/148‑1 (MQ=255)
                                                         gcgcGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTaa                                                                                    >  1:150467/1‑149 (MQ=255)
                                                              gCATGGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACACAAAATCGTCGACGGTTAACAGCACGCCATCGGCACCAGGCACGCTTGCGACTTCATCTAACCTGCGCGCGACACTGGCGTAAGAACCGACTAACGTcc                                                                               <  6:68018/149‑1 (MQ=255)
                                                                 ttGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCAACAGGCACGCTTGAGACTTCATCTAAAATGCGCTCGCCACTGGCGTAAGAACCGACTAACGTCCCCa                                                                            >  8:324105/1‑149 (MQ=255)
                                                                        aGTGGTTTTATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTATCATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGa                                                                      <  1:178602/148‑1 (MQ=255)
                                                                            gTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGtg                                                                  >  8:329071/1‑147 (MQ=255)
                                                                               gAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTAcc                                                              >  6:153374/1‑149 (MQ=255)
                                                                                 aTGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGACCCCATATTGATGTTTAccgc                                                            >  4:197340/1‑149 (MQ=255)
                                                                                 aTGCGCTCGCCGAAGGTTTCGATTCCAGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTAccgc                                                            <  8:212167/149‑1 (MQ=255)
                                                                                  tGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTAccgcc                                                           >  1:146876/1‑149 (MQ=255)
                                                                                  tGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTAccgcc                                                           >  1:151908/1‑149 (MQ=255)
                                                                                  tGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTAccgcc                                                           >  6:124397/1‑149 (MQ=255)
                                                                                       tCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTaa                                                                                               <  2:160277/108‑1 (MQ=255)
                                                                                       tCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTaa                                                                                               >  1:160277/1‑108 (MQ=255)
                                                                                         gCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTgg                                                    <  7:91838/149‑1 (MQ=255)
                                                                                         gCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTgg                                                    <  4:26388/149‑1 (MQ=255)
                                                                                            gAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGAt                                                 <  7:270930/149‑1 (MQ=255)
                                                                                               ggTTTCGATTCCAGACATAAAATCGTAGAAGGTTAACAGCACGCCTTCGGCTCCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGAGTAACGTCCACATATTGATTTTTACCGACGAAGTGGGATCgg                                              <  7:324105/149‑1 (MQ=255)
                                                                                                   tCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCa                                           <  4:37110/148‑1 (MQ=255)
                                                                                                     gATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATc                                         <  4:311282/148‑1 (MQ=255)
                                                                                                     gATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATc                                         <  3:197340/148‑1 (MQ=255)
                                                                                                          ccGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAATGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACg                                    >  6:356778/1‑148 (MQ=255)
                                                                                                           cGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAATGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGaa                                  >  8:369965/1‑149 (MQ=255)
                                                                                                             aCAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACg                                >  7:3363/1‑149 (MQ=255)
                                                                                                             aCAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAAAATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACg                                >  5:181049/1‑149 (MQ=255)
                                                                                                                gAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTg                             <  6:147491/149‑1 (MQ=255)
                                                                                                                  aaaTCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCgg                                              >  4:248001/1‑130 (MQ=255)
                                                                                                                  aaaTCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCgg                                              <  3:248001/130‑1 (MQ=255)
                                                                                                                  aaaTCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGgt                           <  5:39039/149‑1 (MQ=255)
                                                                                                                  aaaTCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGgt                           <  1:249930/149‑1 (MQ=255)
                                                                                                                       gtcgAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCa                       <  4:315278/148‑1 (MQ=255)
                                                                                                                       gtcgAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAg                      <  1:152508/149‑1 (MQ=255)
                                                                                                                          gAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATc                                                >  8:236361/1‑120 (MQ=255)
                                                                                                                          gAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATc                                                <  7:236361/120‑1 (MQ=255)
                                                                                                                                  aCAGCACGCCTTTGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCg            <  3:143361/148‑1 (MQ=255)
                                                                                                                                  aCAGCACGCCTTTGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCg            <  3:148573/148‑1 (MQ=255)
                                                                                                                                  aCAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCgg           >  4:153312/1‑149 (MQ=255)
                                                                                                                                    aGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGt         <  8:3363/149‑1 (MQ=255)
                                                                                                                                      cACGCCTTCGGCACCACGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTAt       <  7:237824/149‑1 (MQ=255)
                                                                                                                                           cTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTc  >  6:336834/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CGGTCGAGCGGTTAAGGTCGTCATCATGCCACCTCCTGAGTCAGCGCAGGTAGATGGGCGCGGCATTGCATCAGTGGTTGAATGCGCTCGCCGAAGGTTTCGATTCCCGACAGAAAATCGTCGAAGGTTAACAGCACGCCTTCGGCACCAGGCACGCTTGCGACTTCATCTAACATGCGCGCGACACTGGCGTAAGAACCGACTAACGTCCCCATATTGATGTTTACCGCCGAAGTGGGATCGGCCATCTGACGAACGTTGGTGTCAGTACCGGAGCGGGTATCTTTC  >  NC_000913/1072840‑1073127

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: ‑