BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,539,774	G→A	A578T (GCC→ACC)	feoB →	ferrous iron transport protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,539,774	0	G	A	100.0%	1022.1 / NA	280	A578T (GCC→ACC)	feoB	ferrous iron transport protein B
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (139/141); total (139/141)

Only 100 of 289 total aligned reads displayed.

GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863

  gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCAt                                                                                       >  2:4113417/1‑93 (MQ=255)

   tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATg                                                                                      >  1:1136906/1‑93 (MQ=255)

      gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCg                                                                                   <  1:32940/93‑1 (MQ=255)

      gtcACCCGGGTGATCGCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCg                                                                                   <  1:3257893/93‑1 (MQ=255)

         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                <  1:2309522/93‑1 (MQ=255)

         aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGaaa                                                                                <  2:3353930/93‑1 (MQ=255)

          gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGCaag                                                                               >  2:1640696/1‑93 (MQ=255)

          gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaag                                                                               >  2:5658522/1‑93 (MQ=255)

           ccGGGTAATCACCACGGCCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaaga                                                                              <  1:5105491/93‑1 (MQ=255)

            cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaa                                                                             <  2:132963/93‑1 (MQ=255)

             gggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaag                                                                            <  2:1987688/93‑1 (MQ=255)

             gggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaag                                                                            <  1:1405822/93‑1 (MQ=255)

              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAaagaagg                                                                           >  2:741916/1‑92 (MQ=255)

              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           <  2:221907/93‑1 (MQ=255)

              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           <  1:1622539/93‑1 (MQ=255)

              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           >  1:2392575/1‑93 (MQ=255)

              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           <  2:4710819/93‑1 (MQ=255)

              ggTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGt                                                                           >  2:5794981/1‑93 (MQ=255)

                   tCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAg                                                                      >  2:5714101/1‑93 (MQ=255)

                    cACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGt                                                                     <  1:2130762/93‑1 (MQ=255)

                     aCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTg                                                                    >  2:2466253/1‑93 (MQ=255)

                     aCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTg                                                                    >  2:1847936/1‑93 (MQ=255)

                        ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGt                                                                 <  2:4940561/93‑1 (MQ=255)

                         cGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTa                                                                >  1:2683577/1‑93 (MQ=255)

                         cGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTa                                                                <  1:5883290/93‑1 (MQ=255)

                         cGGTCTTAAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTa                                                                <  1:194464/93‑1 (MQ=255)

                            tCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGc                                                             <  1:5742890/93‑1 (MQ=255)

                             cTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCt                                                            <  1:4002820/93‑1 (MQ=255)

                               tCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                                                          <  1:2274898/93‑1 (MQ=255)

                                 aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                                                        >  2:1735218/1‑93 (MQ=255)

                                     cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                                    <  2:5141755/93‑1 (MQ=255)

                                     cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                                    >  2:855653/1‑93 (MQ=255)

                                     cAATTGGAGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                                    <  1:5286311/93‑1 (MQ=255)

                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTCCGCTCAACAACct                                                   >  1:1535952/1‑93 (MQ=255)

                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   >  1:678923/1‑93 (MQ=255)

                                      aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct                                                   <  2:115307/93‑1 (MQ=255)

                                       aTTGGTGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                                  >  2:36481/1‑93 (MQ=255)

                                       aTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                                  <  2:4433944/93‑1 (MQ=255)

                                       aTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                                  <  2:4877890/93‑1 (MQ=255)

                                       aTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                                  >  1:3866415/1‑93 (MQ=255)

                                        ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                                 >  1:4797114/1‑93 (MQ=255)

                                        ttGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                                 >  1:5970839/1‑93 (MQ=255)

                                          ggCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac                                               <  1:3558969/93‑1 (MQ=255)

                                           gCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTaca                                              >  2:3974979/1‑93 (MQ=255)

                                             gTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACAcc                                            <  1:268028/93‑1 (MQ=255)

                                               gCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc                                          >  2:3687453/1‑93 (MQ=255)

                                                cATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCa                                         <  2:4659434/93‑1 (MQ=255)

                                                cATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCATCTACACCGCa                                         >  1:280354/1‑93 (MQ=255)

                                                 aTGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                        <  2:4310521/93‑1 (MQ=255)

                                                 aTGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                        >  2:610532/1‑93 (MQ=255)

                                                   gAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaa                                      >  1:1128010/1‑93 (MQ=255)

                                                   gAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaa                                      <  2:3119955/93‑1 (MQ=255)

                                                    aaGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaa                                     >  2:6434679/1‑93 (MQ=255)

                                                    aaGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaa                                     >  2:6263360/1‑93 (MQ=255)

                                                    aaGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaa                                     <  1:5205711/93‑1 (MQ=255)

                                                    aaGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaa                                     >  1:4686279/1‑93 (MQ=255)

                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    >  1:2086167/1‑93 (MQ=255)

                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    <  2:3585586/93‑1 (MQ=255)

                                                     aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaaa                                    >  1:2929927/1‑93 (MQ=255)

                                                      gATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAat                                   >  2:938027/1‑93 (MQ=255)

                                                      gATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAat                                   <  1:749589/93‑1 (MQ=255)

                                                      gATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAat                                   <  1:5938065/93‑1 (MQ=255)

                                                         aaCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATAtt                                <  1:4115885/93‑1 (MQ=255)

                                                          aCTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTc                               >  2:3003916/1‑93 (MQ=255)

                                                           cTGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                              >  2:5982629/1‑93 (MQ=255)

                                                            tGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAg                             <  2:1066614/93‑1 (MQ=255)

                                                            tGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAg                             >  1:3152469/1‑93 (MQ=255)

                                                            tGGCAGGCAACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAg                             >  1:2183646/1‑93 (MQ=255)

                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  1:1728355/1‑93 (MQ=255)

                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  1:4408494/1‑93 (MQ=255)

                                                             ggcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAgg                            >  2:1425206/1‑93 (MQ=255)

                                                              gcaggcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGa                           >  1:6602082/1‑93 (MQ=255)

                                                                  gcaACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACgca                       >  2:2805434/1‑91 (MQ=255)

                                                                   caACGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAag                      <  2:5493054/93‑1 (MQ=255)

                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     <  2:376557/93‑1 (MQ=255)

                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     <  1:5536098/93‑1 (MQ=255)

                                                                    aaCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAaga                     <  1:4227996/93‑1 (MQ=255)

                                                                     aCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag                    <  1:5671099/93‑1 (MQ=255)

                                                                     aCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag                    <  1:2530212/93‑1 (MQ=255)

                                                                     aCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag                    >  1:3702606/1‑93 (MQ=255)

                                                                     aCGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag                    <  1:3938469/93‑1 (MQ=255)

                                                                      cGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt                   <  1:3801499/93‑1 (MQ=255)

                                                                      cGGTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt                   <  2:2722648/93‑1 (MQ=255)

                                                                       ggTTGGCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGtt                  >  1:6316371/1‑93 (MQ=255)

                                                                           ggCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAAt              <  1:3083906/93‑1 (MQ=255)

                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  1:124984/93‑1 (MQ=255)

                                                                            gCCTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATc             <  1:820158/93‑1 (MQ=255)

                                                                             ccTGTTTACAGGTACCTTGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCTTGTAGACCGGAGCACATAGTCAGGAAGCGGAGTTCCATcg            >  1:4723881/1‑92 (MQ=255)

                                                                             ccTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATcc            >  2:5402732/1‑93 (MQ=255)

                                                                             ccTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATcc            >  2:5030588/1‑93 (MQ=255)

                                                                             ccTGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATcc            >  1:653627/1‑93 (MQ=255)

                                                                               tGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          >  1:4552711/1‑93 (MQ=255)

                                                                               tGTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCgg          >  2:2331501/1‑93 (MQ=255)

                                                                                gTTTACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGc         <  1:4949649/93‑1 (MQ=255)

                                                                                   tACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGa      >  2:5608692/1‑93 (MQ=255)

                                                                                   tACAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCAGGCAGa      >  2:3207396/1‑93 (MQ=255)

                                                                                   tACAGGTACCATGGCGAAAGAAGTGGCAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGa      >  1:1027148/1‑93 (MQ=255)

                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  1:5989253/93‑1 (MQ=255)

                                                                                    aCAGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGaa     <  1:6155006/93‑1 (MQ=255)

                                                                                      aGGTACCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAAtt   <  2:3495398/93‑1 (MQ=255)

GCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAGAATTT  >  NC_000913/3539684‑3539863

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