BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,493	T→C	C124C (TGT→TGC)	yidX →	putative replicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,493	0	T	C	98.6%	952.6 / ‑1.2	293	C124C (TGT→TGC)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base T (3/1); new base C (135/154); total (138/155)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 297 total aligned reads displayed.

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTA  >  NC_000913/3873403‑3873584

cGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCAc                                                                                           <  2:3112598/93‑1 (MQ=255)

 ggCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCAcc                                                                                          >  2:2561436/1‑93 (MQ=255)

 ggCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCAcc                                                                                          <  2:2178634/93‑1 (MQ=255)

 ggCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCAcc                                                                                          >  1:4430275/1‑93 (MQ=255)

  gCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCa                                                                                         <  2:3726015/93‑1 (MQ=255)

    acaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAATACTGCACCAgg                                                                                       >  1:4645147/1‑93 (MQ=255)

    acaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAgg                                                                                       <  1:4292965/93‑1 (MQ=255)

    acaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAgg                                                                                       >  1:3585926/1‑93 (MQ=255)

     caATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGt                                                                                      >  2:3028655/1‑93 (MQ=255)

      aaTATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTc                                                                                     <  1:4277186/93‑1 (MQ=255)

      aaTATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTc                                                                                     >  1:3781121/1‑93 (MQ=255)

           cGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTgg                                                                                <  2:3620949/93‑1 (MQ=255)

            gCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGa                                                                               >  2:866096/1‑93 (MQ=255)

             cTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaa                                                                              >  2:120642/1‑93 (MQ=255)

              tttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaa                                                                             <  1:3095495/93‑1 (MQ=255)

              tttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaa                                                                             <  1:288224/93‑1 (MQ=255)

               ttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaaa                                                                            <  2:487320/93‑1 (MQ=255)

                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                        >  1:2026568/1‑93 (MQ=255)

                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                       >  1:621551/1‑93 (MQ=255)

                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                       >  1:339368/1‑93 (MQ=255)

                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTaaaa                                                                     >  2:4346840/1‑93 (MQ=255)

                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                     >  2:2375588/1‑93 (MQ=255)

                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                     <  2:3321109/93‑1 (MQ=255)

                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                     >  1:421935/1‑93 (MQ=255)

                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                     <  2:3542872/93‑1 (MQ=255)

                       actGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAg                                                                    <  2:2325645/93‑1 (MQ=255)

                       actGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAg                                                                    >  2:3448086/1‑93 (MQ=255)

                       actGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGCTTTGACTGTGGTGTTCAACACTGCACCAAGTCACTGGAAAACCTAAAAg                                                                    >  1:552155/1‑93 (MQ=255)

                          ggAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAACACCTAACAGGca                                                                 >  2:3336654/1‑93 (MQ=255)

                          ggAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGca                                                                 >  1:2725453/1‑93 (MQ=255)

                          ggAAGAACGAAGAAGAAAAGGCACTGCCAACGCAGTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGca                                                                 <  1:2297659/93‑1 (MQ=255)

                           gAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGcac                                                                <  2:1353413/93‑1 (MQ=255)

                           gAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGcaa                                                                >  2:1555033/1‑92 (MQ=255)

                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                              >  2:1663572/1‑93 (MQ=255)

                              tAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAAt                                                             <  1:2334993/93‑1 (MQ=255)

                              tAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAAt                                                             >  2:745274/1‑93 (MQ=255)

                               aaCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATc                                                            >  1:4010830/1‑93 (MQ=255)

                                aCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATcc                                                           <  2:1326140/93‑1 (MQ=255)

                                aCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATcc                                                           <  1:1246163/93‑1 (MQ=255)

                                aCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATcc                                                           <  1:3199069/93‑1 (MQ=255)

                                 cGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCa                                                          <  2:555130/93‑1 (MQ=255)

                                 cGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCa                                                          <  1:875780/93‑1 (MQ=255)

                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAt                                                         >  2:3501816/1‑93 (MQ=255)

                                   aagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCATGTCACTGGAAAACCTAAAAGGCACAATCCATg                                                        >  1:3815185/1‑93 (MQ=255)

                                   aagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCGCTGGAAAACCTAAAAGGCACAATCCATg                                                        <  2:3710766/93‑1 (MQ=255)

                                   aagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATg                                                        >  1:2978739/1‑93 (MQ=255)

                                    agaagaaAAGGCACTGGCAACGCATTATGGTTGTTCCTGTGGTGTTCAACACTGCACCAGTTCACTGGAAAACCTAAAAGGCACAATCCATga                                                       <  2:2020177/93‑1 (MQ=255)

                                       agaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaa                                                    >  2:1175102/1‑93 (MQ=255)

                                          aaaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaaa                                                 <  1:3884201/93‑1 (MQ=255)

                                          aaaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaaa                                                 <  1:4120372/93‑1 (MQ=255)

                                            aGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATa                                               >  1:967459/1‑93 (MQ=255)

                                             ggCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaa                                              >  1:1194036/1‑93 (MQ=255)

                                             ggCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaa                                              <  2:2421958/93‑1 (MQ=255)

                                             ggCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaa                                              >  1:2802430/1‑93 (MQ=255)

                                                aCTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaaa                                           <  1:4532293/93‑1 (MQ=255)

                                                aCTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAtaaaaataaaaa                                           <  1:2634155/93‑1 (MQ=255)

                                                aCTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaaa                                           >  1:3155014/1‑93 (MQ=255)

                                                 cTGGCAACGCCTTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                          >  1:723236/1‑93 (MQ=255)

                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                          <  1:3459268/93‑1 (MQ=255)

                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                          >  1:3003455/1‑93 (MQ=255)

                                                  tGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACa                                         <  1:1341318/93‑1 (MQ=255)

                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                      <  2:2818775/93‑1 (MQ=255)

                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                      <  2:4741516/93‑1 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                     <  2:169117/93‑1 (MQ=255)

                                                      aaCGCATTAAGGTTTTGCTTGTGGAGTTCAACACTGCACCAGGTCACTGGAAAGCCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                     <  1:2236373/93‑1 (MQ=255)

                                                       aCGCATTATTGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAt                                    <  1:1493464/93‑1 (MQ=255)

                                                       aCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAt                                    <  1:2564856/93‑1 (MQ=255)

                                                       aCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAt                                    <  1:2491624/93‑1 (MQ=255)

                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                   <  2:2683167/93‑1 (MQ=255)

                                                         gCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTa                                  >  1:4603441/1‑93 (MQ=255)

                                                          cATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTAc                                 <  2:426826/93‑1 (MQ=255)

                                                           aTTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACt                                >  1:181594/1‑93 (MQ=255)

                                                                ggTTTTGCCTGTGGTGTTCAACATTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAg                           >  2:249783/1‑93 (MQ=255)

                                                                ggTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAg                           <  2:115532/93‑1 (MQ=255)

                                                                ggTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAg                           <  2:320025/93‑1 (MQ=255)

                                                                ggTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAg                           <  2:4681092/93‑1 (MQ=255)

                                                                 gTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTAGTCAAAgc                          <  2:2883286/93‑2 (MQ=255)

                                                                 gTTGTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAgg                          <  2:4433919/93‑1 (MQ=255)

                                                                  ttttGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAggg                         >  2:4028646/1‑92 (MQ=255)

                                                                  ttttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGt                         <  2:1237559/93‑1 (MQ=255)

                                                                  ttttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGt                         <  2:337393/93‑1 (MQ=255)

                                                                   tttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTg                        <  2:3069843/93‑1 (MQ=255)

                                                                    ttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGa                       >  2:2407508/1‑93 (MQ=255)

                                                                      gCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATg                     >  2:413277/1‑93 (MQ=255)

                                                                       ccTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATgg                    <  1:1029039/93‑1 (MQ=255)

                                                                       ccTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATgg                    <  2:4539444/93‑1 (MQ=255)

                                                                        cTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGc                   >  1:4730327/1‑93 (MQ=255)

                                                                         tgtgGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCg                  >  2:689208/1‑93 (MQ=255)

                                                                         tgtgGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCg                  >  1:1047227/1‑93 (MQ=255)

                                                                            ggtgTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTc               <  2:2147200/93‑1 (MQ=255)

                                                                            ggtgTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTc               <  2:185778/93‑1 (MQ=255)

                                                                             gtgtTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCt              <  1:100594/93‑1 (MQ=255)

                                                                               gtTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAc            <  2:4558860/93‑1 (MQ=255)

                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc        <  2:1091866/93‑1 (MQ=255)

                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc        <  1:1614723/93‑1 (MQ=255)

                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca      >  1:4334952/1‑93 (MQ=255)

                                                                                      acTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatccat     >  1:3504405/1‑93 (MQ=255)

                                                                                       cTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAtt    <  1:3742831/93‑1 (MQ=255)

                                                                                       cTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAtt    >  1:1143272/1‑93 (MQ=255)

                                                                                         gTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTa  >  1:1734229/1‑93 (MQ=255)

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTA  >  NC_000913/3873403‑3873584

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