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.100.0% 561.4 / NA 122intergenic (‑104/+47)ygaG/gshAhypothetical protein/gamma‑glutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (73/49);  total (73/49)

Only 100 of 130 total aligned reads displayed.

CGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATC  >  NC_000913/2812766‑2812946
                                                                                            |                                                                                        
cGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTc                                                                                          >  1:4138907/1‑93 (MQ=255)
cGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTc                                                                                          >  2:2370937/1‑93 (MQ=255)
 ggTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTca                                                                                         >  2:2690440/1‑92 (MQ=255)
  gTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                                        >  1:1140301/1‑91 (MQ=255)
  gTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                                        >  2:397203/1‑91 (MQ=255)
     aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAgaga                                                                                    >  1:534630/1‑93 (MQ=255)
      tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGAt                                                                                   >  2:474850/1‑93 (MQ=255)
      tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGAt                                                                                   >  1:3385495/1‑93 (MQ=255)
       ttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATg                                                                                  <  2:3961363/93‑1 (MQ=255)
       ttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATg                                                                                  <  2:1707034/93‑1 (MQ=255)
         ttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGaa                                                                                >  2:2777315/1‑93 (MQ=255)
          tttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGaaa                                                                               >  2:1423507/1‑93 (MQ=255)
          tttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGaaa                                                                               <  1:287456/93‑1 (MQ=255)
            tAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtt                                                                             <  2:1183809/93‑1 (MQ=255)
            tAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtt                                                                             <  2:3433822/93‑1 (MQ=255)
             aaaaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAttt                                                                            >  1:390168/1‑93 (MQ=255)
               aaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTg                                                                          >  2:4557770/1‑93 (MQ=255)
                 aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGc                                                                        >  2:1665366/1‑93 (MQ=255)
                      cTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTc                                                                   >  1:2258099/1‑93 (MQ=255)
                        gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                 >  2:2065106/1‑93 (MQ=255)
                        gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                 <  1:2360154/93‑1 (MQ=255)
                         aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                                                >  2:2757020/1‑93 (MQ=255)
                          aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                               <  2:1621809/93‑1 (MQ=255)
                          aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                               <  2:1408438/93‑1 (MQ=255)
                          aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                               >  2:3291862/1‑93 (MQ=255)
                           ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                                              <  1:4525863/93‑1 (MQ=255)
                            tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                             >  2:3144205/1‑93 (MQ=255)
                            tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                             >  1:2120053/1‑93 (MQ=255)
                            tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                             >  1:4039691/1‑93 (MQ=255)
                            tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                             <  2:4039691/93‑1 (MQ=255)
                             ctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTg                                                            >  1:1703166/1‑93 (MQ=255)
                              tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                           >  1:4581582/1‑93 (MQ=255)
                              tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                           >  2:3092416/1‑93 (MQ=255)
                              tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                           >  1:2837963/1‑93 (MQ=255)
                                  ttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                                       >  2:2904589/1‑93 (MQ=255)
                                  ttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                                       <  2:1296138/93‑1 (MQ=255)
                                    cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                                                     <  1:4049148/93‑1 (MQ=255)
                                    cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                                                     >  1:4415192/1‑93 (MQ=255)
                                    cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                                                     >  1:3392446/1‑93 (MQ=255)
                                    cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGACTGGCCtttt                                                     <  2:3143307/93‑1 (MQ=255)
                                     ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                    >  2:626407/1‑93 (MQ=255)
                                      cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCATGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttc                                                   <  2:4376984/93‑1 (MQ=255)
                                       aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                                                  >  1:2947141/1‑93 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  1:1044235/1‑93 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  2:3754478/1‑93 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  2:1267277/1‑93 (MQ=255)
                                         gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGATATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                                                <  2:70920/93‑1 (MQ=255)
                                         gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                                                <  1:353285/93‑1 (MQ=255)
                                           gAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                              >  2:3363842/1‑93 (MQ=255)
                                            aGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTtttcttttct                                             >  2:914676/1‑93 (MQ=255)
                                            aGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTtttcttttct                                             >  1:4365611/1‑93 (MQ=255)
                                              tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGt                                           >  2:327261/1‑93 (MQ=255)
                                               cTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTc                                          >  2:1010586/1‑93 (MQ=255)
                                               cTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTc                                          >  2:3427837/1‑93 (MQ=255)
                                                tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                         <  2:4277568/93‑1 (MQ=255)
                                                tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                         >  2:4046861/1‑93 (MQ=255)
                                                tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                         >  1:344593/1‑93 (MQ=255)
                                                 gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                        >  2:4189856/1‑93 (MQ=255)
                                                 gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                        <  1:1218799/93‑1 (MQ=255)
                                                   gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                                      <  1:1965483/93‑1 (MQ=255)
                                                   gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                                      <  1:205597/93‑1 (MQ=255)
                                                    tataTGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCg                                     <  2:1418301/93‑1 (MQ=255)
                                                     atatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgt                                    <  1:4450336/93‑1 (MQ=255)
                                                     atatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgt                                    >  2:1646109/1‑93 (MQ=255)
                                                     atatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgt                                    <  2:2542425/93‑1 (MQ=255)
                                                      tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCCCTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                   >  1:478227/1‑93 (MQ=255)
                                                      tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                   >  2:49690/1‑93 (MQ=255)
                                                      tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                   >  1:406040/1‑93 (MQ=255)
                                                      tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                   <  1:364432/93‑1 (MQ=255)
                                                      tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGATGCCTTTTTCTTTTCTGTCAGGCgcg                                   >  2:4691690/1‑91 (MQ=255)
                                                        tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCCCTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                                 <  2:3981284/93‑1 (MQ=255)
                                                        tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                                 <  2:3981283/93‑1 (MQ=255)
                                                        tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                                 <  1:756770/93‑1 (MQ=255)
                                                        tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                                 >  2:1843846/1‑93 (MQ=255)
                                                         gAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt                                >  1:896009/1‑93 (MQ=255)
                                                          aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtttt                               <  1:4735172/93‑1 (MQ=255)
                                                           aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttttt                              >  2:4324443/1‑93 (MQ=255)
                                                           aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttttt                              >  1:1912452/1‑93 (MQ=255)
                                                             gACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTcc                            <  1:1488871/93‑1 (MQ=255)
                                                             gACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTcc                            <  2:1830464/93‑1 (MQ=255)
                                                             gACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTcc                            >  1:3894075/1‑93 (MQ=255)
                                                             gACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTcc                            <  2:2357794/93‑1 (MQ=255)
                                                             gACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTcc                            <  2:2042031/93‑1 (MQ=255)
                                                              aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCa                           >  2:2529180/1‑93 (MQ=255)
                                                                  gcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCca                       >  1:1333142/1‑93 (MQ=255)
                                                                     tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCacac                    >  1:1034058/1‑93 (MQ=255)
                                                                     tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTGTCTGTCAGGCGTGTTTTTCCAGCCacac                    >  1:4204874/1‑93 (MQ=255)
                                                                      ttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACAcc                   <  2:1095799/93‑1 (MQ=255)
                                                                         ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCa                >  1:502638/1‑93 (MQ=255)
                                                                           atcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa              <  1:2991531/93‑1 (MQ=255)
                                                                            tcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAAc             >  1:2826094/1‑93 (MQ=255)
                                                                            tcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAAc             <  2:898412/93‑1 (MQ=255)
                                                                             catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACg            >  2:4716049/1‑93 (MQ=255)
                                                                             catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACg            <  2:87731/93‑1 (MQ=255)
                                                                               tcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGt          <  1:4547601/93‑1 (MQ=255)
                                                                                  tcacTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCg       >  1:1829947/4‑93 (MQ=255)
                                                                                      tGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCTTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt   <  1:1408021/93‑1 (MQ=255)
                                                                                      tGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt   <  1:4133147/93‑1 (MQ=255)
                                                                                      tGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt   <  1:4133148/93‑1 (MQ=255)
                                                                                       gAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGGCACACCGCCAACGGTTCGGTATc  >  2:2596923/1‑93 (MQ=255)
                                                                                            |                                                                                        
CGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATC  >  NC_000913/2812766‑2812946

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