Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 444,641 C→T intergenic (‑126/‑42) panE ← / → yajQ 2‑dehydropantoate reductase, NADPH‑specific/phage Phi6 host factor, ATP/GTP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913444,6410CT100.0% 243.9 / NA 77intergenic (‑126/‑42)panE/yajQ2‑dehydropantoate reductase, NADPH‑specific/phage Phi6 host factor, ATP/GTP binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/33);  total (44/33)

TAATTTTCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGACT  >  NC_000913/444506‑444776
                                                                                                                                       |                                                                                                                                       
tAATTTTCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTggg                                                                                                                                      >  2:684820/1‑139 (MQ=255)
              cctcctCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCa                                                                                                                        >  1:419535/1‑139 (MQ=255)
               ctcctcACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCAt                                                                                                                       <  2:419535/139‑1 (MQ=255)
                      cccGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAAt                                                                                                                >  1:1177950/1‑139 (MQ=255)
                         gCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                             <  1:684820/139‑1 (MQ=255)
                               aCTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGgaga                                                                                                       >  1:229004/1‑139 (MQ=255)
                                 tCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAgg                                                                                                     >  2:469069/1‑139 (MQ=255)
                                            aaTAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          <  1:469069/139‑1 (MQ=255)
                                              tAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                        <  2:410449/139‑1 (MQ=255)
                                                        ccTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                        >  2:98356/1‑129 (MQ=255)
                                                        ccTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATc                                                                                        <  1:98356/129‑1 (MQ=255)
                                                            aaTTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     <  2:576935/96‑1 (MQ=255)
                                                            aaTTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     >  1:576935/1‑96 (MQ=255)
                                                                gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          <  2:450299/119‑1 (MQ=255)
                                                                gtgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          >  1:450299/1‑119 (MQ=255)
                                                                gtacctGCTGTGGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGa                                                                      >  2:781613/6‑139 (MQ=255)
                                                                 tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     >  2:871107/1‑91 (MQ=255)
                                                                 tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCa                                                                                                                     <  1:871107/91‑1 (MQ=255)
                                                                 tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGaa                                                                                                   >  2:644878/1‑109 (MQ=255)
                                                                 tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGaa                                                                                                   <  1:644878/109‑1 (MQ=255)
                                                                 tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGaa                                                                     >  2:475761/1‑139 (MQ=255)
                                                                 tgCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGaa                                                                     >  2:259487/1‑139 (MQ=255)
                                                                  gCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAg                                                                    >  2:905020/1‑139 (MQ=255)
                                                                   cATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGt                                                                   >  2:1275813/1‑139 (MQ=255)
                                                                      gCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                             >  2:114602/1‑94 (MQ=255)
                                                                      gCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                             <  1:114602/94‑1 (MQ=255)
                                                                                cTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGaa                                                      <  1:905020/139‑1 (MQ=255)
                                                                                  aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATg                                                                                                                               <  2:992974/64‑1 (MQ=255)
                                                                                  aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATg                                                                                                                               >  1:992974/1‑64 (MQ=255)
                                                                                  aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                           >  1:1006514/1‑100 (MQ=255)
                                                                                  aTGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                           <  2:1006514/100‑1 (MQ=255)
                                                                                    gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGa                                                                                                              <  2:571634/79‑1 (MQ=255)
                                                                                    gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGa                                                                                                              >  1:571634/1‑79 (MQ=255)
                                                                                    gCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCAc                                                  >  1:1319557/1‑139 (MQ=255)
                                                                                         aTTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAAc                                             <  2:1177950/139‑1 (MQ=255)
                                                                                              gggATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGt                                        >  2:1001684/1‑139 (MQ=255)
                                                                                               ggATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGaa                                                      >  1:213530/1‑124 (MQ=255)
                                                                                               ggATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGaa                                                      <  2:213530/124‑1 (MQ=255)
                                                                                                  tACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGAt                                    >  1:483194/1‑139 (MQ=255)
                                                                                                    cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                  >  2:314054/1‑139 (MQ=255)
                                                                                                    cGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATaa                                  >  2:737568/1‑139 (MQ=255)
                                                                                                      tCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCa                                                   <  1:1310657/120‑1 (MQ=255)
                                                                                                      tCCTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCa                                                   >  2:1310657/1‑120 (MQ=255)
                                                                                                        cTGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAAcgcg                              <  1:314054/139‑1 (MQ=255)
                                                                                                         tGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                                              <  2:1207546/90‑1 (MQ=255)
                                                                                                         tGGTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                                              >  1:1207546/1‑90 (MQ=255)
                                                                                                          ggTGCAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAg                            <  2:1319557/139‑1 (MQ=255)
                                                                                                              cAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            >  2:822507/1‑71 (MQ=255)
                                                                                                              cAGGACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            <  1:822507/71‑1 (MQ=255)
                                                                                                                ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAgg                                                        >  2:710826/1‑105 (MQ=255)
                                                                                                                ggACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAgg                                                        <  1:710826/105‑1 (MQ=255)
                                                                                                                 gACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAACAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCcgcg                       <  1:661400/137‑1 (MQ=255)
                                                                                                                 gACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                                              >  1:939588/1‑82 (MQ=255)
                                                                                                                 gACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATAtt                                                                              <  2:939588/82‑1 (MQ=255)
                                                                                                                 gACTGTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCcgcg                       >  2:661400/1‑137 (MQ=255)
                                                                                                                     gTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                           >  2:316686/1‑65 (MQ=255)
                                                                                                                     gTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGcc                                                                                           <  1:316686/65‑1 (MQ=255)
                                                                                                                     gTCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTgg                 >  1:696709/1‑139 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                             >  1:1150555/1‑46 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGaa                                                                                                             <  2:1150555/46‑1 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            <  1:1194857/63‑1 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGc                                                                                            >  2:1194857/1‑63 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          <  2:570540/65‑1 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCa                                                                                          >  1:570540/1‑65 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAAcgc                                           >  1:1179089/1‑112 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAAcgc                                           <  2:1179089/112‑1 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCcgcg                       >  1:732390/1‑132 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCcgcg                       <  2:732390/132‑1 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGa                >  2:417666/1‑139 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGa                >  1:412206/1‑139 (MQ=255)
                                                                                                                      tCGGTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAAGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGa                >  2:1082289/1‑139 (MQ=255)
                                                                                                                         gTTATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTc             <  2:229004/139‑1 (MQ=255)
                                                                                                                           tATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAATATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTc                                                           <  2:802219/91‑1 (MQ=255)
                                                                                                                           tATTTAACTTTGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTc                                                           >  1:802219/1‑91 (MQ=255)
                                                                                                                                    ttGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGAc   <  1:503567/138‑1 (MQ=255)
                                                                                                                                    ttGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGAc   >  2:503567/1‑138 (MQ=255)
                                                                                                                                    ttGTGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGACt  >  2:431425/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TAATTTTCATTGTTCCTCCTCACCCGCAACCACTCCGGGTGTTCAATAAGGCTATCCCTTAATTGTGCATGCTGTTGCGACTATGCACAATTAAGGGATACGTCCTGGTGCAGGACTGTCGGTTATTTAACTTTGCGGGTATTATGCTTCGCATCAAAAATGAAGGGAGAGGAAAAGATGCCATCTTTCGATATTGTCTCTGAAGTTGATCTTCAGGAAGCACGTAACGCGGTCGATAACGCGAGCCGCGAAGTGGAGTCCCGTTTTGACT  >  NC_000913/444506‑444776

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: