BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,684	+A	coding (400/501 nt)	ytfI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,681	1	.	A	99.3%	2522.1 / 4.6	678	coding (397/501 nt)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (4/1); new base A (303/369); total (308/370)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 746 total aligned reads displayed.

GAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑A‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAT  >  NC_000913/4435589‑4435768

 aGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a                                                                                         >  1:3414108/1‑93 (MQ=255)

 aGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCAGTacaac‑a                                                                                         >  1:5825255/1‑93 (MQ=255)

   gTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aa                                                                                      >  2:8020151/1‑93 (MQ=255)

   gTAGTGGATTTGAATACGAAGAAATACCCATAAAAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAAaac‑a‑aa                                                                                      >  1:8475178/1‑93 (MQ=255)

     aGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaac                                                                                    >  1:6247862/1‑91 (MQ=255)

      gTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaaca                                                                                   >  2:5514588/1‑90 (MQ=255)

      gTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaaca                                                                                   <  2:3219439/93‑4 (MQ=255)

       tGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaacag                                                                                  >  2:3717774/1‑89 (MQ=255)

       tGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a‑aaacag                                                                                  >  1:799044/1‑89 (MQ=255)

       tGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACGA‑AACAg                                                                                   >  2:2223375/1‑93 (MQ=255)

            ttGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAgatgat                                                                              >  1:6494158/1‑93 (MQ=255)

            ttGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAgatgat                                                                              >  2:5775605/1‑93 (MQ=255)

             tGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGAtt                                                                             >  1:3543060/1‑93 (MQ=255)

              gAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTg                                                                            >  2:8662448/1‑93 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTa                                                                          <  1:2491214/93‑1 (MQ=255)

                 tACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTAt                                                                         <  2:3055698/93‑1 (MQ=255)

                 tACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTAt                                                                         <  2:755696/93‑1 (MQ=255)

                 tACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTAt                                                                         <  1:4454267/93‑1 (MQ=255)

                    gaagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGcc                                                                      <  2:6918840/93‑1 (MQ=255)

                     aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCt                                                                     >  1:10008187/1‑93 (MQ=255)

                         aaTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGa                                                                 >  1:2126834/1‑93 (MQ=255)

                           tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGata                                                               >  1:2993086/1‑93 (MQ=255)

                           tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGata                                                               >  2:6709708/1‑93 (MQ=255)

                           tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGata                                                               <  2:8602186/93‑1 (MQ=255)

                            aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGatat                                                              <  2:1532709/93‑1 (MQ=255)

                              ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTc                                                            >  1:6891989/1‑93 (MQ=255)

                                 tAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTAc                                                         >  1:974810/1‑93 (MQ=255)

                                   aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAg                                                       >  1:9635553/1‑93 (MQ=255)

                                   aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAg                                                       >  1:8706819/1‑93 (MQ=255)

                                    cAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGa                                                      >  1:9975739/1‑93 (MQ=255)

                                     aaTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGaa                                                     <  2:8474982/93‑1 (MQ=255)

                                       tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGaaaa                                                   <  1:7793932/93‑1 (MQ=255)

                                       tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGaaaa                                                   >  2:5388638/1‑93 (MQ=255)

                                        aaaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGaaaaa                                                  >  1:7491223/1‑93 (MQ=255)

                                         aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGaaaaaa                                                 <  2:448490/93‑1 (MQ=255)

                                         aaTGAATATAACTCATTCATGGATTTCAAAAACCGGGAATTTCCGTTAGAAC‑AGAACAGATGATTGTATGCCTGGGATATTCTACAGaaacac                                                 >  1:6755636/1‑90 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                <  2:7684482/93‑1 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                <  1:5638887/93‑1 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                >  1:8093692/1‑93 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                <  2:2583634/93‑1 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                >  2:569583/1‑93 (MQ=255)

                                           tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               <  2:95951/93‑1 (MQ=255)

                                           tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               >  1:950336/1‑93 (MQ=255)

                                             aaTATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGt                                             >  2:7134639/1‑93 (MQ=255)

                                             aaTATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGt                                             >  2:4094620/1‑93 (MQ=255)

                                             aaTATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGt                                             >  2:1710605/1‑93 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  2:5973096/93‑1 (MQ=255)

                                               tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                           >  1:3791446/1‑93 (MQ=255)

                                               tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                           <  1:2753362/93‑1 (MQ=255)

                                               tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                           <  1:5273626/93‑1 (MQ=255)

                                               tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                           <  1:8872209/93‑1 (MQ=255)

                                                ataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTg                                          <  2:783163/93‑1 (MQ=255)

                                                 taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                                         >  1:1215024/1‑93 (MQ=255)

                                                 taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                                         >  1:19887/1‑93 (MQ=255)

                                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                                        <  1:5965/93‑1 (MQ=255)

                                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                                        <  2:4742153/93‑1 (MQ=255)

                                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                                        <  2:8994607/93‑1 (MQ=255)

                                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                                        >  1:4284428/1‑93 (MQ=255)

                                                   aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                       <  1:3204752/93‑1 (MQ=255)

                                                    cTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                                      <  1:9147126/93‑1 (MQ=255)

                                                    cTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                                      <  1:4011730/93‑1 (MQ=255)

                                                     tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATa                                     <  1:4383689/93‑1 (MQ=255)

                                                      cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAc                                    <  1:2301614/93‑1 (MQ=255)

                                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAcc                                   <  1:9030053/93‑1 (MQ=255)

                                                        ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                  >  2:1193572/1‑93 (MQ=255)

                                                        ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                  >  1:4930423/1‑93 (MQ=255)

                                                        ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                  >  2:8269599/1‑93 (MQ=255)

                                                         tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATCTTCTACCGAAAAAACAGTCTGACATACCGa                                 >  2:6347800/1‑93 (MQ=255)

                                                         tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa                                 >  1:6356907/1‑93 (MQ=255)

                                                           atGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                               <  1:9355193/93‑1 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                              >  1:1375110/1‑93 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                              >  1:5067100/1‑93 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                              >  2:5253426/1‑93 (MQ=255)

                                                              gATTTCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg                            >  2:9820126/1‑93 (MQ=255)

                                                                 ttCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt                         >  1:857444/1‑93 (MQ=255)

                                                                  tCAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa                        <  1:769171/93‑1 (MQ=255)

                                                                   cAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaa                       <  2:483289/93‑1 (MQ=255)

                                                                   cAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaa                       <  2:2499243/93‑1 (MQ=255)

                                                                   cAAAAACCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaa                       <  2:7290494/93‑1 (MQ=255)

                                                                    aaaaaCCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa                      <  1:5378125/93‑1 (MQ=255)

                                                                      aaaCCGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg                    <  1:6144085/93‑1 (MQ=255)

                                                                          cGTGAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtat                <  2:4410270/93‑1 (MQ=255)

                                                                             gAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATAtt             <  1:8267733/93‑1 (MQ=255)

                                                                             gAATTTCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATAtt             <  2:3335696/93‑1 (MQ=255)

                                                                                tttCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc          <  2:1996442/93‑1 (MQ=255)

                                                                                tttCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc          >  2:9428105/1‑93 (MQ=255)

                                                                                 ttCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAcc         <  1:7278429/93‑1 (MQ=255)

                                                                                 ttCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAcc         >  2:2908799/1‑93 (MQ=255)

                                                                                  tCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa        <  2:3610153/93‑1 (MQ=255)

                                                                                  tCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa        <  2:6810380/93‑1 (MQ=255)

                                                                                  tCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa        <  2:2719440/93‑1 (MQ=255)

                                                                                  tCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa        >  1:4629774/1‑93 (MQ=255)

                                                                                  tCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa        <  1:4754342/93‑1 (MQ=255)

                                                                                  tCCGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa        <  2:5064579/93‑1 (MQ=255)

                                                                                   ccGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAg       <  1:1559690/93‑1 (MQ=255)

                                                                                   ccGTTAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAg       <  2:6492198/93‑1 (MQ=255)

                                                                                     gTTAGAGCAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGc     <  1:7052461/93‑1 (MQ=255)

                                                                                      ttAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCg    >  1:1802231/1‑93 (MQ=255)

                                                                                       tAGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGa   <  1:3664997/93‑1 (MQ=255)

                                                                                        aGAACAA‑AACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAt  >  2:9158404/1‑93 (MQ=255)

GAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑A‑AACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAT  >  NC_000913/4435589‑4435768

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