Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,574,513 T→C intergenic (+1216/+184) yhgN → / ← gntK hypothetical protein/gluconokinase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,574,5130TC99.0% 606.1 / ‑3.9 199intergenic (+1216/+184)yhgN/gntKhypothetical protein/gluconokinase 2
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (99/98);  total (100/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 210 total aligned reads displayed.

CGGCAATCTGATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGTGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTAAAG  >  NC_000913/3574426‑3574602
                                                                                       |                                                                                         
cGGCAATCTGATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTc                                                                                      >  2:1399395/1‑93 (MQ=255)
 ggCAATCTGATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTcc                                                                                     >  1:808825/1‑93 (MQ=255)
  gCAATCTGATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTccc                                                                                    >  1:4106937/1‑93 (MQ=255)
  gCAATCTGATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTccc                                                                                    <  1:363208/93‑1 (MQ=255)
        tGATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCt                                                                              >  2:1008734/1‑93 (MQ=255)
        tGATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCt                                                                              >  1:2430700/1‑93 (MQ=255)
         gATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCtt                                                                             <  1:2813300/93‑1 (MQ=255)
          aTCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGCCAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCttt                                                                            <  1:3113560/93‑1 (MQ=255)
          aTCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCttt                                                                            <  1:3702445/93‑1 (MQ=255)
            cGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCttttt                                                                          >  2:881948/1‑93 (MQ=255)
                tGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCAt                                                                      <  1:1813396/93‑1 (MQ=255)
                 tcgcCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCAtc                                                                     <  1:2856123/92‑1 (MQ=255)
                 tcgcCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCAtc                                                                     <  1:3609621/92‑1 (MQ=255)
                  cgcCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCAtcg                                                                    >  2:3297232/1‑93 (MQ=255)
                  cgcCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCAtcg                                                                    <  1:123333/93‑1 (MQ=255)
                    ctggTTTCATCTAAGATCTTGTCAAACATAGCTCCCAGGGCGACAACCATAGCCAGAGAGCCGCGGGCGCCTTCCATCCCTTTCTCCAtcgtc                                                                  >  2:722586/3‑93 (MQ=255)
                    ccGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCAtcgtc                                                                  >  2:853883/1‑93 (MQ=255)
                      ggTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGc                                                                <  1:963581/93‑1 (MQ=255)
                       gTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCt                                                               >  1:4639952/1‑93 (MQ=255)
                       gTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCt                                                               <  2:1385363/93‑1 (MQ=255)
                           cATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGa                                                           <  2:2410842/93‑1 (MQ=255)
                             tGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGAtt                                                         <  2:1278643/93‑1 (MQ=255)
                             tGTAAGATCTTGCCAAACATAGCGCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGAtt                                                         <  2:180503/93‑1 (MQ=255)
                              gTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATAGTCGCTGCGAttt                                                        >  2:3894735/1‑93 (MQ=255)
                               tAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGAtttt                                                       >  2:3666353/1‑93 (MQ=255)
                                aaGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTa                                                      <  1:1375986/93‑1 (MQ=255)
                                  gATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATc                                                    >  2:1321957/1‑93 (MQ=255)
                                   aTCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCg                                                   <  1:1003739/93‑1 (MQ=255)
                                   aTCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCg                                                   <  2:2838285/93‑1 (MQ=255)
                                   aTCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCg                                                   >  2:3823573/1‑93 (MQ=255)
                                    tCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGa                                                  >  1:2204537/1‑93 (MQ=255)
                                    tCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGa                                                  >  1:1459148/1‑93 (MQ=255)
                                     cTTTCCAAACATCGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAg                                                 <  2:518948/93‑1 (MQ=255)
                                      ttGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGc                                                <  2:3724263/93‑1 (MQ=255)
                                          cAAGCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCCCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCa                                            <  2:1675612/93‑1 (MQ=255)
                                          cAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCa                                            <  2:2929295/93‑1 (MQ=255)
                                           aaaCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCAt                                           <  1:4075025/93‑1 (MQ=255)
                                           aaaCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCAt                                           <  1:3045326/93‑1 (MQ=255)
                                            aaCCTAGCTCCCAGGCCAACAACCGCCGCAAGGAAGCCGAGGGCGCCTCCCATCCCTTTTCCCATCGTCGCTGCGATTTTATCGAGCGGCATa                                          <  2:537387/93‑1 (MQ=255)
                                            aaCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATa                                          <  2:122124/93‑1 (MQ=255)
                                            aaCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATa                                          >  2:3266043/1‑93 (MQ=255)
                                            aaCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATa                                          >  1:1154614/1‑93 (MQ=255)
                                             aCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATAc                                         >  1:932462/1‑93 (MQ=255)
                                             aCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATAc                                         <  2:320725/93‑1 (MQ=255)
                                             aCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATAc                                         <  2:2026527/93‑1 (MQ=255)
                                             aCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATAc                                         <  1:1259353/93‑1 (MQ=255)
                                             aCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATAc                                         >  1:191910/1‑93 (MQ=255)
                                             aCATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATAc                                         <  1:3750168/93‑1 (MQ=255)
                                              cATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATAcc                                        >  1:564117/1‑93 (MQ=255)
                                              cATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATAcc                                        >  1:1662225/1‑93 (MQ=255)
                                               aTAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCa                                       <  2:2733694/93‑1 (MQ=255)
                                                 aGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGa                                     <  1:2899629/93‑1 (MQ=255)
                                                   cTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGaaa                                   >  1:3222207/1‑93 (MQ=255)
                                                    tCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGaaaa                                  <  1:1289250/93‑1 (MQ=255)
                                                    tCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGAGCCTCCCATCCCTTTTTCCATCATCGCTGCGATTTTCTTGAGCGGAATACCAGaaaa                                  >  2:3596392/1‑93 (MQ=255)
                                                     cccAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGaaaaa                                 <  2:1992629/93‑1 (MQ=255)
                                                       cAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAg                               <  1:4291737/93‑1 (MQ=255)
                                                       cAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAg                               >  1:1572984/1‑93 (MQ=255)
                                                       cAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAg                               >  2:4357845/1‑93 (MQ=255)
                                                       cAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAg                               >  2:4522487/1‑93 (MQ=255)
                                                       cAGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAg                               <  1:650109/93‑1 (MQ=255)
                                                        aGGGCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAgg                              <  1:438726/93‑1 (MQ=255)
                                                         gggCGACACCCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGc                             <  2:1721842/93‑1 (MQ=255)
                                                         gggCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGc                             >  1:1524832/1‑93 (MQ=255)
                                                           gCGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCa                           <  2:2320808/93‑1 (MQ=255)
                                                            cGACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAg                          <  1:3875592/93‑1 (MQ=255)
                                                             gACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGc                         <  1:2507093/93‑1 (MQ=255)
                                                             gACAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGc                         <  2:507815/93‑1 (MQ=255)
                                                              aCAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGcc                        <  2:1732683/93‑1 (MQ=255)
                                                               cAACCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGccc                       >  2:4487133/1‑93 (MQ=255)
                                                                aaCCACCGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGcccc                      >  2:393144/1‑93 (MQ=255)
                                                                 accaccGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGccccc                     >  1:1956218/1‑93 (MQ=255)
                                                                 accaccGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGccccc                     >  1:365716/1‑93 (MQ=255)
                                                                  ccaccGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCa                    >  2:4256450/1‑93 (MQ=255)
                                                                  ccaccGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCa                    <  2:4046943/93‑1 (MQ=255)
                                                                  ccaccGCCAGGAAGCAGAGGGCGCCTCCCATCCCTTTGTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCa                    <  1:1858197/93‑1 (MQ=255)
                                                                   caccGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCAt                   <  1:3016382/93‑1 (MQ=255)
                                                                    cccGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATg                  <  2:1592654/92‑1 (MQ=255)
                                                                    accGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATTGAGCGGCATACCAGAAAACAGGCCAGCCCCCATg                  >  1:4342548/1‑93 (MQ=255)
                                                                    accGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATg                  <  1:108044/93‑1 (MQ=255)
                                                                    accGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATg                  >  2:1976377/1‑93 (MQ=255)
                                                                    accGCCAGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATg                  >  1:2719180/1‑93 (MQ=255)
                                                                          aGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCCTGGAcacc            >  1:329448/1‑93 (MQ=255)
                                                                          aGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGAcacc            >  2:2679376/1‑93 (MQ=255)
                                                                          aGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGAcacc            >  2:4203117/1‑93 (MQ=255)
                                                                          aGGAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGAcacc            >  2:3092831/1‑93 (MQ=255)
                                                                           ggAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGAcacca           <  1:3794988/93‑1 (MQ=255)
                                                                            gAAGCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGAcaccac          <  2:3430031/93‑1 (MQ=255)
                                                                               gCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCAt       >  1:4485390/1‑93 (MQ=255)
                                                                               gCCGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCAt       >  2:4283764/1‑93 (MQ=255)
                                                                                ccGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGCCACCACCAta      >  1:3661715/1‑92 (MQ=255)
                                                                                ccGAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCAtt      <  1:352361/93‑1 (MQ=255)
                                                                                  gAGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTaa    >  2:220608/1‑93 (MQ=255)
                                                                                   aGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTaaa   <  1:679523/93‑1 (MQ=255)
                                                                                   aGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTaaa   <  1:2472052/93‑1 (MQ=255)
                                                                                   aGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTaaa   >  2:211533/1‑93 (MQ=255)
                                                                                   aGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTaaa   <  1:4455517/93‑1 (MQ=255)
                                                                                   aGGGCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTaaa   >  1:2040444/1‑93 (MQ=255)
                                                                                    gggCGCCTCCCATCCCTTTTTTCATCCTCGCTGCGCTTTTATCGAGCGGCATACCCGAAAAAAGGCCAGCCCCCATGGACACCACCATTaaac  >  2:1778641/1‑92 (MQ=255)
                                                                                    gggCGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTAAAg  >  2:850286/1‑93 (MQ=255)
                                                                                       |                                                                                         
CGGCAATCTGATCGACTGCGCCGGTTTCATGTAAGATCTTGCCAAACATAGCTCCCAGGGCGACAACCACCGCCAGGAAGCCGAGGGTGCCTCCCATCCCTTTTTCCATCGTCGCTGCGATTTTATCGAGCGGCATACCAGAAAAAAGGCCAGCCCCCATGGACACCACCATTAAAG  >  NC_000913/3574426‑3574602

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