BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,845,852	T→C	T61A (ACC→GCC)	selD ←	selenophosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,845,852	0	T	C	99.6%	734.9 / ‑4.9	230	T61A (ACC→GCC)	selD	selenophosphate synthase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (159/70); total (159/71)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 241 total aligned reads displayed.

GCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAT  >  NC_000913/1845763‑1845941

gCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGcc                                                                                        >  1:723574/1‑93 (MQ=255)

   gATGGCGTTAGTCGCCGCAATGCGGCCCAAATCGAAAAGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAtt                                                                                     >  2:6524185/1‑93 (MQ=255)

   gATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAtt                                                                                     >  1:6890740/1‑93 (MQ=255)

   gATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAtt                                                                                     >  2:1643004/1‑93 (MQ=255)

     tGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGc                                                                                   >  2:8615986/1‑93 (MQ=255)

      ggCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGcc                                                                                  <  2:9178585/93‑1 (MQ=255)

      ggCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGcc                                                                                  <  2:7546198/93‑1 (MQ=255)

       gCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGccc                                                                                 >  2:7238044/1‑93 (MQ=255)

        cGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCa                                                                                <  1:4481031/93‑1 (MQ=255)

         gTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAg                                                                               >  1:1979011/1‑93 (MQ=255)

         gTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAg                                                                               >  1:9210632/1‑93 (MQ=255)

         gTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAg                                                                               >  1:10071/1‑93 (MQ=255)

          ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                              <  1:185490/93‑1 (MQ=255)

          ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                              >  1:8926014/1‑93 (MQ=255)

          ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                              >  1:473690/1‑93 (MQ=255)

           tAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGCCGCCATTGACCAGac                                                                             >  2:1814866/1‑92 (MQ=255)

            aGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTAGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATc                                                                            >  2:2668624/1‑93 (MQ=255)

             gTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCg                                                                           >  1:5449820/1‑93 (MQ=255)

               cgccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTa                                                                         <  1:6752100/93‑1 (MQ=255)

               cgccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTa                                                                         <  1:9735317/93‑1 (MQ=255)

                   gcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAcc                                                                     >  2:8311041/1‑93 (MQ=255)

                    cAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCCTAAAGAAGTCGGTGGTACTGATAACGCTGGCGACATTGCCCAGATCGTACAgcg                                                                    >  2:3069844/1‑93 (MQ=255)

                    cAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccg                                                                    >  1:7918526/1‑93 (MQ=255)

                     aaTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgc                                                                   <  1:5885636/93‑1 (MQ=255)

                      aTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgcc                                                                  >  1:572852/1‑93 (MQ=255)

                      aTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgcc                                                                  >  1:6107170/1‑93 (MQ=255)

                       tGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCAATTGCCCAGATCGCACAccgccg                                                                 >  2:7933590/1‑93 (MQ=255)

                        gCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACcgccgc                                                                <  2:7652672/93‑1 (MQ=255)

                          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAt                                                              >  1:1308414/1‑93 (MQ=255)

                          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAt                                                              <  2:5769061/93‑1 (MQ=255)

                          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAt                                                              >  1:945664/1‑93 (MQ=255)

                             cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                           >  2:5628960/1‑93 (MQ=255)

                              aaaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtc                                                          >  2:68354/1‑93 (MQ=255)

                               aaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtcg                                                         >  2:5376209/1‑93 (MQ=255)

                               aaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtcg                                                         >  1:8409916/1‑93 (MQ=255)

                                 aTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTcgcg                                                       <  1:251919/93‑1 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  2:1884343/1‑93 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  2:2534379/1‑93 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  2:9423861/1‑93 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  2:3619075/1‑93 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  1:4357772/1‑93 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    <  1:7199772/93‑1 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGGGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  2:4271593/1‑93 (MQ=255)

                                      aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCCGATCGTACACCGCCGCATCGTCGCGGGttt                                                  >  1:3058043/1‑93 (MQ=255)

                                      aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                  >  2:8394879/1‑93 (MQ=255)

                                      aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                  >  1:1167895/1‑93 (MQ=255)

                                       aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                 >  1:1328064/1‑93 (MQ=255)

                                           ttATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTa                                             >  1:659251/1‑93 (MQ=255)

                                           ttATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTa                                             <  2:5099324/93‑1 (MQ=255)

                                            tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                            <  1:1647766/93‑1 (MQ=255)

                                            tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                            <  2:49535/93‑1 (MQ=255)

                                            tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                            >  2:5931180/1‑93 (MQ=255)

                                             aTCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGGTTCATTAcc                                           >  2:6261249/1‑93 (MQ=255)

                                              tCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAccc                                          >  2:4342386/1‑93 (MQ=255)

                                              tCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAccc                                          >  1:2434870/1‑93 (MQ=255)

                                               cAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCca                                         <  1:7690707/93‑1 (MQ=255)

                                                  cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                      <  2:4020118/93‑1 (MQ=255)

                                                   gATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAg                                     >  1:8779827/1‑93 (MQ=255)

                                                    aTCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGc                                    <  1:179778/93‑1 (MQ=255)

                                                     tCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCa                                   >  2:7101011/1‑93 (MQ=255)

                                                         cATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGATTCATTACCCACAAGCAAAtc                               >  1:2140914/1‑92 (MQ=255)

                                                         cATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGACGCATCGTCGCGGGTTTCATTACCCACAAGCAAAtt                               >  2:10046491/1‑93 (MQ=255)

                                                          aTAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTc                              >  1:705453/1‑93 (MQ=255)

                                                            aaaGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCgg                            >  2:8038378/1‑93 (MQ=255)

                                                             aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGa                           >  1:9534248/1‑93 (MQ=255)

                                                             aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCAGa                           <  2:4567938/93‑1 (MQ=255)

                                                             aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGCTTCATTACCCACAAGCAAATTCGGa                           >  1:8479244/1‑93 (MQ=255)

                                                              agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                          >  1:2623910/1‑93 (MQ=255)

                                                               gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACACGCAAATTCGGATc                         >  1:2236096/1‑93 (MQ=255)

                                                               gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                         >  2:2935179/1‑93 (MQ=255)

                                                                  gTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaac                      <  2:9064690/93‑1 (MQ=255)

                                                                   tCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaaca                     >  1:2319942/1‑93 (MQ=255)

                                                                   tCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaaca                     <  1:825811/93‑1 (MQ=255)

                                                                   tCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaaca                     <  1:995539/93‑1 (MQ=255)

                                                                   tCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATAGTAAACCGCCGCATCGTCGCGGGTATCATTACCCCCAAGCAAATTCCGATcaaca                     >  2:4794947/1‑93 (MQ=255)

                                                                     ggtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACaaa                   >  2:7355887/1‑93 (MQ=255)

                                                                      gtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAAc                  <  2:5466036/93‑1 (MQ=255)

                                                                      gtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAAc                  >  2:4696317/1‑93 (MQ=255)

                                                                       tggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCTGATCAACAAACt                 >  2:6299464/1‑93 (MQ=255)

                                                                        ggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACtt                <  2:325255/93‑1 (MQ=255)

                                                                         gtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc               <  1:4170265/93‑1 (MQ=255)

                                                                         gtACTGATAACGCTGGCGCCATTGCCCAGATAGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc               <  1:6164888/93‑1 (MQ=255)

                                                                          tACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCg              >  2:2092259/1‑93 (MQ=255)

                                                                          tACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCg              >  1:3781666/1‑93 (MQ=255)

                                                                           aCTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGc             >  2:6955415/1‑93 (MQ=255)

                                                                           aCTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGc             <  1:577844/93‑1 (MQ=255)

                                                                            cTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGcc            >  1:4582764/1‑93 (MQ=255)

                                                                             tGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCGGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCt           >  2:7831922/1‑93 (MQ=255)

                                                                             tGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCt           >  1:3153559/1‑93 (MQ=255)

                                                                              gATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTAATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTg          >  1:7376057/1‑93 (MQ=255)

                                                                              gATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGGTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGACTg          >  2:3729169/1‑93 (MQ=255)

                                                                               aTAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGc         <  2:1366204/93‑1 (MQ=255)

                                                                                tAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCt        <  1:5312126/93‑1 (MQ=255)

                                                                                tAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGCTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCt        >  1:6486250/1‑93 (MQ=255)

                                                                                  aCGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCa      <  2:6492891/93‑1 (MQ=255)

                                                                                  aCGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCa      <  2:1171692/93‑1 (MQ=255)

                                                                                   cGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc     >  1:4020300/1‑93 (MQ=255)

                                                                                     cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa   <  1:3373531/93‑1 (MQ=255)

                                                                                      tGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAt  >  2:9511539/1‑93 (MQ=255)

                                                                                      tGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAt  >  1:1096827/1‑93 (MQ=255)

GCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAT  >  NC_000913/1845763‑1845941

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