Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 374,295 C→T intergenic (+190/‑343) mhpE → / → mhpT 4‑hydroxy‑2‑ketovalerate aldolase/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913374,2950CT98.8% 888.5 / ‑2.0 322intergenic (+190/‑343)mhpE/mhpT4‑hydroxy‑2‑ketovalerate aldolase/putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base T (141/175);  total (146/176)
Fisher's exact test for biased strand distribution p-value = 4.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

Only 100 of 337 total aligned reads displayed.

GTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGG  >  NC_000913/374205‑374384
                                                                                          |                                                                                         
gTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGgtgc                                                                                         <  2:3048155/93‑1 (MQ=37)
gTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGgtgc                                                                                         <  2:1035916/93‑1 (MQ=255)
 tAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGcgcg                                                                                        >  1:5732409/1‑93 (MQ=37)
  aGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCgg                                                                                       <  1:7860924/93‑1 (MQ=255)
  aGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCgg                                                                                       <  2:7892213/93‑1 (MQ=37)
   ggTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCggg                                                                                      >  2:6427150/1‑93 (MQ=37)
     tCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGaa                                                                                    <  2:2328150/93‑1 (MQ=25)
        gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                                                                 <  2:7529987/93‑1 (MQ=18)
        gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                                                                 <  1:598540/93‑1 (MQ=18)
          tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                                                               <  1:579215/93‑1 (MQ=21)
          tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                                                               >  2:1959600/1‑93 (MQ=18)
           aaGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                                                              <  1:3326902/93‑1 (MQ=18)
              aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGc                                                                           <  1:2633794/93‑1 (MQ=18)
               cgtgTTAGCTTCGCATCAGGCAGCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGcc                                                                          <  2:746325/90‑1 (MQ=11)
                  gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                                                       >  1:5624759/1‑93 (MQ=18)
                  gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                                                       <  1:6750056/93‑1 (MQ=18)
                   ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                                                      >  1:7389143/1‑93 (MQ=18)
                   ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                                                      >  1:6258680/1‑93 (MQ=18)
                    tAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATa                                                                     <  1:4875180/93‑1 (MQ=18)
                    tAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCAGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATa                                                                     <  2:7283907/93‑1 (MQ=255)
                     aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                                                    <  2:4697046/93‑1 (MQ=18)
                     aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                                                    <  1:863540/93‑1 (MQ=18)
                     aGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                                                    <  2:6982737/93‑1 (MQ=18)
                      gCGTCGCATCAGGCGTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAg                                                                   >  2:6674648/1‑93 (MQ=255)
                      gCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAg                                                                   >  1:1699799/1‑93 (MQ=18)
                      gCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAg                                                                   >  1:1289900/1‑93 (MQ=18)
                        gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGc                                                                           >  2:7240445/1‑83 (MQ=18)
                         tCGCATCAGGCGTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                                <  1:4751172/93‑1 (MQ=255)
                           gCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                              <  1:3389504/93‑1 (MQ=18)
                              tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                           <  2:440169/93‑1 (MQ=18)
                              tCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                           <  1:1506281/93‑1 (MQ=18)
                              tCAGGCATCTGCGCACAACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                           >  2:1899571/1‑93 (MQ=255)
                               cAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAg                                                          <  2:3228336/93‑1 (MQ=18)
                               cAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAg                                                          >  2:5818024/1‑93 (MQ=33)
                                aGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                         <  1:5044798/93‑1 (MQ=18)
                                 ggCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                        >  1:420367/1‑93 (MQ=18)
                                  gCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                       <  1:6109991/93‑1 (MQ=18)
                                  gCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                       <  2:3430731/93‑1 (MQ=18)
                                   cgtCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTc                                                      <  1:4110186/91‑1 (MQ=16)
                                   cATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTc                                                      <  2:5287020/93‑1 (MQ=18)
                                   cATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTc                                                      <  1:6134519/93‑1 (MQ=18)
                                     tCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGc                                                    >  2:4598138/1‑93 (MQ=18)
                                     tCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGc                                                    >  1:1813119/1‑93 (MQ=18)
                                      cTGCGCACGACTGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                   <  2:5716047/93‑1 (MQ=16)
                                       tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt                                                  >  2:6564104/1‑93 (MQ=18)
                                         cgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                <  1:3081305/93‑1 (MQ=18)
                                           cACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAgg                                              <  2:924343/93‑1 (MQ=18)
                                            aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                             <  1:2219050/93‑1 (MQ=18)
                                            aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                             >  1:6845139/1‑93 (MQ=21)
                                            aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                             <  2:5498625/93‑1 (MQ=18)
                                             cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCa                                            <  2:6131636/93‑1 (MQ=18)
                                             cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCa                                            >  1:119367/1‑93 (MQ=18)
                                             cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCa                                            <  1:4566587/93‑1 (MQ=18)
                                             cGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCa                                            <  1:7329017/93‑1 (MQ=18)
                                              gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                           >  1:1065351/1‑93 (MQ=18)
                                               aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                          <  2:513888/93‑1 (MQ=18)
                                                 tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTg                                        >  1:4395657/1‑93 (MQ=18)
                                                  gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgc                                       >  1:6671024/1‑93 (MQ=18)
                                                  gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGGGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgc                                       >  2:7330621/1‑93 (MQ=18)
                                                    cGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTgcgc                                     >  1:5209527/1‑93 (MQ=25)
                                                     ggATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCa                                    <  2:6888705/93‑1 (MQ=37)
                                                        tGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGa                                 <  2:1666397/93‑1 (MQ=37)
                                                        tGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGa                                 <  2:3998456/93‑1 (MQ=37)
                                                         gcggcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAc                                <  1:7570218/93‑1 (MQ=37)
                                                         gcggcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAc                                <  1:7901459/93‑1 (MQ=37)
                                                         gcggcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAc                                <  2:3167017/93‑1 (MQ=37)
                                                            gcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGTATCAGGCATCTGCGCACGACTGc                             >  2:4304485/1‑93 (MQ=255)
                                                            gcgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGc                             >  1:2298187/1‑93 (MQ=37)
                                                             cgTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGcc                            <  2:1855711/93‑1 (MQ=37)
                                                               tGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCgg                          >  2:4223378/1‑93 (MQ=25)
                                                                gAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGa                         <  2:5319104/93‑1 (MQ=25)
                                                                   cGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgc                      <  2:3120477/93‑1 (MQ=18)
                                                                   cGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgc                      <  2:1228883/93‑1 (MQ=21)
                                                                    gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcg                     >  1:4451053/1‑93 (MQ=18)
                                                                    cccTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGTCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcg                     <  2:2116787/92‑1 (MQ=255)
                                                                     ccTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcgg                    >  1:745554/1‑93 (MQ=21)
                                                                     ccTTATCCGGCCTACGGATGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCTAGAGTCGCATCAGGCATCTGCGCACGACTGCCGGATgccg                    >  2:1939858/1‑91 (MQ=11)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGa                         >  1:689681/1‑86 (MQ=18)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcggcg                  >  2:1123036/1‑93 (MQ=37)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgcggcg                  >  2:3615528/1‑93 (MQ=37)
                                                                          tCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGa               >  1:2801853/1‑93 (MQ=255)
                                                                          tCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGa               >  1:3765038/1‑93 (MQ=255)
                                                                          tCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGa               >  1:7102025/1‑93 (MQ=255)
                                                                           ccGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGaa              <  2:2600411/93‑1 (MQ=255)
                                                                           ccGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGaa              <  1:507421/93‑1 (MQ=255)
                                                                            cGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAc             >  2:7728107/1‑93 (MQ=255)
                                                                            cGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAc             <  2:6346084/93‑1 (MQ=255)
                                                                              gCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGc           <  2:1752384/93‑1 (MQ=255)
                                                                               ccTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc          <  1:2931275/93‑1 (MQ=255)
                                                                               ccTACGGATGGTGCGGGAATTTGTAGGCCTGAGAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc          <  1:3913770/93‑1 (MQ=255)
                                                                                cTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCt         <  1:351309/93‑1 (MQ=255)
                                                                                 tACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCtt        <  1:6544725/93‑1 (MQ=255)
                                                                                  aCGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTa       >  2:489480/1‑93 (MQ=255)
                                                                                   cGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAt      <  2:222076/93‑1 (MQ=255)
                                                                                   cGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAt      <  1:1976106/93‑1 (MQ=255)
                                                                                    ggATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATc     <  1:7662126/93‑1 (MQ=255)
                                                                                    ggATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATc     <  1:5406324/93‑1 (MQ=255)
                                                                                     gATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATcc    <  1:6915850/93‑1 (MQ=255)
                                                                                     gATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATcc    >  2:4818448/1‑93 (MQ=255)
                                                                                     gATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATcc    >  1:215615/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGG  >  NC_000913/374205‑374384

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: