BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,822,370	C→T	R77C (CGC→TGC)	spoU →	putative RNA methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,822,370	0	C	T	87.2%	921.6 / 81.6	327	R77C (CGC→TGC)	spoU	putative RNA methylase
Reads supporting (aligned to +/- strand): ref base C (18/24); new base T (96/189); total (114/213)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 327 total aligned reads displayed.

GAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGC  >  NC_000913/3822280‑3822459

gAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGc                                                                                         <  1:956511/93‑1 (MQ=255)

gAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGc                                                                                         <  1:6708970/93‑1 (MQ=255)

gAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGc                                                                                         <  1:3217100/93‑1 (MQ=255)

gAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGc                                                                                         >  1:19592/1‑93 (MQ=255)

gAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGc                                                                                         >  2:4419838/1‑93 (MQ=255)

   gTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCAcc                                                                                      >  2:5281155/1‑93 (MQ=255)

   gTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCAcc                                                                                      <  1:276903/93‑1 (MQ=255)

   gTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCAcc                                                                                      >  2:5963681/1‑93 (MQ=255)

    ttCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCa                                                                                     <  1:2925256/93‑1 (MQ=255)

     tCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCAt                                                                                    <  1:18599/93‑1 (MQ=255)

     tCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCAt                                                                                    >  1:2148747/1‑93 (MQ=255)

     nnaCGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGCAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAt                                                                                    <  2:3586856/91‑1 (MQ=255)

       aCGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTg                                                                                  >  1:1007737/1‑93 (MQ=255)

       aCGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGCACAGGTGAAAACACACTGCACCATTg                                                                                  >  2:1244409/1‑93 (MQ=255)

        cGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTgg                                                                                 <  1:6417530/93‑1 (MQ=255)

         gCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGTACAGGGGAAAACACACTGCACCATTGGc                                                                                >  1:1286747/1‑93 (MQ=255)

           cGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGa                                                                              <  1:6018834/93‑1 (MQ=255)

           cGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGa                                                                              <  2:3332708/93‑1 (MQ=255)

           cGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGa                                                                              <  1:7291361/93‑1 (MQ=255)

            gTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGAt                                                                             <  1:4315495/93‑1 (MQ=255)

            gTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGAt                                                                             >  1:4589406/1‑93 (MQ=255)

            gTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGGACAGGTGAAAACACACTGCACCATTGGCGAt                                                                             >  1:2384632/1‑93 (MQ=255)

             tCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATg                                                                            <  2:3483970/93‑1 (MQ=255)

               tGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACCGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGcc                                                                          >  2:436044/1‑93 (MQ=255)

               tGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGcc                                                                          <  1:1659749/93‑1 (MQ=255)

                   cTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCg                                                                      >  2:419571/1‑93 (MQ=255)

                   cTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCg                                                                      <  1:4905261/93‑1 (MQ=255)

                    tGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGGACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGc                                                                     >  1:1410529/1‑93 (MQ=255)

                     ggTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCt                                                                    <  2:2109328/93‑1 (MQ=255)

                       tAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCa                                                                  <  2:3160934/93‑1 (MQ=255)

                       tAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCa                                                                  <  2:4806650/93‑1 (MQ=255)

                       tAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCa                                                                  <  1:5949835/93‑1 (MQ=255)

                         gCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCAtc                                                                >  2:1668424/1‑93 (MQ=255)

                         gCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCAtc                                                                <  1:6896422/93‑1 (MQ=255)

                         gCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCAtc                                                                >  1:121238/1‑93 (MQ=255)

                          ccGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCAtct                                                               <  2:6321638/93‑1 (MQ=255)

                          ccGCATGCGCACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCAtct                                                               >  1:6327433/1‑93 (MQ=255)

                           cGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCAtctc                                                              <  2:2244468/93‑1 (MQ=255)

                             cATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCaa                                                            >  2:6249614/1‑93 (MQ=255)

                              aTGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCaaa                                                           <  2:7193848/93‑1 (MQ=255)

                              aTGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCaaa                                                           <  1:849683/93‑1 (MQ=255)

                              aTGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCaaa                                                           <  2:6963443/93‑1 (MQ=255)

                              aTGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCaaa                                                           <  2:6458551/93‑1 (MQ=255)

                              aTGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCaaa                                                           <  1:5666803/93‑1 (MQ=255)

                               tGCGCACCATGGCTTCGGTAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAg                                                          <  2:61577/93‑1 (MQ=255)

                               tGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAg                                                          <  1:5419192/93‑1 (MQ=255)

                                  gcACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                                                       >  1:3343965/1‑93 (MQ=255)

                                  gcACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                                                       >  2:2820182/1‑93 (MQ=255)

                                  gcACCATGGCTTCGGCAGCGGCGGGGAGTAACAGCTGGGGACAGGGGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                                                       >  1:6568124/1‑93 (MQ=255)

                                   cACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCa                                                      <  2:2665246/93‑1 (MQ=255)

                                   cACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCa                                                      >  2:5463481/1‑93 (MQ=255)

                                    aCCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAg                                                     >  2:450124/1‑93 (MQ=255)

                                     ccATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAgg                                                    <  1:45168/93‑1 (MQ=255)

                                      cATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAggg                                                   <  2:2354493/93‑1 (MQ=255)

                                      cATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAggg                                                   >  1:5133310/1‑93 (MQ=255)

                                       aTGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGc                                                  <  2:2015147/93‑1 (MQ=255)

                                        tGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCa                                                 <  2:841581/93‑1 (MQ=255)

                                         ggCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCAt                                                >  2:6002141/1‑93 (MQ=255)

                                           cTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGc                                              >  2:3759256/1‑93 (MQ=255)

                                           cTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGCGCATGc                                              <  2:4722135/93‑1 (MQ=255)

                                            ttCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCa                                             >  1:2289999/1‑93 (MQ=255)

                                            ttCGGCAGCGGCGGGTAGTAACAGCTGGGGACAGGTGAAAACACACTGCACCATTGGCGATGGCGTCGCTCATCTCAAAGGCCAGGGCATGCa                                             >  1:5107598/1‑93 (MQ=255)

                                               ggCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAt                                          <  1:2531832/93‑1 (MQ=255)

                                                gCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAtt                                         >  2:3634252/1‑93 (MQ=255)

                                                   gcggcgGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTg                                      <  1:4110262/93‑1 (MQ=255)

                                                   gcggcgGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTg                                      <  1:3958284/93‑1 (MQ=255)

                                                   gcggcgGGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTg                                      <  2:2074606/93‑1 (MQ=255)

                                                     ggcggGTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGc                                    <  2:1311532/93‑1 (MQ=255)

                                                        gggTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAAc                                 >  2:5334213/1‑93 (MQ=255)

                                                         ggTAGTAACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAAcc                                >  1:4626841/1‑93 (MQ=255)

                                                           tagtaACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCa                              <  2:1579506/93‑1 (MQ=255)

                                                             gtaACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATc                            <  2:5791286/93‑1 (MQ=255)

                                                             gtaACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATc                            <  2:3429066/93‑1 (MQ=255)

                                                             gtaACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATc                            <  1:4973267/93‑1 (MQ=255)

                                                              taACAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCt                           <  2:6629356/93‑1 (MQ=255)

                                                               aaCAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCtt                          >  1:5569433/1‑93 (MQ=255)

                                                                aCAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCttt                         <  2:2963804/93‑1 (MQ=255)

                                                                 cAGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTc                        <  1:2860740/93‑1 (MQ=255)

                                                                  aGCTGGGTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCt                       <  1:4064211/93‑1 (MQ=255)

                                                                     tGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGAt                    >  1:7447417/1‑93 (MQ=255)

                                                                      gggTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATa                   >  2:559535/1‑93 (MQ=255)

                                                                       ggTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATaa                  <  1:6155565/93‑1 (MQ=255)

                                                                       ggTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATaa                  <  2:7099801/93‑1 (MQ=255)

                                                                       ggTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATaa                  <  1:366198/93‑1 (MQ=255)

                                                                       ggTACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATaa                  <  2:2444995/93‑1 (MQ=255)

                                                                         tACAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACg                <  1:1494701/93‑1 (MQ=255)

                                                                          aCAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGc               <  1:5528131/93‑1 (MQ=255)

                                                                          aCAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGc               <  1:6245878/93‑1 (MQ=255)

                                                                          aCAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGc               <  2:1230770/93‑1 (MQ=255)

                                                                           cAGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCt              <  1:3435564/93‑1 (MQ=255)

                                                                            aGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTg             <  2:1906323/93‑1 (MQ=255)

                                                                            aGGTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTg             >  2:4114383/1‑93 (MQ=255)

                                                                             ggTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGt            <  2:1616940/93‑1 (MQ=255)

                                                                             ggTGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGt            <  1:7162194/93‑1 (MQ=255)

                                                                              gTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTc           <  2:3746516/93‑1 (MQ=255)

                                                                               tGAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCg          >  2:745562/1‑93 (MQ=255)

                                                                                gAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGa         <  2:1589708/93‑1 (MQ=255)

                                                                                gAAAACACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGa         >  2:4555861/1‑93 (MQ=255)

                                                                                 aaaaCACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGAt        <  1:6983294/93‑1 (MQ=255)

                                                                                 aaaaCACACTGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGAt        <  2:7345956/93‑1 (MQ=255)

GAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGC  >  NC_000913/3822280‑3822459

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