Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,015,142 T→C S189S (AGT→AGC fliI → flagellum‑specific ATP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,015,1420TC94.7% 477.4 / 8.5 152S189S (AGT→AGCfliIflagellum‑specific ATP synthase
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base C (70/74);  total (74/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 159 total aligned reads displayed.

GGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATT  >  NC_000913/2015052‑2015230
                                                                                          |                                                                                        
ggACACCGGCCTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAgcgt                                                                                        <  2:3514710/93‑1 (MQ=255)
  acacCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGtgc                                                                                      <  1:302628/93‑1 (MQ=255)
   cacCGCCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGtgct                                                                                     <  2:4290151/93‑1 (MQ=255)
    acCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGtgctg                                                                                    <  1:5194967/93‑1 (MQ=255)
      cGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTgctgct                                                                                  <  2:5711484/93‑1 (MQ=255)
      cGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTgctgct                                                                                  <  1:1711128/93‑1 (MQ=255)
       ggCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGctgctg                                                                                 >  2:3777397/1‑93 (MQ=255)
        gCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTgg                                                                                >  1:1619298/1‑93 (MQ=255)
          gTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTgggg                                                                              <  1:5977884/93‑1 (MQ=255)
           tGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTGCGGCGTTGGTAAAAGCGTGCTGCTggggg                                                                             >  1:4938375/1‑92 (MQ=255)
           tGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGga                                                                             >  2:1496109/1‑93 (MQ=255)
           tGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGga                                                                             <  1:3934781/93‑1 (MQ=255)
            gcgcCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGgat                                                                            <  2:1146235/93‑1 (MQ=255)
              gcCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGgatga                                                                          <  2:313284/93‑1 (MQ=255)
               cccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGgatgag                                                                         >  2:1665576/1‑92 (MQ=255)
                ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGatgatg                                                                        <  1:5842763/93‑1 (MQ=255)
                ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGatgatg                                                                        <  1:3693954/93‑1 (MQ=255)
                 caatcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATgg                                                                       >  1:3573372/1‑93 (MQ=255)
                  aatcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGc                                                                      >  2:2898174/1‑93 (MQ=255)
                  aatcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGc                                                                      >  2:1296254/1‑93 (MQ=255)
                  aatcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGc                                                                      >  1:132029/1‑93 (MQ=255)
                  aatcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGc                                                                      <  2:400926/93‑1 (MQ=255)
                   atcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATggag                                                                     >  2:1238544/1‑91 (MQ=255)
                   atcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATggaa                                                                     >  2:2284160/1‑91 (MQ=255)
                     caatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACg                                                                   >  1:1971619/1‑93 (MQ=255)
                     caatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACg                                                                   <  1:4168807/93‑1 (MQ=255)
                      aatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGt                                                                  >  1:4284979/1‑93 (MQ=255)
                      aatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGt                                                                  >  1:2265763/1‑93 (MQ=255)
                      aatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGt                                                                  <  1:3041230/93‑1 (MQ=255)
                       atGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGtt                                                                 >  2:2037867/1‑93 (MQ=255)
                        tGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTa                                                                >  1:6130355/1‑93 (MQ=255)
                         gCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTac                                                               <  1:3822595/93‑1 (MQ=255)
                           ccTGCTTACCGTTGGGCGTGGGCAGCGTATGTGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGTAGGATGGCACGTTacac                                                             >  2:5476374/1‑93 (MQ=255)
                           ccTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTacac                                                             >  2:5106825/1‑93 (MQ=255)
                           ccTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTacac                                                             <  2:268385/93‑1 (MQ=255)
                           ccTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTacac                                                             >  2:6128607/1‑93 (MQ=255)
                            cTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACAcc                                                            <  2:5880152/93‑1 (MQ=255)
                            cTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACAcc                                                            <  1:1546303/93‑1 (MQ=255)
                             tGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACAccc                                                           >  1:435594/1‑93 (MQ=255)
                             tGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACAccc                                                           >  2:6007692/1‑93 (MQ=255)
                             tGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACAccc                                                           >  1:3160477/1‑93 (MQ=255)
                              gCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCcg                                                          <  1:2913779/93‑1 (MQ=255)
                                     gTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGAt                                                   >  1:4394707/1‑93 (MQ=255)
                                      ttGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATg                                                  <  2:2491377/93‑1 (MQ=255)
                                          gCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCAt                                              <  1:5846930/93‑1 (MQ=255)
                                            gTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTg                                            <  1:1600194/93‑1 (MQ=255)
                                             ggggCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTCCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGt                                           <  1:3049642/92‑1 (MQ=255)
                                              gggCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTc                                          <  1:2937597/93‑1 (MQ=255)
                                                gCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGt                                        <  2:775099/93‑1 (MQ=255)
                                                   gCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTggg                                     <  2:4637974/93‑1 (MQ=255)
                                                   gCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTggg                                     >  2:5559728/1‑93 (MQ=255)
                                                     gTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTgggct                                   >  2:3747566/1‑91 (MQ=255)
                                                      tATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGttt                                  >  2:1494876/1‑93 (MQ=255)
                                                      tATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGttt                                  >  1:4524654/1‑93 (MQ=255)
                                                      tATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGttt                                  >  1:5705096/1‑93 (MQ=255)
                                                      tATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGttt                                  >  2:881085/1‑93 (MQ=255)
                                                      tATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGttt                                  >  1:3466938/1‑93 (MQ=255)
                                                       aTGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTg                                 >  1:2687551/1‑93 (MQ=255)
                                                        tGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGa                                >  2:4485204/1‑93 (MQ=255)
                                                          gggCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGAtt                              <  2:976702/93‑1 (MQ=255)
                                                          gggCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGAtt                              <  2:3938043/93‑1 (MQ=255)
                                                           ggCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTg                             >  1:1933735/1‑93 (MQ=255)
                                                            gCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTgg                            <  1:2174652/93‑1 (MQ=255)
                                                            gCTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTgg                            <  2:2735504/93‑1 (MQ=255)
                                                             cTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTggg                           >  1:5025248/1‑92 (MQ=255)
                                                             cTGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGt                           <  2:4777840/93‑1 (MQ=255)
                                                              tGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTg                          >  2:4893040/1‑93 (MQ=255)
                                                              tGTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTg                          >  1:4114939/1‑93 (MQ=255)
                                                              agTTTGCCGGGTCCGGCGTTGGTAAAAGCGTGCAGCTGGGGATGAAGGCACGTTACACCCACGACGATGTCATTGTCGTGGGTTTGATTGGTg                          <  1:4763245/92‑1 (MQ=255)
                                                                tttGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGaa                        <  1:2750778/93‑1 (MQ=255)
                                                                  tGCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACg                      <  1:4444199/93‑1 (MQ=255)
                                                                   gCCGGGTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGt                     <  2:3069802/93‑1 (MQ=255)
                                                                    ccGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTg                    >  1:3400250/1‑93 (MQ=255)
                                                                      gggTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTggg                  <  1:1119804/93‑1 (MQ=255)
                                                                      gggTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTggg                  <  1:2526044/93‑1 (MQ=255)
                                                                      gggTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTggg                  <  1:5223687/93‑1 (MQ=255)
                                                                       ggTCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGgc                 >  2:5824431/1‑93 (MQ=255)
                                                                         tCCGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGgcgc               >  2:5265982/1‑93 (MQ=255)
                                                                         tCCGGCGTTGGTAAAAGCGGGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGgcgc               >  2:3956405/1‑93 (MQ=255)
                                                                          ccGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGcgcg              >  2:5536771/1‑93 (MQ=255)
                                                                          ccGGCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGcgcg              >  1:748841/1‑93 (MQ=255)
                                                                            ggCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGaa            <  2:4545354/93‑1 (MQ=255)
                                                                            ggCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGaa            >  1:4832718/1‑93 (MQ=255)
                                                                             gCGTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAg           >  1:3299906/1‑93 (MQ=255)
                                                                              ggTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGt          <  2:4297667/92‑1 (MQ=255)
                                                                               gTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAgaa         >  1:92488/1‑91 (MQ=255)
                                                                               gTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTa         >  2:352549/1‑93 (MQ=255)
                                                                               gTTGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTa         >  1:4013704/1‑93 (MQ=255)
                                                                                ttGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaa        <  1:4884048/93‑1 (MQ=255)
                                                                                ttGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaa        <  2:4707815/93‑1 (MQ=255)
                                                                                ttGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaa        >  2:2133435/1‑93 (MQ=255)
                                                                                ttGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaa        >  2:2524669/1‑93 (MQ=255)
                                                                                ttGGTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaa        <  1:1663343/93‑1 (MQ=255)
                                                                                  ggTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaaaa      <  2:4949783/93‑1 (MQ=255)
                                                                                  ggTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaaaa      <  2:1793291/93‑1 (MQ=255)
                                                                                   gTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAg     <  2:1229894/93‑1 (MQ=255)
                                                                                   gTAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAg     <  2:360769/93‑1 (MQ=255)
                                                                                    tAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGa    <  1:3200847/93‑1 (MQ=255)
                                                                                    tAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGa    <  1:1775332/93‑1 (MQ=255)
                                                                                    tAAAAGCGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGa    >  1:1094831/1‑93 (MQ=255)
                                                                                          |                                                                                        
GGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATT  >  NC_000913/2015052‑2015230

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