BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,734	A→G	noncoding (409/2906 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,734	0	A	G	100.0%	600.8 / NA	224	noncoding (409/2906 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (127/97); total (127/97)

Only 100 of 235 total aligned reads displayed.

GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAG  >  NC_000913/3941644‑3941819

gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                                                     >  1:4021374/1‑93 (MQ=25)

 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGAGACGTGGTATCCTGCATGAATATggggggg                                                                                    >  2:3860834/1‑93 (MQ=11)

 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATggggggg                                                                                    >  1:1003410/1‑93 (MQ=25)

    cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGAcc                                                                                 >  1:1707448/1‑93 (MQ=255)

     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCa                                                                                >  1:2548844/1‑93 (MQ=255)

     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCa                                                                                <  2:3406201/93‑1 (MQ=255)

     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCa                                                                                <  1:1490614/93‑1 (MQ=255)

     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGGATCCTGTCTGAATATGGGGGGGACCa                                                                                >  2:3490481/1‑93 (MQ=255)

       ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCAtc                                                                              <  2:390760/93‑1 (MQ=255)

            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCa                                                                         >  2:2093441/1‑93 (MQ=255)

             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCaa                                                                        >  1:4567562/1‑93 (MQ=255)

                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGc                                                                     <  2:1722557/93‑1 (MQ=255)

                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGc                                                                     >  1:3080370/1‑93 (MQ=255)

                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGc                                                                     <  2:2577126/93‑1 (MQ=255)

                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGc                                                                     <  1:4047154/93‑1 (MQ=255)

                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGc                                                                     >  2:1143109/1‑93 (MQ=255)

                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGc                                                                     >  2:2897219/1‑93 (MQ=255)

                cccGTACACAAAAATGCACATATTGGGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGCCCATCCTCCAAGGc                                                                     <  1:695895/93‑1 (MQ=255)

                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCTTCCAAGGCt                                                                    >  1:3427945/1‑93 (MQ=255)

                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTa                                                                   >  1:4505561/1‑93 (MQ=255)

                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTaaa                                                                 <  2:32736/93‑1 (MQ=255)

                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTaaa                                                                 >  1:910937/1‑93 (MQ=255)

                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTaaa                                                                 >  2:441797/1‑93 (MQ=255)

                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAAt                                                                <  1:3517312/93‑1 (MQ=255)

                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATa                                                               <  1:1676546/93‑1 (MQ=255)

                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATAc                                                              >  1:3223730/1‑93 (MQ=255)

                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACt                                                             >  2:2348306/1‑93 (MQ=255)

                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACt                                                             >  1:1909763/1‑93 (MQ=255)

                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTc                                                            >  2:3101983/1‑93 (MQ=255)

                          aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTcc                                                           >  1:980858/1‑93 (MQ=255)

                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctg                                                         <  2:1452451/93‑1 (MQ=255)

                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctga                                                        >  1:1352202/1‑93 (MQ=255)

                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgac                                                       >  1:1500385/1‑93 (MQ=255)

                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgac                                                       <  2:3733759/93‑1 (MQ=255)

                               gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATGCTCctgact                                                      <  2:243271/93‑1 (MQ=255)

                               gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgact                                                      >  2:3995417/1‑93 (MQ=255)

                               gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgact                                                      >  2:1202123/1‑93 (MQ=255)

                               gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgact                                                      >  1:1166976/1‑93 (MQ=255)

                               gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgact                                                      >  1:884789/1‑93 (MQ=255)

                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgactga                                                    >  2:138890/1‑93 (MQ=255)

                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCctgactga                                                    <  1:502925/93‑1 (MQ=255)

                                   atatTGTGAGCTCGATTAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                  <  1:4548445/93‑1 (MQ=255)

                                   atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                  <  2:1903075/93‑1 (MQ=255)

                                     atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGa                                                <  2:2260943/93‑1 (MQ=21)

                                     atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGa                                                <  1:1148889/93‑1 (MQ=25)

                                     atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGa                                                <  1:534246/93‑1 (MQ=21)

                                        gtgAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                             >  1:663194/1‑93 (MQ=18)

                                        gtgAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                             <  1:806901/93‑1 (MQ=18)

                                        gtgAGCTCGATGAGTAGAGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAg                                             <  1:4069199/93‑1 (MQ=16)

                                          gAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGCCTGACCGATAGTg                                           <  2:4603961/93‑1 (MQ=12)

                                            gCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGaa                                         >  1:2206846/1‑93 (MQ=18)

                                            gCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGaa                                         >  2:1768828/1‑93 (MQ=18)

                                               cGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCa                                      <  1:2888911/93‑1 (MQ=17)

                                               cGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCa                                      <  2:3145299/93‑1 (MQ=17)

                                                gATGAGTAGGGCGGGCCACGTTGTAGCCTGCCTGAATCTGGGGGGGACCATCCTCCAACGCTAAATACTCCTCACTGTCCGAGAGTTAACCAg                                     >  2:3200620/1‑93 (MQ=11)

                                                gATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAg                                     <  1:566689/93‑1 (MQ=18)

                                                 aTGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGCCCGATAGTGAACCAGt                                    >  2:4278974/1‑93 (MQ=16)

                                                  tGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTa                                   <  2:4654058/93‑1 (MQ=17)

                                                   gAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTAc                                  <  2:3959558/93‑1 (MQ=18)

                                                   gAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTAc                                  >  2:2191333/1‑93 (MQ=18)

                                                     gcaGGGCGGGACACGTGGTATCCTGTCTGGATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCg                                >  2:4367807/3‑93 (MQ=255)

                                                     gTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCg                                >  2:2526284/1‑93 (MQ=18)

                                                     gTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCg                                >  2:1798732/1‑93 (MQ=18)

                                                      gaGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGt                               <  2:1754531/92‑1 (MQ=18)

                                                       aGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCTATAGTGAACCAGTACCGTg                              <  2:4527321/93‑1 (MQ=16)

                                                         ggCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAg                            >  1:2301786/1‑93 (MQ=18)

                                                          gCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAgg                           >  2:2813868/1‑93 (MQ=18)

                                                          gCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAgg                           >  1:3600631/1‑93 (MQ=18)

                                                           cGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAggg                          <  1:292530/93‑1 (MQ=18)

                                                              gACACGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGaaa                       <  1:2691867/93‑1 (MQ=18)

                                                                cacGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAgg                     <  2:1125961/93‑1 (MQ=18)

                                                                cacGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAgg                     <  1:2038734/93‑1 (MQ=18)

                                                                  cGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCg                   >  2:2461268/1‑93 (MQ=17)

                                                                   gTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCCAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGa                  <  1:3199617/93‑1 (MQ=255)

                                                                    tGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaa                 >  1:2820089/1‑93 (MQ=17)

                                                                    tGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaa                 >  2:848180/1‑93 (MQ=17)

                                                                     ggTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaaa                >  2:4383271/1‑93 (MQ=18)

                                                                      gTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaaaa               <  1:2375540/93‑1 (MQ=18)

                                                                      gTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaaaa               >  2:826202/1‑93 (MQ=17)

                                                                      gTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaaaa               >  1:2014500/1‑93 (MQ=17)

                                                                      gTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGaaaa               >  2:4147148/1‑93 (MQ=18)

                                                                       tATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAg              <  1:4286733/93‑1 (MQ=18)

                                                                       tATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAg              >  2:2634184/1‑93 (MQ=17)

                                                                        aTCCTGTCTGAATATGGGGGGGGCCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGa             >  1:4181236/1‑93 (MQ=12)

                                                                        aTCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGa             >  1:2130539/1‑93 (MQ=18)

                                                                         tCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGaa            <  1:4612528/93‑1 (MQ=18)

                                                                         tCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGaa            <  1:446439/93‑1 (MQ=18)

                                                                         tCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGaa            <  1:1771247/93‑1 (MQ=18)

                                                                          ccTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAc           >  1:3150633/1‑93 (MQ=18)

                                                                             gTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAcccc        <  1:4268552/93‑1 (MQ=18)

                                                                             gTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAAcccc        >  1:1246478/1‑93 (MQ=17)

                                                                              tCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCg       >  2:2768418/1‑93 (MQ=17)

                                                                              tCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCg       >  2:1492321/1‑93 (MQ=17)

                                                                              tCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCg       >  2:2890814/1‑93 (MQ=17)

                                                                               cTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCgg      <  1:3277409/93‑1 (MQ=18)

                                                                                tGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGc     <  1:136358/93‑1 (MQ=18)

                                                                                 gAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCg    <  2:4116586/93‑1 (MQ=18)

                                                                                  aaTATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGa   <  2:1640162/93‑1 (MQ=18)

                                                                                   atatGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg  >  1:1636381/1‑93 (MQ=14)

                                                                                   atatGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAg  >  2:1607591/1‑93 (MQ=14)

GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAG  >  NC_000913/3941644‑3941819

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