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	99.8%	1609.2 / ‑5.6	502	C124C (TGT→TGC)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (234/266); total (235/267)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 507 total aligned reads displayed.

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATT  >  NC_000913/3873403‑3873582

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

 ggCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCTTTATGGTTTTGCCTGTGGTGTTCAACACTGCAcc                                                                                        <  1:3035214/93‑1 (MQ=255)

  gCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCa                                                                                       >  1:5538222/1‑93 (MQ=255)

     caATATTGGTTTACTACTACTGGAATAACGAAGCAGAAAACGGCCTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGt                                                                                    <  1:3844357/93‑1 (MQ=255)

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

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

       atatCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCa                                                                                  >  2:134626/1‑93 (MQ=255)

       atatCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCa                                                                                  <  1:2291162/93‑1 (MQ=255)

        tatCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCAc                                                                                 >  1:4380158/1‑93 (MQ=255)

        tatCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCAc                                                                                 <  1:1871549/93‑1 (MQ=255)

          tCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTg                                                                               <  2:4305138/93‑1 (MQ=255)

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

                tactactacTGGAAGAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAAc                                                                         <  1:3736124/93‑1 (MQ=255)

                  ctactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCt                                                                       <  1:2617603/93‑1 (MQ=255)

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

                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                     <  2:4466406/93‑1 (MQ=255)

                     ctactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaa                                                                    <  2:7931227/93‑1 (MQ=255)

                      tgctGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                   <  1:1982863/91‑1 (MQ=255)

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

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

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

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

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

                        ctGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAgg                                                                 >  2:3941799/1‑93 (MQ=255)

                        ctGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAgg                                                                 <  2:3443186/93‑1 (MQ=255)

                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                                <  2:5902349/93‑1 (MQ=255)

                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                                <  1:2426226/93‑1 (MQ=255)

                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                                >  1:5751014/1‑93 (MQ=255)

                          ggAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGca                                                               <  1:160436/93‑1 (MQ=255)

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

                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                            <  2:4033523/93‑1 (MQ=255)

                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                            <  2:7860651/93‑1 (MQ=255)

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

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

                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAt                                                       <  1:860278/93‑1 (MQ=255)

                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAt                                                       <  2:6419006/93‑1 (MQ=255)

                                    agaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATga                                                     >  1:4152499/1‑93 (MQ=255)

                                    agaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATga                                                     <  1:5466754/93‑1 (MQ=255)

                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                    >  2:5220884/1‑93 (MQ=255)

                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                    >  1:7869357/1‑93 (MQ=255)

                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                    >  1:2813919/1‑93 (MQ=255)

                                      aagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgaga                                                   >  2:1065326/1‑93 (MQ=255)

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

                                        gaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaa                                                 >  2:4719937/1‑93 (MQ=255)

                                        gaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaa                                                 >  1:6149110/1‑93 (MQ=255)

                                        gaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaa                                                 <  1:659942/93‑1 (MQ=255)

                                        aaaaaGGCACTGGCAACGGATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaa                                                 <  1:3809805/92‑1 (MQ=255)

                                         aaaaGGCACTGGCAACGCATTATGGTTTTGCTTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaa                                                >  2:6812205/1‑93 (MQ=255)

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

                                          aaaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaaa                                               >  2:5549477/1‑93 (MQ=255)

                                               cACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaa                                          >  1:2720823/1‑93 (MQ=255)

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

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

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

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

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

                                                   ggCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAACCCTAAAAGGCACAATCCATGAGAAAAATAAAAACAt                                      <  1:5080130/93‑1 (MQ=255)

                                                   ggCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACAt                                      <  2:2352335/93‑1 (MQ=255)

                                                    gCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATg                                     <  1:168994/93‑1 (MQ=255)

                                                     cAACACATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                    <  2:2947348/93‑1 (MQ=255)

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

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   >  2:3542372/1‑93 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   >  1:6354993/1‑93 (MQ=255)

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

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

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                   >  1:4517440/1‑93 (MQ=255)

                                                       aCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAt                                  >  2:7305694/1‑93 (MQ=255)

                                                       aCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAt                                  >  1:5298634/1‑93 (MQ=255)

                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                 >  2:2220331/1‑93 (MQ=255)

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

                                                           aTTAAGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACt                              >  1:1280199/1‑93 (MQ=255)

                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                             >  1:5899675/1‑93 (MQ=255)

                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                             >  1:5231408/1‑93 (MQ=255)

                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                             >  1:1107856/1‑93 (MQ=255)

                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                             >  2:2112593/1‑93 (MQ=255)

                                                             tATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCa                            <  2:1247177/93‑1 (MQ=255)

                                                               tGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCGCAATCCATGAGAAAAATAAAAACATGGATTACTCaaa                          <  2:706202/93‑1 (MQ=255)

                                                               tGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaaa                          <  2:3277426/93‑1 (MQ=255)

                                                               tGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaaa                          >  1:7117304/1‑93 (MQ=255)

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

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

                                                                 gTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAgg                        <  2:5113493/93‑1 (MQ=255)

                                                                 gTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAgg                        <  2:7698025/93‑1 (MQ=255)

                                                                  ttttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGt                       >  1:6543474/1‑93 (MQ=255)

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

                                                                   tttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTg                      >  1:6970956/1‑93 (MQ=255)

                                                                    ttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGa                     <  2:5011586/93‑1 (MQ=255)

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

                                                                          gtggtgTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGt               <  1:2199881/93‑1 (MQ=255)

                                                                             gtgtTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATGAAAACATGGATTACTCAAAGGTGATGGCGTTCt            >  1:6232099/1‑93 (MQ=255)

                                                                             gtgtTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCt            >  1:3211170/1‑93 (MQ=255)

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

                                                                                ttCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAcc         <  1:467806/93‑1 (MQ=255)

                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACACGGATTACTCAAAGGTGATGGCGTTCTAccatc      <  1:1407748/93‑1 (MQ=255)

                                                                                    acacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcc     >  2:197164/1‑93 (MQ=255)

                                                                                    acacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcc     <  2:3202809/93‑1 (MQ=255)

                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca    <  1:3411882/93‑1 (MQ=255)

                                                                                      acTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatccat   <  2:633783/93‑1 (MQ=255)

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

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

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATT  >  NC_000913/3873403‑3873582

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: `‑`