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% 834.3 / ‑5.1 254Y221C (TAC→TGC) ushAUDP‑sugar hydrolase (5'‑nucleotidase)
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (146/107);  total (146/108)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 258 total aligned reads displayed.

AGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTG  >  NC_000913/504709‑504888
                                                                                          |                                                                                         
 gCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGa                                                                                        <  1:5753349/93‑1 (MQ=255)
  cccGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGAt                                                                                       <  2:6545668/93‑1 (MQ=255)
  cccGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGAt                                                                                       >  1:4556279/1‑93 (MQ=255)
    cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     <  2:1386984/93‑1 (MQ=255)
    cgccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATaa                                                                                     >  2:5983317/1‑93 (MQ=255)
     gccgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAAt                                                                                    <  2:5061429/93‑1 (MQ=255)
       cgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATgg                                                                                  >  2:1584148/1‑93 (MQ=255)
       cgATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATgg                                                                                  >  2:5694975/1‑93 (MQ=255)
         aTGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTg                                                                                <  2:1054527/93‑1 (MQ=255)
         aTGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTg                                                                                >  2:615657/1‑93 (MQ=255)
           gaagcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAg                                                                              >  1:1824614/1‑93 (MQ=255)
           gaagcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAg                                                                              >  1:736865/1‑93 (MQ=255)
           gaagcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGAGGCATTGCGATAATGGTGAg                                                                              >  2:5880100/1‑93 (MQ=255)
             agcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCa                                                                            >  2:5415/1‑93 (MQ=255)
             agcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCa                                                                            <  2:4286632/93‑1 (MQ=255)
              gcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCAc                                                                           >  1:3833757/1‑93 (MQ=255)
              gcgaagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCAc                                                                           >  1:1986658/1‑93 (MQ=255)
                 aagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGc                                                                        <  1:1326557/93‑1 (MQ=255)
                 aagcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGc                                                                        <  1:332870/93‑1 (MQ=255)
                   gcTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGctc                                                                      <  2:4670173/93‑1 (MQ=255)
                          aTTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCa                                                               >  1:662212/1‑93 (MQ=255)
                          aTTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCa                                                               >  1:2650813/1‑93 (MQ=255)
                           ttCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCAc                                                              >  1:4715835/1‑93 (MQ=255)
                           ttCAGGAGCGGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCAc                                                              <  1:1481214/93‑1 (MQ=255)
                             cAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCg                                                            <  1:6331124/93‑1 (MQ=255)
                              aGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCgg                                                           >  1:1052798/1‑93 (MQ=255)
                              aGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCgg                                                           <  1:2121120/93‑1 (MQ=255)
                               ggAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCagg                                                          >  1:3169101/1‑93 (MQ=255)
                                gAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGc                                                         <  1:100551/93‑1 (MQ=255)
                                  gCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGa                                                       >  2:1551901/1‑93 (MQ=255)
                                  gCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGa                                                       <  2:5542704/93‑1 (MQ=255)
                                   cTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGAt                                                      >  1:2192592/1‑93 (MQ=255)
                                   cTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGAt                                                      <  1:302107/93‑1 (MQ=255)
                                       aaCAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTgg                                                  >  2:2591347/1‑93 (MQ=255)
                                        acagacagAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGTACGGCTCTAACGCACCGGGCGATGTGga                                                 >  2:6373322/1‑93 (MQ=255)
                                        acagacagAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGga                                                 <  2:5793352/93‑1 (MQ=255)
                                        acagacagAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGga                                                 >  1:3794483/1‑93 (MQ=255)
                                          agacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGgaga                                               >  2:4998470/1‑93 (MQ=255)
                                          agacagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGgaga                                               >  2:177586/1‑93 (MQ=255)
                                            acagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATg                                             >  1:2447579/1‑93 (MQ=255)
                                            acagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATg                                             <  1:1656755/93‑1 (MQ=255)
                                             cagaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATgg                                            <  1:2229026/93‑1 (MQ=255)
                                              agaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGc                                           >  2:199387/1‑93 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  2:6053778/1‑93 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  1:487965/1‑93 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  1:5386709/1‑93 (MQ=255)
                                               gaAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          >  1:6710462/1‑93 (MQ=255)
                                               gaAAAGCCAGACATGATTATCGCGGCGACCCATATGGGGCATTGCGAGAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCa                                          <  2:1409460/93‑1 (MQ=255)
                                                aaaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAc                                         <  1:3652306/93‑1 (MQ=255)
                                                 aaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcg                                        <  1:6636711/93‑1 (MQ=255)
                                                 aaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcg                                        <  2:5290167/93‑1 (MQ=255)
                                                 aaaGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcg                                        >  1:4611763/1‑93 (MQ=255)
                                                   aGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcg                                      <  1:11449/93‑1 (MQ=255)
                                                   aGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcg                                      >  2:2347614/1‑93 (MQ=255)
                                                   aGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCAcgcg                                      <  1:1360466/93‑1 (MQ=255)
                                                      cAGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACgcgcgc                                   <  1:5925002/93‑1 (MQ=255)
                                                       aGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGct                                  <  2:1311948/93‑1 (MQ=255)
                                                       aGACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGct                                  >  2:4494789/1‑93 (MQ=255)
                                                        gACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGGTGTGGAGATGGCACGCGCGctg                                 >  2:4279568/1‑93 (MQ=255)
                                                        gACATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctg                                 >  1:2559501/1‑93 (MQ=255)
                                                         aCATTATTATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgc                                <  2:4114709/93‑1 (MQ=255)
                                                           attattATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcct                              >  2:4561565/1‑93 (MQ=255)
                                                           attattATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcct                              >  2:2948224/1‑93 (MQ=255)
                                                           attattATCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcct                              >  2:3831592/1‑93 (MQ=255)
                                                            ttattaTCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcctg                             >  1:5095690/1‑93 (MQ=255)
                                                            ttattaTCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcctg                             >  2:4410570/1‑93 (MQ=255)
                                                             tattatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGctgcctgc                            <  2:1947850/93‑1 (MQ=255)
                                                              attatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCtgcctgcc                           >  2:4833117/1‑93 (MQ=255)
                                                               ttatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCg                          <  1:5637347/93‑1 (MQ=255)
                                                               ttatCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCg                          >  2:522490/1‑93 (MQ=255)
                                                                 atCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGa                        <  1:1623961/93‑1 (MQ=255)
                                                                  tCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGAt                       >  1:4410264/1‑93 (MQ=255)
                                                                  tCGCGGCGACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGAt                       >  2:2694938/1‑93 (MQ=255)
                                                                  tCGCGGCGACCCATATGGGGCATTGCGATAACGGTGAGCACAGCTCTAACGAACCGGGCGATGTGGAGATGGCACGCGCGATGCATGACGGac                       >  1:2886500/1‑92 (MQ=255)
                                                                      ggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCt                   <  2:4057356/93‑1 (MQ=255)
                                                                      ggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCt                   >  2:6118006/1‑93 (MQ=255)
                                                                      ggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCt                   <  1:1774142/93‑1 (MQ=255)
                                                                      ggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCt                   <  2:2032378/93‑1 (MQ=255)
                                                                      ggcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCt                   <  2:190576/93‑1 (MQ=255)
                                                                       gcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTg                  <  2:2473744/93‑1 (MQ=255)
                                                                       gcgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTg                  >  1:1095978/1‑93 (MQ=255)
                                                                        cgACCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTgg                 <  1:1889061/93‑1 (MQ=255)
                                                                          aCCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCg               <  2:3086839/93‑1 (MQ=255)
                                                                          aCCCATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCg               <  1:5804156/93‑1 (MQ=255)
                                                                           cccATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCga              >  2:3522048/1‑93 (MQ=255)
                                                                            ccATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgat             >  1:2761052/1‑93 (MQ=255)
                                                                             cATATGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatg            >  1:6074051/1‑93 (MQ=255)
                                                                              atatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatga           <  2:1260717/93‑1 (MQ=255)
                                                                              atatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatga           <  1:3103857/93‑1 (MQ=255)
                                                                              atatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatga           >  2:2535503/1‑93 (MQ=255)
                                                                               tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTTGCgatgat          <  2:752783/93‑1 (MQ=255)
                                                                               tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          <  1:5620452/93‑1 (MQ=255)
                                                                               tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          <  1:1748102/93‑1 (MQ=255)
                                                                               tatGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCgatgat          >  1:3286071/1‑93 (MQ=255)
                                                                                atGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtc         >  1:3888445/1‑93 (MQ=255)
                                                                                 tGGGGCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcg        >  2:3063326/1‑93 (MQ=255)
                                                                                  ggggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgt       <  2:3603861/93‑1 (MQ=255)
                                                                                  ggggCATTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGAtcgt       >  1:4328550/1‑93 (MQ=255)
                                                                                       aTTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCggtg  >  2:6281630/1‑93 (MQ=255)
                                                                                       aTTGCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATAGTCggcg  >  1:4025556/1‑91 (MQ=255)
                                                                                          |                                                                                         
AGCCCGCCGATGAAGCGAAGCTGGTGATTCAGGAGCTGCAACAGACAGAAAAGCCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTG  >  NC_000913/504709‑504888

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