Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,812,859 Δ1 bp intergenic (‑105/+46) 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,8580C.99.4% 792.7 / ‑2.3 169intergenic (‑104/+47)ygaG/gshAhypothetical protein/gamma‑glutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (103/65);  minor base T (1/0);  total (104/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

Only 100 of 183 total aligned reads displayed.

CGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATC  >  NC_000913/2812766‑2812946
                                                                                            |                                                                                        
cGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTc                                                                                          <  1:513256/93‑1 (MQ=255)
cGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTc                                                                                          <  2:2369996/93‑1 (MQ=255)
cGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTc                                                                                          <  2:492906/93‑1 (MQ=255)
  gTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                                        >  2:1517178/1‑91 (MQ=255)
  gTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                                        <  2:3362623/93‑3 (MQ=255)
    aaTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcagag                                                                                      >  2:5079583/1‑89 (MQ=255)
    aaTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcagag                                                                                      <  2:4277795/93‑5 (MQ=255)
     aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAgaga                                                                                    >  1:5025136/1‑93 (MQ=255)
     aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAgaga                                                                                    <  2:3583475/93‑1 (MQ=255)
     aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGCCTGAGTATATGAAAGACCCGCATTTGTTATCATGATCCCTGAATTCAgaga                                                                                    >  1:5843772/1‑93 (MQ=255)
     aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGAGCGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCACTGAATTCAgaga                                                                                    >  1:5768862/1‑93 (MQ=255)
      tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGAt                                                                                   <  1:5977092/93‑1 (MQ=255)
       ttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATg                                                                                  >  1:790610/1‑93 (MQ=255)
       ttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATg                                                                                  >  2:1273598/1‑93 (MQ=255)
       ttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATg                                                                                  >  1:5769604/1‑93 (MQ=255)
        tttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGa                                                                                 <  2:4825197/93‑1 (MQ=255)
          tttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGaaa                                                                               <  2:2348080/93‑1 (MQ=255)
           tgaaaaaTTTGATGAAGTCTTTCTTCCCAGGCAAGTCCGAGTATAGGAAAGCCGCGCAGTTGTTGTCATCATCCCTGAATTCAGAGATGAAAt                                                                              <  1:1670954/91‑1 (MQ=255)
            tAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtt                                                                             <  1:3241392/93‑1 (MQ=255)
             aaaaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAttt                                                                            >  2:2220593/1‑93 (MQ=255)
              aaaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtttt                                                                           <  1:2191022/93‑1 (MQ=255)
                 aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGc                                                                        >  1:3908291/1‑93 (MQ=255)
                  ttttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGcc                                                                       >  2:6001670/1‑93 (MQ=255)
                  ttttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGCATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGcc                                                                       <  2:5816922/93‑1 (MQ=255)
                   tttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCa                                                                      <  2:570478/93‑1 (MQ=255)
                    ttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCAc                                                                     <  2:3160524/93‑1 (MQ=255)
                      cTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTc                                                                   <  2:1318671/93‑1 (MQ=255)
                       tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                                  >  2:4310920/1‑93 (MQ=255)
                       tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                                  >  1:4113854/1‑93 (MQ=255)
                        gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                 >  2:2786771/1‑93 (MQ=255)
                        gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                 >  1:332555/1‑93 (MQ=255)
                         aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                                                >  2:821680/1‑93 (MQ=255)
                         aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                                                >  2:6197600/1‑93 (MQ=255)
                           ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                                              >  2:1483473/1‑93 (MQ=255)
                           ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                                              >  2:3900829/1‑93 (MQ=255)
                            tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAgg                                                             >  2:1146842/1‑92 (MQ=255)
                            tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                             <  1:357062/93‑1 (MQ=255)
                             ctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTg                                                            >  1:3401812/1‑93 (MQ=255)
                              tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                           <  1:1138275/93‑1 (MQ=255)
                               ttcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGc                                                          >  2:744062/1‑93 (MQ=255)
                                tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCTCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGcc                                                         <  2:6384278/93‑1 (MQ=255)
                                  ttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                                       >  1:845365/1‑93 (MQ=255)
                                   tcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttt                                                      >  1:4016880/1‑93 (MQ=255)
                                     ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                    >  2:5719394/1‑93 (MQ=255)
                                     ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                    >  1:6346224/1‑93 (MQ=255)
                                      cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttc                                                   >  1:1171142/1‑93 (MQ=255)
                                      cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttc                                                   <  1:4589702/93‑1 (MQ=255)
                                      cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttc                                                   >  1:5857755/1‑93 (MQ=255)
                                        ggCGAGTCTTAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 <  1:1040761/93‑1 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  1:5726304/1‑93 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  1:5615211/1‑93 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  2:5371245/1‑93 (MQ=255)
                                          cGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctttt                                               >  1:2057740/1‑93 (MQ=255)
                                           gAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                              <  2:274342/93‑1 (MQ=255)
                                              tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGt                                           <  2:1142286/93‑1 (MQ=255)
                                                 gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                        >  1:5019799/1‑93 (MQ=255)
                                                 gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                        >  1:3262754/1‑93 (MQ=255)
                                                 gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                        >  2:6226983/1‑93 (MQ=255)
                                                   gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                                      >  2:1702581/1‑93 (MQ=255)
                                                   gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                                      <  1:152722/93‑1 (MQ=255)
                                                   gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                                      >  1:4322072/1‑93 (MQ=255)
                                                    tataTGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCg                                     <  1:4561497/93‑1 (MQ=255)
                                                     atatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgt                                    >  1:3048662/1‑93 (MQ=255)
                                                      tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                   >  1:1660274/1‑93 (MQ=255)
                                                       atGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtgt                                  <  2:3262344/93‑1 (MQ=255)
                                                        tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                                 <  1:3508892/93‑1 (MQ=255)
                                                          aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtttt                               >  1:1260724/1‑93 (MQ=255)
                                                          aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtttt                               >  1:1814077/1‑93 (MQ=255)
                                                           aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttttt                              >  2:557186/1‑93 (MQ=255)
                                                           aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttttt                              >  2:3819175/1‑93 (MQ=255)
                                                            aGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc                             >  2:4471743/1‑93 (MQ=255)
                                                            aGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc                             >  2:3556364/1‑93 (MQ=255)
                                                             gACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTcc                            >  2:5743618/1‑93 (MQ=255)
                                                              aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCa                           <  1:5564191/93‑1 (MQ=255)
                                                                 cgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGcc                        <  1:6459/93‑1 (MQ=255)
                                                                   cATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac                      >  2:141456/1‑93 (MQ=255)
                                                                    aTTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcaca                     >  1:4780706/1‑93 (MQ=255)
                                                                    aTTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcaca                     <  2:681771/93‑1 (MQ=255)
                                                                    aTTTGTTAACATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcaca                     >  1:6290970/1‑93 (MQ=255)
                                                                        gTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGc                 >  2:2436413/1‑93 (MQ=255)
                                                                        gTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGAGTGTTTTTCCAGCCACACCGc                 >  1:2105695/1‑93 (MQ=255)
                                                                           atcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa              <  2:5989134/93‑1 (MQ=255)
                                                                           atcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa              >  2:4798601/1‑93 (MQ=255)
                                                                           atcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa              <  1:518893/93‑1 (MQ=255)
                                                                            tcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAAc             >  1:2908723/1‑93 (MQ=255)
                                                                            tcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAAc             <  1:4008588/93‑1 (MQ=255)
                                                                            tcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAAc             <  1:4856248/93‑1 (MQ=255)
                                                                             catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACg            <  1:5125843/93‑1 (MQ=255)
                                                                               tcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGt          >  1:2784924/1‑93 (MQ=255)
                                                                                catcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt         >  2:3428074/1‑93 (MQ=255)
                                                                                 atcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc        >  2:2923444/1‑93 (MQ=255)
                                                                                 atcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc        >  1:3852106/1‑93 (MQ=255)
                                                                                  tcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCg       >  1:5286966/1‑93 (MQ=255)
                                                                                  tcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCg       <  1:2743974/93‑1 (MQ=255)
                                                                                  tcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCg       <  1:5039358/93‑1 (MQ=255)
                                                                                      tGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt   >  1:2571786/1‑93 (MQ=255)
                                                                                      tGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt   >  1:3715267/1‑93 (MQ=255)
                                                                                      tGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt   >  1:4914178/1‑93 (MQ=255)
                                                                                       gAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATc  <  2:3468460/93‑1 (MQ=255)
                                                                                       gAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATc  >  2:3047187/1‑93 (MQ=255)
                                                                                            |                                                                                        
CGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATC  >  NC_000913/2812766‑2812946

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