BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,241,465	T→C	D397G (GAC→GGC)	malF ←	MalF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,241,465	0	T	C	99.7%	1068.1 / ‑5.6	340	D397G (GAC→GGC)	malF	MalF
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (175/164); total (175/165)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 357 total aligned reads displayed.

GTCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGTCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTGCGGGCG  >  NC_000913/4241375‑4241554

 tCAGCGGTTTAATCAGCAGCGGCAGCGTAATTTTAAAGAAGTTCTTGAACTGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCg                                                                                        >  2:6076420/1‑93 (MQ=255)

 tCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGTCATCTATTGCGGAGGCTACATACAAATCGCCCg                                                                                        >  1:7040062/1‑93 (MQ=255)

 tCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCg                                                                                        >  2:3645269/1‑93 (MQ=255)

 tCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCg                                                                                        >  1:4338045/1‑93 (MQ=255)

   aGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGa                                                                                      >  1:27550/1‑93 (MQ=255)

   aGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCACTGCTGAGGCTTCATACAAATCGCCCGGa                                                                                      >  1:5996804/1‑93 (MQ=255)

   aGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAATTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAACTCGCCCGGa                                                                                      >  2:6211560/1‑93 (MQ=255)

    gCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGaa                                                                                     <  1:6832717/93‑1 (MQ=255)

     cGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAAt                                                                                    >  1:6244559/1‑93 (MQ=255)

      ggTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATc                                                                                   >  1:1361029/1‑93 (MQ=255)

      ggTTTAATCAGCAGCGGCAGCCTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATc                                                                                   >  1:1388089/1‑93 (MQ=255)

       gTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCg                                                                                  >  1:7481631/1‑93 (MQ=255)

        tttAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGc                                                                                 <  2:5812670/93‑1 (MQ=255)

        tttAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGc                                                                                 <  1:2924713/93‑1 (MQ=255)

        tttAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCAATGCTGAGGCTTCATACAAATCGCCCGGAATCGc                                                                                 >  1:7654799/1‑93 (MQ=255)

         ttAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCt                                                                                <  2:4165354/93‑1 (MQ=255)

           aaTCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCttt                                                                              >  2:4675672/1‑93 (MQ=255)

            aTCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTc                                                                             <  1:4753569/93‑1 (MQ=255)

             tCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCa                                                                            <  1:1455326/93‑1 (MQ=255)

             tCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCa                                                                            <  2:1719609/93‑1 (MQ=255)

                gcagcggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCa                                                                         >  2:3783118/1‑93 (MQ=255)

                gcagcggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCa                                                                         >  1:4393580/1‑93 (MQ=255)

                gcagcggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATAGCCCGGAATCGCTTTCAGCa                                                                         >  2:3706932/1‑93 (MQ=255)

                     ggcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGccc                                                                    <  1:368234/93‑1 (MQ=255)

                      gcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCa                                                                   >  1:3565658/1‑93 (MQ=255)

                      gcagcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCa                                                                   >  2:7584985/1‑93 (MQ=255)

                         gcgTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGc                                                                >  1:6925738/1‑93 (MQ=255)

                           gTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCag                                                              >  1:1911989/1‑93 (MQ=255)

                            tAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCaga                                                             >  2:3086059/1‑93 (MQ=255)

                            tAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCaga                                                             <  1:7892771/93‑1 (MQ=255)

                             aaTCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCagag                                                            >  1:5360010/1‑93 (MQ=255)

                             aaTCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCagag                                                            >  2:2149004/1‑93 (MQ=255)

                              aTCTTAAAGAAGTTCTGGAACGGACGTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGGTTTCAGCAAGCCCATGCAGAgg                                                           >  1:4985495/1‑93 (MQ=255)

                               tCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGa                                                          >  1:3199152/1‑93 (MQ=255)

                               tCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGa                                                          <  2:6780046/93‑1 (MQ=255)

                                cTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGat                                                         <  1:1954840/93‑1 (MQ=255)

                                   aaaGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGatcat                                                      <  1:1428372/93‑1 (MQ=255)

                                   aaaGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTCTCAGCAAGCCCATGCAGAGGatcat                                                      <  2:449870/93‑1 (MQ=255)

                                    aagaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGatcatc                                                     >  2:3965467/1‑93 (MQ=255)

                                    aagaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGatcatc                                                     >  2:5165152/1‑93 (MQ=255)

                                     agaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGAtcatca                                                    >  1:6155550/1‑93 (MQ=255)

                                      gaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATcatcat                                                   >  1:3760421/1‑93 (MQ=255)

                                      gaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATcatcat                                                   <  1:3315414/93‑1 (MQ=255)

                                      gaagTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATcatcat                                                   >  1:2963974/1‑93 (MQ=255)

                                        agTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGt                                                 >  1:920311/1‑93 (MQ=255)

                                        agTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGt                                                 >  2:3467974/1‑93 (MQ=255)

                                         gTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTa                                                <  1:1020575/93‑1 (MQ=255)

                                         attCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTa                                                <  1:2549377/92‑1 (MQ=255)

                                          ttCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTAc                                               >  2:2261305/1‑93 (MQ=255)

                                              ggAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGAt                                           >  2:1469385/1‑93 (MQ=255)

                                               gAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATa                                          <  1:5142642/93‑1 (MQ=255)

                                               gAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATa                                          <  1:7806285/93‑1 (MQ=255)

                                                aaCGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATaa                                         >  1:5860033/1‑93 (MQ=255)

                                                 aCGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAAc                                        >  1:1163419/1‑93 (MQ=255)

                                                   ggACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAAccc                                      <  1:3696617/93‑1 (MQ=255)

                                                   ggACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAAccc                                      <  2:642223/93‑1 (MQ=255)

                                                    gACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACcca                                     <  2:1182374/93‑1 (MQ=255)

                                                    gACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACcca                                     >  1:4695045/1‑93 (MQ=255)

                                                    gACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACcca                                     >  1:552347/1‑93 (MQ=255)

                                                    gACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACcca                                     >  2:2667012/1‑93 (MQ=255)

                                                    gACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACcca                                     <  1:6582630/93‑1 (MQ=255)

                                                      ccTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagc                                   >  2:6676373/1‑93 (MQ=255)

                                                      ccTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagc                                   <  1:2096434/93‑1 (MQ=255)

                                                      ccTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagc                                   <  2:5573357/93‑1 (MQ=255)

                                                       cTGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagcc                                  >  2:6922330/1‑93 (MQ=255)

                                                        tGCGCCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagcca                                 <  2:2646922/93‑1 (MQ=255)

                                                         gcgcCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagccag                                <  1:2291667/93‑1 (MQ=255)

                                                         gcgcCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagccag                                <  2:3199791/93‑1 (MQ=255)

                                                         gcgcCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACccagccag                                <  2:57405/93‑1 (MQ=255)

                                                           gcCATCCATTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGt                              <  2:2395666/93‑1 (MQ=255)

                                                            ccatccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTa                             <  1:5816234/93‑1 (MQ=255)

                                                            ccatccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTa                             <  2:6175939/93‑1 (MQ=255)

                                                             catccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTAt                            >  1:3887284/1‑93 (MQ=255)

                                                              atccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTAtt                           >  1:3348826/1‑93 (MQ=255)

                                                              atccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTAtt                           <  1:717549/93‑1 (MQ=255)

                                                              atccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTAtt                           >  1:422200/1‑93 (MQ=255)

                                                               tccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTAtta                          >  2:5350795/1‑92 (MQ=255)

                                                               tccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTg                          >  2:484366/1‑93 (MQ=255)

                                                                ccatTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGa                         <  1:5641165/93‑1 (MQ=255)

                                                                ccatTGCTGAGGCTTCATACAAACCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGa                         <  2:4295805/93‑1 (MQ=255)

                                                                 catTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGAc                        >  1:16957/1‑93 (MQ=255)

                                                                 catTGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGAc                        <  2:2070952/93‑1 (MQ=255)

                                                                    tGCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGAt                     <  2:7747557/93‑1 (MQ=255)

                                                                     gCTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATa                    <  2:7243667/93‑1 (MQ=255)

                                                                      cTGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATaa                   <  2:1530315/93‑1 (MQ=255)

                                                                       tGAGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAAt                  >  2:7182427/1‑93 (MQ=255)

                                                                         acgCTTCATACAAAGCGCCCGGAATCTCTGTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTa                <  1:6400703/91‑1 (MQ=255)

                                                                         aGGCTTCATACACATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCCGGTATTGACGATAATTa                <  2:7030901/93‑1 (MQ=255)

                                                                         aGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTa                <  1:4024611/93‑1 (MQ=255)

                                                                         aGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTa                <  2:7868283/93‑1 (MQ=255)

                                                                         aGGCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTa                <  1:50555/93‑1 (MQ=255)

                                                                           gCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGc              >  2:6537031/1‑93 (MQ=255)

                                                                           gCTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGc              >  1:3500526/1‑93 (MQ=255)

                                                                            cTTCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCa             <  2:5034987/93‑1 (MQ=255)

                                                                             ttCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCAt            >  2:6361675/1‑93 (MQ=255)

                                                                             ttCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCAt            <  2:434577/93‑1 (MQ=255)

                                                                              tCATACAAATCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATc           <  2:27550/93‑1 (MQ=255)

                                                                                    aaaTCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgcgg     <  1:4587301/93‑1 (MQ=255)

                                                                                     aaTCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgcggg    <  2:3000711/93‑1 (MQ=255)

                                                                                       tCGCCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTgcgggcg  >  1:4126610/1‑93 (MQ=255)

GTCAGCGGTTTAATCAGCAGCGGCAGCGTAATCTTAAAGAAGTTCTGGAACGGACCTGCGCCATCCATTGCTGAGGCTTCATACAAATCGTCCGGAATCGCTTTCAGCAAGCCCATGCAGAGGATCATCATGTACGGATAACCCAGCCAGGTATTGACGATAATTAGCATCGTGCGGGCG  >  NC_000913/4241375‑4241554

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