Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,810,758 G→A 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,7580GA93.3% 279.6 / 6.7 90intergenic (+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/5);  new base A (44/40);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGA  >  NC_000913/2810663‑2810848
                                                                                               |                                                                                          
tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCa                                                                                        <  2:1830531/100‑1 (MQ=255)
  aTGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATa                                                                                      >  1:275003/1‑100 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAc                                                                                     >  2:91779/1‑100 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAc                                                                                     >  1:1759641/1‑100 (MQ=255)
     gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAc                                                                                     >  2:1645168/1‑98 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATg                                                                              <  2:1467723/100‑1 (MQ=255)
           atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATgg                                                                             >  2:1910775/1‑100 (MQ=255)
             ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGAt                                                                           <  1:1520557/100‑1 (MQ=255)
             ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGAt                                                                           <  2:513633/100‑1 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATa                                                                          >  2:1797199/1‑100 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATa                                                                          <  1:1829482/100‑1 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATa                                                                          <  2:1015984/100‑1 (MQ=255)
                aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGt                                                                        >  1:826611/1‑100 (MQ=255)
                 cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGtt                                                                       >  1:1301465/1‑100 (MQ=255)
                 cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGtt                                                                       <  2:68763/100‑1 (MQ=255)
                  ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTc                                                                      >  1:461266/1‑100 (MQ=255)
                  ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTc                                                                      <  2:100848/100‑1 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCg                                                                     >  2:1391185/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCg                                                                     >  2:608244/1‑100 (MQ=255)
                      aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGtt                                                                   <  2:1289345/99‑1 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTa                                                                 >  2:495930/1‑100 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTa                                                                 <  2:1519528/100‑1 (MQ=255)
                          aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGc                                                              <  1:571056/100‑1 (MQ=255)
                          aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGc                                                              <  2:704055/100‑1 (MQ=255)
                            tccgatcttAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGccc                                                            <  2:540215/92‑1 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGccc                                                            >  2:1744787/1‑100 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGccc                                                            <  1:1083539/100‑1 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGccc                                                            >  2:801202/1‑100 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGccc                                                            <  1:1083525/100‑1 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGccc                                                            >  2:1746065/1‑100 (MQ=255)
                             aattttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCa                                                           >  1:1374924/3‑100 (MQ=255)
                              tttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCAt                                                          <  2:1479259/100‑1 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                        >  2:563907/1‑100 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                        <  2:1813228/100‑1 (MQ=255)
                                 ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                       <  1:1744801/100‑1 (MQ=255)
                                 ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                       <  1:1746079/100‑1 (MQ=255)
                                  tCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGaa                                                      >  1:960617/1‑100 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAGATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                     >  1:1884700/1‑100 (MQ=255)
                                    tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCagatcgga                                                    >  1:540216/1‑93 (MQ=255)
                                    tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                    >  1:660623/1‑100 (MQ=255)
                                    tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                    <  2:2005167/100‑1 (MQ=255)
                                     aaaaCAAGACAGTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAACCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                   >  2:1079276/1‑100 (MQ=255)
                                       aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                 >  1:1027277/1‑100 (MQ=255)
                                       aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                 <  1:824780/100‑1 (MQ=255)
                                        aCAATACATGTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                <  2:916935/100‑1 (MQ=255)
                                         cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                               <  1:801205/100‑1 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                          >  2:1757958/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                          >  1:1687728/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaa                                           >  2:1568194/1‑99 (MQ=255)
                                               aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                         >  2:1000789/1‑100 (MQ=255)
                                                tttACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAtt                                        >  2:1061158/1‑100 (MQ=255)
                                                 ttACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAttt                                       <  2:368899/100‑1 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                      >  2:1211740/1‑100 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                      <  1:1727515/100‑1 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                         <  2:521518/97‑1 (MQ=255)
                                                   aCTTTATTTGTCAGTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                     >  2:1265881/1‑100 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                    <  1:608246/100‑1 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                   >  1:1576866/1‑100 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACACCCATTGAACAAATGCTAAAATTTcgcg                                   <  2:309990/100‑1 (MQ=255)
                                                       tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                 >  2:652842/1‑100 (MQ=255)
                                                       tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                 >  2:70015/1‑100 (MQ=255)
                                                         tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                               <  1:495931/100‑1 (MQ=255)
                                                           tGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                             >  2:1405580/1‑100 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                            >  2:1687537/1‑100 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                          >  1:584602/1‑100 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                          >  1:1187201/1‑100 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGGTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                          >  1:549560/1‑100 (MQ=255)
                                                               aCTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCa                         <  1:70015/100‑1 (MQ=255)
                                                                 tGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACg                       >  1:436937/1‑100 (MQ=255)
                                                                  gTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                      >  2:1259201/1‑100 (MQ=255)
                                                                  gTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                      >  2:503117/1‑100 (MQ=255)
                                                                    cGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                    >  1:323621/1‑100 (MQ=255)
                                                                     gTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                   <  1:1391194/100‑1 (MQ=255)
                                                                          tataTCGGCTTAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTc              <  1:1266270/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc             <  2:1968524/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc             <  1:1548840/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc             <  1:1757972/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc             <  2:1301458/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc             <  1:1000794/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt            <  1:1259207/100‑1 (MQ=255)
                                                                              tCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTa          <  2:458791/100‑1 (MQ=255)
                                                                              tCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTa          >  2:238776/1‑100 (MQ=255)
                                                                               cGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt         <  1:1784189/100‑1 (MQ=255)
                                                                               cGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt         <  1:1154850/100‑1 (MQ=255)
                                                                                 gCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa       >  1:161384/1‑100 (MQ=255)
                                                                                 gCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa       <  2:584600/100‑1 (MQ=255)
                                                                                    gAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga    <  2:323620/100‑1 (MQ=255)
                                                                                    gAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga    >  2:798053/1‑100 (MQ=255)
                                                                                     aaaTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgag   <  1:91780/100‑1 (MQ=255)
                                                                                      aaTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga  <  2:232051/100‑1 (MQ=255)
                                                                                               |                                                                                          
TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGA  >  NC_000913/2810663‑2810848

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