BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,675,076	+G	coding (920/2076 nt)	yhjG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,675,072	1	.	G	82.1%	777.2 / 170.0	279	coding (924/2076 nt)	yhjG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (19/31); new base G (128/101); total (147/132)
Fisher's exact test for biased strand distribution p-value = 2.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 317 total aligned reads displayed.

TCACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGACACCGCC  >  NC_000913/3674980‑3675158

tCACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc                                                                                         <  1:3145454/93‑1 (MQ=255)

 cACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGTCGACCATCCGTTTCAAAcggc‑g                                                                                       <  2:2708528/93‑1 (MQ=255)

  gccAATTCGCCCATTAAAACCACGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑gg                                                                                      <  1:7432768/92‑1 (MQ=255)

  aCCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑gg                                                                                      >  2:3694032/1‑93 (MQ=255)

  aCCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑gg                                                                                      <  1:3382102/93‑1 (MQ=255)

  aCCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑gg                                                                                      >  1:217830/1‑93 (MQ=255)

  aCCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAAcggc‑gg                                                                                      <  1:7019865/93‑1 (MQ=255)

   ccAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑ggg                                                                                     >  1:3768977/1‑93 (MQ=255)

   ccAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑ggg                                                                                     <  2:7325825/93‑1 (MQ=255)

   ccAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑ggg                                                                                     <  2:6891145/93‑1 (MQ=255)

    cAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gggg                                                                                    <  1:189000/93‑1 (MQ=255)

      aTTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gggggt                                                                                  >  2:651418/1‑91 (MQ=255)

       ttCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gggggta                                                                                 >  2:2482403/1‑90 (MQ=255)

       ttCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gggggta                                                                                 <  2:6270401/93‑4 (MQ=255)

       ttCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gggggta                                                                                 <  1:2478478/93‑4 (MQ=255)

       ttCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑gggggga                                                                                 >  1:1709236/1‑90 (MQ=255)

          gCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCg                                                                               >  1:6969880/1‑93 (MQ=255)

            ccATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCggg                                                                             >  2:5809726/1‑93 (MQ=255)

            ccATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCggg                                                                             >  2:1723457/1‑93 (MQ=255)

            ccATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCggg                                                                             <  2:3750560/93‑1 (MQ=255)

            ccATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCggg                                                                             >  2:2584582/1‑93 (MQ=255)

            ccATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCggc                                                                             <  2:7409120/93‑2 (MQ=255)

             cATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGgc                                                                            >  2:2013330/1‑93 (MQ=255)

              aTTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGgca                                                                           >  1:1145137/1‑93 (MQ=255)

              aTTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCATTTCAAACGGCGGGGGTATCGGgca                                                                           >  1:4746939/1‑93 (MQ=255)

               ttAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGgcag                                                                          >  1:6367494/1‑93 (MQ=255)

               ttAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGgcag                                                                          >  2:5543568/1‑93 (MQ=255)

               ttAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGgcag                                                                          >  1:1452441/1‑93 (MQ=255)

                tAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGgcagc                                                                         <  2:4837879/93‑1 (MQ=255)

                 aaaaCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGcagcag                                                                       >  1:2979077/1‑93 (MQ=255)

                  aaaCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCAGGcagcag                                                                       >  1:4188391/1‑93 (MQ=255)

                   aaCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGa                                                                      >  1:2200223/1‑93 (MQ=255)

                   aaCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGaa                                                                     <  1:4489736/93‑1 (MQ=255)

                    cccGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGaa                                                                     <  1:6235643/92‑1 (MQ=255)

                    aCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGaa                                                                     <  1:2045891/93‑1 (MQ=255)

                    aCCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAAc                                                                    <  1:92773/93‑1 (MQ=255)

                       gcgATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGc                                                                  >  2:7206967/1‑93 (MQ=255)

                       gcgATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGc                                                                  >  1:6452951/1‑93 (MQ=255)

                        cgATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGcc                                                                 <  2:5446155/93‑1 (MQ=255)

                        cgATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGcc                                                                 <  1:4130455/93‑1 (MQ=255)

                         gATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGccc                                                                >  1:96173/1‑93 (MQ=255)

                         gATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGccc                                                                >  2:3378593/1‑93 (MQ=255)

                            aaTCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTc                                                             <  1:5698792/93‑1 (MQ=255)

                             aTCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCa                                                            <  2:3485731/93‑1 (MQ=255)

                                aaaGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGGATCGGGCAGCAGAACGCCCGTCAGTTc                                                        >  1:1006386/1‑93 (MQ=255)

                                 aaGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGCCCGTCAGTTc                                                        >  1:6834058/1‑93 (MQ=255)

                                   gACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCAt                                                      <  2:5213292/93‑1 (MQ=255)

                                   gACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATa                                                     <  2:846965/93‑1 (MQ=255)

                                      cgacgaTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATaga                                                   >  2:3584693/1‑93 (MQ=255)

                                       gacgaTTTTTCAGTGTCGATTTTCGCTACCAGCCAACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATagag                                                  >  2:4674011/1‑93 (MQ=255)

                                        acgaTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGGAGCAGAAAGCCCGTCAGTTCATagaggt                                                >  1:4842156/1‑91 (MQ=255)

                                         cgaTTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGAt                                                <  1:4273236/93‑1 (MQ=255)

                                            tttttCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCAc                                             <  2:3016677/93‑1 (MQ=255)

                                            tttttCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCAc                                             >  1:2381333/1‑93 (MQ=255)

                                             ttttCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCAcc                                            >  1:7118366/1‑93 (MQ=255)

                                                 cAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGt                                        <  2:2299188/93‑1 (MQ=255)

                                                  aGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGTCCGTCAGTTCATAGAGATCACCCAGTg                                       >  1:6771157/1‑93 (MQ=255)

                                                   gtgtCGATTTTCGCTACCAGCCGGCCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGa                                      <  1:4923226/93‑1 (MQ=255)

                                                   gtgtCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGa                                      >  1:6533286/1‑93 (MQ=255)

                                                    tgtCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAAt                                    >  2:3895053/1‑93 (MQ=255)

                                                     gtCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGGAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAAt                                    >  1:5351751/1‑93 (MQ=255)

                                                      tCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATc                                   <  1:1921955/93‑1 (MQ=255)

                                                      tCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATc                                   <  2:7054648/93‑1 (MQ=255)

                                                       cGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCg                                  <  2:2649316/93‑1 (MQ=255)

                                                          ttttCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCa                               >  2:6789768/1‑93 (MQ=255)

                                                           tttCGCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAg                              >  1:3579547/1‑93 (MQ=255)

                                                               gCTACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGaaaa                          >  1:3532351/1‑93 (MQ=255)

                                                                cTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAAt                         >  2:5791211/1‑93 (MQ=255)

                                                                cTACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAAt                         <  1:5578782/93‑1 (MQ=255)

                                                                 tACCAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAAtt                        >  1:3687167/1‑93 (MQ=255)

                                                                 tACCAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAAtt                        >  2:3860722/1‑93 (MQ=255)

                                                                   ccAGCCGACCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTa                      >  1:4406244/1‑93 (MQ=255)

                                                                    cAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCTCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTaa                     <  1:3737935/93‑1 (MQ=255)

                                                                    cAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTaa                     <  1:6032844/93‑1 (MQ=255)

                                                                    cAGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTaa                     <  1:820730/93‑1 (MQ=255)

                                                                     aGCCGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAg                    <  2:5391005/93‑1 (MQ=255)

                                                                      gCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGcc                  <  2:2393025/93‑1 (MQ=255)

                                                                        cGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCg                 >  2:181783/1‑93 (MQ=255)

                                                                        cGACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCg                 <  1:547661/93‑1 (MQ=255)

                                                                         gACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGt                <  2:2798993/93‑1 (MQ=255)

                                                                         gACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGt                >  2:2857772/1‑93 (MQ=255)

                                                                         gACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGt                >  2:3225987/1‑93 (MQ=255)

                                                                         gACCATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGt                >  2:2906076/1‑93 (MQ=255)

                                                                          aCCATCCGTTTCAAACGGCGGGGGGATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAAATTAAGCCGTa               >  1:5742792/1‑93 (MQ=255)

                                                                           ccATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaa              >  1:4250024/1‑93 (MQ=255)

                                                                           ccATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaaa             <  2:161808/93‑1 (MQ=255)

                                                                            cATCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTaaa             >  2:4523131/1‑93 (MQ=255)

                                                                            cATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAAt            <  2:5027231/93‑1 (MQ=255)

                                                                             aTCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATc           <  1:5112802/93‑1 (MQ=255)

                                                                             aTCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATc           <  2:6103561/93‑1 (MQ=255)

                                                                              tCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATc           >  1:5933937/1‑93 (MQ=255)

                                                                              tCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATc           >  1:284424/1‑93 (MQ=255)

                                                                              tCCGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATc           <  1:4560553/93‑1 (MQ=255)

                                                                                cGTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGa         <  2:4331763/93‑1 (MQ=255)

                                                                                 gTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGac        <  2:5060895/93‑1 (MQ=255)

                                                                                 gTTTCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGac        <  1:2539068/93‑1 (MQ=255)

                                                                                   ttCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGacac      <  2:661014/93‑1 (MQ=255)

                                                                                   ttCAAACGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGacac      <  2:1302111/93‑1 (MQ=255)

                                                                                     cAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGACACCGc   <  2:1140068/93‑1 (MQ=255)

                                                                                       aaCGGCGGGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGACACCGcc  >  2:1476292/1‑93 (MQ=255)

TCACCAATTCGCCCATTAAAACCGCGATAATCAAAGACCGACGATTTTTCAGTGTCGATTTTCGCTACCAGCCGACCATCCGTTTCAAACGGC‑GGGGTATCGGGCAGCAGAACGCCCGTCAGTTCATAGAGATCACCCAGTGAATCGCCAGAAAATTTAAGCCGTAAATCGACACCGCC  >  NC_000913/3674980‑3675158

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