Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,812,848 T→C intergenic (‑94/+57) ygaG ← / ← gshA hypothetical protein/gamma‑glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,812,8480TC98.5% 805.4 / ‑0.5 260intergenic (‑94/+57)ygaG/gshAhypothetical protein/gamma‑glutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (168/88);  total (171/89)
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.83e-01

Only 100 of 287 total aligned reads displayed.

GCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGT  >  NC_000913/2812758‑2812938
                                                                                          |                                                                                           
gCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCAccc                                                                                          >  2:2440775/1‑93 (MQ=255)
  cACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCt                                                                                        >  1:2708577/1‑93 (MQ=255)
    ccTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGa                                                                                      >  1:6913047/1‑93 (MQ=255)
    ccTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGa                                                                                      >  2:7022070/1‑93 (MQ=255)
    ccTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGa                                                                                      >  2:4853722/1‑93 (MQ=255)
    ccTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGa                                                                                      >  1:7698010/1‑93 (MQ=255)
     cTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGaa                                                                                     <  1:5221014/93‑1 (MQ=255)
     cTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGaa                                                                                     >  2:1227086/1‑93 (MQ=255)
     cTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGaa                                                                                     >  2:671778/1‑93 (MQ=255)
         ggTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTca                                                                                 >  2:6418043/1‑92 (MQ=255)
         ggTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTca                                                                                 >  2:6203938/1‑92 (MQ=255)
           tAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTcaga                                                                               <  2:7093192/93‑4 (MQ=255)
             aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAgaga                                                                            >  1:3364421/1‑93 (MQ=255)
             aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAgaga                                                                            <  1:984715/93‑1 (MQ=255)
              tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGAt                                                                           >  1:2473960/1‑93 (MQ=255)
              tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGAt                                                                           >  1:784975/1‑93 (MQ=255)
              tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCCGAATTCAGAGAt                                                                           >  2:4266389/1‑93 (MQ=255)
               ttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATg                                                                          <  2:6003017/93‑1 (MQ=255)
                tttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGa                                                                         >  2:6083092/1‑93 (MQ=255)
                 ttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGaa                                                                        >  1:4534284/1‑93 (MQ=255)
                 ttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGaa                                                                        <  1:1286532/93‑1 (MQ=255)
                  tttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGaaa                                                                       <  2:1570614/93‑1 (MQ=255)
                  tttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGaaa                                                                       >  1:1576336/1‑93 (MQ=255)
                   ttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAt                                                                      >  1:6735453/1‑93 (MQ=255)
                   ttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAAt                                                                      <  1:4975956/93‑1 (MQ=255)
                        aaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCCCCCCTGAATTCAGAGATGAAATTTTgg                                                                 <  2:358059/93‑1 (MQ=255)
                           tttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCa                                                              >  1:7279834/1‑93 (MQ=255)
                           tttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCa                                                              <  2:1594359/93‑1 (MQ=255)
                            ttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCAc                                                             >  1:6283866/1‑93 (MQ=255)
                            ttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCAc                                                             >  1:2927324/1‑93 (MQ=255)
                            ttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCAc                                                             <  1:6465960/93‑1 (MQ=255)
                              cTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTc                                                           >  1:5975589/1‑93 (MQ=255)
                               tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCCCCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                          <  1:4218741/93‑1 (MQ=255)
                               tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                          >  1:2613306/1‑93 (MQ=255)
                               tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                          >  2:3709063/1‑93 (MQ=255)
                                gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                         >  2:6289175/1‑93 (MQ=255)
                                 aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                                        >  1:1529978/1‑93 (MQ=255)
                                 aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                                        >  2:1699506/1‑93 (MQ=255)
                                  aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                       >  1:356724/1‑93 (MQ=255)
                                   ctctTTCTTCCCAGGCGAGTCTGGGAATATGAAAGACGCGCATTTGTTATCATCCCCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                                      >  2:1961770/1‑93 (MQ=255)
                                   ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                       <  2:5756941/93‑1 (MQ=255)
                                   ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                                      <  1:3826121/93‑1 (MQ=255)
                                    tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                     >  1:4715102/1‑93 (MQ=255)
                                     ctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTg                                                    <  2:4698222/93‑1 (MQ=255)
                                       ttcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGc                                                  >  1:1720907/1‑93 (MQ=255)
                                        tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGcc                                                 <  1:6949666/93‑1 (MQ=255)
                                        tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGcc                                                 <  2:316897/93‑1 (MQ=255)
                                         cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCt                                                <  1:4023571/93‑1 (MQ=255)
                                          ttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                               >  2:2891715/1‑93 (MQ=255)
                                          ttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                               >  2:3382207/1‑93 (MQ=255)
                                             ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGCAATTTTGGCCCCTCCCGCGCGGCCttttt                                            >  1:67165/1‑93 (MQ=255)
                                              cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttc                                           >  1:1172859/1‑93 (MQ=255)
                                               aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                                          >  2:267102/1‑93 (MQ=255)
                                               aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                                          <  2:6638557/93‑1 (MQ=255)
                                               aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                                          <  2:919628/93‑1 (MQ=255)
                                               aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                                          >  1:1913730/1‑93 (MQ=255)
                                                 gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                                        >  2:1441864/1‑93 (MQ=255)
                                                   gAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                      >  2:5349759/1‑93 (MQ=255)
                                                   gAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                      <  1:211418/93‑1 (MQ=255)
                                                    aGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTtttcttttct                                     >  1:6615888/1‑93 (MQ=255)
                                                     gTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTg                                    >  2:3343378/1‑93 (MQ=255)
                                                       cTGAGTATATGAAAGACGCGCATTTGTTATCATCGCCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTTGCCTTTTTCTTTTCTGTc                                  <  2:3924085/93‑1 (MQ=255)
                                                        tGAGTATATGAAAGACGCGCATTTGTTATCATCTCCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                 <  2:3017303/93‑1 (MQ=255)
                                                        tGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                 <  2:1845281/93‑1 (MQ=255)
                                                        tGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                 >  2:784686/1‑93 (MQ=255)
                                                        tGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                 <  2:2291951/93‑1 (MQ=255)
                                                         gAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                <  2:2176512/93‑1 (MQ=255)
                                                         gAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                <  2:4168615/93‑1 (MQ=255)
                                                           gTATATGAAAGACGCGCATCTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                              <  1:3972411/93‑1 (MQ=255)
                                                             atatGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgt                            >  1:1465596/1‑93 (MQ=255)
                                                              tatGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGGCACTCACGAGCGGCCTTTTTCTTTTATGTCAGGCgtg                           >  1:2978822/1‑93 (MQ=255)
                                                              tatGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                           >  1:7498563/1‑93 (MQ=255)
                                                              tatGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                           >  2:4934223/1‑93 (MQ=255)
                                                               atGAAAGACGCGCTTTTGTTATCATCTCCCCTGAAGTCAGAGATGAAGTTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtgt                          <  1:581375/93‑1 (MQ=255)
                                                               atGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtgt                          >  2:6493807/1‑93 (MQ=255)
                                                               atGAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtgt                          >  2:2600635/1‑93 (MQ=255)
                                                                 gAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt                        <  1:5652849/93‑1 (MQ=255)
                                                                 gAAAGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt                        >  1:5381151/1‑93 (MQ=255)
                                                                  aaaGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtttt                       <  2:893825/93‑1 (MQ=255)
                                                                    aGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc                     >  1:5557828/1‑93 (MQ=255)
                                                                    aGACGCGCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc                     >  1:843049/1‑93 (MQ=255)
                                                                       cgcgCATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAg                  >  2:253151/1‑93 (MQ=255)
                                                                        gcgcATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGc                 >  2:7926780/1‑93 (MQ=255)
                                                                          gcATTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCca               >  2:3499740/1‑93 (MQ=255)
                                                                           cATTTGTTATCATCATCCCTGAATTCGGAGATGAAATTTTGTCAACTCACGAGTGGCCTTTTGCTTTTCTGTCAGACGTGGTTTTccacccac              >  1:4320388/1‑93 (MQ=255)
                                                                            aTTTGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcaca             >  1:325341/1‑93 (MQ=255)
                                                                             tttGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                 >  1:1097721/1‑72 (MQ=255)
                                                                             tttGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                 <  2:1097721/72‑1 (MQ=255)
                                                                             tttGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCacac            >  1:856956/1‑93 (MQ=255)
                                                                               tGTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTTCAGCCACACCg          >  2:4314674/1‑93 (MQ=255)
                                                                                gTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGc         >  2:6873555/1‑93 (MQ=255)
                                                                                gTTATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGc         >  1:6361159/1‑93 (MQ=255)
                                                                                  tATCATCACCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaa       >  2:7581230/1‑93 (MQ=255)
                                                                                   atcatcaccCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa      >  2:4201543/1‑93 (MQ=255)
                                                                                     catcaccCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACCCCGCAAACg    <  1:1343990/93‑1 (MQ=255)
                                                                                      atcaccCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACgg   >  1:1615973/1‑93 (MQ=255)
                                                                                      atcaccCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACgg   >  2:1166492/1‑93 (MQ=255)
                                                                                       tcaccCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGt  >  1:2829305/1‑93 (MQ=255)
                                                                                       tcaccCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGt  <  2:7309618/93‑1 (MQ=255)
                                                                                       tcaccCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGt  >  2:5015596/1‑93 (MQ=255)
                                                                                          |                                                                                           
GCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGT  >  NC_000913/2812758‑2812938

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: