Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,085,695 T→C K581R (AAG→AGG)  frvR ← putative frv operon regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,085,6950TC97.1% 524.7 / 3.3 171K581R (AAG→AGG) frvRputative frv operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (91/75);  total (95/76)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 184 total aligned reads displayed.

CGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTA  >  NC_000913/4085604‑4085784
                                                                                           |                                                                                         
cGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtg                                                                                         <  1:2116469/93‑1 (MQ=255)
 gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                                        >  1:3868847/1‑93 (MQ=255)
  aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgtg                                                                                       >  1:3712563/1‑93 (MQ=255)
    aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGaa                                                                                     <  1:669731/93‑1 (MQ=255)
     tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAg                                                                                    >  2:1475152/1‑93 (MQ=255)
      ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGt                                                                                   >  1:2771917/1‑93 (MQ=255)
       gTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTa                                                                                  >  2:1194333/1‑93 (MQ=255)
         tgatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAAc                                                                                <  1:3591515/93‑1 (MQ=255)
          gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACt                                                                               >  1:3893832/1‑93 (MQ=255)
          gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACt                                                                               >  1:904822/1‑93 (MQ=255)
           atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                                              >  1:2793366/1‑93 (MQ=255)
            tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTcc                                                                             >  1:2198266/1‑93 (MQ=255)
            tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTcc                                                                             >  2:3997439/1‑93 (MQ=255)
             gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCa                                                                            >  1:234726/1‑93 (MQ=255)
             gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCa                                                                            <  1:3039732/93‑1 (MQ=255)
              aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCAt                                                                           >  1:4504326/1‑93 (MQ=255)
              aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                                           >  2:356396/1‑93 (MQ=255)
              aCTATAATAGCCTGCGTTAAACGTTTTTGACGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                                           <  2:4636712/93‑1 (MQ=255)
               cTATAATAGCCTGCGTTACACGTTTTTGATGATAAATGCAGTGTCGAATGCGGCGCATCCGACACAGTTAACCCCCTGTGAAGTAACTCCATa                                                                          <  2:2916547/93‑1 (MQ=255)
               cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                                          <  2:2825244/93‑1 (MQ=255)
                  taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                                       >  1:3475180/1‑93 (MQ=255)
                  taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                                       >  1:4417822/1‑93 (MQ=255)
                     taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGCAGTAACTCCATAAATGcc                                                                    >  1:1663606/1‑93 (MQ=255)
                     taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGcc                                                                    <  2:3512575/93‑1 (MQ=255)
                     taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGGTCACCCCCTGTGAAGTAACTCCATAAATGcc                                                                    >  1:3858286/1‑93 (MQ=255)
                      aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCt                                                                   <  1:319084/93‑1 (MQ=255)
                       gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                                                  >  2:3110518/1‑93 (MQ=255)
                       gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                                                  >  1:1082908/1‑93 (MQ=255)
                         cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCat                                                                >  1:389532/1‑93 (MQ=255)
                          tGCGTTAAACGTTTTTGATGATAAATGAAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                                               >  2:2741211/1‑93 (MQ=255)
                           gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                                              >  2:2093630/1‑93 (MQ=255)
                           gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                                              <  2:3191941/93‑1 (MQ=255)
                           gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                                              >  2:4554499/1‑93 (MQ=255)
                            cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                                             >  1:4657868/1‑93 (MQ=255)
                             gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                                            >  1:477744/1‑93 (MQ=255)
                             gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATAGcc                                                            >  2:732929/1‑93 (MQ=255)
                              ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                           <  2:717588/93‑1 (MQ=255)
                              ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                           >  2:1692742/1‑93 (MQ=255)
                               tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                                          <  2:4439393/93‑1 (MQ=255)
                               tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                                          >  2:819229/1‑93 (MQ=255)
                                 aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTgct                                                        >  2:2890834/1‑91 (MQ=255)
                                   cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                                      <  1:3680915/93‑1 (MQ=255)
                                    gTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                                                     >  1:3630531/1‑93 (MQ=255)
                                     tttttGATGATAAATGCAGGGTCGGATGCGGCGCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                                    <  1:4158185/93‑1 (MQ=255)
                                      ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                                                   >  2:1343460/1‑93 (MQ=255)
                                         tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                                <  1:463390/93‑1 (MQ=255)
                                          gatgatAAATGCAGTGTCGGATGCGGTGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCNTCATATCCTGTTAACCCGGCa                                               >  1:635305/1‑93 (MQ=255)
                                          gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATCAATGCCTCATATCCTGTTAACCCGGCa                                               >  1:1327472/1‑93 (MQ=255)
                                           atgatAACTGCAGTGTCGGATGCTGCGCATCCGACACAGTGCACCCCCTGTGAAGTAACTCCATAAATGCATCATATCCTGTTAACCCGGCaa                                              <  2:235329/93‑1 (MQ=255)
                                           atgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa                                              >  1:2443884/1‑93 (MQ=255)
                                            tgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                                             >  1:3760758/1‑93 (MQ=255)
                                              atAAATGCAGTGTGGGATGCGGCGCATCCAACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                                           <  1:1544196/93‑1 (MQ=255)
                                               tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                                          >  2:4715245/1‑93 (MQ=255)
                                                 aaTGCAGCGTCGGATGCGTCGGATCCGACACAGTTCACGCCCTGTGAAGTGACTCCATAAATGACTCGTATCCGGTTAACCCGGCAATAAtct                                        >  2:3186238/1‑93 (MQ=255)
                                                  aTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                                       >  1:4373590/1‑93 (MQ=255)
                                                   tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                                      >  1:3332017/1‑93 (MQ=255)
                                                    gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                                     >  1:2072757/1‑93 (MQ=255)
                                                     cAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgc                                    <  1:2814919/93‑1 (MQ=255)
                                                      aGTGTCGGGTGCGGCGCATCCGACACAGGTCACCCCCTGTGAAGGAACTCCATAAATGCCTCATATCCTGTTTAACCCGCAATAATCTccgcc                                   >  2:1876007/1‑93 (MQ=255)
                                                        tgtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                                 >  1:2916392/1‑93 (MQ=255)
                                                         gtgggATGCGGCGCATCCGACGCAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                <  1:2263315/90‑1 (MQ=255)
                                                         gtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                >  2:1081663/1‑93 (MQ=255)
                                                          tCGGATGGTGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                               <  1:898794/93‑1 (MQ=255)
                                                          tCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCggct                               >  2:1823263/1‑91 (MQ=255)
                                                           tggctGCGGCGCATACCACACAGTTCACCCGCTGTGCAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                              <  1:2216499/89‑1 (MQ=255)
                                                           cGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                              >  1:3188780/1‑93 (MQ=255)
                                                           cGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                              >  2:3316711/1‑93 (MQ=255)
                                                            ggctGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc                             <  1:2477385/90‑1 (MQ=255)
                                                             gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGcc                            <  1:202186/93‑2 (MQ=255)
                                                             gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt                            >  2:1398651/1‑93 (MQ=255)
                                                             gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt                            >  1:4018969/1‑93 (MQ=255)
                                                             gATGCGGCGCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt                            <  1:4196061/93‑1 (MQ=255)
                                                               tGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa                          >  2:411085/1‑93 (MQ=255)
                                                               tGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGACGGATGCTGa                          <  2:3057309/93‑1 (MQ=255)
                                                                gcggcgCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac                         <  1:3629859/93‑1 (MQ=255)
                                                                    gcgcaTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGAcaca                      >  2:4182818/5‑93 (MQ=255)
                                                                    gcgcaTCCGACACAGTTCACCCCCTGCGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGAcaca                      >  1:1609027/5‑93 (MQ=255)
                                                                     cgcaTCAGACACAGTTCCCCCCCTGTGAAGTAACGCCATAACTGCCTCATATCCTGTTAGCCCGGCAATAATCTCCGCTGGATGCTGACACaa                     <  2:2654533/90‑1 (MQ=255)
                                                                      gcaTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                    <  1:583657/91‑1 (MQ=255)
                                                                      gcaTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                    <  1:4508780/91‑1 (MQ=255)
                                                                       caTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   <  1:1320149/92‑1 (MQ=255)
                                                                       caTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   <  1:750317/92‑1 (MQ=255)
                                                                        aTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat                  <  2:595816/93‑1 (MQ=255)
                                                                         tCCGCCACATTTCACCACCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTATCCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                 <  1:4320234/93‑1 (MQ=255)
                                                                           cGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt               <  2:2478327/93‑1 (MQ=255)
                                                                            gACACAGTGCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt              <  1:2640176/93‑1 (MQ=255)
                                                                             acacAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTg             >  2:3744060/1‑93 (MQ=255)
                                                                               acaGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGcc           <  1:927984/93‑1 (MQ=255)
                                                                                caGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCa          <  1:791467/93‑1 (MQ=255)
                                                                                caGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCa          <  2:4620269/93‑1 (MQ=255)
                                                                                 aGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAg         <  1:1450000/93‑1 (MQ=255)
                                                                                   ttCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt       <  2:2204136/93‑1 (MQ=255)
                                                                                    tCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATa      <  1:1035421/93‑1 (MQ=255)
                                                                                    tCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATa      <  2:2927581/93‑1 (MQ=255)
                                                                                      aCCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGc    >  2:425178/1‑93 (MQ=255)
                                                                                      aCCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGc    <  2:1052528/93‑1 (MQ=255)
                                                                                       ccccTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATCGCt   >  2:1580269/1‑93 (MQ=255)
                                                                                       ccccCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCt   <  2:3048362/93‑1 (MQ=255)
                                                                                        ccccTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTa  >  1:1278641/1‑93 (MQ=255)
                                                                                        ccccTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTa  >  1:2559959/1‑93 (MQ=255)
                                                                                           |                                                                                         
CGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTA  >  NC_000913/4085604‑4085784

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