Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 652,748 A→G T431A (ACC→GCC)  citA → putative sensor‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913652,7480AG99.6% 770.8 / ‑5.1 242T431A (ACC→GCC) citAputative sensor‑type protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (102/138);  total (103/139)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

Only 100 of 248 total aligned reads displayed.

GGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAATTA  >  NC_000913/652658‑652837
                                                                                          |                                                                                         
ggTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGcc                                                                                         >  2:6103488/1‑93 (MQ=255)
     attGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGTATAGCGCCGAGtt                                                                                    <  1:1823685/91‑1 (MQ=255)
     aGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGtt                                                                                    <  2:5584416/93‑1 (MQ=255)
     aGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGtt                                                                                    >  2:29191/1‑93 (MQ=255)
     aGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGtt                                                                                    <  1:3187808/93‑1 (MQ=255)
     aGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGtt                                                                                    <  1:5647563/93‑1 (MQ=255)
      gTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGttt                                                                                   <  1:3975881/93‑1 (MQ=255)
      gTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGttt                                                                                   >  1:5473804/1‑93 (MQ=255)
        gCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGc                                                                                 <  2:5245037/93‑1 (MQ=255)
              cgcCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCAt                                                                           <  1:2674964/93‑1 (MQ=255)
               gcCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCAtt                                                                          <  1:2200094/93‑1 (MQ=255)
                cccGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATtg                                                                         <  1:4225238/93‑1 (MQ=255)
                cccGGGAACTGGAGTTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATtg                                                                         <  2:859391/93‑1 (MQ=255)
                 ccGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATtgt                                                                        <  1:1493291/93‑1 (MQ=255)
                 ccGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATtgt                                                                        <  1:492225/93‑1 (MQ=255)
                  cGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGGCGAGTTTGCAGCCATtgtg                                                                       >  1:2312772/1‑93 (MQ=255)
                  cGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATtgtg                                                                       <  1:872415/93‑1 (MQ=255)
                   gggAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTgg                                                                      <  2:2043256/93‑1 (MQ=255)
                   gggAACTGGGGCTAAAAATGATCATTGTCCACGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTgg                                                                      <  2:961719/93‑1 (MQ=255)
                    ggAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTggg                                                                     >  1:695498/1‑93 (MQ=255)
                     gAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGc                                                                    >  2:6492363/1‑93 (MQ=255)
                       aCTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCaa                                                                  >  1:5065434/1‑93 (MQ=255)
                            ggCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTAc                                                             <  2:4261450/93‑1 (MQ=255)
                            ggCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTAc                                                             >  1:3922209/1‑93 (MQ=255)
                             gCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACt                                                            <  1:1290768/93‑1 (MQ=255)
                               tAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTg                                                          >  1:4103135/1‑93 (MQ=255)
                               tAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTg                                                          >  1:4948483/1‑93 (MQ=255)
                                 aaaaTGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGCCTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGAt                                                        <  2:4322326/93‑1 (MQ=255)
                                 aaaaTGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGAt                                                        >  2:5755640/1‑93 (MQ=255)
                                  aaaTGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATa                                                       >  2:2973936/1‑93 (MQ=255)
                                  aaaTGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATa                                                       <  1:576625/93‑1 (MQ=255)
                                   aaTGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATaa                                                      >  2:5459108/1‑93 (MQ=255)
                                   aaTGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATaa                                                      <  2:1623169/93‑1 (MQ=255)
                                    aTGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAAc                                                     <  1:5948460/93‑1 (MQ=255)
                                    aTGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAAc                                                     >  1:3414984/1‑93 (MQ=255)
                                     tGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACg                                                    <  1:3524055/93‑1 (MQ=255)
                                     tGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACg                                                    >  1:4020313/1‑93 (MQ=255)
                                     tGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACg                                                    >  1:4781594/1‑93 (MQ=255)
                                       gtCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGcc                                                  <  2:4899215/92‑1 (MQ=255)
                                       aTCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGcc                                                  <  1:613959/93‑1 (MQ=255)
                                       aTCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGcc                                                  >  2:4203871/1‑93 (MQ=255)
                                        tCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCt                                                 <  1:4653780/93‑1 (MQ=255)
                                        tCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCt                                                 >  1:375664/1‑93 (MQ=255)
                                        tCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCt                                                 >  1:727774/1‑93 (MQ=255)
                                        tCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCt                                                 <  1:6389403/93‑1 (MQ=255)
                                        tCATTGTCCCCGGTACCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCt                                                 >  1:1285334/1‑93 (MQ=255)
                                          aTTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTc                                               <  1:2860586/93‑1 (MQ=255)
                                             gcccccGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTGCTTGATAACGCCTTCGaa                                            <  2:468504/91‑1 (MQ=255)
                                             gTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGaa                                            <  2:6217932/93‑1 (MQ=255)
                                                cccGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAagcc                                         >  1:2054688/1‑93 (MQ=255)
                                                cacGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAagcc                                         <  2:344671/91‑1 (MQ=255)
                                                 ccGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAagcca                                        <  1:6142775/93‑1 (MQ=255)
                                                 ccGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAagcca                                        <  2:2187670/93‑1 (MQ=255)
                                                  cGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAagccag                                       >  2:1958583/1‑93 (MQ=255)
                                                  cGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAagccag                                       >  2:751513/1‑93 (MQ=255)
                                                    gTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAagccagcc                                     >  2:4653827/1‑93 (MQ=255)
                                                      aGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTg                                   >  2:72458/1‑93 (MQ=255)
                                                         cAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGt                                <  1:5237097/93‑1 (MQ=255)
                                                         cAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGt                                >  1:4543026/1‑93 (MQ=255)
                                                         cAGCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGt                                >  1:4758239/1‑93 (MQ=255)
                                                           gCTTTCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAg                              >  2:3700353/1‑93 (MQ=255)
                                                             tttCGCAACTGCCGCGAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCg                            <  2:5314332/93‑1 (MQ=255)
                                                             tttCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCg                            <  1:3023055/93‑1 (MQ=255)
                                                             tttCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCg                            >  2:2716546/1‑93 (MQ=255)
                                                             tttCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCg                            <  2:329533/93‑1 (MQ=255)
                                                              ttCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGa                           <  1:5278962/93‑1 (MQ=255)
                                                              ttCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGa                           <  1:5976639/93‑1 (MQ=255)
                                                              ttCGCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGa                           >  1:3767328/1‑93 (MQ=255)
                                                                 gCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGa                        <  2:1198911/93‑1 (MQ=255)
                                                                 gCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGa                        <  2:4823736/93‑1 (MQ=255)
                                                                 gCAACTGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGa                        >  2:3517509/1‑93 (MQ=255)
                                                                   aaCTGCCGCCAGGACTGGATAGCGCCGAGTTGGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAg                      <  1:1713296/93‑1 (MQ=255)
                                                                      tGCCGCCAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGaa                   >  2:1459469/1‑93 (MQ=255)
                                                                       gccgccAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGaaa                  <  2:4067740/93‑1 (MQ=255)
                                                                       gccgccAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGaaa                  <  2:1523403/93‑1 (MQ=255)
                                                                        ccgccAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAAc                 <  2:6433421/93‑1 (MQ=255)
                                                                          gccAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACaa               <  2:5913527/93‑1 (MQ=255)
                                                                           ccAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAg              <  1:3715482/93‑1 (MQ=255)
                                                                           ccAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAg              >  2:1851146/1‑93 (MQ=255)
                                                                            cAGGACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGa             >  2:4074770/1‑93 (MQ=255)
                                                                              ggACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATc           >  1:3821297/1‑93 (MQ=255)
                                                                              ggACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATc           >  1:1196277/1‑93 (MQ=255)
                                                                              ggACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATc           <  2:231291/93‑1 (MQ=255)
                                                                              ggACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATc           >  2:211282/1‑93 (MQ=255)
                                                                               gACTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCg          <  2:6159673/93‑1 (MQ=255)
                                                                                aCTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGt         <  1:2793834/93‑1 (MQ=255)
                                                                                aCTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGt         <  1:1696750/93‑1 (MQ=255)
                                                                                 cTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGtt        >  2:4938437/1‑93 (MQ=255)
                                                                                 cTGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGtt        <  1:3061504/93‑1 (MQ=255)
                                                                                  tGGATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTg       >  2:4737639/1‑93 (MQ=255)
                                                                                   ggATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGa      >  1:6249599/1‑93 (MQ=255)
                                                                                   ggATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGa      >  1:5484380/1‑93 (MQ=255)
                                                                                    gATAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGaa     <  2:1233215/93‑1 (MQ=255)
                                                                                     aTAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAat    <  1:2650654/93‑1 (MQ=255)
                                                                                     aTAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAat    <  2:2463375/93‑1 (MQ=255)
                                                                                     aTAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAat    <  1:5207589/93‑1 (MQ=255)
                                                                                     aTAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAat    <  2:257863/93‑1 (MQ=255)
                                                                                      tAGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAatt   <  2:6535741/93‑1 (MQ=255)
                                                                                       aGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAatta  <  1:1534510/93‑1 (MQ=255)
                                                                                       aGCGCCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAatta  >  1:3297199/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGCGTAGCGATGAAGGAAACAAGATCGTTGAATTA  >  NC_000913/652658‑652837

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: