BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,209	A→G	T30A (ACT→GCT)	yidX →	putative replicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,209	0	A	G	85.0%	716.6 / 114.8	306	T30A (ACT→GCT)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base A (25/21); new base G (146/114); total (171/135)
Fisher's exact test for biased strand distribution p-value = 8.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 307 total aligned reads displayed.

GCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGAT  >  NC_000913/3873118‑3873301

gcTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCAc                                                                                             >  1:4851972/1‑93 (MQ=255)

 cTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCt                                                                                            >  1:1687501/1‑93 (MQ=255)

  tGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTg                                                                                           >  2:6082993/1‑93 (MQ=255)

   gATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTgg                                                                                          >  2:4710566/1‑93 (MQ=255)

   gATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCttg                                                                                          >  2:6566734/1‑91 (MQ=255)

   gATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTgg                                                                                          <  1:6760912/93‑1 (MQ=255)

   gATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTgg                                                                                          <  1:1924160/93‑1 (MQ=255)

   gATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTgg                                                                                          >  2:1476618/1‑93 (MQ=255)

    aTGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTggg                                                                                         <  2:4647907/93‑1 (MQ=255)

    aTGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTggg                                                                                         >  1:4078403/1‑93 (MQ=255)

    aTGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTggg                                                                                         >  2:4043179/1‑93 (MQ=255)

      gCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCt                                                                                       <  1:6476224/93‑1 (MQ=255)

          tCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCg                                                                                   >  2:5336406/1‑93 (MQ=255)

           cAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGc                                                                                  <  2:6026560/93‑1 (MQ=255)

            aGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCa                                                                                 <  1:1503688/93‑1 (MQ=255)

            aGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCa                                                                                 <  2:1549614/93‑1 (MQ=255)

             ggCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAg                                                                                <  2:4322399/93‑1 (MQ=255)

             ggCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAg                                                                                <  2:4994/93‑1 (MQ=255)

             ggCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAg                                                                                <  1:3314160/93‑1 (MQ=255)

              gCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAgg                                                                               >  2:26950/1‑93 (MQ=255)

                 gTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCaa                                                                            <  2:3833173/93‑1 (MQ=255)

                  tATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCaaa                                                                           >  1:4711167/1‑93 (MQ=255)

                  tATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCaaa                                                                           >  2:7018787/1‑93 (MQ=255)

                    tctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAg                                                                         <  1:3044273/93‑1 (MQ=255)

                     ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGa                                                                        <  1:3726295/93‑1 (MQ=255)

                     ctcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGa                                                                        >  2:4229124/1‑93 (MQ=255)

                      tcGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGAt                                                                       >  1:3531401/1‑93 (MQ=255)

                       cGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATa                                                                      <  1:5264868/93‑1 (MQ=255)

                         tATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGt                                                                    >  1:1361432/1‑93 (MQ=255)

                         tATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGt                                                                    >  1:121794/1‑93 (MQ=255)

                         tATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCACACACCATCGCTGGGCTATCGCAGGCAAAAGATAGt                                                                    >  2:4790989/1‑93 (MQ=255)

                          aTGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTa                                                                   >  2:1292792/1‑93 (MQ=255)

                            tctttGGTTTCCCATACCGCAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAAt                                                                 <  1:1428909/89‑1 (MQ=255)

                              tctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATgg                                                               >  2:1826816/1‑93 (MQ=255)

                               ctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGa                                                              <  2:7095692/93‑1 (MQ=255)

                               ctGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGa                                                              >  1:2159806/1‑93 (MQ=255)

                                tGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGaa                                                             >  2:265168/1‑93 (MQ=255)

                                tGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGaa                                                             >  1:6369231/1‑93 (MQ=255)

                                tGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGaa                                                             <  2:2443396/93‑1 (MQ=255)

                                   tttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaa                                                          >  1:7497359/1‑93 (MQ=255)

                                   tttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaa                                                          >  2:4907446/1‑93 (MQ=255)

                                   tttCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaa                                                          >  1:7005389/1‑93 (MQ=255)

                                     tCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaaaa                                                        <  2:4586278/93‑1 (MQ=255)

                                       ccATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAgg                                                      >  1:855258/1‑93 (MQ=255)

                                        cATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGc                                                     <  2:1213683/93‑1 (MQ=255)

                                        cATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTAGGCTATCGCAGGCAACAGATAGTAATGGAACAAAAGGc                                                     >  2:6990782/1‑93 (MQ=255)

                                          tACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTa                                                   >  2:668044/1‑93 (MQ=255)

                                          tACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTa                                                   >  2:2333285/1‑93 (MQ=255)

                                           aCCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTAt                                                  <  2:6796126/93‑1 (MQ=255)

                                            ccGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATg                                                 <  2:5841291/93‑1 (MQ=255)

                                            ccGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATg                                                 <  1:2544242/93‑1 (MQ=255)

                                            ccGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATg                                                 >  2:3583446/1‑93 (MQ=255)

                                             cGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGt                                                <  1:5672898/93‑1 (MQ=255)

                                              gCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtt                                               >  1:5509865/1‑93 (MQ=255)

                                              gCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtt                                               <  2:1636289/93‑1 (MQ=255)

                                              gCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtt                                               >  1:5747835/1‑93 (MQ=255)

                                              gCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtt                                               >  1:5473200/1‑93 (MQ=255)

                                                aaaGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtttt                                             >  2:955868/1‑93 (MQ=255)

                                                aaaGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGtttt                                             >  1:7339671/1‑93 (MQ=255)

                                                 aaGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGttttt                                            >  2:3799565/1‑93 (MQ=255)

                                                 aaGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGttttt                                            >  1:23621/1‑93 (MQ=255)

                                                  aGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTg                                           <  2:2229405/93‑1 (MQ=255)

                                                    ggTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTa                                         >  2:5515959/1‑93 (MQ=255)

                                                    ggTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTa                                         >  1:4581992/1‑93 (MQ=255)

                                                     gTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAg                                        >  2:5230889/1‑93 (MQ=255)

                                                      tAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAgg                                       >  2:3229109/1‑93 (MQ=255)

                                                      tAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAgg                                       >  1:3711712/1‑93 (MQ=255)

                                                       aGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAggg                                      >  1:2808368/1‑93 (MQ=255)

                                                        gTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAgggg                                     >  2:7238122/1‑93 (MQ=255)

                                                        gTTCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAgggg                                     <  2:5548210/93‑1 (MQ=255)

                                                         ttCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGa                                    >  1:3180215/1‑93 (MQ=255)

                                                         ttCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGa                                    >  2:3049455/1‑93 (MQ=255)

                                                         ttCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGa                                    >  1:4370370/1‑93 (MQ=255)

                                                          tCAGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGaa                                   >  2:5911302/1‑93 (MQ=255)

                                                            aGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATc                                 <  2:673749/93‑1 (MQ=255)

                                                            aGGAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATc                                 >  1:1124559/1‑93 (MQ=255)

                                                             ggAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCg                                >  1:134598/1‑93 (MQ=255)

                                                             ggAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCg                                >  2:7437239/1‑93 (MQ=255)

                                                              gAAAGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGc                               <  2:6478576/93‑1 (MQ=255)

                                                                aaGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTg                             <  1:7145052/93‑1 (MQ=255)

                                                                 aGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTgg                            >  1:1774201/1‑93 (MQ=255)

                                                                 aGTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTgg                            >  2:1136884/1‑93 (MQ=255)

                                                                  gTATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGa                           <  2:2040489/93‑1 (MQ=255)

                                                                   tATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGAt                          <  2:3346912/93‑1 (MQ=255)

                                                                   tATCAATCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGAt                          <  1:4434197/93‑1 (MQ=255)

                                                                       aaTCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTAcc                      <  1:1015618/93‑1 (MQ=255)

                                                                       aaTCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTAcc                      >  2:6321856/1‑93 (MQ=255)

                                                                        aTCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCt                     <  1:6614750/93‑1 (MQ=255)

                                                                        aTCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCt                     <  1:2064586/93‑1 (MQ=255)

                                                                         tCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCtt                    >  1:5875652/1‑93 (MQ=255)

                                                                         tCGCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCtt                    <  2:2808050/93‑1 (MQ=255)

                                                                           gCAGTCAGACACCATCGCTGGGCTATCGCAGGCAAAAGAGAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTAt                  >  1:6219519/1‑93 (MQ=255)

                                                                             aGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCa                >  1:6057180/1‑93 (MQ=255)

                                                                                 aGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGa            <  1:3797999/93‑1 (MQ=255)

                                                                                 aGACACCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGa            <  2:558804/93‑1 (MQ=255)

                                                                                  gACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGAt           >  2:576212/1‑93 (MQ=255)

                                                                                    cacCATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATgg         <  2:685271/93‑1 (MQ=255)

                                                                                       cATCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGc      <  2:6656153/93‑1 (MQ=255)

                                                                                        aTCGCTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGcc     >  1:3479416/1‑93 (MQ=255)

                                                                                        aTCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGcc     <  1:2611726/93‑1 (MQ=255)

GCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGAT  >  NC_000913/3873118‑3873301

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