Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,007,730 T→C W316R (TGG→CGG)  ygeY → putative deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,007,7300TC99.7% 1253.5 / NA 361W316R (TGG→CGG) ygeYputative deacetylase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (169/191);  total (170/191)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 370 total aligned reads displayed.

CAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTGGA  >  NC_000913/3007640‑3007820
                                                                                          |                                                                                          
 aGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGa                                                                                         <  1:100398/93‑1 (MQ=255)
 aGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGa                                                                                         >  1:8313686/1‑93 (MQ=255)
  gAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGaa                                                                                        <  2:4960656/93‑1 (MQ=255)
   aaaGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGaaa                                                                                       >  1:5671456/1‑93 (MQ=255)
    aaGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAg                                                                                      <  2:1495386/93‑1 (MQ=255)
      gCTAACGCTGTTGTTTCTATGTACAACTACGACCGGCCGTCCTGGACTGGCCTGGTTTGCCCAACCGAATGCTACTTCCCGACCCGGAAAGTg                                                                                    <  2:2291802/93‑1 (MQ=255)
       gtAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTgg                                                                                   <  1:6609564/92‑1 (MQ=255)
        tAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGga                                                                                  >  1:2032250/1‑93 (MQ=255)
        tAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGGTTACCCAACCGAATGCTACTTCCCGACCCGGCAAGTGga                                                                                  >  1:4877740/1‑93 (MQ=255)
         aaCGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGgaa                                                                                 <  2:6924909/93‑1 (MQ=255)
          aCGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGgaag                                                                                >  1:4764675/1‑93 (MQ=255)
          aCGCGGTTGTTTCTATGTACAAATACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAACGCTACTTTCCGCCCCGGATAGCTgaag                                                                                >  2:2451785/1‑93 (MQ=255)
           cgcgGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGgaaga                                                                               <  1:342047/93‑1 (MQ=255)
             cgGTTGGTTCTATGTACAACTACGACCGTCCGTCCTTGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGaagaag                                                                             <  1:7370350/93‑1 (MQ=255)
              ggTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAagaaga                                                                            >  1:2210696/1‑93 (MQ=255)
              ggTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAagaaga                                                                            >  1:4548084/1‑93 (MQ=255)
              ggTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAagaaga                                                                            <  2:63090/93‑1 (MQ=255)
                 tgttTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCa                                                                         <  1:1498897/93‑1 (MQ=255)
                       tATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCAc                                                                   >  1:5852271/1‑93 (MQ=255)
                       tATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCAc                                                                   >  2:1665717/1‑93 (MQ=255)
                        aTGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCAcg                                                                  >  2:8307987/1‑92 (MQ=255)
                        aTGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCAcc                                                                  <  2:1950588/93‑1 (MQ=255)
                         tGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCg                                                                 >  1:6892629/1‑93 (MQ=255)
                         tGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCg                                                                 <  2:8549993/93‑1 (MQ=255)
                          gTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGt                                                                <  2:2210497/93‑1 (MQ=255)
                          gTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGt                                                                <  2:5764111/93‑1 (MQ=255)
                          gTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGt                                                                <  2:5209819/93‑1 (MQ=255)
                          gTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAAAACACTACAACGt                                                                >  1:7183872/1‑93 (MQ=255)
                             cAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTaa                                                             <  2:1693363/93‑1 (MQ=255)
                               aCTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTaaac                                                           >  2:6329060/1‑92 (MQ=255)
                                 tACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGCAGAAGATCACTTCACCGTTAAAGCa                                                         >  1:1265441/1‑93 (MQ=255)
                                 tACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCa                                                         >  2:762201/1‑93 (MQ=255)
                                  aCGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCAc                                                        >  2:4549386/1‑93 (MQ=255)
                                   cGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACt                                                       <  2:1078371/93‑1 (MQ=255)
                                     aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                     >  1:5293845/1‑93 (MQ=255)
                                     aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                     >  2:3036851/1‑93 (MQ=255)
                                     aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                     >  1:3548467/1‑93 (MQ=255)
                                     aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                     >  1:1861054/1‑93 (MQ=255)
                                     aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                     >  2:8148686/1‑93 (MQ=255)
                                     aCCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTgg                                                     >  1:4953489/1‑93 (MQ=255)
                                      ccgtccgtCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGt                                                    <  2:3158752/93‑1 (MQ=255)
                                         tccgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGaa                                                 >  1:6469031/1‑93 (MQ=255)
                                          ccgtccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAAt                                                >  2:1195862/1‑93 (MQ=255)
                                             tccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGcc                                             <  2:5576977/93‑1 (MQ=255)
                                             tccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGcc                                             <  2:1256671/93‑1 (MQ=255)
                                              ccTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCt                                            <  1:3696362/93‑1 (MQ=255)
                                               cTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTa                                           <  2:2737096/93‑1 (MQ=255)
                                                tGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTAc                                          <  2:3365898/93‑1 (MQ=255)
                                                  gACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGa                                        >  1:60076/1‑93 (MQ=255)
                                                   gcTGGCCTGGTTTACCCAACCGATTGCTACTACCCGACCCGGAAAGTGGAAGAAGAGCACTTCACCGTTAAAGCACTGGTGAATGCCTACGaa                                       <  2:4756068/92‑1 (MQ=255)
                                                    cTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAg                                      <  2:5982448/93‑1 (MQ=255)
                                                    cTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAg                                      >  2:2639410/1‑93 (MQ=255)
                                                    cTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAg                                      <  1:2330742/93‑1 (MQ=255)
                                                    cTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAg                                      <  2:2092995/93‑1 (MQ=255)
                                                     tGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAgg                                     >  1:7898335/1‑93 (MQ=255)
                                                     tGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAgg                                     >  1:4242579/1‑93 (MQ=255)
                                                     tGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAgg                                     <  1:1456917/93‑1 (MQ=255)
                                                      ggCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGCAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAggg                                    >  1:7793172/1‑93 (MQ=255)
                                                       gCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGc                                   >  1:8823364/1‑93 (MQ=255)
                                                       gCCTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGc                                   <  2:554226/93‑1 (MQ=255)
                                                        ccTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCt                                  <  1:519818/93‑1 (MQ=255)
                                                         cTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTg                                 <  2:3964635/93‑1 (MQ=255)
                                                         cTGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTg                                 >  2:8423709/1‑93 (MQ=255)
                                                          tGGTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGt                                >  1:770277/1‑93 (MQ=255)
                                                           tgTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGtt                               <  1:3213805/92‑1 (MQ=255)
                                                           ggTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGtt                               <  2:4834625/93‑1 (MQ=255)
                                                            gTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGttt                              <  2:4151111/93‑1 (MQ=255)
                                                            gTTTACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGttt                              <  1:3327346/93‑1 (MQ=255)
                                                             tttACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTg                             <  2:3531879/93‑1 (MQ=255)
                                                             tttACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTg                             <  1:2995530/93‑1 (MQ=255)
                                                              ttACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTgg                            <  2:7171224/93‑1 (MQ=255)
                                                              ttACCCAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTgg                            <  1:7016392/93‑1 (MQ=255)
                                                                  ccAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCaaa                        <  2:7403011/93‑1 (MQ=255)
                                                                  ccAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCaaa                        <  1:3617051/93‑1 (MQ=255)
                                                                   cAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAg                       <  1:7996983/93‑1 (MQ=255)
                                                                   cAACCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAg                       >  2:6126886/1‑93 (MQ=255)
                                                                     aCCGAATGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAgcg                     >  1:8452557/1‑93 (MQ=255)
                                                                           tGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgtt               <  2:6680044/93‑1 (MQ=255)
                                                                           tGCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgtt               <  2:6513855/93‑1 (MQ=255)
                                                                            gCTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgttg              >  2:5849523/1‑93 (MQ=255)
                                                                             cTACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgttgt             <  2:8665908/93‑1 (MQ=255)
                                                                              tACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgttgtt            <  1:6175978/93‑1 (MQ=255)
                                                                              tACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgttgtt            <  2:8094791/93‑1 (MQ=255)
                                                                              tACTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGgttgtt            <  2:3572987/93‑1 (MQ=255)
                                                                               aCTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGttgttg           >  1:1108979/1‑93 (MQ=255)
                                                                               aCTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGttgttg           >  1:5411985/1‑93 (MQ=255)
                                                                                cTTCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGGGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGa          >  1:6087632/1‑93 (MQ=255)
                                                                                 ttCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGAt         >  2:2372008/1‑93 (MQ=255)
                                                                                  tCCCGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATa        >  2:3673719/1‑93 (MQ=255)
                                                                                   cccGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATaa       <  1:156553/93‑1 (MQ=255)
                                                                                    ccGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAg      >  2:534356/1‑93 (MQ=255)
                                                                                    ccGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAg      <  1:4326750/93‑1 (MQ=255)
                                                                                    ccGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAg      >  1:3926073/1‑93 (MQ=255)
                                                                                    ccGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAg      <  1:2622901/93‑1 (MQ=255)
                                                                                     cGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATacgg     >  2:883441/1‑90 (MQ=255)
                                                                                     cGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGt     <  1:281007/93‑1 (MQ=255)
                                                                                     cGACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGt     <  1:249063/93‑1 (MQ=255)
                                                                                      gACCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTg    >  1:4020921/1‑93 (MQ=255)
                                                                                       aCCCGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGATGTTGCTAAGTgg   >  2:765968/1‑93 (MQ=255)
                                                                                        tcTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGGTGATTAGTGGa  >  2:7835426/2‑93 (MQ=255)
                                                                                          |                                                                                          
CAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTGGA  >  NC_000913/3007640‑3007820

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