BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,627	Δ1 bp	coding (506/672 nt)	yidX →	putative replicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,625	0	G	.	99.7%	2872.0 / ‑1.9	631	coding (504/672 nt)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base G (0/2); new base . (328/301); total (328/303)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 695 total aligned reads displayed.

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

  aaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCgga                                                                                        >  2:6606566/1‑92 (MQ=255)

  aaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCgga                                                                                        <  2:7010657/93‑2 (MQ=255)

  aaaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCgga                                                                                        >  1:3348798/1‑92 (MQ=255)

  aaaTAAAAACATGGATTACTCAAAGGGGATGGCGTTCGACCATCCATTTAAAGAGCGATGTTATGAATACTATTCACCCAGAGGCATTCCgga                                                                                        >  2:8981299/1‑92 (MQ=255)

   aaTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggat                                                                                       >  2:1643650/1‑91 (MQ=255)

    aTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatg                                                                                      <  2:7350673/93‑4 (MQ=255)

    aTAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatg                                                                                      <  2:952754/93‑4 (MQ=255)

     tAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatgg                                                                                     >  2:3628894/1‑89 (MQ=255)

     tAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatgg                                                                                     <  1:9118381/93‑5 (MQ=255)

     tAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCggatgg                                                                                     <  1:2625544/93‑5 (MQ=255)

      aaaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgt                                                                                   >  1:926038/1‑93 (MQ=255)

       aaaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtg                                                                                  <  1:5520053/93‑1 (MQ=255)

        aaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtgt                                                                                 <  1:5887538/93‑1 (MQ=255)

        aaaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGgtgt                                                                                 >  2:817973/1‑93 (MQ=255)

         aaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGtt                                                                                >  2:1946028/1‑93 (MQ=255)

         aaCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGtt                                                                                <  2:5032298/93‑1 (MQ=255)

          aCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGttt                                                                               >  2:4817165/1‑93 (MQ=255)

          aCATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGttt                                                                               <  1:4555349/93‑1 (MQ=255)

           cATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTc                                                                              <  2:2250010/93‑1 (MQ=255)

           cATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTc                                                                              >  2:5617507/1‑93 (MQ=255)

            aTGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGtttac                                                                             >  1:9650232/1‑91 (MQ=255)

             tGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCg                                                                            >  2:3442690/1‑93 (MQ=255)

                aTTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcag                                                                         >  2:64601/1‑93 (MQ=255)

                  tACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcagca                                                                       >  1:8209162/1‑93 (MQ=255)

                  tACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCgcagca                                                                       <  1:10021421/93‑1 (MQ=255)

                   aCTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCAt                                                                      <  1:1528039/93‑1 (MQ=255)

                     tCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTa                                                                    >  2:4041090/1‑93 (MQ=255)

                     tCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTa                                                                    >  2:1448201/1‑93 (MQ=255)

                      cAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTAc                                                                   >  2:9614347/1‑93 (MQ=255)

                      cAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTAc                                                                   <  2:7230323/93‑1 (MQ=255)

                        aaGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTg                                                                 >  2:9093908/1‑93 (MQ=255)

                         aGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGc                                                                <  2:7577130/93‑1 (MQ=255)

                           gTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCa                                                              <  2:8598677/93‑1 (MQ=255)

                            tGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAg                                                             <  2:2710421/93‑1 (MQ=255)

                            tGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAg                                                             <  2:5873823/93‑1 (MQ=255)

                             gATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGACGCATTCC‑GGATGGTGTTTGCGCAGCATTACTGCCAgg                                                            >  2:1790800/1‑92 (MQ=255)

                               tGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGa                                                          <  1:1681520/93‑1 (MQ=255)

                                 gCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACt                                                        <  1:6612216/93‑1 (MQ=255)

                                  cGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTg                                                       <  1:4321813/93‑1 (MQ=255)

                                   gTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACCGt                                                      >  2:9738029/1‑93 (MQ=255)

                                    ttCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGtt                                                     >  1:8684289/1‑93 (MQ=255)

                                    ttCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGtt                                                     >  1:7418884/1‑93 (MQ=255)

                                    ttCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGtt                                                     <  1:7035080/93‑1 (MQ=255)

                                     tCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGttc                                                    >  1:2181723/1‑92 (MQ=255)

                                     tCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTa                                                    >  1:6549836/1‑93 (MQ=255)

                                       tACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACg                                                  >  1:250160/1‑93 (MQ=255)

                                       tACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACg                                                  >  2:3302845/1‑93 (MQ=255)

                                         ccatccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCt                                                <  2:1358966/93‑1 (MQ=255)

                                           atccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTgg                                              >  2:7934258/1‑93 (MQ=255)

                                            tccatTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGa                                             >  1:4120077/1‑93 (MQ=255)

                                              catTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAca                                           <  1:4952007/93‑1 (MQ=255)

                                              catTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAca                                           >  1:6413115/1‑93 (MQ=255)

                                               atTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcat                                          <  2:306874/93‑1 (MQ=255)

                                                tttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatc                                         >  2:435384/1‑93 (MQ=255)

                                                 ttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatca                                        <  2:3948279/93‑1 (MQ=255)

                                                 ttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatca                                        >  1:7829405/1‑93 (MQ=255)

                                                 ttAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAcatca                                        >  1:950188/1‑93 (MQ=255)

                                                   aaaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCAtt                                      <  2:10026594/93‑1 (MQ=255)

                                                   aaaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCAtt                                      >  1:985158/1‑93 (MQ=255)

                                                   aaaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCAtt                                      >  2:2812956/1‑93 (MQ=255)

                                                    aaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTa                                     >  2:6319595/1‑93 (MQ=255)

                                                    aaGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTa                                     >  1:5978768/1‑93 (MQ=255)

                                                      gTGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACt                                   >  1:3586083/1‑93 (MQ=255)

                                                       tGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACtc                                  >  1:7631241/1‑92 (MQ=255)

                                                       tGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTg                                  >  1:3862027/1‑93 (MQ=255)

                                                       tGCGATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTg                                  >  2:9349735/1‑93 (MQ=255)

                                                          gATTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGGTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAc                               >  1:1114331/1‑93 (MQ=255)

                                                           aTTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAcc                              >  1:6306139/1‑93 (MQ=255)

                                                           aTTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAcc                              >  1:3511953/1‑93 (MQ=255)

                                                           aTTTTATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCAcc                              <  2:5227556/93‑1 (MQ=255)

                                                             tttATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGc                            >  1:2417248/1‑93 (MQ=255)

                                                               tATGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTg                          >  1:8500641/1‑93 (MQ=255)

                                                               tATGAATACTATTCACAGAGAGGCATT‑CGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTg                          <  2:540295/93‑1 (MQ=255)

                                                                aTGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGc                         >  2:8318976/1‑93 (MQ=255)

                                                                aTGAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGc                         <  1:318081/93‑1 (MQ=255)

                                                                  gAATACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCaa                       >  2:4241510/1‑93 (MQ=255)

                                                                   aaTACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAt                      <  2:7657050/93‑1 (MQ=255)

                                                                   aaTACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAt                      >  2:2081466/1‑93 (MQ=255)

                                                                    aTACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAtt                     <  1:5261973/93‑1 (MQ=255)

                                                                    aTACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAtt                     >  2:4464746/1‑93 (MQ=255)

                                                                     tACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAttt                    >  1:9361352/1‑93 (MQ=255)

                                                                     tACTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAAttt                    >  2:1093699/1‑93 (MQ=255)

                                                                      aCTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTc                   >  1:6967850/1‑93 (MQ=255)

                                                                      aCTATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTc                   <  1:371060/93‑1 (MQ=255)

                                                                        tATTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTg                 <  1:1261474/93‑1 (MQ=255)

                                                                         aTTCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGCCATCCTTACTCCACCGCTGCACTTTCTgt                >  1:2056294/1‑92 (MQ=255)

                                                                          ttCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGt               <  1:1324499/93‑1 (MQ=255)

                                                                           tCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGtt              <  2:6728457/93‑1 (MQ=255)

                                                                           tCACCCAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGtt              >  1:6819749/1‑93 (MQ=255)

                                                                              cccAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGta           >  1:4756543/1‑93 (MQ=255)

                                                                                cAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtata         >  1:7957273/1‑93 (MQ=255)

                                                                                cAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtata         <  2:8878493/93‑1 (MQ=255)

                                                                                cAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGtata         <  1:3035493/93‑1 (MQ=255)

                                                                                cAGAGGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGAAATCATTACTGCACGGCTGCAATTTCTGGTTGtata         <  1:1813953/93‑1 (MQ=255)

                                                                                 agagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTatat        >  1:7229246/1‑93 (MQ=255)

                                                                                  gagGCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCCGCAATTTCTGGTTGTATATg       >  2:4661287/1‑93 (MQ=255)

                                                                                   agGCATTCC‑GGATGGTGTTTCCGCAGCATTTCTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGc      >  2:2401508/1‑93 (MQ=255)

                                                                                    ggCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCa     <  2:3216332/93‑1 (MQ=255)

                                                                                    ggCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGCTGTATATGCa     <  2:2643037/93‑1 (MQ=255)

                                                                                     gCATTCC‑GGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAg    <  2:6194672/93‑1 (MQ=255)

AAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCATTTAAAGTGCGATTTTATGAATACTATTCACCCAGAGGCATTCCGGGATGGTGTTTCCGCAGCATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTA  >  NC_000913/3873533‑3873713

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