Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,717 T→C 11.3% intergenic (+145/‑295) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,7170TC11.3% 107.7 / 4.2 62intergenic (+145/‑295)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (17/31);  new base C (5/3);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCG  >  NC_000913/349571‑349865
                                                                                                                                                   |                                                                                                                                                    
aaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                                                                                                                    <  2:173377/149‑1 (MQ=255)
       gCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATAGCGTCGCATCCGACAACCGCTGCCTGCTGCGCCGCTGACGGGAGTTCTCa                                                                                                                                             >  1:513168/1‑149 (MQ=17)
          cGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGc                                                                                                                                          >  2:396219/1‑149 (MQ=255)
            ttGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCt                                                                                                                                        >  1:636062/1‑149 (MQ=255)
            tggATTGTCTTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCt                                                                                                                                        <  1:342434/147‑1 (MQ=255)
              gATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTAc                                                                                                                                      >  1:464463/1‑149 (MQ=255)
                 tGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                   <  1:252803/149‑1 (MQ=255)
                     ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGc                                                                                                                               <  2:497533/149‑1 (MQ=255)
                      tCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGAGGCGCCGCTGACGGGGCTTATCAGGCCTACGGGGGGCa                                                                                                                              >  2:154300/1‑149 (MQ=255)
                      gcGACAGCCGATGCCTGCTGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCTTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCa                                                                                                                              <  2:371614/148‑1 (MQ=255)
                        gACAGCCGATGCCTGATGCTGCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACt                                                                                                                            <  2:589785/149‑1 (MQ=255)
                           aGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATTCGCCGCTGACGGGACTTACCAGGCCTACGGGGCGCACTGaa                                                                                                                         >  1:279736/1‑149 (MQ=255)
                           aGCCGACGCCTTATGCGTCGCTGACGCCACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGCAAACTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                                                                         <  2:483041/149‑1 (MQ=17)
                              cGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGGTC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTg                                                                                                                      <  1:317707/149‑1 (MQ=255)
                                 tggcggATGCGTTGCTGAAGCGACTTATCAGGCCTACGGGGTGCGAAGGTACAAT‑GGAAGGAAAGGACCGAGGGCCTCGCATCCGCCATTCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                                                                   <  2:517228/144‑1 (MQ=9)
                                 tGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                                                                   <  1:472051/149‑1 (MQ=255)
                                 ggCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                                                                   <  1:344599/148‑1 (MQ=255)
                                        tGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                                                                            >  1:132203/1‑149 (MQ=255)
                                                       aCTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                             <  2:284837/149‑1 (MQ=40)
                                                       aCTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                             <  2:369791/149‑1 (MQ=39)
                                                       aCTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                             <  2:369811/149‑1 (MQ=255)
                                                       aCTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                             <  1:603150/149‑1 (MQ=40)
                                                       aCTTAGCAGGCGTACGAGGTGC‑ACTGGACTGTAGGTAGGGTAAGCCGGATGCCGTCGGCACCGACAACCGATGCCTAATGCTCCGCTTACGCTACTTATCAAGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                             <  2:174074/149‑1 (MQ=11)
                                                         ttATCAGGCCTACGAGGTGC‑ACTGAACTTTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGc                                                                                           >  1:390922/1‑149 (MQ=255)
                                                         ttATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCCTCCGACAACCCATACCTGATGCGCCGCTAACGCAACTTATCAGGCCCACGCGGCGCACTGAACTGTAGGTCGGATAAGACCCGAAGcgccgc                                                                                           >  2:94788/1‑149 (MQ=11)
                                                         ttATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGGGCACTGCACTGTAGGACGGCTAAGACGCATAGCGTCGc                                                                                           >  2:215816/1‑149 (MQ=255)
                                                           aTCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                   >  1:576850/1‑107 (MQ=255)
                                                            tCAGGCCTTCGAGGTGC‑ACTGAACGGTAGGTCGGATAAGACGGATGGCGTCGCTTCCGACAACCGATGCCTGATGCGCCGCTGACGTGATTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAATACGCATAGCGTCGCATc                                                                                        <  2:542183/149‑1 (MQ=14)
                                                              aGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGGCTTATCAGGCCTACGGGGGGAACTGAACTGTATGGCGGATAAGACGCATAGGGGCGCGTCCg                                                                                      >  2:466524/1‑149 (MQ=11)
                                                               ggCCTACGATGGGC‑ACTGAAATGTGGGTCGGAAAAGCCGCCTGGCGTCGGATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                                                                     <  2:376005/149‑1 (MQ=11)
                                                                gCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                    <  1:373711/149‑1 (MQ=255)
                                                                gCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGGCGCAGCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                    >  2:551824/1‑149 (MQ=255)
                                                                    gcGAGGGGC‑ACGGGACGGGAGGTCGGATAAGACGGATGCCGGCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                                                                <  2:40126/148‑1 (MQ=14)
                                                                         gTGC‑AGTGTACTGCAGGTCGGATAAGACGGATGGCGTCGTATTCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                           <  2:215620/149‑1 (MQ=11)
                                                                                tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCACCGCCTGATg                                                                     >  2:683764/1‑149 (MQ=255)
                                                                                      ggaGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                                                               <  1:156358/147‑1 (MQ=37)
                                                                                      gTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                                                               >  1:471058/1‑149 (MQ=38)
                                                                                          gccGGATAAGACGGAGGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAc                                                           <  1:187720/147‑1 (MQ=37)
                                                                                          gTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGCTGcgccggggag                                                           >  2:148599/1‑144 (MQ=255)
                                                                                                    ccgcacggCATCCCATCCGACAACCGATGCCTGACGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGCAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                 <  2:45042/143‑1 (MQ=11)
                                                                                                    aCGGATGGCGTTGCATCCGACAACCGATGCCTGGTGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                 <  2:464463/149‑1 (MQ=12)
                                                                                                      ggATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGTCGACGCGACTTATCAg                                               >  1:626269/1‑149 (MQ=255)
                                                                                                              tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgg                                       >  1:731791/1‑149 (MQ=32)
                                                                                                               cGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                      <  1:704176/149‑1 (MQ=32)
                                                                                                               cGCATCCGACAACCGATGCCAGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                      <  1:639133/149‑1 (MQ=17)
                                                                                                                  aTCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTg                                   <  2:721824/149‑1 (MQ=32)
                                                                                                                   tCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGc                                  <  1:296431/149‑1 (MQ=32)
                                                                                                                      gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt                               <  2:366203/149‑1 (MQ=21)
                                                                                                                      gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGAGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt                               >  1:503324/1‑149 (MQ=25)
                                                                                                                         aaCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGACCTGTAGGTCGGATAAGAAGCAAAGCGTCGCATCCGACAACCGCTGACAGGTGCGAGGCTGACGGGACTGATAGGGCATAAGGAGTGGCCTGaa                            >  2:121080/1‑149 (MQ=11)
                                                                                                                         aaCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGTCAACCGATGCCTGATGAGCCGCGGACGCGACTTATCAGGCACACGGGGTGCagcaca                            >  2:443928/1‑144 (MQ=11)
                                                                                                                         aaCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                            <  1:331767/149‑1 (MQ=35)
                                                                                                                         aaCCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                            <  1:207822/149‑1 (MQ=17)
                                                                                                                            cGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGGGCACTGAACtc                         >  2:253035/1‑148 (MQ=25)
                                                                                                                             gATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGt                        <  2:129072/149‑1 (MQ=17)
                                                                                                                              aTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTAATCAGGCCTACGGGGTGCACTGAGCGGTa                       >  1:564500/1‑149 (MQ=255)
                                                                                                                                  cTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCCTCCAGCAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGCGCCCTGAACTGTAGGTc                   >  1:536159/1‑149 (MQ=11)
                                                                                                                                   tGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATTAGGCCTACGGGGTGCACTGAACTGTAGGTCg                  >  1:179059/1‑149 (MQ=21)
                                                                                                                                   tGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                  <  2:658080/149‑1 (MQ=35)
                                                                                                                                   tGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                  <  1:725011/149‑1 (MQ=35)
                                                                                                                                      tccgccgcGGACGTGACTGACTGGGCCTACGGGGTGCACTGAACTGGAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt               <  2:13404/147‑1 (MQ=11)
                                                                                                                                                  gTGACTTATCAGGCCTACGGGGAGCACTGAACTGTAGGGCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCCCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGc   >  2:723389/1‑149 (MQ=16)
                                                                                                                                                   tGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg  >  1:306082/1‑149 (MQ=32)
                                                                                                                                                   |                                                                                                                                                    
AAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGC‑ACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCG  >  NC_000913/349571‑349865

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑