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,8980AG99.6% 847.8 / ‑4.6 260Q378R (CAG→CGG) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (123/136);  minor base T (1/0);  total (124/136)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

Only 100 of 277 total aligned reads displayed.

GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987
                                                                                          |                                                                                         
gcgcAGTTTGGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCACGTCTGTTAATGCTCGATATCCGCGCGTACTGGAAGATTGATCggg                                                                                         >  2:6726796/1‑93 (MQ=255)
gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         <  1:1089599/93‑1 (MQ=255)
gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         >  2:1193878/1‑93 (MQ=255)
 cgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGt                                                                                        <  1:3147182/93‑1 (MQ=255)
  gcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTg                                                                                       <  1:2025938/93‑1 (MQ=255)
   cAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGc                                                                                      >  1:1551541/1‑93 (MQ=255)
      tgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTg                                                                                   >  2:4627647/1‑93 (MQ=255)
         tGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAgg                                                                                <  1:729789/93‑1 (MQ=255)
            gCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGgcgc                                                                             >  1:1676796/1‑93 (MQ=255)
                    gTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAg                                                                     >  2:743969/1‑93 (MQ=255)
                     tGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGt                                                                    <  2:6149523/93‑1 (MQ=255)
                      gAGCGAGTGGGGTTGGGACGAAGCTCGTCTTGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGtt                                                                   <  1:6371509/93‑1 (MQ=255)
                         cGAGTGGGGTTTGGACGAGGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCt                                                                <  2:1726665/93‑1 (MQ=255)
                         cGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCt                                                                <  2:6330631/93‑1 (MQ=255)
                         cGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCt                                                                >  1:5039590/1‑93 (MQ=255)
                          gAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTc                                                               <  1:5122866/93‑1 (MQ=255)
                          gAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTc                                                               <  2:4455825/93‑1 (MQ=255)
                            gTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             <  2:5140300/93‑1 (MQ=255)
                                  ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       <  2:6438473/93‑1 (MQ=255)
                                  ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       <  1:690439/93‑1 (MQ=255)
                                  ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       <  1:3474014/93‑1 (MQ=255)
                                   tGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatg                                                      >  2:2488920/1‑93 (MQ=255)
                                    gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     <  1:5885962/93‑1 (MQ=255)
                                    gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     >  1:1041089/1‑93 (MQ=255)
                                    gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     <  1:3618224/93‑1 (MQ=255)
                                      gACGAAGCGCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGatgatg                                                   >  1:6535744/1‑93 (MQ=255)
                                       aCGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGTCGCACAACGAGTTGCTCGACGCGATGATGc                                                  >  2:1286862/1‑93 (MQ=255)
                                       aCGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGc                                                  <  1:6001218/93‑1 (MQ=255)
                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  1:3425535/93‑1 (MQ=255)
                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 >  1:2076282/1‑93 (MQ=255)
                                         gAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCag                                                >  2:4820675/1‑93 (MQ=255)
                                          aaGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCaga                                               >  2:565475/1‑93 (MQ=255)
                                          aaGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCaga                                               <  2:6052301/93‑1 (MQ=255)
                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              >  2:3651711/1‑93 (MQ=255)
                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              <  2:3274360/93‑1 (MQ=255)
                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              >  1:3858496/1‑93 (MQ=255)
                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              <  1:1364752/93‑1 (MQ=255)
                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              <  1:4369546/93‑1 (MQ=255)
                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              >  1:569240/1‑93 (MQ=255)
                                            gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             <  2:819986/93‑1 (MQ=255)
                                            gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             >  2:6227870/1‑93 (MQ=255)
                                            gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             >  2:5817042/1‑93 (MQ=255)
                                              tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTa                                           >  1:5712015/1‑93 (MQ=255)
                                              tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTa                                           >  2:1007563/1‑93 (MQ=255)
                                                   tGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTa                                      <  2:5485896/93‑1 (MQ=255)
                                                    ggCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTaa                                     <  1:3942790/93‑1 (MQ=255)
                                                     gCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAAt                                    >  1:540578/1‑93 (MQ=255)
                                                     gCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAAt                                    <  2:1458140/93‑1 (MQ=255)
                                                      cAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATc                                   <  1:6054631/93‑1 (MQ=255)
                                                      cAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATc                                   <  1:6217122/93‑1 (MQ=255)
                                                       aaTGCTCGATAACCGCGCTAACTGGAATATTGATCGGGTGCGTGAGGCGCACAACGAGTTCCTCGACGCGATGATGCAGAGCTACCGTAATCt                                  <  1:5586025/93‑1 (MQ=255)
                                                       aaTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCt                                  <  2:3588519/93‑1 (MQ=255)
                                                       aaTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCt                                  <  2:3354368/93‑1 (MQ=255)
                                                        aTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTg                                 >  2:3152661/1‑93 (MQ=255)
                                                         tGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGa                                >  1:2947899/1‑93 (MQ=255)
                                                          gCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGAt                               <  1:4476083/93‑1 (MQ=255)
                                                          gCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGAt                               <  1:4933935/93‑1 (MQ=255)
                                                           cTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              <  1:6689114/93‑1 (MQ=255)
                                                           cTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              >  1:2516866/1‑93 (MQ=255)
                                                              gATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGc                           >  1:1860842/1‑93 (MQ=255)
                                                              gATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGc                           <  2:1785845/93‑1 (MQ=255)
                                                               aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          >  1:5796019/1‑93 (MQ=255)
                                                               aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          <  2:3896825/93‑1 (MQ=255)
                                                               aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          <  2:4472821/93‑1 (MQ=255)
                                                                tAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCtt                         <  1:2996869/93‑1 (MQ=255)
                                                                tAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCtt                         <  1:5161261/93‑1 (MQ=255)
                                                                tAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCtt                         <  1:4174000/93‑1 (MQ=255)
                                                                tAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCtt                         >  2:3010683/1‑93 (MQ=255)
                                                                 aaCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCttt                        >  1:3730107/1‑93 (MQ=255)
                                                                        cTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTc                 <  2:6502161/93‑1 (MQ=255)
                                                                         tAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCg                >  2:4547806/1‑93 (MQ=255)
                                                                         tAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCg                >  1:805640/1‑93 (MQ=255)
                                                                          aaCTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGc               >  1:2438708/1‑93 (MQ=255)
                                                                          aaCTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGc               >  1:1234031/1‑93 (MQ=255)
                                                                          aaCTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGc               <  2:2622690/93‑1 (MQ=255)
                                                                            cTGGAAGATTGAGCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCaa             <  2:4831905/93‑1 (MQ=255)
                                                                             tGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAAt            >  2:1918437/1‑93 (MQ=255)
                                                                             tGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAAt            >  1:1490460/1‑93 (MQ=255)
                                                                              ggAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGTTGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATa           >  2:1880444/1‑93 (MQ=255)
                                                                              ggAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATa           >  2:1953665/1‑93 (MQ=255)
                                                                              ggAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATa           >  2:2845501/1‑93 (MQ=255)
                                                                                aaGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAc         <  1:2362574/93‑1 (MQ=255)
                                                                                aaGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAc         >  2:2292434/1‑93 (MQ=255)
                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        >  2:2918653/1‑93 (MQ=255)
                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        <  1:4571038/93‑1 (MQ=255)
                                                                                 aGATTGATCGGGTGCGTGAGGCGCAAAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        <  2:5669617/93‑1 (MQ=255)
                                                                                  gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       >  2:1572724/1‑93 (MQ=255)
                                                                                  gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       <  2:6690630/93‑1 (MQ=255)
                                                                                  gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       <  2:4309520/93‑1 (MQ=255)
                                                                                   aTTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCtt      >  2:4874400/1‑93 (MQ=255)
                                                                                   aTTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCtt      <  1:4651980/93‑1 (MQ=255)
                                                                                   aTTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCtt      <  2:2016100/93‑1 (MQ=255)
                                                                                    ttGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCCTa     >  1:4729839/1‑93 (MQ=255)
                                                                                     tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    <  1:5006818/93‑1 (MQ=255)
                                                                                     tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    <  1:2774055/93‑1 (MQ=255)
                                                                                     tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    <  1:2903919/93‑1 (MQ=255)
                                                                                      gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   >  1:2616672/1‑93 (MQ=255)
                                                                                      gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   >  2:6600202/1‑93 (MQ=255)
                                                                                      gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   >  2:3037903/1‑93 (MQ=255)
                                                                                       aTCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCg  >  2:2228914/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987

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