Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 373,351 T→C L87P (CTG→CCG)  mhpE → 4‑hydroxy‑2‑ketovalerate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913373,3510TC92.0% 740.5 / 45.6 276L87P (CTG→CCG) mhpE4‑hydroxy‑2‑ketovalerate aldolase
Reads supporting (aligned to +/- strand):  ref base T (9/13);  new base C (113/141);  total (122/154)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

Only 100 of 303 total aligned reads displayed.

CTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGAAGCTG  >  NC_000913/373259‑373440
                                                                                            |                                                                                         
  ttAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTtgccgc                                                                                         <  1:1171951/93‑1 (MQ=255)
   tAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTgccgcc                                                                                        >  1:1498690/1‑93 (MQ=255)
    aaCTATGGTTTCGGCGCACATAGCGACTTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTGCCa                                                                                       <  1:4304428/93‑1 (MQ=255)
    aaCTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCa                                                                                       <  1:591784/93‑1 (MQ=255)
    aaCTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCa                                                                                       <  1:8569755/93‑1 (MQ=255)
     aCTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAg                                                                                      <  1:3798980/93‑1 (MQ=255)
       tATGGTTTCGGCGCACATAGCGACCTTGAATGGAATGAAGCGGCGGCGGATGTGGAGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGa                                                                                    >  2:8655473/1‑93 (MQ=255)
        aTGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGaa                                                                                   <  2:5344670/93‑1 (MQ=255)
         tGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAAt                                                                                  >  2:7678862/1‑93 (MQ=255)
         tGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAAt                                                                                  >  2:2957752/1‑93 (MQ=255)
           gTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGGAtgg                                                                                >  1:9163600/1‑91 (MQ=255)
            tttCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCgg                                                                               <  2:8787590/93‑1 (MQ=255)
             ttCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCTAAAATCGCGACGTTGTTGCCGCCAGGAATCGGc                                                                              <  2:671520/93‑1 (MQ=255)
             ttCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGc                                                                              <  1:3014644/93‑1 (MQ=255)
              tCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGcc                                                                             >  2:3379693/1‑92 (MQ=255)
              tCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCa                                                                             <  1:7792154/93‑1 (MQ=255)
                 gcgcACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTa                                                                          >  2:889102/1‑93 (MQ=255)
                  cgcACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTAt                                                                         <  2:3115812/93‑1 (MQ=255)
                  cgcACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTAt                                                                         >  2:5008683/1‑93 (MQ=255)
                   gcACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTACCAGGAATCGGCACTAtt                                                                        >  1:1366404/1‑93 (MQ=255)
                   gcACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTAtt                                                                        <  2:1305604/93‑1 (MQ=255)
                    cacaTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCGAAAATCGGGAGTGTGTTGCCGCCAGGAATGGGCCTTATTc                                                                       >  1:5371239/1‑93 (MQ=255)
                    cacaTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTGCCAGGAATCGGCACTATTc                                                                       <  2:5212705/93‑1 (MQ=255)
                    cacaTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTc                                                                       <  2:6388854/93‑1 (MQ=255)
                     ccaTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCa                                                                      <  2:5662571/92‑1 (MQ=255)
                     acaTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCTACGTTGTTGCCGCCAGGAATCGGCACTATTCa                                                                      <  1:1863556/93‑1 (MQ=255)
                     acaTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCa                                                                      <  2:1809915/93‑1 (MQ=255)
                      caTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCAc                                                                     >  2:1170794/1‑93 (MQ=255)
                       aTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACg                                                                    >  2:6551789/1‑93 (MQ=255)
                       aTAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACg                                                                    >  1:6630518/1‑93 (MQ=255)
                        tAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGa                                                                   >  2:8483654/1‑93 (MQ=255)
                         aGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGAt                                                                  <  1:9959501/93‑1 (MQ=255)
                         aGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGAt                                                                  >  1:2585263/1‑93 (MQ=255)
                          gCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATc                                                                 <  1:92450/93‑1 (MQ=255)
                          gCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATc                                                                 >  2:4503624/1‑93 (MQ=255)
                            gACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTg                                                               <  1:4484411/93‑1 (MQ=255)
                            gACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTg                                                               >  1:9658131/1‑93 (MQ=255)
                              ccTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGaa                                                             <  1:3973430/93‑1 (MQ=255)
                              ccTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGaa                                                             >  2:9681182/1‑93 (MQ=255)
                               cTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGaaa                                                            >  1:8940125/1‑93 (MQ=255)
                                ttGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGaaaa                                                           >  2:3436547/1‑93 (MQ=255)
                                 tGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTGCCAGGAATCGGCACTATTCACGATCTGaaaaa                                                          <  1:3568857/93‑1 (MQ=255)
                                 tGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGaaaaa                                                          <  2:3181198/93‑1 (MQ=255)
                                  gAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTGCCAGGAATCGGCACTATTCACGATCTGAAAAAt                                                         >  2:5110494/1‑93 (MQ=255)
                                  gAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAAt                                                         >  2:1549601/1‑93 (MQ=255)
                                   aaTGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATg                                                        >  1:3426648/1‑93 (MQ=255)
                                     tGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCCAAATCGCGACGTTGTTGCTGCCAGGAATCGGCACTATTCACGATCTGAAAAATGcc                                                      <  1:1072782/93‑1 (MQ=255)
                                     tGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGcc                                                      <  2:1593348/93‑1 (MQ=255)
                                      ggATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCt                                                     <  1:6050477/93‑1 (MQ=255)
                                      ggATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCt                                                     <  2:3673155/93‑1 (MQ=255)
                                       gATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAACTCGCGACGTTGTTGCCGCCAGGATTCGGCACTATTCACGATCTGAAAAATGCCTg                                                    <  2:257620/93‑1 (MQ=255)
                                       gATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTg                                                    <  2:7291580/93‑1 (MQ=255)
                                       gATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTg                                                    >  2:6189311/1‑93 (MQ=255)
                                       gATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTg                                                    <  1:5164368/93‑1 (MQ=255)
                                         ttGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGc                                                  <  1:2275439/93‑1 (MQ=255)
                                          tGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCa                                                 <  2:9861911/93‑1 (MQ=255)
                                           taaGCGGCCGCAGATGTTGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAg                                                <  1:8110217/92‑1 (MQ=255)
                                           gAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAg                                                >  2:5538779/1‑93 (MQ=255)
                                            aaGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAgg                                               >  1:363225/1‑93 (MQ=255)
                                             aGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGc                                              <  2:4939818/93‑1 (MQ=255)
                                              gcggcggcgGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCt                                             <  2:2917706/93‑1 (MQ=255)
                                                 gcggcggATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGgc                                          >  2:4584031/1‑93 (MQ=255)
                                                    gcggATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGgcgcg                                       >  1:2119287/1‑93 (MQ=255)
                                                      ggATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGcgcgcg                                     <  2:2755963/93‑1 (MQ=255)
                                                       gATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCgg                                    >  2:9341639/1‑93 (MQ=255)
                                                       gATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCgg                                    <  2:6958427/93‑1 (MQ=255)
                                                       gATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCgg                                    >  1:5603986/1‑93 (MQ=255)
                                                       gATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCgg                                    >  2:5896155/1‑93 (MQ=255)
                                                         tgtgGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCgggg                                  >  1:8287001/1‑92 (MQ=255)
                                                         tgtgGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGggt                                  <  1:5004277/93‑1 (MQ=255)
                                                         tgtgGTGAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGggt                                  <  1:6853918/93‑1 (MQ=255)
                                                          gtggtgAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGggtg                                 <  1:7943746/93‑1 (MQ=255)
                                                          gtggtgAAGCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGggtg                                 <  1:5213944/93‑1 (MQ=255)
                                                                  gCACGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCgtgt                         <  2:7358576/93‑1 (MQ=255)
                                                                     cGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGc                      >  1:4242899/1‑93 (MQ=255)
                                                                     cGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGc                      >  2:2643400/1‑93 (MQ=255)
                                                                     cGCCAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGc                      <  1:1712967/93‑1 (MQ=255)
                                                                       ccAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCaa                    <  1:6191552/93‑1 (MQ=255)
                                                                       ccAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCaa                    <  1:391089/93‑1 (MQ=255)
                                                                       ccAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCaa                    <  1:5781787/93‑1 (MQ=255)
                                                                       ccAAAATCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCaa                    >  1:735892/1‑93 (MQ=255)
                                                                        ctaaaTCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAAc                   <  1:9481301/91‑1 (MQ=255)
                                                                          aaaTCGCGACGTTGTTGCCGCTAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCGACGc                 >  2:2199522/1‑93 (MQ=255)
                                                                           aaTCGCGACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCa                >  1:7694604/1‑93 (MQ=255)
                                                                            aTCGCGCCGTTCTTGCCGCCAGAAATCGGCACTATTCGCGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGTGTGGTTCGTGTGGCAACGCAc               <  2:353001/93‑1 (MQ=255)
                                                                               gcgACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGt            >  2:8118556/1‑93 (MQ=255)
                                                                                cgACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTa           <  1:2442089/93‑1 (MQ=255)
                                                                                 gACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTAc          <  1:8347112/93‑1 (MQ=255)
                                                                                 gACGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTAc          <  2:210427/93‑1 (MQ=255)
                                                                                  aCGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTAcc         >  2:392467/1‑93 (MQ=255)
                                                                                   cGTTGTTGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCg        >  1:3402019/1‑93 (MQ=255)
                                                                                    gttgttGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGa       >  2:4289888/1‑93 (MQ=255)
                                                                                    gttgttGCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGa       >  2:9581362/1‑93 (MQ=255)
                                                                                     ttgttgCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGaa      >  1:3755501/1‑93 (MQ=255)
                                                                                      tgttgCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGAAg     <  2:7120562/93‑1 (MQ=255)
                                                                                      tgttgCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGAAg     <  2:7905538/93‑1 (MQ=255)
                                                                                      tgttgCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGAAg     <  1:3247467/93‑1 (MQ=255)
                                                                                      tgttgCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGAAg     >  2:4776531/1‑93 (MQ=255)
                                                                                      tgttgCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGAAg     >  2:4069782/1‑93 (MQ=255)
                                                                                        ttgCCGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGAAGCt   >  1:7452664/1‑93 (MQ=255)
                                                                                            |                                                                                         
CTTTAACTATGGTTTCGGCGCACATAGCGACCTTGAATGGATTGAAGCGGCGGCGGATGTGGTGAAGCACGCCAAAATCGCGACGTTGTTGCTGCCAGGAATCGGCACTATTCACGATCTGAAAAATGCCTGGCAGGCTGGCGCGCGGGTGGTTCGTGTGGCAACGCACTGTACCGAAGCTG  >  NC_000913/373259‑373440

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