Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 504,799 A→G Y221C (TAC→TGC)  ushA → UDP‑sugar hydrolase (5'‑nucleotidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913504,7990AG99.6% 832.2 / ‑5.1 264Y221C (TAC→TGC) ushAUDP‑sugar hydrolase (5'‑nucleotidase)
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (137/126);  total (138/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

Only 100 of 270 total aligned reads displayed.

GCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTG  >  NC_000913/504710‑504888
                                                                                         |                                                                                         
gCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGa                                                                                        >  2:7436795/1‑93 (MQ=255)
  ccgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATa                                                                                      <  1:950658/93‑1 (MQ=255)
  ccgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATa                                                                                      >  2:7005250/1‑93 (MQ=255)
  ccgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATa                                                                                      >  1:7014595/1‑93 (MQ=255)
  ccgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATa                                                                                      <  1:2376267/93‑1 (MQ=255)
   cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     <  2:7285589/93‑1 (MQ=255)
   cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     <  2:3698705/93‑1 (MQ=255)
   cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     <  1:3428090/93‑1 (MQ=255)
   cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     >  2:5371911/1‑93 (MQ=255)
   cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     <  2:1565070/93‑1 (MQ=255)
       gATGAAACGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGAGCATTGCGATAATGGt                                                                                 <  2:7025102/93‑1 (MQ=255)
              cgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACg                                                                          >  1:7349048/1‑93 (MQ=255)
              cgaagcTGGTGATTAAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACg                                                                          <  2:4360274/93‑1 (MQ=255)
               gaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACgg                                                                         >  2:1331691/1‑93 (MQ=255)
               gaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACgg                                                                         >  1:6287512/1‑93 (MQ=255)
                 agcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGct                                                                       <  2:2194347/93‑1 (MQ=255)
                   cTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGctct                                                                     <  1:2781440/93‑1 (MQ=255)
                    tGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTa                                                                    >  2:354616/1‑93 (MQ=255)
                    tGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTa                                                                    >  2:4241290/1‑93 (MQ=255)
                     ggTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTaa                                                                   <  2:5695932/93‑1 (MQ=255)
                     ggTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTaa                                                                   <  1:2622235/93‑1 (MQ=255)
                     ggTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTaa                                                                   <  1:7160652/93‑1 (MQ=255)
                     ggTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTaa                                                                   <  2:512903/93‑1 (MQ=255)
                       tGATTCAGGAGATGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACg                                                                 >  1:7786228/1‑93 (MQ=255)
                        gATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGc                                                                >  1:3347518/1‑93 (MQ=255)
                         aTTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCa                                                               >  2:6850659/1‑93 (MQ=255)
                         aTTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCa                                                               <  1:4666556/93‑1 (MQ=255)
                         aTTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCa                                                               <  1:1568339/93‑1 (MQ=255)
                             aGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCAccag                                                           >  1:3836290/1‑91 (MQ=255)
                             aGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCgg                                                           >  1:1757906/1‑93 (MQ=255)
                                 gCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGa                                                       <  2:7914953/93‑1 (MQ=255)
                                  cTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGAt                                                      <  1:3332529/93‑1 (MQ=255)
                                   tGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGAtg                                                     >  2:5022254/1‑93 (MQ=255)
                                   tGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGAtg                                                     <  2:1195486/93‑1 (MQ=255)
                                      aaCAGACAGAAACGCCAGACATTATTATCGCGGCGCCCCATATGGGGCGTTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTgg                                                  <  1:4675070/93‑1 (MQ=255)
                                      aaCAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTgg                                                  >  2:3330135/1‑93 (MQ=255)
                                       acagacagAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGga                                                 <  2:4082198/93‑1 (MQ=255)
                                         agacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGgaga                                               <  1:540995/93‑1 (MQ=255)
                                         agacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGgaga                                               <  2:1352276/93‑1 (MQ=255)
                                          gacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGAt                                              >  1:1799230/1‑93 (MQ=255)
                                           acagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATg                                             >  2:1041496/1‑93 (MQ=255)
                                            cagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATgg                                            <  1:5142953/93‑1 (MQ=255)
                                             agaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGc                                           >  1:3809470/1‑93 (MQ=255)
                                              gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  1:3033995/1‑93 (MQ=255)
                                              gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  1:6454700/1‑93 (MQ=255)
                                              gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  1:7588195/1‑93 (MQ=255)
                                               aaaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAc                                         >  2:5955759/1‑93 (MQ=255)
                                               aaaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAc                                         <  1:770872/93‑1 (MQ=255)
                                               aaaaGCCAGACATTATTATCGCGGCCACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAc                                         <  2:6501476/93‑1 (MQ=255)
                                                aaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcg                                        >  2:6317732/1‑93 (MQ=255)
                                                  aGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcg                                      >  2:3155709/1‑93 (MQ=255)
                                                  aGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcg                                      >  1:1022010/1‑93 (MQ=255)
                                                   gCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcgc                                     <  2:3280987/93‑1 (MQ=255)
                                                     cAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACgcgcgc                                   >  1:4180959/1‑93 (MQ=255)
                                                     cAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACgcgcgc                                   >  2:7869495/1‑93 (MQ=255)
                                                     cAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACgcgcgc                                   <  1:6489841/93‑1 (MQ=255)
                                                      aGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGgaga                                               <  1:5806708/80‑1 (MQ=255)
                                                      aGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACgcgcgcg                                  >  1:1408028/1‑92 (MQ=255)
                                                      aGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGct                                  >  1:6613564/1‑93 (MQ=255)
                                                      aGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGct                                  >  1:6085287/1‑93 (MQ=255)
                                                       gACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctg                                 >  1:1601954/1‑93 (MQ=255)
                                                        aCATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgc                                >  2:5185064/1‑93 (MQ=255)
                                                        aCATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgc                                <  1:2685318/93‑1 (MQ=255)
                                                        aCATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgc                                <  1:706047/93‑1 (MQ=255)
                                                         cATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcc                               <  2:5894517/93‑1 (MQ=255)
                                                         cATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcc                               <  1:1474267/93‑1 (MQ=255)
                                                         cATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcc                               <  1:5646295/93‑1 (MQ=255)
                                                         cATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcc                               >  1:4414072/1‑93 (MQ=255)
                                                         cATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcc                               <  1:7683320/93‑1 (MQ=255)
                                                          attattATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcct                              <  2:1665452/93‑1 (MQ=255)
                                                           ttattaTCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcctg                             >  2:5553318/1‑93 (MQ=255)
                                                            tattatCGCGGCGACCCATATGGGGCACTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcctgc                            <  2:2328798/93‑1 (MQ=255)
                                                             attatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCtgcctgcc                           <  1:4750254/93‑1 (MQ=255)
                                                             attatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCtgcctgcc                           >  1:6915958/1‑93 (MQ=255)
                                                                 tCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGAt                       <  2:4509147/93‑1 (MQ=255)
                                                                  cgcgGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATc                      <  1:4291986/93‑1 (MQ=255)
                                                                  cgcgGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATc                      <  2:3000768/93‑1 (MQ=255)
                                                                  cgcgGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATc                      >  2:282833/1‑93 (MQ=255)
                                                                   gcggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCg                     >  1:2167624/1‑93 (MQ=255)
                                                                        gACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGc                <  2:4906499/93‑1 (MQ=255)
                                                                        gACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGc                <  1:5282080/93‑1 (MQ=255)
                                                                        gACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGc                >  2:1745333/1‑93 (MQ=255)
                                                                         aCCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCg               >  1:2049587/1‑93 (MQ=255)
                                                                           ccATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgat             <  2:4770088/93‑1 (MQ=255)
                                                                             atatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatga           >  1:3896528/1‑93 (MQ=255)
                                                                             atatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatga           <  2:3936211/93‑1 (MQ=255)
                                                                              tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          <  1:3569718/93‑1 (MQ=255)
                                                                              tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          >  1:223595/1‑93 (MQ=255)
                                                                              tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          >  2:5873398/1‑93 (MQ=255)
                                                                              tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          >  2:2590423/1‑93 (MQ=255)
                                                                              tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          >  2:6602803/1‑93 (MQ=255)
                                                                               atGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtc         >  1:1157056/1‑93 (MQ=255)
                                                                                tGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcg        >  2:254028/1‑93 (MQ=255)
                                                                                 ggggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgt       <  2:6708828/93‑1 (MQ=255)
                                                                                 ggggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgt       >  2:6112262/1‑93 (MQ=255)
                                                                                 ggggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgt       <  1:882710/93‑1 (MQ=255)
                                                                                  gggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGCTGAtcgtc      >  2:818188/1‑93 (MQ=255)
                                                                                   ggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgtcg     <  2:7654278/93‑1 (MQ=255)
                                                                                   ggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgtcg     <  1:7452902/93‑1 (MQ=255)
                                                                                      aTTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCggtg  >  2:3806194/1‑93 (MQ=255)
                                                                                         |                                                                                         
GCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTG  >  NC_000913/504710‑504888

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