BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,689	C→T	D135D (GAC→GAT)	ytfI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,689	0	C	T	99.9%	2700.0 / ‑6.1	726	D135D (GAC→GAT)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (327/398); total (328/398)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 740 total aligned reads displayed.

TTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTG  >  NC_000913/4435601‑4435778

 tGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGAtt                                                                                       >  2:3117735/1‑93 (MQ=255)

  gAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTg                                                                                      >  1:1315008/1‑93 (MQ=255)

  gAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTg                                                                                      >  2:8662447/1‑93 (MQ=255)

  gAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTg                                                                                      >  1:8185823/1‑93 (MQ=255)

    aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTa                                                                                    >  1:1990264/1‑93 (MQ=255)

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

    aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTa                                                                                    <  1:7312818/93‑1 (MQ=255)

    aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTa                                                                                    <  2:8533430/93‑1 (MQ=255)

      aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATg                                                                                  <  2:2076598/93‑1 (MQ=255)

      aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATg                                                                                  >  1:6708999/1‑93 (MQ=255)

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

         aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCt                                                                               <  1:3846703/93‑1 (MQ=255)

          agaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTg                                                                              >  2:2660357/1‑93 (MQ=255)

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

                aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGatat                                                                        >  1:8787019/1‑93 (MQ=255)

                aCCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGatat                                                                        <  1:4723028/93‑1 (MQ=255)

                  ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGACTTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTc                                                                      >  1:1529658/1‑93 (MQ=255)

                   cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCt                                                                     >  2:2565844/1‑93 (MQ=255)

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

                      aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACa                                                                  >  1:9073827/1‑93 (MQ=255)

                       aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAg                                                                 >  1:8322807/1‑93 (MQ=255)

                       aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAg                                                                 >  1:7025225/1‑93 (MQ=255)

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

                        cAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGa                                                                <  1:2295917/93‑1 (MQ=255)

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

                        cAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGa                                                                <  1:5068530/93‑1 (MQ=255)

                         aaTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaa                                                               >  2:5713799/1‑93 (MQ=255)

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

                           tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaa                                                             >  1:6231396/1‑93 (MQ=255)

                            aaaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaaa                                                            <  1:5308660/93‑1 (MQ=255)

                            aaaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaaa                                                            >  1:3672541/1‑93 (MQ=255)

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

                            aaaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaaa                                                            >  1:7697799/1‑93 (MQ=255)

                            aaaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaaa                                                            >  2:8246021/1‑93 (MQ=255)

                             aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaaaa                                                           <  1:5759713/93‑1 (MQ=255)

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

                               tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                                         <  1:6117423/93‑1 (MQ=255)

                                 aaTATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGt                                                       <  1:7993430/93‑1 (MQ=255)

                                  atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                                      <  2:7036847/93‑1 (MQ=255)

                                   tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                                     >  2:6630503/1‑93 (MQ=255)

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

                                     taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                                                   <  2:6790993/93‑1 (MQ=255)

                                     taACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTATGa                                                   >  2:8408982/1‑93 (MQ=255)

                                     taACTCATTCATGGATTTCAAAAACCGTGAATCTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGa                                                   >  2:7373786/1‑93 (MQ=255)

                                       aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                                 >  1:7876216/1‑93 (MQ=255)

                                       aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                                 <  1:4045623/93‑1 (MQ=255)

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

                                       aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                                 <  2:946262/93‑1 (MQ=255)

                                       aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                                 <  2:4624055/93‑1 (MQ=255)

                                       aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                                 >  2:2835628/1‑93 (MQ=255)

                                         tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATa                                               <  1:9191092/93‑1 (MQ=255)

                                         tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATa                                               >  1:356313/1‑93 (MQ=255)

                                         tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATa                                               >  2:5560039/1‑93 (MQ=255)

                                            ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                            >  2:396141/1‑93 (MQ=255)

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

                                             tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa                                           >  1:6356906/1‑93 (MQ=255)

                                             tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa                                           >  1:4317185/1‑93 (MQ=255)

                                              catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa                                          >  1:9492636/1‑93 (MQ=255)

                                              catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa                                          >  1:3445269/1‑93 (MQ=255)

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

                                               atGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                                         <  1:4925425/93‑1 (MQ=255)

                                                tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                                        >  1:8689416/1‑93 (MQ=255)

                                                tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                                        <  1:6839002/93‑1 (MQ=255)

                                                tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                                        >  1:5468043/1‑93 (MQ=255)

                                                 ggATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg                                       <  1:6487002/93‑1 (MQ=255)

                                                  gATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg                                      <  2:5816529/93‑1 (MQ=255)

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

                                                   aTTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa                                     <  2:305308/93‑1 (MQ=255)

                                                    tttCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt                                    <  2:9656748/93‑1 (MQ=255)

                                                      tCAAAAACCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa                                  >  1:5948219/1‑93 (MQ=255)

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

                                                        aaaaaCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa                                <  2:3267190/93‑1 (MQ=255)

                                                         aaaaCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAg                               <  2:6057421/93‑1 (MQ=255)

                                                          aaaCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg                              <  2:6397901/93‑1 (MQ=255)

                                                          aaaCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg                              <  2:3373708/93‑1 (MQ=255)

                                                           acccGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                             <  2:1722689/91‑1 (MQ=255)

                                                            aCCGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGtt                            >  2:9082946/1‑93 (MQ=255)

                                                             ccGTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTta                           <  2:7819753/93‑1 (MQ=255)

                                                               gTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtata                         >  2:4778345/1‑93 (MQ=255)

                                                               gTGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtata                         >  2:7954538/1‑93 (MQ=255)

                                                                tGAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTatat                        <  2:3664345/93‑1 (MQ=255)

                                                                ggAATTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTatat                        <  2:7601019/92‑1 (MQ=255)

                                                                   aTTTCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCa                     <  1:439771/93‑1 (MQ=255)

                                                                     ttCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAcc                   <  1:1713982/93‑1 (MQ=255)

                                                                      tCCGTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa                  <  1:4439611/93‑1 (MQ=255)

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

                                                                         gTTAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGc               >  2:4431565/1‑93 (MQ=255)

                                                                          ttAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCg              >  2:5046229/1‑93 (MQ=255)

                                                                           tAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGa             >  1:8524895/1‑93 (MQ=255)

                                                                           tAGAACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGa             <  2:6850046/93‑1 (MQ=255)

                                                                                caaaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGt        <  2:9306525/91‑1 (MQ=255)

                                                                                  aaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGa      <  1:8891897/93‑1 (MQ=255)

                                                                                  aaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGa      <  2:2263248/93‑1 (MQ=255)

                                                                                  aaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGa      <  2:5245257/93‑1 (MQ=255)

                                                                                  aaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGa      <  1:120660/93‑1 (MQ=255)

                                                                                  aaacaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCCAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGa      <  1:2144572/93‑1 (MQ=255)

                                                                                    acaGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACg    <  2:8840155/93‑1 (MQ=255)

                                                                                     caGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACgt   >  2:4285548/1‑93 (MQ=255)

                                                                                     caGATGATTGTATGCCTGGGATAGTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACgt   <  2:7190878/93‑1 (MQ=255)

                                                                                      aGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACgtg  >  1:3027751/1‑93 (MQ=255)

TTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGATTGGTGACGTG  >  NC_000913/4435601‑4435778

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