BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,014,968	G→A	M296I (ATG→ATA)	b2878 →	putative oxidoreductase Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,014,968	0	G	A	100.0%	1090.0 / NA	282	M296I (ATG→ATA)	b2878	putative oxidoreductase Fe‑S subunit
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (154/128); total (154/128)

Only 100 of 285 total aligned reads displayed.

TGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATGCTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAATAAAGG  >  NC_000913/3014878‑3015057

tGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACt                                                                                         <  1:2344215/93‑1 (MQ=255)

tGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACt                                                                                         >  2:725073/1‑93 (MQ=255)

  agaTCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTgg                                                                                       >  2:5859307/1‑93 (MQ=255)

  agaTCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTgg                                                                                       >  2:3230396/1‑93 (MQ=255)

  agaTCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTgg                                                                                       >  2:2802370/1‑93 (MQ=255)

   gaTCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGa                                                                                      >  2:4888374/1‑93 (MQ=255)

   gaTCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGa                                                                                      >  2:2560188/1‑93 (MQ=255)

    aTCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGaa                                                                                     <  2:2244749/93‑1 (MQ=255)

    aTCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGaa                                                                                     <  2:3933961/93‑1 (MQ=255)

    aTCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGaa                                                                                     <  2:6201167/93‑1 (MQ=255)

    aTCCTCGATGACTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGaa                                                                                     >  2:1844210/1‑93 (MQ=255)

       cTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGc                                                                                  <  1:5869116/93‑1 (MQ=255)

       cTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGc                                                                                  >  2:6215028/1‑93 (MQ=255)

        tCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGcc                                                                                 >  2:2894870/1‑93 (MQ=255)

         cGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCt                                                                                <  1:6034716/93‑1 (MQ=255)

           aTGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGa                                                                              >  2:4257179/1‑93 (MQ=255)

           aTGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGa                                                                              >  2:3065279/1‑93 (MQ=255)

           aTGTCTGCGGGTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGa                                                                              >  1:5550008/1‑93 (MQ=255)

            tGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGAt                                                                             >  2:1546180/1‑93 (MQ=255)

            tGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGAt                                                                             <  2:5962807/93‑1 (MQ=255)

             gTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATg                                                                            <  2:2749816/93‑1 (MQ=255)

             gTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATg                                                                            >  1:4225086/1‑93 (MQ=255)

             gTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATg                                                                            >  2:3724878/1‑93 (MQ=255)

                  cGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACt                                                                       <  1:4542432/93‑1 (MQ=255)

                   ggTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTg                                                                      <  1:923633/93‑1 (MQ=255)

                   ggTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATAGCACTg                                                                      >  1:6395088/1‑93 (MQ=255)

                       tCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCaa                                                                  <  2:4452897/93‑1 (MQ=255)

                        cGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCaaa                                                                 >  1:1531217/1‑93 (MQ=255)

                         ggTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCaaaa                                                                <  2:3309316/93‑1 (MQ=255)

                           tttcATAAGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGa                                                              <  1:4355941/90‑1 (MQ=255)

                           ttACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGa                                                              >  1:4079497/1‑93 (MQ=255)

                             aCATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGAAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGaaa                                                            <  2:6364586/93‑1 (MQ=255)

                              cATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGaaaa                                                           >  2:5327861/1‑93 (MQ=255)

                               aTAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGaaaaa                                                          <  1:5366569/93‑1 (MQ=255)

                                tAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAAt                                                         >  1:1211232/1‑93 (MQ=255)

                                tAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAAt                                                         <  2:2303596/93‑1 (MQ=255)

                                tAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAAt                                                         <  2:1943611/93‑1 (MQ=255)

                                 aGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATc                                                        >  1:1358153/1‑93 (MQ=255)

                                    cTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACt                                                     >  2:1681483/1‑93 (MQ=255)

                                    cTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACt                                                     <  1:1935748/93‑1 (MQ=255)

                                      tAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCg                                                   >  2:2430935/1‑93 (MQ=255)

                                        aaaGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGc                                                 <  1:5968812/93‑1 (MQ=255)

                                        aaaGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGc                                                 <  2:2861484/93‑1 (MQ=255)

                                        aaaGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGc                                                 >  2:76866/1‑93 (MQ=255)

                                         aagaagAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCt                                                <  2:2089937/93‑1 (MQ=255)

                                         aagaagAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCt                                                >  2:2158259/1‑93 (MQ=255)

                                          agaagaGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCtt                                               <  1:5458238/93‑1 (MQ=255)

                                          agaagaGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCtt                                               <  1:450406/93‑1 (MQ=255)

                                          agaagaGTCATTTGATCACGACCTCAAGCTGACCCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCtt                                               <  1:3930537/93‑1 (MQ=255)

                                           gaagaGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCttt                                              >  2:5090425/1‑93 (MQ=255)

                                           gaagaGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCttt                                              >  1:4907909/1‑93 (MQ=255)

                                           gaagaGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCttt                                              >  1:2189570/1‑93 (MQ=255)

                                            aagaGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTg                                             <  1:3588157/93‑1 (MQ=255)

                                             agagTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTgg                                            >  2:5388394/1‑93 (MQ=255)

                                              gagTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGc                                           >  1:2981146/1‑93 (MQ=255)

                                               agTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCg                                          <  1:4383664/93‑1 (MQ=255)

                                                gTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGt                                         >  2:6059392/1‑93 (MQ=255)

                                                   aTTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTaaa                                      >  1:5632892/1‑93 (MQ=255)

                                                    tttGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTaaaa                                     >  1:113593/1‑93 (MQ=255)

                                                    tttGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTCGGCGTaaaa                                     >  1:1014422/1‑93 (MQ=255)

                                                      tGATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACt                                   <  1:5247269/93‑1 (MQ=255)

                                                       gATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTg                                  <  1:2628878/93‑1 (MQ=255)

                                                       gATCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTg                                  <  1:2841382/93‑1 (MQ=255)

                                                        aTCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGa                                 >  1:2014501/1‑93 (MQ=255)

                                                        aTCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGa                                 >  2:1288275/1‑93 (MQ=255)

                                                         tCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGAc                                <  1:4493047/93‑1 (MQ=255)

                                                         tCACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGAc                                <  1:3207478/93‑1 (MQ=255)

                                                          cACGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACt                               >  1:658962/1‑93 (MQ=255)

                                                           aCGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTa                              >  1:898847/1‑93 (MQ=255)

                                                            cGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTaa                             >  2:1592320/1‑93 (MQ=255)

                                                            cGACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTaa                             >  1:5716405/1‑93 (MQ=255)

                                                             gACCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAac                            >  2:2685151/1‑93 (MQ=255)

                                                              aCCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAaca                           >  2:6295239/1‑93 (MQ=255)

                                                              aCCTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAaca                           >  2:1102510/1‑93 (MQ=255)

                                                               ccTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAacac                          <  2:3651141/93‑1 (MQ=255)

                                                               ccTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAacac                          <  2:2019252/93‑1 (MQ=255)

                                                                cTCAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACAct                         >  2:3629711/1‑93 (MQ=255)

                                                                  cAAGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACActct                       <  1:1731646/93‑1 (MQ=255)

                                                                   aaGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACtctc                      >  1:1816799/1‑93 (MQ=255)

                                                                    aGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCg                     >  1:3761767/1‑93 (MQ=255)

                                                                    aGCTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGAATAACACTCTCg                     <  2:1642610/93‑1 (MQ=255)

                                                                      cTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGc                   >  1:2657548/1‑93 (MQ=255)

                                                                      cTGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGc                   <  2:4398172/93‑1 (MQ=255)

                                                                       tGACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCa                  >  1:2449625/1‑93 (MQ=255)

                                                                        gACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCAc                 >  2:5934320/1‑93 (MQ=255)

                                                                        gACGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCAc                 >  2:1790152/1‑93 (MQ=255)

                                                                         aCGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCAcc                <  2:5028579/93‑1 (MQ=255)

                                                                          cGCAAGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCa               >  1:4718223/1‑93 (MQ=255)

                                                                          cGCAAGCACTGGAAATACTGGAACGCCTAATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCa               <  1:6106225/93‑1 (MQ=255)

                                                                             aaGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATca            >  2:870513/1‑93 (MQ=255)

                                                                             aaGCACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATca            <  1:2512810/93‑1 (MQ=255)

                                                                             aaGCACTGCAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATca            >  2:1234668/1‑93 (MQ=255)

                                                                                cACTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATcaaca         <  2:3420351/93‑1 (MQ=255)

                                                                                  cTGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAAt       <  2:1380139/93‑1 (MQ=255)

                                                                                   tGGAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAATa      <  2:620630/93‑1 (MQ=255)

                                                                                    ggAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCTCTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAATaa     <  2:4593718/93‑1 (MQ=255)

                                                                                     gAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAATaaa    <  2:1028790/93‑1 (MQ=255)

                                                                                     gAAATACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAATaaa    <  2:4962141/93‑1 (MQ=255)

                                                                                       aaTACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAATAAAgg  <  2:554478/93‑1 (MQ=255)

                                                                                       aaTACTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAATAAAgg  >  2:4662459/1‑93 (MQ=255)

TGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATGCTGGAACGCCTGATGGCACTGGCAAAAGAAAAATCACTCGGCTTTGGCGTAAAACTGACTAACACTCTCGGCACCATCAACAATAAAGG  >  NC_000913/3014878‑3015057

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