BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	577,031	T→C	Y66H (TAT→CAT)	ybcS →	lysozyme‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	577,031	0	T	C	100.0%	892.4 / NA	272	Y66H (TAT→CAT)	ybcS	lysozyme‑like protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (134/138); total (134/138)

Only 100 of 281 total aligned reads displayed.

GGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCG  >  NC_000913/576941‑577120

ggTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGcat                                                                                         <  2:550481/93‑1 (MQ=255)

ggTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGcat                                                                                         <  1:5360135/93‑1 (MQ=255)

ggTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGcat                                                                                         >  2:5887430/1‑93 (MQ=255)

   gtCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATAcc                                                                                      <  1:2307198/93‑1 (MQ=255)

   gtCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATAcc                                                                                      <  2:3009683/93‑1 (MQ=255)

    tCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCa                                                                                     >  2:5631208/1‑93 (MQ=255)

    tCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCa                                                                                     >  1:4856506/1‑93 (MQ=255)

      aGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCaaa                                                                                   <  1:18930/93‑1 (MQ=255)

      aGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCaaa                                                                                   <  2:2934869/93‑1 (MQ=255)

      aGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCaaa                                                                                   >  1:2708052/1‑93 (MQ=255)

       gCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAg                                                                                  <  2:3846194/93‑1 (MQ=255)

         tACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCa                                                                                >  2:4104148/1‑93 (MQ=255)

         tACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCa                                                                                <  1:4504623/93‑1 (MQ=255)

         tACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCa                                                                                <  2:1672759/93‑1 (MQ=255)

         tACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCa                                                                                <  2:127177/93‑1 (MQ=255)

           cataccatacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGa                                                                              <  1:4370868/93‑1 (MQ=255)

             taccatacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAAt                                                                            >  2:5319404/1‑93 (MQ=255)

             taccatacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAAt                                                                            <  1:3884980/93‑1 (MQ=255)

             taccatacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAAt                                                                            <  1:3080644/93‑1 (MQ=255)

             taccatacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAAt                                                                            >  2:5733119/1‑93 (MQ=255)

             taccatacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAAt                                                                            >  1:2511662/1‑93 (MQ=255)

                 atacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCaa                                                                        >  1:4982232/1‑93 (MQ=255)

                 atacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCaa                                                                        >  1:5402696/1‑93 (MQ=255)

                  tacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCaaa                                                                       >  1:2780122/1‑93 (MQ=255)

                  tacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCaaa                                                                       >  1:2447956/1‑93 (MQ=255)

                  tacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCaaa                                                                       <  2:6171662/93‑1 (MQ=255)

                  tacAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCaaa                                                                       >  1:6284873/1‑93 (MQ=255)

                   acAAAGATATTGTTGGTGTATGCACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAg                                                                      <  2:2116121/93‑1 (MQ=255)

                    cAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGc                                                                     <  2:4512853/93‑1 (MQ=255)

                    cAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGc                                                                     >  2:5599678/1‑93 (MQ=255)

                         atatTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCActct                                                                <  1:1553480/93‑1 (MQ=255)

                         atatTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCActct                                                                >  1:5440798/1‑93 (MQ=255)

                          tatTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCtt                                                               <  2:2532507/93‑1 (MQ=255)

                            ttgttgGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGa                                                             <  1:1601744/93‑1 (MQ=255)

                              gttgGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAAt                                                           <  2:2418711/93‑1 (MQ=255)

                              gttgGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAAt                                                           >  2:6127267/1‑93 (MQ=255)

                                 ggTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATaaa                                                        >  2:743869/1‑93 (MQ=255)

                                 ggTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATaaa                                                        <  2:4240813/93‑1 (MQ=255)

                                  gtgtATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAg                                                       >  1:5116316/1‑93 (MQ=255)

                                    gtATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGAc                                                     >  1:2657563/1‑93 (MQ=255)

                                    gtATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGAc                                                     <  2:691987/93‑1 (MQ=255)

                                      aTGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCt                                                   >  2:3190558/1‑93 (MQ=255)

                                      aTGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCt                                                   >  2:2398663/1‑93 (MQ=255)

                                       tGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCtt                                                  <  2:1398721/93‑1 (MQ=255)

                                         gACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGc                                                >  1:564031/1‑93 (MQ=255)

                                         gACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGc                                                >  1:941347/1‑93 (MQ=255)

                                          aCTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGcc                                               >  1:5278761/1‑93 (MQ=255)

                                          aCTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGcc                                               <  2:4869731/93‑1 (MQ=255)

                                          aCTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGcc                                               >  2:1616247/1‑93 (MQ=255)

                                           cTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCa                                              >  1:1273890/1‑93 (MQ=255)

                                            tGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCAc                                             >  2:2617467/1‑93 (MQ=255)

                                             gTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACt                                            >  2:885664/1‑93 (MQ=255)

                                             gTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACt                                            >  1:5911021/1‑93 (MQ=255)

                                              tATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTg                                           >  2:4021508/1‑93 (MQ=255)

                                              tATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGAACTTGCAACTg                                           >  2:2034871/1‑93 (MQ=255)

                                               aTGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGt                                          >  2:1460005/1‑93 (MQ=255)

                                                 gTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCg                                        >  2:4975560/1‑93 (MQ=255)

                                                   cACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGcc                                      <  1:4620453/93‑1 (MQ=255)

                                                    aCGGACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCa                                     <  1:3683475/93‑1 (MQ=255)

                                                      ggACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGa                                   <  2:4177058/93‑1 (MQ=255)

                                                      ggACACACCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGACACTGTCGCCAGa                                   <  2:337069/93‑1 (MQ=255)

                                                        acacacCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACa                                 >  2:1812690/1‑93 (MQ=255)

                                                           cacCGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAAt                              <  2:6016701/93‑1 (MQ=255)

                                                             ccGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTa                            <  2:5036964/93‑1 (MQ=255)

                                                             ccGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTa                            >  1:4018105/1‑93 (MQ=255)

                                                              cGGAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTaa                           >  1:5561451/1‑93 (MQ=255)

                                                               ggAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAAc                          >  1:5493903/1‑93 (MQ=255)

                                                                gAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAAcc                         >  2:4325631/1‑93 (MQ=255)

                                                                gAAAAGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAAcc                         <  1:338546/93‑1 (MQ=255)

                                                                 aaaaGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAAccc                        <  2:1723520/93‑1 (MQ=255)

                                                                   aaGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGt                      >  1:6017063/1‑93 (MQ=255)

                                                                   aaGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGt                      >  1:4106311/1‑93 (MQ=255)

                                                                    aGACATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGta                     >  2:345456/1‑93 (MQ=255)

                                                                      aCATCATGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGtata                   <  1:3179681/93‑1 (MQ=255)

                                                                       catcatGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTatat                  >  1:716892/1‑93 (MQ=255)

                                                                       catcatGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTatat                  >  1:4234752/1‑93 (MQ=255)

                                                                       catcatGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTatat                  <  1:852065/93‑1 (MQ=255)

                                                                        atcatGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATc                 >  1:1823062/1‑93 (MQ=255)

                                                                          catGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCaa               >  2:5090522/1‑93 (MQ=255)

                                                                           atGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCaaa              <  1:3035227/93‑1 (MQ=255)

                                                                           atGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCaaa              >  2:1145280/1‑93 (MQ=255)

                                                                           atGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCaaa              >  2:1291651/1‑93 (MQ=255)

                                                                            tGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAg             >  1:1002090/1‑93 (MQ=255)

                                                                            tGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAg             <  1:6248415/93‑1 (MQ=255)

                                                                            tGCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAg             <  1:5861853/93‑1 (MQ=255)

                                                                             gCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGt            <  1:2008270/93‑1 (MQ=255)

                                                                             gCTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGt            >  2:6096124/1‑93 (MQ=255)

                                                                              cTCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTc           >  2:4695661/1‑93 (MQ=255)

                                                                               tCGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCg          >  1:5775313/1‑93 (MQ=255)

                                                                                cGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGa         <  1:5049790/93‑1 (MQ=255)

                                                                                cGGTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGa         >  1:5959173/1‑93 (MQ=255)

                                                                                 ggTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGat        <  2:3727209/93‑1 (MQ=255)

                                                                                  gTAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGata       >  1:2232862/1‑93 (MQ=255)

                                                                                   tAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGatat      <  2:1328293/93‑1 (MQ=255)

                                                                                   tAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGatat      >  1:151476/1‑93 (MQ=255)

                                                                                   tAAAACGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGatat      >  2:4271966/1‑93 (MQ=255)

                                                                                     aaaCGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATAc    >  2:2672051/1‑93 (MQ=255)

                                                                                      aaCGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATAcc   >  1:3136676/1‑93 (MQ=255)

                                                                                      aaCGCATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATAcc   <  1:3922501/93‑1 (MQ=255)

                                                                                       aCGCATACCAAAGCCGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCg  <  1:3607820/93‑1 (MQ=255)

GGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCG  >  NC_000913/576941‑577120

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