Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,098,512 C→A 66.5% R259R (CGC→CGA rhaS → transcriptional activator of rhaBAD and rhaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,098,5120CA66.5% 149.0 / 135.5 144R259R (CGC→CGArhaStranscriptional activator of rhaBAD and rhaT
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base A (46/48);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 8.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 146 total aligned reads displayed.

GAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTG  >  NC_000913/4098413‑4098607
                                                                                                   |                                                                                               
gaAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGc                                                                                                 <  1:1664353/100‑1 (MQ=255)
  aaGCCCGACATCTTCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACg                                                                                               <  1:1101209/100‑1 (MQ=255)
  aaGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACg                                                                                               <  1:2259047/100‑1 (MQ=255)
  aaGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACg                                                                                               <  1:1100314/100‑1 (MQ=255)
    gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCgag                                                                                             <  1:3111353/100‑1 (MQ=255)
    gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                             <  2:2082566/100‑1 (MQ=255)
    gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                             <  2:2335291/100‑1 (MQ=255)
      ccGACATCTGCTACGCCACAGCTAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACg                                                                                               <  1:1293817/96‑1 (MQ=255)
      ccGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGagag                                                                                           <  2:893852/100‑1 (MQ=255)
      ccGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGagag                                                                                           >  2:673130/1‑100 (MQ=255)
       cgggatgtGCTACGCCCCAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGagag                                                                                           <  1:3169903/92‑1 (MQ=255)
             cTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACt                                                                                    >  1:2531074/1‑100 (MQ=255)
             cTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACt                                                                                    <  2:1476345/100‑1 (MQ=255)
              tGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTg                                                                                   <  1:177271/100‑1 (MQ=255)
               gCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTgg                                                                                  <  2:24774/100‑1 (MQ=255)
                cTACGCCACAGCGAGGCCAGTGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAATTGGt                                                                                 >  2:1105784/1‑100 (MQ=255)
                cTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGCGTTTAACt                                                                                    >  1:2839742/1‑97 (MQ=255)
                 tACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTc                                                                                <  1:293747/100‑1 (MQ=255)
                  aCGCCACATCGAGGCAAGCTTTACTGACATCGCCTATCCCTGTAGATTCAGCGATCTGAACCACTTTACGACTCTATTTCGAAGAAAGTTTAACTgg                                                                                  <  1:1029169/97‑1 (MQ=255)
                  aCGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCa                                                                               <  1:1183593/100‑1 (MQ=255)
                  aCGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCa                                                                               <  2:1184545/100‑1 (MQ=255)
                   cGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCAc                                                                              <  1:535371/100‑1 (MQ=255)
                     ccACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcg                                                                            >  2:434771/1‑100 (MQ=255)
                       acaGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACcgcg                                                                          <  2:3392365/100‑1 (MQ=255)
                       acaGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcgcg                                                                          <  2:247625/100‑1 (MQ=255)
                        caGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGt                                                                         <  2:3083019/100‑1 (MQ=255)
                          gCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCCCCTTTCGACGCGCTCTCGACGAGAGTTTAACTCGTCACCGCGTGa                                                                       >  1:102841/1‑100 (MQ=255)
                           cGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGat                                                                      <  2:3217467/100‑1 (MQ=255)
                           cGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGat                                                                      >  2:3135332/1‑100 (MQ=255)
                              ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                   >  1:2245705/1‑100 (MQ=255)
                              ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                   >  2:821709/1‑100 (MQ=255)
                              ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                   >  1:1999141/1‑100 (MQ=255)
                              ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                   >  2:1646452/1‑100 (MQ=255)
                              ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                   >  2:1927611/1‑100 (MQ=255)
                               gCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTTACCGCGTGATATTc                                                                  >  1:1568698/1‑100 (MQ=255)
                                ccAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGCGTTTAACTGGTCACCGCGTGATATTCg                                                                 >  1:3036693/1‑100 (MQ=255)
                                 cAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTTTTTTTTTTTTTAACTGGTCACCGCGTGATAtt                                                                   >  1:3301742/1‑97 (MQ=255)
                                   gCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTCAACTGGTCACCGCGTGATATTCGCCa                                                              <  1:3062543/100‑1 (MQ=255)
                                   gCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCa                                                              <  2:2042965/100‑1 (MQ=255)
                                      ttACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAggg                                                           >  2:2729397/1‑100 (MQ=255)
                                      ttACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAggg                                                           <  1:2738107/100‑1 (MQ=255)
                                       tACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGa                                                          >  2:3171349/1‑100 (MQ=255)
                                         cTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACg                                                        <  1:2350570/100‑1 (MQ=255)
                                          tGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACgg                                                       >  2:1234904/1‑100 (MQ=255)
                                           gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                      >  1:339133/1‑100 (MQ=255)
                                           gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                      >  2:533664/1‑100 (MQ=255)
                                           gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                      >  2:1589432/1‑100 (MQ=255)
                                             cATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCAAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGAt                                                    >  1:2424488/1‑100 (MQ=255)
                                             cATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGAt                                                    >  1:872622/1‑100 (MQ=255)
                                             cATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGAt                                                    >  2:2674150/1‑100 (MQ=255)
                                                cGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                                 >  2:1993162/1‑100 (MQ=255)
                                                cGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                                 >  2:1993045/1‑100 (MQ=255)
                                                 gCCTCTCGCTGTGGATTCAGCGACAGTAACCTCGTTGCGACGCTGTATCTACTTGAGTTTACTTGTTCACcgcg                                                                          <  1:3273640/74‑1 (MQ=255)
                                                 gCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCt                                                >  2:1521431/1‑100 (MQ=255)
                                                     aacGCTGTGGATTGAGCGACAGTAACCACTTATCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCt                                            >  2:2284946/3‑100 (MQ=255)
                                                      tCGCTGTGGATTCAGCGACAGTAGCCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTg                                           <  1:2878522/100‑1 (MQ=255)
                                                      tCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTg                                           <  2:1550422/100‑1 (MQ=255)
                                                      tCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTg                                           <  1:1004606/100‑1 (MQ=255)
                                                       cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                          <  1:2141700/100‑1 (MQ=255)
                                                        gCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCa                                         <  2:291731/100‑1 (MQ=255)
                                                        gCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCa                                         >  2:651512/1‑100 (MQ=255)
                                                         cTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCaa                                        <  2:1487922/100‑1 (MQ=255)
                                                           gtgGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATa                                      >  1:3391695/1‑100 (MQ=255)
                                                           gaatccacAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATa                                      >  1:1634685/8‑100 (MQ=255)
                                                             ggATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAc                                    <  1:1918487/100‑1 (MQ=255)
                                                             ggATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAc                                    <  1:1917951/100‑1 (MQ=255)
                                                             ggATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACAGGATGGCTTTCTGCAATAAc                                    >  2:811003/1‑100 (MQ=255)
                                                              gATTCGGCGACAGTAACCACTTTTCGACGTTTTTTCGCCGGTAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcg                                   <  1:919553/100‑1 (MQ=255)
                                                              gATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcg                                   >  1:2124019/1‑100 (MQ=255)
                                                              gATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcg                                   >  2:180188/1‑100 (MQ=255)
                                                               aTTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACAGGATGGCTTTCTGCAATAAcgc                                  >  1:1812508/1‑100 (MQ=255)
                                                                ttCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgcg                                 <  2:2763428/100‑1 (MQ=255)
                                                                ttCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgcg                                 <  2:2568420/100‑1 (MQ=255)
                                                                ttCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgag                                 >  1:895731/1‑98 (MQ=255)
                                                                   aGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                              >  1:1329684/1‑100 (MQ=255)
                                                                   aGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                              >  1:1455553/1‑100 (MQ=255)
                                                                   aGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                              >  2:1270346/1‑100 (MQ=255)
                                                                   aGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                              >  2:1723815/1‑100 (MQ=255)
                                                                    gCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtc                             >  2:1929556/1‑100 (MQ=255)
                                                                    gCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtc                             >  1:3025525/1‑100 (MQ=255)
                                                                    gCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtc                             <  1:1152664/100‑1 (MQ=255)
                                                                    gCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtc                             <  1:1927598/100‑1 (MQ=255)
                                                                     cGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtct                            >  2:853530/1‑100 (MQ=255)
                                                                     cGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                              <  2:1287189/98‑1 (MQ=255)
                                                                      gACAGTAACCACTTTTCGACGCTTTTTTTTTTTTTTTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtctt                           >  1:2001086/1‑100 (MQ=255)
                                                                        cAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttct                         <  1:240045/100‑1 (MQ=255)
                                                                         aGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTtctc                        <  1:831871/100‑1 (MQ=255)
                                                                         aGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTtctc                        <  1:2965175/100‑1 (MQ=255)
                                                                          gTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                       <  2:2218038/100‑1 (MQ=255)
                                                                          gTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                       >  2:1084922/1‑100 (MQ=255)
                                                                            aaCCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACAGGATGGCTTTCTGCAATAACGCGAATCTTCTCAAc                     <  1:811000/100‑1 (MQ=255)
                                                                                 cTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAtt                <  2:2919473/100‑1 (MQ=255)
                                                                                    ttCGCCGCTTTTTCGCCGAGAGTTTATCTCGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGt             <  1:841697/100‑1 (MQ=255)
                                                                                     tCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTa            <  1:760778/100‑1 (MQ=255)
                                                                                      cGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc           >  1:2496351/1‑100 (MQ=255)
                                                                                      cGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTAc           <  1:2655362/100‑1 (MQ=255)
                                                                                       gACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACg          >  2:715007/1‑100 (MQ=255)
                                                                                         cGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGcc        >  1:581699/1‑100 (MQ=255)
                                                                                             ttttCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCatat    >  2:1609307/1‑100 (MQ=255)
                                                                                               ttCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTg  <  2:3275144/100‑1 (MQ=255)
                                                                                                   |                                                                                               
GAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTG  >  NC_000913/4098413‑4098607

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