Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,333,842 C→T G98G (GGG→GGA murA ← UDP‑N‑glucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,8420CT99.5% 681.0 / NA 192G98G (GGG→GGAmurAUDP‑N‑glucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (93/98);  minor base A (1/0);  total (94/98)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 199 total aligned reads displayed.

GCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAG  >  NC_000913/3333752‑3333930
                                                                                          |                                                                                        
gCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGtcc                                                                                        >  1:3483237/1‑93 (MQ=255)
  cAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCa                                                                                      >  1:1324946/1‑93 (MQ=255)
   aGAAATGTGTAGATCAACCGGACGGGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAg                                                                                     <  2:2537908/93‑1 (MQ=255)
   aGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAg                                                                                     <  1:2304741/93‑1 (MQ=255)
       aTGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGcc                                                                                 <  2:2293239/93‑1 (MQ=255)
         gtgtAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCa                                                                               <  1:2186636/93‑1 (MQ=255)
         gtgtAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCa                                                                               >  2:4953803/1‑93 (MQ=255)
           gtAGATCAGCCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAATGCGCGCTACCAGCGGTCCCAGCGCCCAGa                                                                             <  1:4173436/93‑1 (MQ=255)
           gtAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGa                                                                             <  1:2450969/93‑1 (MQ=255)
           gtAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGa                                                                             <  1:1061637/93‑1 (MQ=255)
            tAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCACCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGAt                                                                            <  2:4739202/93‑1 (MQ=255)
              gagcAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAg                                                                          <  2:2337531/90‑1 (MQ=255)
              gATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAg                                                                          >  1:1168948/1‑93 (MQ=255)
              gATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAg                                                                          >  2:5957271/1‑93 (MQ=255)
               aTCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGa                                                                         >  2:3925520/1‑93 (MQ=255)
               aTCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGa                                                                         <  1:5009833/93‑1 (MQ=255)
               aTCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGa                                                                         <  2:922931/93‑1 (MQ=255)
                 cAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAg                                                                       <  1:3931773/93‑1 (MQ=255)
                 cAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAg                                                                       >  2:5941732/1‑93 (MQ=255)
                  aaCCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGc                                                                      <  2:3503322/93‑1 (MQ=255)
                    ccGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCAc                                                                    >  1:4693633/1‑93 (MQ=255)
                      ggACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGc                                                                  <  1:5377492/93‑1 (MQ=255)
                       gACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGCAGCACGCa                                                                 >  1:4833838/1‑93 (MQ=255)
                       gACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCa                                                                 >  1:5499116/1‑93 (MQ=255)
                         cgcgCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATg                                                               >  1:3392057/1‑93 (MQ=255)
                         cgcgCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATg                                                               <  2:5653266/93‑1 (MQ=255)
                           cgcACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGt                                                             <  2:2193800/93‑1 (MQ=255)
                           cgcACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGt                                                             <  2:550592/93‑1 (MQ=255)
                              aCCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGtttt                                                          <  1:5400264/93‑1 (MQ=255)
                              aCCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGcttt                                                          >  1:3624343/1‑93 (MQ=255)
                               ccGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTa                                                         >  1:3313028/1‑93 (MQ=255)
                                cGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTaa                                                        >  2:1212041/1‑93 (MQ=255)
                                cGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTaa                                                        >  1:4383773/1‑93 (MQ=255)
                                cGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTaa                                                        <  1:2502805/93‑1 (MQ=255)
                                cGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTaa                                                        >  1:1670380/1‑93 (MQ=255)
                                 gATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAAc                                                       >  1:360307/1‑93 (MQ=255)
                                  aTCGTACAACCGCCAGGGAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAAca                                                      >  1:1969780/1‑92 (MQ=255)
                                   tCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTTACCa                                                     >  2:30992/1‑93 (MQ=255)
                                      tACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGAt                                                  >  1:5583409/1‑93 (MQ=255)
                                         aaCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGt                                               <  1:259704/93‑1 (MQ=255)
                                         aaCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGt                                               >  1:2930770/1‑93 (MQ=255)
                                          aCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTa                                              >  1:1350603/1‑93 (MQ=255)
                                           ccGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTaa                                             <  2:2627973/93‑1 (MQ=255)
                                           ccGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTaa                                             <  1:1140708/93‑1 (MQ=255)
                                            cGCCAGGTAGTGAAACTTGCGCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAg                                            <  2:4827455/93‑1 (MQ=255)
                                            cGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAg                                            <  1:4020975/93‑1 (MQ=255)
                                             gCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAgg                                           <  1:5894747/93‑1 (MQ=255)
                                              ccAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGt                                          <  1:5784706/93‑1 (MQ=255)
                                              ccAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGt                                          <  1:5835299/93‑1 (MQ=255)
                                                cggTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTgc                                        <  1:3075883/92‑1 (MQ=255)
                                                aGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTgc                                        >  2:619269/1‑93 (MQ=255)
                                                 ggTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTgcg                                       >  2:1651638/1‑93 (MQ=255)
                                                 ggTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTgcg                                       >  1:3166900/1‑93 (MQ=255)
                                                  gTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTgcgc                                      <  2:6252573/93‑1 (MQ=255)
                                                    aGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAg                                    >  1:881784/1‑93 (MQ=255)
                                                      tGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGaa                                  >  1:3498078/1‑93 (MQ=255)
                                                      tGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGaa                                  <  1:2159374/93‑1 (MQ=255)
                                                       gAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAAt                                 <  1:5230580/93‑1 (MQ=255)
                                                        aaaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATa                                <  1:2401782/93‑1 (MQ=255)
                                                         aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAc                               >  1:1026192/1‑93 (MQ=255)
                                                         aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGAAGAATAc                               >  1:3027370/1‑93 (MQ=255)
                                                          aCTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACa                              <  2:5116637/93‑1 (MQ=255)
                                                          aCTTGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACa                              <  1:1860109/93‑1 (MQ=255)
                                                            ttGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACAtt                            <  2:2850365/93‑1 (MQ=255)
                                                            ttGCCCCTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACAtt                            <  2:4079195/93‑1 (MQ=255)
                                                               ccccTGCCCAAAGCGCGCCACCAGCAGTCCCAGCGCCAAGATAGAAGCACGCATGGTTTTACCCAGATAGAAAGGGACGCAGAACACATAAAc                         <  2:2328625/93‑1 (MQ=255)
                                                               ccccTGACCAAAGCGCGCTCCCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAAc                         <  2:621034/93‑1 (MQ=255)
                                                               ccccTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAAc                         <  1:6254175/93‑1 (MQ=255)
                                                                cccTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAAct                        >  2:4553960/1‑92 (MQ=255)
                                                                cccTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACg                        >  2:2498784/1‑93 (MQ=255)
                                                                cccTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACg                        >  2:5712653/1‑93 (MQ=255)
                                                                 ccTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGt                       >  1:2720451/1‑93 (MQ=255)
                                                                 ccTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGt                       >  2:5315417/1‑93 (MQ=255)
                                                                 ccTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGAGCGCAGAATACATTAACGt                       >  1:3813791/1‑93 (MQ=255)
                                                                  cTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTc                      >  1:4058567/1‑93 (MQ=255)
                                                                  cTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTc                      >  1:4159410/1‑93 (MQ=255)
                                                                  cTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTc                      >  2:4217715/1‑93 (MQ=255)
                                                                  cTGACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGCTTTAACCAGATCGTAAGGTGCGCAGACTACATTAACCTc                      >  2:2053068/1‑93 (MQ=255)
                                                                  cTGACCAAAGCCCCCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTc                      <  2:5798657/93‑1 (MQ=255)
                                                                    gACCAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcgc                    <  1:2360712/93‑1 (MQ=255)
                                                                      ccAAAGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCgg                  <  1:5663528/93‑1 (MQ=255)
                                                                         aaGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCa               >  1:2629841/1‑93 (MQ=255)
                                                                         aaGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCa               <  1:210046/93‑1 (MQ=255)
                                                                          aGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCAt              <  1:4214525/93‑1 (MQ=255)
                                                                          aGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCAt              <  2:385110/93‑1 (MQ=255)
                                                                          aGCGCGCTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCAt              >  2:2984860/1‑93 (MQ=255)
                                                                           gcgcgcTACCAGCGGTCCCAGCGCCCAGATAGAAGAACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATc             >  2:5415054/1‑93 (MQ=255)
                                                                           gcgcgcTACCAGCGGTCCCAGCGCACAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATc             <  2:1343138/93‑1 (MQ=255)
                                                                             gcgcTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCaa           <  1:4490935/93‑1 (MQ=255)
                                                                               gcTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAata         <  1:1458251/93‑1 (MQ=255)
                                                                               gcTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAata         >  1:4642171/1‑93 (MQ=255)
                                                                                cTACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAatat        <  2:3139929/93‑1 (MQ=255)
                                                                                 tACCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATg       <  2:4396802/93‑1 (MQ=255)
                                                                                  aCCAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGc      >  2:4239198/1‑93 (MQ=255)
                                                                                   ccAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGca     >  2:3012480/1‑93 (MQ=255)
                                                                                    cAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcac    <  2:3963602/93‑1 (MQ=255)
                                                                                    cAGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcac    <  1:1942664/93‑1 (MQ=255)
                                                                                     aGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcaca   >  1:1584577/1‑93 (MQ=255)
                                                                                     aGCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcaca   >  2:1810450/1‑93 (MQ=255)
                                                                                      gCGGTCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAg  >  1:6076993/1‑93 (MQ=255)
                                                                                          |                                                                                        
GCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAG  >  NC_000913/3333752‑3333930

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