Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,230,246 T→C 5.5% H47R (CAT→CGT)  dusC ← tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,230,2460TC 5.5% 226.8 / 3.2 74H47R (CAT→CGT) dusCtRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base C (2/2);  total (38/36)
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.71e-01

GAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCAGTAACACACGCATATCATCAC  >  NC_000913/2230104‑2230393
                                                                                                                                              |                                                                                                                                                   
gagaTCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAg                                                                                                                                                <  8:180430/148‑1 (MQ=255)
  gaTCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCCGCGTACCAGATGGTGTCCGGCCGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACt                                                                                                                                             >  6:23780/1‑149 (MQ=255)
     ccACGCCCCAGGACCCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACttt                                                                                                                                           >  2:258103/1‑148 (MQ=255)
      cACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTa                                                                                                                                         >  5:30639/1‑149 (MQ=255)
          ccccAGGACCCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCgg                                                                                                                                     <  4:232256/149‑1 (MQ=255)
           cccAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGACGAAAGACTTTTACCgg                                                                                                                                     >  3:316330/1‑148 (MQ=255)
           cccAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGACGAAAGACTTTTACCgg                                                                                                                                     >  3:311234/1‑148 (MQ=255)
             cAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGgcag                                                                                                                                  >  6:195724/1‑149 (MQ=255)
             cAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGgcag                                                                                                                                  <  6:32745/149‑1 (MQ=255)
                    cTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCTCGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGACGAAAGACTTTTACCGGCAGCAGTTGa                                                                                                                           <  4:311234/149‑1 (MQ=255)
                    cTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGACGAAAGACTTTTACCGGCAGCAGTTGa                                                                                                                           <  4:316330/149‑1 (MQ=255)
                       aCTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATc                                                                                                                         <  3:178708/148‑1 (MQ=255)
                        cTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAATACTTTTACCGGCAGCAGTTGATcca                                                                                                                       >  6:186051/1‑149 (MQ=255)
                        cTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATcca                                                                                                                       <  4:140312/149‑1 (MQ=255)
                                   gggCGGCGTTCTCTGCCAGCCATTGTGGGAACTGCCCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCGCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGAc                                                                                                             >  6:288155/1‑148 (MQ=255)
                                     gcggcgTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACaaa                                                                                                          >  7:325657/1‑149 (MQ=255)
                                       ggcgTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACt                                                                                                        >  7:370330/1‑149 (MQ=255)
                                             ctctGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTg                                                                                                   <  5:186051/148‑1 (MQ=255)
                                                tGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGa                                                                                                  <  1:198157/146‑1 (MQ=255)
                                                  ccagccaTTGTTGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTACGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACATTTACCGTCAGCAGTTCATCCACCACGCGGACAAACTCGGTGATGca                                                                                              >  2:311278/1‑148 (MQ=255)
                                                     gccaTTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGca                                                                                              >  5:108861/1‑145 (MQ=255)
                                                     gccaTTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTAAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCCCCACGCGGACAAACTCGGTGATgac                                                                                              >  6:53461/1‑143 (MQ=255)
                                                        aTTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTg                                                                                                                            <  6:103196/112‑1 (MQ=255)
                                                        aTTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGAACAATATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGAAGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCg                                                                                       >  1:313238/1‑149 (MQ=255)
                                                         ttGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCg                                                                                       >  7:347794/1‑148 (MQ=255)
                                                           gtgGGAACTGACCTAACAACTGCACGCGCACCATCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACATTCTCGGTGATGCACAGATCGTa                                                                                     <  4:192987/148‑1 (MQ=255)
                                                           gtgGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCCGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTa                                                                                     <  5:23780/148‑1 (MQ=255)
                                                                aaCTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGATCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGt                                                                                >  2:213292/1‑148 (MQ=255)
                                                                   tGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCCCCACGCGGACAAACTCGGTGATGCACAGCTCGTAGTCGTTaaa                                                                            >  1:312807/1‑148 (MQ=255)
                                                                    gACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGAt                                                                                         >  7:321350/1‑135 (MQ=255)
                                                                    gACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGAt                                                                                         <  8:321350/135‑1 (MQ=255)
                                                                     aCCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACtt                                                                          >  8:1865/1‑149 (MQ=255)
                                                                     aCCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACtt                                                                          >  2:130869/1‑149 (MQ=255)
                                                                      ccTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGATCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTc                                                                         <  1:213292/149‑1 (MQ=255)
                                                                        tAACAACTGCACGCGCCCCAGCGTACCAGATTGTGTCCGGCTGCCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGCCTTTTACCGGCAGCAGTTTATCCACCACTCTGACAAAATCGGTGATGCACAGATCGTAGTCGTTAACtgagg                                                                       >  4:62474/1‑145 (MQ=255)
                                                                         aacaacTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATc                                                                                        >  8:79312/1‑131 (MQ=255)
                                                                         aacaacTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATc                                                                                        <  7:79312/131‑1 (MQ=255)
                                                                         aacaacTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGAGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGt                                                                      >  4:208620/1‑149 (MQ=255)
                                                                          acaacTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGt                                                                      <  7:1865/148‑1 (MQ=255)
                                                                          acaacTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGt                                                                      <  6:30639/148‑1 (MQ=255)
                                                                          acaacTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCA‑GGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCa                                                                    >  2:66418/1‑149 (MQ=255)
                                                                            aacTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGt                                                                      >  5:194495/1‑146 (MQ=255)
                                                                             acTTCACGCTCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGTGCAAATGAGATGAAAGACTTTTAACGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGc                                                                  <  1:311278/149‑1 (MQ=255)
                                                                                 cACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAAt                                                               <  8:347794/148‑1 (MQ=255)
                                                                                 cACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCA‑GGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTc                                                             <  1:66418/149‑1 (MQ=255)
                                                                                    gcgcACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGGCAGCAATTCAc                                                           <  1:258103/149‑1 (MQ=255)
                                                                                         ccAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCAcc                                                      <  5:195724/149‑1 (MQ=255)
                                                                                            gCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCaga                                                   <  3:208620/149‑1 (MQ=255)
                                                                                            gCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCa                                                     >  5:69389/1‑147 (MQ=255)
                                                                                                 ccAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTc                                              <  6:108861/149‑1 (MQ=255)
                                                                                                  cAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTc                                              >  1:165535/1‑148 (MQ=255)
                                                                                                    gATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAg                                           <  1:130869/149‑1 (MQ=255)
                                                                                                    gATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAg                                           >  6:312453/1‑149 (MQ=255)
                                                                                                    gATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGCCAAATGCGATGAAAGACTTATACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAg                                           >  1:103532/1‑149 (MQ=255)
                                                                                                          gtCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGAAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCActcc                                     >  7:298319/1‑149 (MQ=255)
                                                                                                           tCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCActccc                                    >  4:6562/1‑149 (MQ=255)
                                                                                                            ccGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTc                                                                                  >  8:222965/1‑102 (MQ=255)
                                                                                                            ccGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTc                                                                                  <  7:222965/102‑1 (MQ=255)
                                                                                                                   gCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCg                             <  5:288155/148‑1 (MQ=255)
                                                                                                                      ttttGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGc                         >  8:328029/1‑149 (MQ=255)
                                                                                                                         tGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCa                       <  5:53461/148‑1 (MQ=255)
                                                                                                                         tGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCa                       <  5:296765/148‑1 (MQ=255)
                                                                                                                         tGTAGCTCAGGGCAAATGCGATGAAAGACTTTAACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTTGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCa                       <  4:42378/148‑1 (MQ=255)
                                                                                                                          gTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCACCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACATATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCATAGAGTCAAGCACTCCGTCCATCGGTGCCAg                      <  4:3015/148‑1 (MQ=255)
                                                                                                                             gCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCAGTaa                   <  5:312453/148‑1 (MQ=255)
                                                                                                                                 aGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCAGTAacaca               <  2:165535/148‑1 (MQ=255)
                                                                                                                                  gggCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGc                                                                                               >  2:76705/1‑67 (MQ=255)
                                                                                                                                  gggCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGc                                                                                               <  1:76705/67‑1 (MQ=255)
                                                                                                                                      aaaTGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGAAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCAGTAACACACGCata         >  5:250049/1‑149 (MQ=255)
                                                                                                                                       aaTGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCAGTAACACACGCatat        >  7:312443/1‑149 (MQ=255)
                                                                                                                                       aaTGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGGCAAGCACTCCCTCCATCGGTGCCAGTAACACACGCatat        <  8:325657/149‑1 (MQ=255)
                                                                                                                                        aTGCGATGAAAGACTTTTACCGGCAGCAGTTGATCAACCACGCGGACAAACTCGGTGATGCACAGAGAGTAGTAGTTCACTTCTGTCAGCAAGTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTTCCATTAACACACGCatat        >  5:297508/1‑148 (MQ=255)
                                                                                                                                           cGATGAAAGACTTTTACCGGCAGCAGTTGATCCCCCACGCGGACAAACTCGGTGATGCACAGCTCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCAGTAACACACGCATatcat     >  5:114657/1‑148 (MQ=255)
                                                                                                                                           cGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCAGTAACACACGCATatcatc    <  3:7006/149‑1 (MQ=255)
                                                                                                                                             ttGAAAGACGTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGTTGCACAGATCGTAGGCGTTACCTTCGTTCAGCAATTCACGCCCCAGAGAGTCAAGCACTCCCTCCTTCGGTGCCAGTAACACACGCATTTCATCAc  <  2:313238/148‑1 (MQ=255)
                                                                                                                                              tGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCAc                                                       >  7:272182/1‑95 (MQ=255)
                                                                                                                                              tGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCAc                                                       <  8:272182/95‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GAGATCCACGCCCCAGGAACCTAACTCCACCGCACGGGCGGCGTTCTCTGCCAGCCATTGTGGGAACTGACCTAACAACTGCACGCGCACCAGCGTACCAGATGGTGTCCGGCTGGCGTTTTGTAGCTCAGGGCAAATGCGATGAAAGACTTTTACCGGCAGCAGTTGATCCACCACGCGGACAAACTCGGTGATGCACAGATCGTAGTCGTTAACTTCGGTCAGCAATTCACGCACCAGAGAGTCAAGCACTCCCTCCATCGGTGCCAGTAACACACGCATATCATCAC  >  NC_000913/2230104‑2230393

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