Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,894,054 G→T A62D (GCC→GAC)  yieG ← putative membrane/transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,894,0540GT98.0% 673.1 / 3.2 205A62D (GCC→GAC) yieGputative membrane/transport protein
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base T (112/88);  total (113/92)
Fisher's exact test for biased strand distribution p-value = 3.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 218 total aligned reads displayed.

ATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATC  >  NC_000913/3893964‑3894143
                                                                                          |                                                                                         
 tAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcgc                                                                                        >  1:1737630/1‑93 (MQ=255)
  aGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCt                                                                                       >  2:1223311/1‑93 (MQ=255)
  aGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCt                                                                                       >  1:1416518/1‑93 (MQ=255)
  aGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCt                                                                                       <  1:3552575/93‑1 (MQ=255)
   gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                                                                                      <  2:373433/93‑1 (MQ=255)
   gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                                                                                      >  2:3547351/1‑93 (MQ=255)
   gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                                                                                      <  1:4314292/93‑1 (MQ=255)
     gggTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGt                                                                                    >  1:4721062/1‑93 (MQ=255)
     gggTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGt                                                                                    >  1:2148790/1‑93 (MQ=255)
      ggTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTa                                                                                   <  1:3774567/93‑1 (MQ=255)
      ggTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTa                                                                                   >  1:4554722/1‑93 (MQ=255)
          ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                                                                               <  1:4082712/93‑1 (MQ=255)
          ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGGATcct                                                                               >  2:472503/1‑92 (MQ=255)
            aGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATg                                                                             >  2:3422835/1‑93 (MQ=255)
              tGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGcc                                                                           <  2:1569629/93‑1 (MQ=255)
              tGCAACTGACAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGcc                                                                           <  2:1318170/93‑1 (MQ=255)
               gCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCa                                                                          <  2:3682854/93‑1 (MQ=255)
                 aaCTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGc                                                                        <  2:4585439/93‑1 (MQ=255)
                       cAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGcc                                                                  <  2:2460951/93‑1 (MQ=255)
                         ggTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCaa                                                                >  2:223400/1‑93 (MQ=255)
                         ggTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCaa                                                                <  2:2449015/93‑1 (MQ=255)
                         ggTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCaa                                                                <  1:1810960/93‑1 (MQ=255)
                          gTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAg                                                               >  2:3188948/1‑93 (MQ=255)
                           ttggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGa                                                              <  1:2669393/90‑1 (MQ=255)
                           ttggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGa                                                              >  2:667317/4‑93 (MQ=255)
                            tggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGaa                                                             <  1:501471/91‑1 (MQ=255)
                            tggCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGaa                                                             <  2:777676/91‑1 (MQ=255)
                              gCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAAtt                                                           <  1:4138399/93‑1 (MQ=255)
                              gCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAAtt                                                           <  2:1821345/93‑1 (MQ=255)
                               cAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAAttt                                                          <  1:1802888/93‑1 (MQ=255)
                                aaaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTg                                                         >  1:763786/1‑93 (MQ=255)
                                 aaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGc                                                        <  1:3031581/93‑1 (MQ=255)
                                 aaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGc                                                        <  2:662570/93‑1 (MQ=255)
                                 aaCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGc                                                        >  1:3334169/1‑93 (MQ=255)
                                  aCAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                                       <  2:2664887/93‑1 (MQ=255)
                                   cAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGATGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCgg                                                      <  2:4553270/93‑1 (MQ=255)
                                     gTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGt                                                    >  1:611147/1‑93 (MQ=255)
                                     gTCCCATCATAATACCGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGt                                                    >  2:228386/1‑93 (MQ=255)
                                      tCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGtt                                                   >  2:1191711/1‑93 (MQ=255)
                                       cccATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTa                                                  >  1:1226848/1‑93 (MQ=255)
                                        ccATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTaa                                                 >  2:1983961/1‑93 (MQ=255)
                                        ccATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTaa                                                 >  2:2125231/1‑93 (MQ=255)
                                        ccATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTaa                                                 >  2:2817642/1‑93 (MQ=255)
                                         catcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAAc                                                <  1:644570/93‑1 (MQ=255)
                                         caccatAATCATGACGAATGCAGCGATCAGACAGGTAGTGACCAAGACGTCGCTCGTATCAATGCCAACACCGCCACGAATTTGCGCGTCAAc                                                <  2:2095125/90‑1 (MQ=255)
                                          atcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACa                                               <  1:387462/93‑1 (MQ=255)
                                          atcatAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCAACCAACGCAAGGAATTGGCGGGTTACCa                                               >  2:2679203/1‑93 (MQ=255)
                                             ataataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaa                                            >  2:3101342/1‑93 (MQ=255)
                                             ataataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaa                                            >  1:525350/1‑93 (MQ=255)
                                              taataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaaa                                           <  1:3744021/93‑1 (MQ=255)
                                              taataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaaa                                           <  2:1648961/93‑1 (MQ=255)
                                              taataCTGCCGAATGCAACGATCAGATAGGTAGTGACGAAGACGTCGCGGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaaa                                           <  2:4050051/93‑1 (MQ=255)
                                               aataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAAc                                          >  1:1646357/1‑93 (MQ=255)
                                                ataCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACg                                         >  1:1328715/1‑93 (MQ=255)
                                                  aCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAACACGAt                                       >  1:2228442/1‑93 (MQ=255)
                                                    tGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGt                                     >  2:3022175/1‑93 (MQ=255)
                                                      ccGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTaa                                   <  1:3978590/93‑1 (MQ=255)
                                                         aaTGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAAcc                                <  2:2605714/93‑1 (MQ=255)
                                                          aggCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCTGCAACGCCAAGAATTTGCGGGTTAATTAAAACGATGTAAACCa                               >  2:4222437/3‑93 (MQ=255)
                                                          aTGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCa                               >  2:2050151/1‑93 (MQ=255)
                                                          aTGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCa                               <  1:3656470/93‑1 (MQ=255)
                                                           tGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                              >  1:1258830/1‑93 (MQ=255)
                                                            gCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATc                             <  2:3660859/93‑1 (MQ=255)
                                                             cAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCg                            >  1:3268617/1‑93 (MQ=255)
                                                              aGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGt                           >  1:2995854/1‑93 (MQ=255)
                                                              aGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGt                           >  2:1009878/1‑93 (MQ=255)
                                                               gCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTc                          <  2:2496340/93‑1 (MQ=255)
                                                               gCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTc                          >  1:727690/1‑93 (MQ=255)
                                                                cGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCa                         <  2:102927/93‑1 (MQ=255)
                                                                cGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCa                         <  2:4527680/93‑1 (MQ=255)
                                                                   tCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGa                      >  2:2511651/1‑93 (MQ=255)
                                                                    cAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGaa                     <  1:4267569/93‑1 (MQ=255)
                                                                    cAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGaa                     <  2:709977/93‑1 (MQ=255)
                                                                     aGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAg                    <  1:1006421/93‑1 (MQ=255)
                                                                     aGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAg                    >  1:2456150/1‑93 (MQ=255)
                                                                      gACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAATGATGTAAACCATCGTCAGTAAgg                   >  2:3241436/1‑93 (MQ=255)
                                                                      gACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAgg                   >  1:3210085/1‑93 (MQ=255)
                                                                      gACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAgg                   >  2:3772044/1‑93 (MQ=255)
                                                                       gcAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggt                  <  2:1246002/92‑1 (MQ=255)
                                                                       aCAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggt                  <  1:3451026/93‑1 (MQ=255)
                                                                        cAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtg                 >  2:1724845/1‑93 (MQ=255)
                                                                        cAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtg                 <  1:1191849/93‑1 (MQ=255)
                                                                          ggTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtggt               <  1:2705744/93‑1 (MQ=255)
                                                                          ggTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtggt               >  2:3678648/1‑93 (MQ=255)
                                                                          ggTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtgga               >  2:2952410/1‑92 (MQ=255)
                                                                           gTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTa              >  2:3837561/1‑93 (MQ=255)
                                                                           gTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTa              <  2:3979643/93‑1 (MQ=255)
                                                                             aGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTaaa            >  2:1498929/1‑93 (MQ=255)
                                                                             aGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTaaa            >  1:4353311/1‑93 (MQ=255)
                                                                              gTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTaaaa           >  2:4667244/1‑93 (MQ=255)
                                                                                gACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAAcc         <  1:4529613/93‑1 (MQ=255)
                                                                                   gAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAACACGATGTAAACCATCGTCAAGAAGGTGGTAAAACCggg      >  2:4255466/1‑92 (MQ=255)
                                                                                   gAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGc      <  2:4367043/93‑1 (MQ=255)
                                                                                   gAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGc      >  2:550510/1‑93 (MQ=255)
                                                                                    aaGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCg     >  1:3675825/1‑93 (MQ=255)
                                                                                     aGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAACCCATCGTCAGGAAGGTGGTAAAACCGGCGa    <  2:2832125/93‑1 (MQ=255)
                                                                                     aGACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGa    >  1:2364096/1‑93 (MQ=255)
                                                                                      gACGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGAt   <  1:4503034/93‑1 (MQ=255)
                                                                                       aCGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATc  <  2:44517/93‑1 (MQ=255)
                                                                                       aCGTCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATc  <  1:1146655/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGGTAAAACCGGCGATC  >  NC_000913/3893964‑3894143

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: