BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,211,972	C→G	T289S (ACC→AGC)	dnaG →	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,211,972	0	C	G	99.1%	406.5 / ‑4.0	115	T289S (ACC→AGC)	dnaG	DNA primase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base G (64/50); total (65/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 122 total aligned reads displayed.

GAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAACCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTGGCGAG  >  NC_000913/3211875‑3212067

gAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAagca                                                                                               >  1:89887/1‑100 (MQ=255)

  aCCCAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAagcacc                                                                                             >  1:737514/1‑100 (MQ=255)

  aCCCAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAagcacc                                                                                             >  1:536443/1‑100 (MQ=255)

  aCCCAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAagcacc                                                                                             <  2:2982064/100‑1 (MQ=255)

  aCCCAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAagcacc                                                                                             <  2:2427201/100‑1 (MQ=255)

  aCCCAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAagcac                                                                                              >  2:408354/1‑99 (MQ=255)

     cAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCAccgcc                                                                                          >  1:1328960/1‑100 (MQ=255)

     cAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCAccgcc                                                                                          <  2:1328947/100‑1 (MQ=255)

     cAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCAccgcc                                                                                          >  2:1454408/1‑100 (MQ=255)

         cGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATc                                                                                      <  2:416771/100‑1 (MQ=255)

         cGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATc                                                                                      <  1:1631148/100‑1 (MQ=255)

          gTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATca                                                                                     <  1:74157/100‑1 (MQ=255)

             tGCTTGTGGTTGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACAt                                                                                  <  1:2537178/100‑1 (MQ=255)

             tGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACAt                                                                                  >  1:396730/1‑100 (MQ=255)

                ttGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACa                                                                               <  1:1204765/100‑1 (MQ=255)

                 tgtgGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTCGCGTCGATAGGTACGTCAAGCACCGCCGATCACATACaa                                                                              >  2:2128243/1‑100 (MQ=255)

                   tgGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACt                                                                            >  1:326083/1‑100 (MQ=255)

                    ggTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtg                                                                           <  1:71726/100‑1 (MQ=255)

                    ggTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtg                                                                           >  1:1499491/1‑100 (MQ=255)

                    ggTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtg                                                                           <  1:1254059/100‑1 (MQ=255)

                    ggTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtg                                                                           >  2:1059214/1‑100 (MQ=255)

                     gTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtgt                                                                          >  1:841210/1‑100 (MQ=255)

                     gTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtg                                                                           >  2:1405694/1‑99 (MQ=255)

                      tCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtgtt                                                                         >  2:1892890/1‑100 (MQ=255)

                      tCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtgtt                                                                         <  1:351351/100‑1 (MQ=255)

                       cGAAGGCTATATGGACGTGGTGGCGCTGGCGCGATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtgttg                                                                        >  1:2876241/1‑100 (MQ=255)

                       cGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtgttg                                                                        >  1:2105954/1‑100 (MQ=255)

                       cGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtgttg                                                                        >  1:2102148/1‑100 (MQ=255)

                       cGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACtgttg                                                                        <  1:611500/100‑1 (MQ=255)

                          aGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACCACtgttgt                                                                       >  2:1202890/1‑98 (MQ=255)

                           ggCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTcc                                                                    <  2:2102130/100‑1 (MQ=255)

                            gCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCcg                                                                   <  2:339844/100‑1 (MQ=255)

                              tataTGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCcgcg                                                                 >  2:428745/1‑100 (MQ=255)

                               atatGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCcgcg                                                                 >  2:2955019/1‑99 (MQ=255)

                               atatGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACCGTTGTTCcgcgc                                                                >  2:1245041/1‑100 (MQ=255)

                                tatGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCcgcgcg                                                               >  1:453283/1‑100 (MQ=255)

                                tatGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCcgcgcg                                                               <  2:1159040/100‑1 (MQ=255)

                                  tGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGAc                                                             <  1:1810266/100‑1 (MQ=255)

                                   ggACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGAcc                                                            >  2:2006661/1‑100 (MQ=255)

                                   ggACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGAcc                                                            >  2:2438226/1‑100 (MQ=255)

                                   ggACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGAc                                                             >  1:1161876/1‑99 (MQ=255)

                                       gtggtgGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACcaaca                                                        <  2:295796/100‑1 (MQ=255)

                                        tggtggCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACcaacaa                                                       >  2:3019962/1‑100 (MQ=255)

                                        tggtggCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACcaacaa                                                       >  2:154349/1‑100 (MQ=255)

                                        tggtggCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACcaacaa                                                       <  1:365795/100‑1 (MQ=255)

                                        tggtggCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACcaacaa                                                       <  2:829059/100‑1 (MQ=255)

                                         ggtggCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCATTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAAt                                                      >  2:1050128/1‑100 (MQ=255)

                                           tggcgctggcgcAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGt                                                    <  1:1322276/100‑1 (MQ=255)

                                              cgctggcgcAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCAt                                                 >  2:1842536/1‑100 (MQ=255)

                                               gctggcgcAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCAtt                                                >  2:1559725/1‑100 (MQ=255)

                                                  ggcgcAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGc                                             >  2:453440/1‑100 (MQ=255)

                                                    cgcAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTg                                           <  1:2438248/100‑1 (MQ=255)

                                                    cgcAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTg                                           <  1:876586/100‑1 (MQ=255)

                                                      cAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGtt                                         >  1:2608841/1‑100 (MQ=255)

                                                       aaTACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTa                                        <  1:2491194/100‑1 (MQ=255)

                                                       aaTACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTa                                        <  1:1487888/100‑1 (MQ=255)

                                                             gCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACgg                                  >  2:1026623/1‑100 (MQ=255)

                                                              cATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGc                                 >  2:2028130/1‑100 (MQ=255)

                                                              cATTAATTACGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGc                                 <  2:560098/100‑1 (MQ=255)

                                                               attaattaCGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCg                                <  1:99626/100‑1 (MQ=255)

                                                               attaattaCGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCg                                <  2:640891/100‑1 (MQ=255)

                                                                ttaattaCGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGa                               <  2:89886/100‑1 (MQ=255)

                                                                 taattaCGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGAc                              >  1:139741/1‑100 (MQ=255)

                                                                  tcttccgatcttGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGAcc                             <  2:2833959/90‑1 (MQ=255)

                                                                  aattaCGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGAcc                             >  2:2627779/1‑100 (MQ=255)

                                                                  aattaCGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGAcc                             <  1:248229/100‑1 (MQ=255)

                                                                   attaCGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCg                            >  2:2369094/1‑100 (MQ=255)

                                                                       cGCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCa                        >  2:1700187/1‑100 (MQ=255)

                                                                       cGCCGTTGCGTCGTTAGGTACGTCAACCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCa                        >  1:2365710/1‑100 (MQ=255)

                                                                        gCCGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAg                       >  1:2843116/1‑100 (MQ=255)

                                                                         ccGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAgg                      >  1:1765367/1‑100 (MQ=255)

                                                                         ccGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAgg                      >  1:2773127/1‑100 (MQ=255)

                                                                         ccGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAgg                      >  2:259345/1‑100 (MQ=255)

                                                                          cGTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGc                     <  2:53318/100‑1 (MQ=255)

                                                                           gTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGcc                    >  2:2213846/1‑100 (MQ=255)

                                                                           gTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGcc                    <  1:1542928/100‑1 (MQ=255)

                                                                           gTTGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGc                     <  2:52239/99‑1 (MQ=255)

                                                                            ttGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCcg                   >  1:2148659/1‑100 (MQ=255)

                                                                            ttGCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCcg                   >  1:1085443/1‑100 (MQ=255)

                                                                              gCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCcgcg                 <  2:1016/100‑1 (MQ=255)

                                                                              gCGTCGTTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCcgcg                 <  2:2725168/100‑1 (MQ=255)

                                                                               cgtcgtTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACCATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGa                >  2:753950/1‑100 (MQ=255)

                                                                                gtcgtTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACTATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGAt               <  1:1180146/100‑1 (MQ=255)

                                                                                  cgtTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATgc             >  2:1031659/1‑100 (MQ=255)

                                                                                  cgtTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATgc             >  1:732857/1‑100 (MQ=255)

                                                                                   gtTAGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATgcc            >  1:552030/1‑100 (MQ=255)

                                                                                      aGGTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATgccgcc         >  1:609564/1‑100 (MQ=255)

                                                                                       ggTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCt        >  1:1930435/1‑100 (MQ=255)

                                                                                       ggTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCt        >  1:298724/1‑100 (MQ=255)

                                                                                        gTACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTg       <  2:799742/100‑1 (MQ=255)

                                                                                         tACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTgg      <  1:1026633/100‑1 (MQ=255)

                                                                                         tACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTgg      >  1:1898278/1‑100 (MQ=255)

                                                                                         tACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTgg      >  2:2157839/1‑100 (MQ=255)

                                                                                         tACGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTgg      >  2:2637625/1‑100 (MQ=255)

                                                                                          aCGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACTGTGCAGGCCGCGATGCCGCCTGGc     <  2:2851000/100‑1 (MQ=255)

                                                                                           cGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTGGCg    <  1:901856/100‑1 (MQ=255)

                                                                                           cGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTGGCg    <  1:1006492/100‑1 (MQ=255)

                                                                                           cGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTGGCg    >  1:1862560/1‑100 (MQ=255)

                                                                                           cGTCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTGGCg    <  2:2818616/100‑1 (MQ=255)

                                                                                             tCAAGCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTGGCGAg  >  2:420879/1‑100 (MQ=255)

GAACCCAATCGTCTGCTTGTGGTCGAAGGCTATATGGACGTGGTGGCGCTGGCGCAATACGGCATTAATTACGCCGTTGCGTCGTTAGGTACGTCAACCACCGCCGATCACATACAACTGTTGTTCCGCGCGACCAACAATGTCATTTGCTGTTATGACGGCGACCGTGCAGGCCGCGATGCCGCCTGGCGAG  >  NC_000913/3211875‑3212067

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