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.6% 1249.8 / ‑3.8 272intergenic (‑104/+47)ygaG/gshAhypothetical protein/gamma‑glutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (162/109);  total (163/109)
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.92e-01

Only 100 of 283 total aligned reads displayed.

CGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATC  >  NC_000913/2812766‑2812946
                                                                                            |                                                                                        
 ggTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTca                                                                                         >  2:348814/1‑92 (MQ=255)
  nnaaTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                                        <  2:1444748/91‑3 (MQ=255)
  gTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                                        >  1:1588894/1‑91 (MQ=255)
   tAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcaga                                                                                       <  2:7049510/93‑4 (MQ=255)
   tAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcaga                                                                                       >  2:2563815/1‑90 (MQ=255)
     aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAgaga                                                                                    <  2:6502094/93‑1 (MQ=255)
     aTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAgaga                                                                                    <  1:408979/93‑1 (MQ=255)
         ttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATCCAGAGATGaa                                                                                >  1:3150959/1‑93 (MQ=255)
          tttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGaaa                                                                               >  1:2536829/1‑93 (MQ=255)
            tAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtt                                                                             >  1:1155653/1‑93 (MQ=255)
             aaaaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAttt                                                                            <  1:7023307/93‑1 (MQ=255)
              aaaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtttt                                                                           >  1:2381519/1‑93 (MQ=255)
              aaaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtttt                                                                           <  1:4317639/93‑1 (MQ=255)
                aaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTgg                                                                         <  2:1071854/93‑1 (MQ=255)
                 aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGc                                                                        <  1:677160/93‑1 (MQ=255)
                 aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGc                                                                        >  1:5119841/1‑93 (MQ=255)
                 aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGc                                                                        >  2:5735066/1‑93 (MQ=255)
                  ttttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGcc                                                                       >  1:6103090/1‑93 (MQ=255)
                       tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                                  >  2:4240693/1‑93 (MQ=255)
                        gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                 >  2:2423318/1‑93 (MQ=255)
                        gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                 >  1:6288296/1‑93 (MQ=255)
                        gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                 >  2:1836389/1‑93 (MQ=255)
                        gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                 >  2:2493626/1‑93 (MQ=255)
                         aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                                                >  2:4174379/1‑93 (MQ=255)
                          aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                               <  2:2285148/93‑1 (MQ=255)
                          aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                               <  2:3274623/93‑1 (MQ=255)
                          aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                               <  1:4761808/93‑1 (MQ=255)
                           ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                                              >  1:6363457/1‑93 (MQ=255)
                            tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                             <  1:3415325/93‑1 (MQ=255)
                            tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                                             >  1:2064036/1‑93 (MQ=255)
                              tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                           >  1:4404658/1‑93 (MQ=255)
                              tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                           >  1:4219899/1‑93 (MQ=255)
                               ttcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGc                                                          >  1:462350/1‑93 (MQ=255)
                                tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGcc                                                         <  2:4464725/93‑1 (MQ=255)
                                tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGcc                                                         >  2:2916716/1‑93 (MQ=255)
                                tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGcc                                                         >  2:2209977/1‑93 (MQ=255)
                                 cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCt                                                        >  1:7428429/1‑93 (MQ=255)
                                  ttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                                       >  2:4676512/1‑93 (MQ=255)
                                   tcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttt                                                      <  1:134450/93‑1 (MQ=255)
                                    cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                                                     >  1:5339767/1‑93 (MQ=255)
                                     ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                    <  2:6961020/93‑1 (MQ=255)
                                     ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                    >  2:5643759/1‑93 (MQ=255)
                                      cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttc                                                   >  2:7223039/1‑93 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 <  1:5602656/93‑1 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 <  1:2837224/93‑1 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  2:6298809/1‑93 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  2:731726/1‑93 (MQ=255)
                                        ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                 >  2:5849386/1‑93 (MQ=255)
                                         gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                                                >  2:3804077/1‑93 (MQ=255)
                                           gAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                              <  2:504757/93‑1 (MQ=255)
                                           gAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                              >  2:5396126/1‑93 (MQ=255)
                                             gTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTg                                            >  2:5381792/1‑93 (MQ=255)
                                             gTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTg                                            <  1:1031867/93‑1 (MQ=255)
                                                tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                         <  1:180256/93‑1 (MQ=255)
                                                tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                         >  2:1007952/1‑93 (MQ=255)
                                                tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                         <  1:2880084/93‑1 (MQ=255)
                                                tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                                         >  1:3903943/1‑93 (MQ=255)
                                                 gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                        >  1:3258583/1‑93 (MQ=255)
                                                 gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                                        <  2:6001295/93‑1 (MQ=255)
                                                  aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAgg                                       <  2:1595901/93‑1 (MQ=255)
                                                  aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAgg                                       >  1:5916028/1‑93 (MQ=255)
                                                  aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAgg                                       <  2:4780173/93‑1 (MQ=255)
                                                  aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAgg                                       >  1:6899809/1‑93 (MQ=255)
                                                   gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                                      <  1:5379370/93‑1 (MQ=255)
                                                      tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                   >  2:6632987/1‑93 (MQ=255)
                                                      tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                   >  1:1370901/1‑93 (MQ=255)
                                                        tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                                 >  1:7052401/1‑93 (MQ=255)
                                                        tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                                 >  1:7249370/1‑93 (MQ=255)
                                                        tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                                 >  1:2469529/1‑93 (MQ=255)
                                                         gAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt                                >  2:186822/1‑93 (MQ=255)
                                                          aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtttt                               >  2:2149707/1‑93 (MQ=255)
                                                           aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttttt                              <  1:1705892/93‑1 (MQ=255)
                                                           aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttttt                              <  1:6507985/93‑1 (MQ=255)
                                                            aGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc                             <  1:1480473/93‑1 (MQ=255)
                                                             gACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTcc                            <  1:953788/93‑1 (MQ=255)
                                                              aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCa                           >  1:2540548/1‑93 (MQ=255)
                                                                gcgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGc                         >  1:501478/1‑93 (MQ=255)
                                                                  gcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCca                       >  1:790359/1‑93 (MQ=255)
                                                                   cATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac                      >  2:1281394/1‑93 (MQ=255)
                                                                   cATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac                      >  1:2968846/1‑93 (MQ=255)
                                                                   cATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac                      >  2:4943485/1‑93 (MQ=255)
                                                                    aTTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcaca                     >  1:6654138/1‑93 (MQ=255)
                                                                     tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCacac                    <  2:4040680/93‑1 (MQ=255)
                                                                     tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCacac                    <  2:3731850/93‑1 (MQ=255)
                                                                        gTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGc                 >  2:7003243/1‑93 (MQ=255)
                                                                          tATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaa               >  2:2738615/1‑93 (MQ=255)
                                                                           atcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa              <  2:6103090/93‑1 (MQ=255)
                                                                           atcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa              >  1:5851838/1‑93 (MQ=255)
                                                                            tcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAAc             >  1:3928688/1‑93 (MQ=255)
                                                                             catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACg            <  2:4683009/93‑1 (MQ=255)
                                                                             catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACg            <  1:5397733/93‑1 (MQ=255)
                                                                                catcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt         >  2:5552832/1‑93 (MQ=255)
                                                                                catcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGtt         <  2:76406/93‑1 (MQ=255)
                                                                                  tcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCg       <  2:6372100/93‑1 (MQ=255)
                                                                                  tcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCg       >  1:4938111/1‑93 (MQ=255)
                                                                                   cccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCgg      >  1:3694864/1‑93 (MQ=255)
                                                                                    ccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACAACGCAAACGGTTCGGt     >  1:3764747/1‑93 (MQ=255)
                                                                                     cTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTa    >  1:4804802/1‑93 (MQ=255)
                                                                                       gAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATc  <  1:3530580/93‑1 (MQ=255)
                                                                                       gAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATc  >  2:837632/1‑93 (MQ=255)
                                                                                            |                                                                                        
CGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATC  >  NC_000913/2812766‑2812946

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