Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,989,898 A→G Q378R (CAG→CGG)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,989,8980AG100.0% 681.6 / NA 212Q378R (CAG→CGG) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (94/118);  total (94/118)

Only 100 of 226 total aligned reads displayed.

GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987
                                                                                          |                                                                                         
gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         >  2:3386584/1‑93 (MQ=255)
gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         <  1:4691442/93‑1 (MQ=255)
gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         <  2:2655719/93‑1 (MQ=255)
gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         <  2:4944449/93‑1 (MQ=255)
gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         >  2:5139628/1‑93 (MQ=255)
  gcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTg                                                                                       >  1:3790604/1‑93 (MQ=255)
  gcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTg                                                                                       <  2:1816661/93‑1 (MQ=255)
  gcAGTGTTGAGCCAGTTAGTGAGCGAGGGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTg                                                                                       >  2:2947455/1‑93 (MQ=255)
   cAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGc                                                                                      >  2:2136527/1‑93 (MQ=255)
   cAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGc                                                                                      <  2:3572240/93‑1 (MQ=255)
     gtgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGt                                                                                    <  2:1084823/93‑1 (MQ=255)
      tgtTGAGCCAGCTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTg                                                                                   >  2:5863708/1‑93 (MQ=255)
       gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  <  1:1123422/93‑1 (MQ=255)
       gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  <  2:408178/93‑1 (MQ=255)
       gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  <  1:348762/93‑1 (MQ=255)
        ttGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAg                                                                                 <  1:4238395/93‑1 (MQ=255)
        ttGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAg                                                                                 <  1:5391678/93‑1 (MQ=255)
          gAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTTCGTGAGgc                                                                               <  2:6156779/93‑1 (MQ=255)
            gCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAGCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGgcgc                                                                             <  1:6112415/93‑1 (MQ=255)
              cAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcac                                                                           >  2:5203808/1‑93 (MQ=255)
               aGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcaca                                                                          <  1:4377613/93‑1 (MQ=255)
                gTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACaa                                                                         >  2:2517222/1‑93 (MQ=255)
                gTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACaa                                                                         >  1:5547411/1‑93 (MQ=255)
                gTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACaa                                                                         >  2:843972/1‑93 (MQ=255)
                 ttAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAAc                                                                        <  2:4512538/93‑1 (MQ=255)
                 ttAGTGAGCGAGCGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcaccac                                                                        >  1:550894/1‑93 (MQ=255)
                    gTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAg                                                                     <  2:4393656/93‑1 (MQ=255)
                     tGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGt                                                                    <  2:432127/93‑1 (MQ=255)
                        gCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGc                                                                 <  2:5675130/93‑1 (MQ=255)
                           aGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCg                                                              >  1:924568/1‑93 (MQ=255)
                            atGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             >  1:5943708/2‑93 (MQ=255)
                              ggggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcg                                                           <  2:5591759/93‑1 (MQ=255)
                               gggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcgc                                                          <  1:5588434/93‑1 (MQ=255)
                                 gTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCga                                                        <  2:1915165/93‑1 (MQ=255)
                                 gTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCga                                                        <  1:2946014/93‑1 (MQ=255)
                                  ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       >  1:3659089/1‑93 (MQ=255)
                                   tGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatg                                                      <  1:5111195/93‑1 (MQ=255)
                                    gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     <  2:1293653/93‑1 (MQ=255)
                                    gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     <  1:1360329/93‑1 (MQ=255)
                                     ggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatgat                                                    <  1:2507402/93‑1 (MQ=255)
                                      gACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGatgatg                                                   <  2:4107912/93‑1 (MQ=255)
                                      gACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGatgatg                                                   <  1:5644868/93‑1 (MQ=255)
                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  2:2854273/93‑1 (MQ=255)
                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  2:3041020/93‑1 (MQ=255)
                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  1:1045416/93‑1 (MQ=255)
                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 >  1:3096443/1‑93 (MQ=255)
                                           aGCTCGTCTGGCAATGCTCGCTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGAGCACAACGAGTTGCTCGACGCGATGATGCagag                                              >  2:2701899/1‑93 (MQ=255)
                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              >  2:1446755/1‑93 (MQ=255)
                                              tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTa                                           >  2:5051272/1‑93 (MQ=255)
                                              tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTa                                           >  1:1448231/1‑93 (MQ=255)
                                              tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTa                                           >  2:4710243/1‑93 (MQ=255)
                                              tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTa                                           <  1:1661432/93‑1 (MQ=255)
                                               cGTCTGCCAATGCTCGATAACCGCGCTAACTGGAAGATGGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTAc                                          <  1:1106513/93‑1 (MQ=255)
                                                gTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTAcc                                         >  2:5779439/1‑93 (MQ=255)
                                                 tCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCg                                        >  1:3754296/1‑93 (MQ=255)
                                                  cTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGt                                       >  1:5126519/1‑93 (MQ=255)
                                                  cTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGt                                       <  2:1437937/93‑1 (MQ=255)
                                                   tGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCCGAGCTACCGTa                                      >  2:5473921/1‑93 (MQ=255)
                                                   tGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTa                                      >  1:3020249/1‑93 (MQ=255)
                                                     gCAATGCTCGACAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAAt                                    <  1:975888/93‑1 (MQ=255)
                                                       aaTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCt                                  >  1:4077908/1‑93 (MQ=255)
                                                        aTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTg                                 <  1:4026349/93‑1 (MQ=255)
                                                        aTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTg                                 <  1:774394/93‑1 (MQ=255)
                                                         tGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGa                                >  1:2039488/1‑93 (MQ=255)
                                                          gCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGAt                               <  2:883921/93‑1 (MQ=255)
                                                           cTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              >  1:178059/1‑93 (MQ=255)
                                                            tCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATcc                             >  1:1734167/1‑93 (MQ=255)
                                                            tCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATcc                             <  1:1859119/93‑1 (MQ=255)
                                                            tCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATcc                             >  1:1580383/1‑93 (MQ=255)
                                                            tCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATcc                             <  2:721190/93‑1 (MQ=255)
                                                               aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          >  1:6287885/1‑93 (MQ=255)
                                                               aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          >  2:3980434/1‑93 (MQ=255)
                                                                tAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCtt                         <  2:5406758/93‑1 (MQ=255)
                                                                 aaCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCttt                        <  2:6039789/93‑1 (MQ=255)
                                                                   ccGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgc                      <  1:5237215/93‑1 (MQ=255)
                                                                   ccGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgc                      <  1:3261504/93‑1 (MQ=255)
                                                                    cgcgCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcg                     <  1:2717152/93‑1 (MQ=255)
                                                                    cgcgCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcg                     <  1:2353667/93‑1 (MQ=255)
                                                                     gcgcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcgc                    <  2:1511470/93‑1 (MQ=255)
                                                                     gcgcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcgc                    >  2:1583310/1‑93 (MQ=255)
                                                                     gcgcTAACTGGAAGATTGATCGGGTGCGTGAGGAGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcgc                    >  2:3783961/1‑93 (MQ=255)
                                                                      cgcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGcgcg                   <  2:316678/93‑1 (MQ=255)
                                                                       gcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGt                  >  1:5758033/1‑93 (MQ=255)
                                                                       gcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGt                  >  1:4601707/1‑93 (MQ=255)
                                                                        cTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTc                 >  1:4203138/1‑93 (MQ=255)
                                                                        cTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTc                 <  1:395183/93‑1 (MQ=255)
                                                                        cTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTc                 <  2:4234265/93‑1 (MQ=255)
                                                                            cTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCaa             <  2:5017268/93‑1 (MQ=255)
                                                                              ggAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATa           >  1:3413614/1‑93 (MQ=255)
                                                                               gAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATaa          <  2:2326725/93‑1 (MQ=255)
                                                                               gAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGCAATCTGATCCGCTTTGCGCGTCGCAATaa          >  2:4802944/1‑93 (MQ=255)
                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        <  2:2176534/93‑1 (MQ=255)
                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        <  1:4955426/93‑1 (MQ=255)
                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        >  1:1023349/1‑93 (MQ=255)
                                                                                   aTTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCtt      <  1:4303832/93‑1 (MQ=255)
                                                                                   aTTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCtt      <  1:4746691/93‑1 (MQ=255)
                                                                                    ttGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTAGCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTa     <  2:2890333/93‑1 (MQ=255)
                                                                                    ttGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTa     >  2:2933436/1‑93 (MQ=255)
                                                                                     tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    <  1:821421/93‑1 (MQ=255)
                                                                                      gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTATCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   <  1:4510204/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987

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