BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,656,923	T→C	P277P (CCT→CCC)	b1587 →	putative oxidoreductase major subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,656,923	0	T	C	100.0%	717.8 / NA	218	P277P (CCT→CCC)	b1587	putative oxidoreductase major subunit
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (88/130); total (88/130)

Only 100 of 235 total aligned reads displayed.

AAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAA  >  NC_000913/1656833‑1657012

aaaaGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCAt                                                                                         >  2:1842788/1‑93 (MQ=255)

 aaaGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCAtt                                                                                        >  1:2636034/1‑93 (MQ=255)

 aaaGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCAtt                                                                                        >  1:3876672/1‑93 (MQ=255)

      gcgcgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGccc                                                                                   <  2:2702211/93‑1 (MQ=255)

       cgcgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCg                                                                                  >  1:19703/1‑93 (MQ=255)

        gcgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCgg                                                                                 <  2:5955693/93‑1 (MQ=255)

          gagaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGc                                                                               >  2:2369819/1‑93 (MQ=255)

          gagaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGc                                                                               <  1:2150849/93‑1 (MQ=255)

          gagaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGc                                                                               >  1:4101137/1‑93 (MQ=255)

          gagaAATCGAACACCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGc                                                                               <  1:18685/93‑1 (MQ=255)

           agaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCa                                                                              >  2:111156/1‑93 (MQ=255)

            gaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCAc                                                                             >  1:3109771/1‑93 (MQ=255)

             aaaTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATAGCTCCCCATTCGCCCGGGCAcc                                                                            >  1:5899167/1‑93 (MQ=255)

              aaTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCg                                                                           >  1:2913107/1‑93 (MQ=255)

              aaTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCg                                                                           >  2:4277354/1‑93 (MQ=255)

               accaaaCGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGa                                                                          <  1:4002187/89‑1 (MQ=255)

               aTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGa                                                                          <  2:132685/93‑1 (MQ=255)

               aTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGa                                                                          >  2:473226/1‑93 (MQ=255)

                 cGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATg                                                                        >  1:1037063/1‑93 (MQ=255)

                   aaCGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGcc                                                                      >  1:5068629/1‑93 (MQ=255)

                    aCGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCcg                                                                     >  2:5878266/1‑93 (MQ=255)

                    aCGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGATGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCcg                                                                     >  1:4021269/1‑93 (MQ=255)

                     cGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCcgc                                                                    >  1:5169127/1‑93 (MQ=255)

                       ccAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCgcgc                                                                  <  1:2930860/93‑1 (MQ=255)

                        cAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCt                                                                 >  2:4922160/1‑93 (MQ=255)

                          aaaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTgg                                                               <  2:923400/93‑1 (MQ=255)

                           aaTTATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGt                                                              <  2:5082178/93‑1 (MQ=255)

                           aaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGTCGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGt                                                              <  2:3554533/93‑1 (MQ=255)

                            aTGATTGTTATCGACCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGtt                                                             <  2:6213494/93‑1 (MQ=255)

                              gATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGc                                                           <  2:2256405/93‑1 (MQ=255)

                                 tGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCggg                                                        <  1:2717373/93‑1 (MQ=255)

                                  gTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGGCCGGGCACCGATGCCGCGCGCGCTGCgggg                                                       >  2:3497179/1‑92 (MQ=255)

                                   ttATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCTGGTa                                                      <  1:5527795/93‑1 (MQ=255)

                                   ttATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTa                                                      <  2:4064103/93‑1 (MQ=255)

                                    tATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTAt                                                     <  2:2754742/93‑1 (MQ=255)

                                     aTCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTAtt                                                    >  1:4843242/1‑93 (MQ=255)

                                      tCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTAttc                                                   >  1:922691/1‑92 (MQ=255)

                                      tCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTg                                                   <  1:6675552/93‑1 (MQ=255)

                                        gATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGcc                                                 <  1:2421796/93‑1 (MQ=255)

                                         gtCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCt                                                <  2:216817/92‑1 (MQ=255)

                                          tCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTg                                               >  1:5313163/1‑93 (MQ=255)

                                          tCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTg                                               >  2:2991241/1‑93 (MQ=255)

                                           ccGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTgg                                              <  1:3930926/93‑1 (MQ=255)

                                            cgcgATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTggg                                             <  1:5771756/93‑1 (MQ=255)

                                               gATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTAt                                          <  1:3574908/93‑1 (MQ=255)

                                               gATATACCGATACGGCTGCCGGTCGGGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTAt                                          <  1:3369626/93‑1 (MQ=255)

                                                atatACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTAtt                                         <  2:5851619/93‑1 (MQ=255)

                                                 tataCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTg                                        >  1:1745365/1‑93 (MQ=255)

                                                  ataCCGATACGGCTGCCGGTCGTGAAGACGAATGGTTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGa                                       <  1:86985/93‑1 (MQ=255)

                                                  ataCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGa                                       >  1:4220349/1‑93 (MQ=255)

                                                    aCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGAtt                                     <  2:1438325/93‑1 (MQ=255)

                                                     ccGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTa                                    <  2:3508568/93‑1 (MQ=255)

                                                     ccGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTa                                    >  1:1385370/1‑93 (MQ=255)

                                                     ccGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTa                                    >  2:5898395/1‑93 (MQ=255)

                                                     ccGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTa                                    >  1:3958569/1‑93 (MQ=255)

                                                      cGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTaa                                   <  2:2503738/93‑1 (MQ=255)

                                                        aTACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACg                                 <  2:426375/93‑1 (MQ=255)

                                                        aTACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACg                                 >  1:2069209/1‑93 (MQ=255)

                                                         tACGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGa                                <  1:6458025/93‑1 (MQ=255)

                                                          aCGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGaa                               <  1:6246662/93‑1 (MQ=255)

                                                           cGGCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGaaa                              >  2:5324553/1‑93 (MQ=255)

                                                            ggCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGaaaa                             >  2:221859/1‑93 (MQ=255)

                                                            ggCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGaaaa                             <  1:2487510/93‑1 (MQ=255)

                                                            ggCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGaaaa                             <  2:4620894/93‑1 (MQ=255)

                                                             gCTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAt                            >  2:5767055/1‑93 (MQ=255)

                                                             gCTGCCGGTCGTGAAGACGAATGCCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAt                            <  1:3333850/93‑1 (MQ=255)

                                                              cTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtc                           <  2:5170644/93‑1 (MQ=255)

                                                              cTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtc                           >  2:5571569/1‑93 (MQ=255)

                                                              cTGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtc                           <  2:3694286/93‑1 (MQ=255)

                                                               tGCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtct                          >  1:1502795/1‑93 (MQ=255)

                                                               tGCCGGTCGTGAAGACGAATGGCTCCCCAGTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtct                          <  2:3163025/93‑1 (MQ=255)

                                                                gCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCATCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtctc                         <  1:1196762/93‑1 (MQ=255)

                                                                gCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtctc                         <  1:6490427/93‑1 (MQ=255)

                                                                gCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtctc                         <  2:2357160/93‑1 (MQ=255)

                                                                gCCGGTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtctc                         <  1:3594040/93‑1 (MQ=255)

                                                                    gTCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTg                     >  2:3799863/1‑93 (MQ=255)

                                                                     tCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGa                    <  1:5207525/93‑1 (MQ=255)

                                                                     tCGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGa                    <  1:5127006/93‑1 (MQ=255)

                                                                      cGTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGAt                   <  2:864411/93‑1 (MQ=255)

                                                                       gTGAAGACGAATGGCTCTCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATc                  <  1:2220266/93‑1 (MQ=255)

                                                                       gTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATc                  <  2:2059214/93‑1 (MQ=255)

                                                                       gTGAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGGATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATc                  >  1:3180578/1‑93 (MQ=255)

                                                                         gAAGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCaa                <  1:4478875/93‑1 (MQ=255)

                                                                          aaGACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAAc               <  2:3954901/93‑1 (MQ=255)

                                                                            gACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCt             <  2:4447976/93‑1 (MQ=255)

                                                                            gACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCt             <  2:6506960/93‑1 (MQ=255)

                                                                            gACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCt             <  1:3474617/93‑1 (MQ=255)

                                                                            gACGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCt             <  1:1574410/93‑1 (MQ=255)

                                                                              cGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCttt           <  1:2562100/93‑1 (MQ=255)

                                                                              cGAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCttt           <  2:6329652/93‑1 (MQ=255)

                                                                              cGAATGGCACCCCATTCGCGCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCttt           <  1:165740/93‑1 (MQ=255)

                                                                               gAATGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCtttt          >  2:6088953/1‑93 (MQ=255)

                                                                                agtGGCTCCGCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtc         <  1:4566129/91‑1 (MQ=255)

                                                                                  tGGCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtctc       <  2:6332686/93‑1 (MQ=255)

                                                                                    gCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGa     <  2:171348/93‑1 (MQ=255)

                                                                                    gCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGa     >  2:5652579/1‑93 (MQ=255)

                                                                                    gCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGa     <  1:2781994/93‑1 (MQ=255)

                                                                                    gCTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGa     <  2:3288286/93‑1 (MQ=255)

                                                                                     cTCCCCATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGAt    >  1:3221612/1‑93 (MQ=255)

                                                                                      cccccATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATa   <  2:4596097/92‑1 (MQ=255)

AAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAA  >  NC_000913/1656833‑1657012

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