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,7990AG98.9% 559.7 / ‑3.2 175Y221C (TAC→TGC) ushAUDP‑sugar hydrolase (5'‑nucleotidase)
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (110/63);  total (112/63)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 180 total aligned reads displayed.

AGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTG  >  NC_000913/504709‑504888
                                                                                          |                                                                                         
aGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCg                                                                                         >  2:2313706/1‑93 (MQ=255)
 gccggccgATGAAGCGAAGCTGGGAATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGa                                                                                        <  2:2596623/93‑1 (MQ=255)
   ccgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCGGACATTATTATCGCGCCGACCCATATGGGGCATTGCGATa                                                                                      >  2:4384995/1‑93 (MQ=255)
    cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     <  2:953007/93‑1 (MQ=255)
    cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     >  2:363188/1‑93 (MQ=255)
    cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     <  2:3932984/93‑1 (MQ=255)
     gccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAAt                                                                                    >  1:1407633/1‑93 (MQ=255)
      ccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATg                                                                                   >  1:1018536/1‑93 (MQ=255)
        gATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGt                                                                                 <  2:2678866/93‑1 (MQ=255)
        gATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGt                                                                                 <  1:540969/93‑1 (MQ=255)
         aTGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTg                                                                                >  2:3711535/1‑93 (MQ=255)
          tGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGa                                                                               >  2:3699899/1‑93 (MQ=255)
          tGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGa                                                                               >  2:1612136/1‑93 (MQ=255)
          tGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGa                                                                               >  1:1062011/1‑93 (MQ=255)
             agcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCa                                                                            <  1:4060979/93‑1 (MQ=255)
             agcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCa                                                                            >  1:3036108/1‑93 (MQ=255)
             agcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCa                                                                            >  1:3041224/1‑93 (MQ=255)
             agcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCa                                                                            <  1:1771166/93‑1 (MQ=255)
             agcgaagcTGGGGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCa                                                                            >  2:1918427/1‑93 (MQ=255)
              gcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCAc                                                                           <  2:2147954/93‑1 (MQ=255)
              gcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGAGGCATTGCGATAATGGTGAGCAc                                                                           >  2:2596010/1‑93 (MQ=255)
               cgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACg                                                                          <  2:875481/93‑1 (MQ=255)
                gaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACgg                                                                         <  2:4455369/93‑1 (MQ=255)
                 aagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGc                                                                        <  2:2528572/93‑1 (MQ=255)
                 aagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGc                                                                        >  2:3136949/1‑93 (MQ=255)
                   gcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGctc                                                                      >  2:2096877/1‑93 (MQ=255)
                   gcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGctc                                                                      <  1:819945/93‑1 (MQ=255)
                   gcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGctc                                                                      <  1:99002/93‑1 (MQ=255)
                   gcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGctc                                                                      >  2:1028276/1‑93 (MQ=255)
                    cTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGctct                                                                     <  2:2493730/93‑1 (MQ=255)
                      ggTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTaa                                                                   <  1:3796422/93‑1 (MQ=255)
                       gTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAAc                                                                  >  2:1398591/1‑93 (MQ=255)
                                      cAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCAGATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCTCCGGGCGAtgtg                                                   <  2:3705653/93‑1 (MQ=255)
                                       aaCAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTgg                                                  >  2:4504686/1‑93 (MQ=255)
                                       aaCAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTgg                                                  >  1:3557867/1‑93 (MQ=255)
                                        acagacagAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGga                                                 >  2:1969177/1‑93 (MQ=255)
                                        acagacagAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGga                                                 <  1:4504884/93‑1 (MQ=255)
                                          agacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGgaga                                               <  2:2267937/93‑1 (MQ=255)
                                          agacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGgaga                                               <  1:3291498/93‑1 (MQ=255)
                                           gacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGAt                                              >  2:2789110/1‑93 (MQ=255)
                                           gacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGAt                                              <  2:1771545/93‑1 (MQ=255)
                                           gacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGAt                                              >  1:993952/1‑93 (MQ=255)
                                            acagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATg                                             <  2:1433498/93‑1 (MQ=255)
                                            acagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATg                                             <  2:4091228/93‑1 (MQ=255)
                                              agaAAAGCTAGACATTATTATCGCGGCGACCGAAATGGGCCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGc                                           <  2:4706989/93‑1 (MQ=255)
                                              agaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGc                                           >  1:4743081/1‑93 (MQ=255)
                                              agaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGc                                           >  1:3562361/1‑93 (MQ=255)
                                              agaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGc                                           >  1:1985163/1‑93 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  1:4706045/1‑93 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  2:3640899/1‑93 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          <  1:604819/93‑1 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGAGGGCa                                          >  1:1302011/1‑93 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGGGGCGATGGcc                                          >  1:3037177/1‑92 (MQ=255)
                                                aaaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAc                                         >  1:2617695/1‑93 (MQ=255)
                                                 aaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcg                                        <  1:972990/93‑1 (MQ=255)
                                                   aGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcg                                      >  2:2944074/1‑93 (MQ=255)
                                                    accAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcgc                                     <  1:4012586/92‑1 (MQ=255)
                                                      cAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACgcgcgc                                   <  1:1937639/93‑1 (MQ=255)
                                                       aGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGct                                  <  1:2974409/93‑1 (MQ=255)
                                                       aGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGct                                  >  1:2631804/1‑93 (MQ=255)
                                                        gACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctg                                 >  2:1055351/1‑93 (MQ=255)
                                                         aCATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgc                                >  2:3952574/1‑93 (MQ=255)
                                                         aCATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgc                                >  1:1676075/1‑93 (MQ=255)
                                                          aattattATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcc                               <  2:1235036/92‑1 (MQ=255)
                                                           attattATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcct                              >  1:2164989/1‑93 (MQ=255)
                                                           attattATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcct                              >  2:3481600/1‑93 (MQ=255)
                                                             tattatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctcactgc                            >  2:846048/1‑93 (MQ=255)
                                                              attatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCtgcctgcc                           >  2:523898/1‑93 (MQ=255)
                                                               ttatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCg                          >  2:3363544/1‑93 (MQ=255)
                                                               ttatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCg                          >  2:1388836/1‑93 (MQ=255)
                                                                 atCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGa                        >  1:1024081/1‑93 (MQ=255)
                                                                 atCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGa                        <  2:774213/93‑1 (MQ=255)
                                                                 atCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGa                        >  2:3511915/1‑93 (MQ=255)
                                                                  tCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGAt                       >  2:3526061/1‑93 (MQ=255)
                                                                  tCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGAt                       >  2:4234286/1‑93 (MQ=255)
                                                                  tCGCGGCGACCCATATGGGGCATTGCGATAATGGCGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGAt                       >  2:474566/1‑93 (MQ=255)
                                                                   cgcgGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATc                      <  1:76777/93‑1 (MQ=255)
                                                                    gcggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCg                     >  1:2352816/1‑93 (MQ=255)
                                                                    gcggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCg                     >  1:2895141/1‑93 (MQ=255)
                                                                     cggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGc                    >  2:288707/1‑93 (MQ=255)
                                                                      ggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCt                   >  1:4642823/1‑93 (MQ=255)
                                                                      ggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCAGGCTCGCt                   >  1:3627305/1‑93 (MQ=255)
                                                                      ggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCGCCCTGCCGGATCgtg                   >  2:3891451/1‑91 (MQ=255)
                                                                         gACCCATATGGGGCATTGCGATAATGTTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTgtc                <  1:1143635/93‑3 (MQ=255)
                                                                           cccATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCga              >  2:2709145/1‑93 (MQ=255)
                                                                            ccATATGGTGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgat             <  1:2708957/93‑1 (MQ=255)
                                                                            ccATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCggt             >  1:3094497/1‑91 (MQ=255)
                                                                            ccATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgat             >  1:3251894/1‑93 (MQ=255)
                                                                            ccATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgat             >  1:1586125/1‑93 (MQ=255)
                                                                            ccATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCNGCCTGCCGGATCGCTGGCgan             >  1:4122482/1‑92 (MQ=255)
                                                                             cATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatg            >  1:2207042/1‑93 (MQ=255)
                                                                             cATATGGGCCATTGCGATAATGGTGAGCACGGCTCTAACGCGCCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatg            <  1:3139219/93‑1 (MQ=255)
                                                                              atatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTTGCgatga           >  2:803647/1‑93 (MQ=255)
                                                                               tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          >  1:1868196/1‑93 (MQ=255)
                                                                               tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          <  1:716661/93‑1 (MQ=255)
                                                                               tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          >  2:3815778/1‑93 (MQ=255)
                                                                                   gggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgtc      <  2:1564507/93‑1 (MQ=255)
                                                                                       aTTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCggtg  <  2:1085358/93‑1 (MQ=255)
                                                                                       aTTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCggtg  >  2:1723667/1‑93 (MQ=255)
                                                                                       aTTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCggtg  <  2:2783244/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTG  >  NC_000913/504709‑504888

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: