Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,400 +GAAA coding (483/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3981.A100.0% 280.4 / NA 77E161K (GAA→AAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (33/44);  total (33/44)
*NC_0009134,008,3982.A98.7% 280.3 / ‑3.4 77coding (481/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (33/43);  total (33/44)
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.73e-01
*NC_0009134,008,3983.G98.7% 242.5 / ‑3.5 77coding (481/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (33/43);  total (33/44)
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.16e-01
*NC_0009134,008,3984.A98.7% 276.9 / ‑3.4 77coding (481/516 nt)b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (33/43);  total (33/44)
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.35e-01

CGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGC‑‑‑G‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008310‑4008457
                                                                                               ||||                                                                 
cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aa                                                                 <  1:4432852/93‑1 (MQ=255)
 gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aag                                                                >  1:1895079/1‑92 (MQ=255)
 gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aag                                                                >  2:6599853/1‑92 (MQ=255)
  aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTTAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aaga                                                               >  2:3194996/1‑91 (MQ=255)
  aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aaga                                                               >  2:6472005/1‑91 (MQ=255)
   tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGTTTCAAAAAGCCATGCGCTGGCAGCGGTGGC‑‑‑G‑‑‑‑aagaa                                                              >  1:6384297/1‑90 (MQ=255)
   tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aagaa                                                              >  1:6660115/1‑90 (MQ=255)
   tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aagaa                                                              <  2:1931685/93‑4 (MQ=255)
    cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑G‑‑‑‑aagaaa                                                             <  1:4598119/93‑5 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                <  2:93170/93‑1 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                >  2:6713870/1‑93 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                >  2:1132533/1‑93 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAGGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                <  1:2981338/93‑1 (MQ=255)
                                    tCTGGAAGTGATGGCGGGCGGCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                <  1:6109615/93‑1 (MQ=255)
                                     cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                               >  2:6652401/1‑92 (MQ=255)
                                     cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                               >  1:990386/1‑92 (MQ=255)
                                     cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                               <  1:1673509/93‑2 (MQ=255)
                                     cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                               <  1:1567329/93‑2 (MQ=255)
                                       ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                             >  2:3952663/1‑90 (MQ=255)
                                       ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                             >  2:5544515/1‑90 (MQ=255)
                                       ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                             >  2:6385580/1‑90 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  2:4370028/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  2:4097734/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  2:5650687/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  2:6498533/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  1:832008/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  1:1500944/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                            <  1:5372490/93‑5 (MQ=255)
                                        gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggg                            >  2:1014114/1‑89 (MQ=255)
                                              aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                      <  1:6058621/93‑1 (MQ=255)
                                              aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                      <  1:3213125/93‑1 (MQ=255)
                                              aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                      <  2:3791411/93‑1 (MQ=255)
                                               tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                     <  2:1844340/93‑1 (MQ=255)
                                                ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                    <  2:5203299/93‑1 (MQ=255)
                                                ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                    <  2:3067848/93‑1 (MQ=255)
                                                 gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                   >  1:5841435/1‑93 (MQ=255)
                                                 gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                   >  1:4191617/1‑93 (MQ=255)
                                                  cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCTAAGA‑‑‑AACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                  <  2:912291/93‑1 (MQ=255)
                                                  cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                  >  1:2192480/1‑93 (MQ=255)
                                                  cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                  <  1:31479/93‑1 (MQ=255)
                                                  cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                  <  2:4998234/93‑1 (MQ=255)
                                                   gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 <  1:47385/93‑1 (MQ=255)
                                                   gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 <  2:3841035/93‑1 (MQ=255)
                                                   gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 >  2:806217/1‑93 (MQ=255)
                                                     gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg               >  2:274444/1‑93 (MQ=255)
                                                     gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg               >  2:2993890/1‑93 (MQ=255)
                                                      cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga              <  2:3073690/93‑1 (MQ=255)
                                                       ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             <  1:5688350/93‑1 (MQ=255)
                                                       ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             >  2:5179197/1‑93 (MQ=255)
                                                       ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             <  2:2367591/93‑1 (MQ=255)
                                                       ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             <  2:1368824/93‑1 (MQ=255)
                                                        gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            >  1:2108765/1‑93 (MQ=255)
                                                        gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            <  2:1665924/93‑1 (MQ=255)
                                                         cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc           >  2:777217/1‑93 (MQ=255)
                                                         cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc           <  1:4782314/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGGATTGCGTTTGGTGTCGGGATGAACGACGcc          >  2:6280598/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          >  1:3137376/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          <  2:3562144/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          >  2:3784278/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          >  2:6297100/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          >  2:2100640/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          <  1:4708802/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          >  2:4869242/1‑93 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGa        <  1:6379746/89‑1 (MQ=255)
                                                             ttcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGa        <  1:5667309/89‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa      <  2:1507336/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa      <  2:5979161/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa      <  2:3475704/91‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa      <  1:2934883/91‑1 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt     <  2:5637763/92‑1 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt     <  2:3991366/92‑1 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg    >  1:4506931/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg    <  2:6010772/93‑1 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg    >  1:2205579/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg    >  1:721326/1‑93 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc   <  1:3895288/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc   <  1:6289377/93‑1 (MQ=255)
                                                                  aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc   <  2:5981706/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  1:3236866/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  2:6691263/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  1:296548/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  2:6020789/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  <  2:3264646/93‑1 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  2:2939531/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  2:2767946/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGAAGCCGAAATGcg  >  1:1083384/1‑92 (MQ=255)
                                                                                               ||||                                                                 
CGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGC‑‑‑G‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008310‑4008457

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