BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,026,063	A→G	D129G (GAC→GGC)	ygfR →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,063	0	A	G	96.8%	577.5 / 1.8	187	D129G (GAC→GGC)	ygfR	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (2/4); new base G (108/72); total (110/77)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

CTACTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGACTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGTA  >  NC_000913/3025973‑3026152

 tACTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTc                                                                                        <  2:2071369/93‑1 (MQ=255)

  aCTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCa                                                                                       >  1:3014865/1‑93 (MQ=255)

   cTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAg                                                                                      <  1:1434819/93‑1 (MQ=255)

    tGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGGCAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGt                                                                                     <  2:1888875/93‑1 (MQ=255)

      gCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTaa                                                                                   >  1:1048241/1‑93 (MQ=255)

       cACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAg                                                                                  >  1:136107/1‑93 (MQ=255)

       cACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAg                                                                                  >  1:630877/1‑93 (MQ=255)

         ccATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGtt                                                                                >  2:928772/1‑93 (MQ=255)

         ccATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGACCTGACCAGTGGCTCAGTAAGtt                                                                                >  2:3224866/1‑93 (MQ=255)

          cATCGGTGCCGTCGCCGGTCAGGTGAATTTGTTGGATAAAGACAACTAGACCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTc                                                                               <  2:2588172/93‑1 (MQ=255)

          cATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTc                                                                               >  1:1186353/1‑93 (MQ=255)

          cATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTc                                                                               <  1:1559910/93‑1 (MQ=255)

           aTCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGTCCTGACCAGTGGCTCAGTAAGTTCa                                                                              >  1:344818/1‑93 (MQ=255)

           aTCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCa                                                                              >  1:4095/1‑93 (MQ=255)

             ccGTGCCGTCGCCGTTCAGGCTAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAat                                                                            <  2:2149414/93‑1 (MQ=255)

             ccGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAat                                                                            <  2:2149415/93‑1 (MQ=255)

              cGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAata                                                                           <  2:1631573/93‑1 (MQ=255)

              cGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAata                                                                           >  1:3547654/1‑93 (MQ=255)

              cGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAata                                                                           >  2:1191346/1‑93 (MQ=255)

              cGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAata                                                                           <  1:4522186/93‑1 (MQ=255)

                tGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATAtt                                                                         >  1:2558479/1‑93 (MQ=255)

                 gCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtc                                                                        <  1:2620403/93‑1 (MQ=255)

                  ccGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtct                                                                       >  2:2365867/1‑93 (MQ=255)

                   cGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtctc                                                                      <  1:2093487/93‑1 (MQ=255)

                   cGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATtctc                                                                      <  2:1824652/93‑1 (MQ=255)

                   cGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCGCAGTAAGTTCAATATtctc                                                                      <  2:3595978/93‑1 (MQ=255)

                    gTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTcc                                                                     >  2:3150924/1‑93 (MQ=255)

                      cGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCgg                                                                   >  2:3686174/1‑93 (MQ=255)

                       gCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGc                                                                  <  2:1138947/93‑1 (MQ=255)

                        ccGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGcc                                                                 <  1:2499101/93‑1 (MQ=255)

                        ccGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGcc                                                                 <  2:3891044/93‑1 (MQ=255)

                         cGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCt                                                                <  2:320790/93‑1 (MQ=255)

                         cGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCt                                                                <  2:859677/93‑1 (MQ=255)

                          ggTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtg                                                               <  1:1344550/93‑1 (MQ=255)

                          ggTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtg                                                               <  2:202842/93‑1 (MQ=255)

                           gTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtgg                                                              >  1:3492357/1‑93 (MQ=255)

                            gcAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtggt                                                             <  1:974379/92‑1 (MQ=255)

                             cAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtggtg                                                            >  2:13932/1‑93 (MQ=255)

                              aGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtggtgg                                                           >  2:4647709/1‑93 (MQ=255)

                              aGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCtggtgg                                                           <  2:1061539/93‑1 (MQ=255)

                               ggCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTggg                                                          >  2:3147938/1‑93 (MQ=255)

                               ggCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTggg                                                          >  2:1415788/1‑93 (MQ=255)

                               ggCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTggg                                                          <  1:4705226/93‑1 (MQ=255)

                               ggCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTggg                                                          <  1:2338178/93‑1 (MQ=255)

                                   aaTTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCa                                                      >  1:2993505/1‑93 (MQ=255)

                                    aTTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCCGCCTGGGGGGCGCAg                                                     >  1:2277014/1‑93 (MQ=255)

                                     tttGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgc                                                    <  1:1768481/93‑1 (MQ=255)

                                     tttGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgc                                                    <  2:36167/93‑1 (MQ=255)

                                      ttgttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgcg                                                   <  2:792862/93‑1 (MQ=255)

                                       tgttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAgcgc                                                  <  1:3978250/93‑1 (MQ=255)

                                         ttgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGccc                                                <  1:3031970/93‑1 (MQ=255)

                                          tgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCg                                               >  1:2910939/1‑93 (MQ=255)

                                          tgGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCg                                               <  1:3891553/93‑1 (MQ=255)

                                                  gACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGtt                                       <  2:4541963/93‑1 (MQ=255)

                                                  gACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGtt                                       >  1:4160176/1‑93 (MQ=255)

                                                     aaCCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTtat                                    >  1:2710831/1‑93 (MQ=255)

                                                     aaCCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTtat                                    <  2:3566271/93‑1 (MQ=255)

                                                      aCCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTtata                                   >  1:1043354/1‑93 (MQ=255)

                                                        cAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATc                                 >  2:2647626/1‑93 (MQ=255)

                                                        cAGATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATc                                 >  2:2539534/1‑93 (MQ=255)

                                                          gATCATCAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGa                               >  1:2431687/1‑93 (MQ=255)

                                                           atcatcAACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGaa                              >  1:1887577/1‑93 (MQ=255)

                                                             catcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATc                            >  1:2221834/1‑93 (MQ=255)

                                                             aatcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATc                            <  2:4330437/92‑1 (MQ=255)

                                                               tcaACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAg                          <  1:1713206/93‑1 (MQ=255)

                                                                caACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGc                         >  2:4490560/1‑93 (MQ=255)

                                                                caACGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGc                         >  2:3754739/1‑93 (MQ=255)

                                                                 aaCGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCg                        >  2:3821970/1‑93 (MQ=255)

                                                                 aaCGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCg                        >  1:2172102/1‑93 (MQ=255)

                                                                 aaCGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCg                        >  1:4642812/1‑93 (MQ=255)

                                                                 aaCGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCg                        >  1:3695174/1‑93 (MQ=255)

                                                                 aaCGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCg                        >  2:3108764/1‑93 (MQ=255)

                                                                  aCGGCGGCAAAGCCCTGACCTGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCgg                       >  1:3942192/1‑93 (MQ=255)

                                                                  aCGGCGGCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCgg                       >  1:3599087/1‑93 (MQ=255)

                                                                   cggcggCAAAGCCCTTGCCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGc                      <  2:2618014/93‑1 (MQ=255)

                                                                   cggcggCAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGc                      >  2:4404728/1‑93 (MQ=255)

                                                                    ggcggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCa                     <  1:3670145/93‑1 (MQ=255)

                                                                    ggcggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCa                     <  2:508350/93‑1 (MQ=255)

                                                                     gcggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAg                    >  1:3940308/1‑93 (MQ=255)

                                                                     gcggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAg                    >  2:401726/1‑93 (MQ=255)

                                                                     gcggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTACGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAg                    >  2:672743/1‑93 (MQ=255)

                                                                       ggcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGa                  <  1:2960268/93‑1 (MQ=255)

                                                                        gcAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGaa                 <  2:3274235/93‑1 (MQ=255)

                                                                         cAAAGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAc                >  2:2330176/1‑93 (MQ=255)

                                                                           aaGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccg              >  2:2411126/1‑93 (MQ=255)

                                                                            aGCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgc             >  1:1608044/1‑93 (MQ=255)

                                                                             gCCCTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgcc            >  2:167288/1‑93 (MQ=255)

                                                                               ccTGACCAGTGTCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgccgc          <  2:1265568/93‑1 (MQ=255)

                                                                               ccTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgccgc          <  2:4095/93‑1 (MQ=255)

                                                                                cTGACCCGTGACTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGCCTCATCGGCAGGACccgccccc         >  2:3003997/1‑93 (MQ=255)

                                                                                cTGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAacgccgcg         >  2:1267587/1‑92 (MQ=255)

                                                                                 tGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgccgccg        >  2:3416150/1‑93 (MQ=255)

                                                                                 tGACCAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAAccgccgccg        >  2:1533263/1‑93 (MQ=255)

                                                                                  gACCAGTGACTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCgg       <  2:731455/93‑1 (MQ=255)

                                                                                    ccAGTGGCTCAGTAAGTTCAATATTCTGCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGACGCCgggg     >  2:1605612/1‑93 (MQ=255)

                                                                                     cAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCggggg    <  2:3222740/93‑1 (MQ=255)

                                                                                     cAGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCggggg    <  1:2106847/93‑1 (MQ=255)

                                                                                      aGTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGt   <  1:2316227/93‑1 (MQ=255)

                                                                                       gTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGTGGGTa  <  1:796690/93‑1 (MQ=255)

                                                                                       gTGGCTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGTa  <  1:366276/93‑1 (MQ=255)

CTACTGGCACCATCCGTGCCGTCGCCGGTCAGGCGAATTTGTTGGATAAAGACAACCAGATCATCAACGGCGGCAAAGCCCTGACCAGTGACTCAGTAAGTTCAATATTCTCCGGCCTGGTGGGCGCAGCGCCCGCAGCGGTTTATATCGAATCAGCGGCAGGAACCGCCGCCGGGGGTA  >  NC_000913/3025973‑3026152

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