Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,758 G→A 100% intergenic (+79/‑12) ygaH → / → mprA putative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,7580GA98.8% 530.1 / ‑3.2 161intergenic (+79/‑12)ygaH/mprAputative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (68/89);  total (71/90)
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.25e-01

Only 100 of 161 total aligned reads displayed.

TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTC  >  NC_000913/2810663‑2810854
                                                                                               |                                                                                                
tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCa                                                                                              <  2:1264153/100‑1 (MQ=255)
tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCa                                                                                              <  1:3189244/100‑1 (MQ=255)
 gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCAt                                                                                             >  2:1610739/1‑100 (MQ=255)
  aTGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATa                                                                                            <  1:3883979/100‑1 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTTTTTGTCACTGTCGTTACTATATCGTCTGAAATTAATAAAGTCATAcc                                                                                          >  2:734163/1‑100 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAcc                                                                                          <  1:3033294/100‑1 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAcc                                                                                          >  1:245033/1‑100 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAcc                                                                                          <  2:3973147/100‑1 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAcc                                                                                          <  1:3739311/100‑1 (MQ=255)
     gCGATTATATAACCCACAAGAATCATTTTTTTTTTTTTTTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAccc                                                                                         >  2:1831531/1‑100 (MQ=255)
     gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAccc                                                                                         <  1:1832472/100‑1 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCa                                                                                        <  2:4039868/100‑1 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCa                                                                                        >  2:87575/1‑100 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCa                                                                                        >  2:3702715/1‑100 (MQ=255)
         ttATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAAt                                                                                     <  1:3638047/100‑1 (MQ=255)
         ttATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAAt                                                                                     <  1:1961394/100‑1 (MQ=255)
         ttATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAAt                                                                                     <  1:1143186/100‑1 (MQ=255)
           atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATgg                                                                                   >  2:1204949/1‑100 (MQ=255)
            tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGa                                                                                  >  1:74128/1‑100 (MQ=255)
            tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGa                                                                                  >  2:4174828/1‑100 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATa                                                                                >  2:2703497/1‑100 (MQ=255)
               aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAg                                                                               <  2:2712558/100‑1 (MQ=255)
                 cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGtt                                                                             <  2:510969/100‑1 (MQ=255)
                 cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGtt                                                                             <  1:1765614/100‑1 (MQ=255)
                  ccACAAGAATCATTTTTCTAAAACGATACCTTTACTGTATTTGTAACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGt                                                                              <  2:3562115/98‑1 (MQ=255)
                  ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCTGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTc                                                                            >  2:3931845/1‑100 (MQ=255)
                  ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTc                                                                            <  1:2042046/100‑1 (MQ=255)
                  ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTc                                                                            <  1:839005/100‑1 (MQ=255)
                  ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTc                                                                            <  2:2294021/100‑1 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCg                                                                           >  2:1532066/1‑100 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTa                                                                       <  2:196002/100‑1 (MQ=255)
                        gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGtt                                                                         <  2:3473214/97‑1 (MQ=255)
                        gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTAc                                                                      >  2:1262764/1‑100 (MQ=255)
                        gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTAc                                                                      <  1:266722/100‑1 (MQ=255)
                        gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTAc                                                                      >  2:3256393/1‑100 (MQ=255)
                        gAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTAc                                                                      >  2:3361450/1‑100 (MQ=255)
                         aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACg                                                                     <  2:3310611/100‑1 (MQ=255)
                          aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGc                                                                    >  1:1368283/1‑100 (MQ=255)
                          aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGc                                                                    <  1:3450762/100‑1 (MQ=255)
                          aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGc                                                                    <  2:1085085/100‑1 (MQ=255)
                           tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGcc                                                                   <  1:3045700/100‑1 (MQ=255)
                           tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGcc                                                                   <  2:3014380/100‑1 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                              <  1:1000015/100‑1 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTATCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                              <  2:2765046/100‑1 (MQ=255)
                                 ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                             <  1:1173684/100‑1 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                           >  2:2702066/1‑100 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                           >  2:1812120/1‑100 (MQ=255)
                                     aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAACTGGATAGTTCGTTTACACCTATTGAACa                                                          >  2:2624693/1‑99 (MQ=255)
                                       aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                       >  1:1269709/1‑100 (MQ=255)
                                       aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                       <  2:2420332/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      <  2:817048/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      <  2:1989065/100‑1 (MQ=255)
                                         cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                                     <  1:2959780/100‑1 (MQ=255)
                                           aTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTa                                                   >  1:1184071/1‑100 (MQ=255)
                                           aTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTa                                                   >  2:375895/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                >  1:2226878/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                >  1:404056/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                >  1:1287863/1‑100 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                            >  1:2495898/1‑100 (MQ=255)
                                                   aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                           <  2:1881006/100‑1 (MQ=255)
                                                    cTTTATTGGTCCCTGTCGTTACTATATCGGCTGAAATTAATAACGACATACACAAAACAATATTTATTTTATGCCACTTGAACAAATGATAAGAttt                                             >  2:2179389/1‑97 (MQ=255)
                                                         tgtGTCACTGTCGTTGCTATATCGGCTGCTATTAATAAAGTAATACCCAAATTGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                                     <  2:2773125/98‑1 (MQ=255)
                                                           tGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                   >  1:1996385/1‑100 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                                  >  1:943742/1‑100 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGACTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                                  >  2:1135587/1‑100 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTATGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                                  >  2:3033783/1‑98 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                <  1:2702082/100‑1 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                <  1:1641358/100‑1 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                <  1:911029/100‑1 (MQ=255)
                                                               aCTGTCATTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcta                               <  2:1686779/100‑3 (MQ=255)
                                                                cTGTCGTACATATCTCGGCTGCAATTAATAAAGTCATACCCAAATGGATAGTTCGTTCACGCCCATTGAACAAATGCTAAAAATTCGGGCCAGCCGCCa                               <  2:352713/99‑1 (MQ=255)
                                                                 tGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAcn                             <  2:3014609/100‑2 (MQ=255)
                                                                    cGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                          <  2:832856/100‑1 (MQ=255)
                                                                     gTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                         >  1:1547601/1‑100 (MQ=255)
                                                                        aCTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAttt                      <  1:843348/100‑1 (MQ=255)
                                                                        aCTATATCGGCTGAAATTAATAAAGACATACCAAAATGGATAGTTCGTTTACGCCCATTGAACAGAAGAAAAAATTTCGCGCCAACCACCACGAAGAtt                       >  2:1173108/1‑99 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                   <  2:645859/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                   <  1:103723/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                   <  1:2598735/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                   <  1:3474099/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                  <  1:3652267/100‑1 (MQ=255)
                                                                             atCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt                 >  2:3341942/1‑100 (MQ=255)
                                                                             atCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt                 >  2:3640660/1‑100 (MQ=255)
                                                                               cGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt               >  2:2126890/1‑100 (MQ=255)
                                                                               cGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt               <  1:448878/100‑1 (MQ=255)
                                                                                ggCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc              >  1:2043600/1‑100 (MQ=255)
                                                                                 gCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa             <  1:3256416/100‑1 (MQ=255)
                                                                                 gCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa             >  1:2428252/1‑100 (MQ=255)
                                                                                 cttGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa             <  2:728669/98‑1 (MQ=255)
                                                                                  cTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg            <  2:3830439/100‑1 (MQ=255)
                                                                                   tGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCaag           >  1:728673/1‑98 (MQ=255)
                                                                                   tGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAgg           <  2:3650976/100‑1 (MQ=255)
                                                                                    gAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga          <  1:1166059/100‑1 (MQ=255)
                                                                                     aaaTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgag         <  2:598415/100‑1 (MQ=255)
                                                                                       aTTAATAAAGTCATACCCAAATGGATAGTTCGTTTATGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATGTTCCTTATCAGGAGAt       >  2:952702/1‑100 (MQ=255)
                                                                                       aTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaggt       >  1:3172048/1‑98 (MQ=255)
                                                                                       aTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt       >  1:1788526/1‑100 (MQ=255)
                                                                                       aTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt       <  1:2855844/100‑1 (MQ=255)
                                                                                       aTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt       >  2:3181332/1‑100 (MQ=255)
                                                                                            taaGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc  >  1:1149071/3‑100 (MQ=255)
                                                                                               |                                                                                                
TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTC  >  NC_000913/2810663‑2810854

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