BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,828,674	T→C	T197T (ACT→ACC)	yicH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,828,674	0	T	C	99.3%	1383.1 / ‑2.4	436	T197T (ACT→ACC)	yicH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (257/176); total (259/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 450 total aligned reads displayed.

GTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCAT  >  NC_000913/3828584‑3828763

gTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCaa                                                                                         >  2:1467535/1‑93 (MQ=255)

gTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCaa                                                                                         >  2:413544/1‑93 (MQ=255)

 tttAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAAc                                                                                        <  2:6346405/93‑1 (MQ=255)

   tAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAAAGATCGCGTTACGCTGACCAACCt                                                                                      >  1:6376614/1‑93 (MQ=255)

    aGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTg                                                                                     <  1:952024/93‑1 (MQ=255)

      tGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTggg                                                                                   <  1:2935668/93‑1 (MQ=255)

         cGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCACCCTGGGTGc                                                                                >  1:2171936/1‑93 (MQ=255)

         cGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGc                                                                                <  1:5128904/93‑1 (MQ=255)

          ggATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGcc                                                                               >  1:1833273/1‑93 (MQ=255)

          ggATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGcc                                                                               >  1:7285108/1‑93 (MQ=255)

          ggATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGcc                                                                               >  2:1586408/1‑93 (MQ=255)

          ggATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGcc                                                                               >  2:4551096/1‑93 (MQ=255)

           gATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCg                                                                              >  2:2632762/1‑93 (MQ=255)

           gATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACCATCGCGTTACGCTGACCAACCTGGGTGCCg                                                                              >  2:2127103/1‑93 (MQ=255)

            aTCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGa                                                                             >  1:5918844/1‑93 (MQ=255)

            aTCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGa                                                                             <  2:5738658/93‑1 (MQ=255)

             tCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGAc                                                                            >  1:1671294/1‑93 (MQ=255)

             tCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGAc                                                                            >  1:4479147/1‑93 (MQ=255)

                cTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATc                                                                         <  1:5477158/93‑1 (MQ=255)

                   tCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGcc                                                                      <  2:3185755/93‑1 (MQ=255)

                      cTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCcgc                                                                   <  1:6124600/93‑1 (MQ=255)

                       tCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCcgcg                                                                  >  1:4569183/1‑93 (MQ=255)

                       tCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCcgcg                                                                  >  2:3890010/1‑93 (MQ=255)

                         aaTGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCggg                                                                <  2:2240786/93‑1 (MQ=255)

                         aaTGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCggg                                                                <  1:1986180/93‑1 (MQ=255)

                         aaTGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCggg                                                                >  1:5300474/1‑93 (MQ=255)

                         aaTGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCggg                                                                >  2:2576196/1‑93 (MQ=255)

                           tgatgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGAc                                                              >  2:2882447/1‑93 (MQ=255)

                           tgatgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGAc                                                              >  1:1495122/1‑93 (MQ=255)

                           tgatgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGAc                                                              >  1:5937478/1‑93 (MQ=255)

                           taatgTACCTACCACCATTGTACTGATTGAAGGCAATATTAATAACGAAAACGATCCGCTGACCAACCAGGGTGCAAACATCGCCCGCGGGAc                                                              <  2:7009207/91‑1 (MQ=255)

                            gatgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACa                                                             <  1:2833163/93‑1 (MQ=255)

                            gatgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACa                                                             >  1:6052567/1‑93 (MQ=255)

                             atgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACAt                                                            >  1:3651766/1‑93 (MQ=255)

                             atgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACAt                                                            <  1:6295173/93‑1 (MQ=255)

                              tgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACAtt                                                           <  1:4270660/93‑1 (MQ=255)

                              tgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACAtt                                                           >  1:353892/1‑93 (MQ=255)

                              tgTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGAAAtt                                                           >  2:3028167/1‑93 (MQ=255)

                                ttCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTac                                                         >  2:6619631/1‑92 (MQ=255)

                                ttCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTaa                                                         <  1:4210050/93‑1 (MQ=255)

                                 tCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAAc                                                        >  1:5678303/1‑93 (MQ=255)

                                 tCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAAc                                                        <  2:7494323/93‑1 (MQ=255)

                                 tCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAAc                                                        <  2:6191616/93‑1 (MQ=255)

                                 tCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAAc                                                        <  2:597465/93‑1 (MQ=255)

                                 tCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAAc                                                        >  2:1916974/1‑93 (MQ=255)

                                   cTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCg                                                      <  2:5514569/93‑1 (MQ=255)

                                   cTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCg                                                      >  2:1489729/1‑93 (MQ=255)

                                     gCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGa                                                    >  2:523016/1‑93 (MQ=255)

                                       caccaATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGaaa                                                  <  1:3931835/93‑1 (MQ=255)

                                       caccaATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGaaa                                                  <  1:7116386/93‑1 (MQ=255)

                                        accaATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAAc                                                 >  1:925503/1‑93 (MQ=255)

                                          caATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATTGCCCGCGGGACATTAACCGGAAAcgc                                               >  2:7164905/1‑93 (MQ=255)

                                          caATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAAcgc                                               >  2:3809688/1‑93 (MQ=255)

                                           aaTGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAAcgcg                                              >  2:1070255/1‑93 (MQ=255)

                                             tGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCa                                            >  2:2367649/1‑93 (MQ=255)

                                             tGTACTGATTGAAGGCAGCATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCa                                            >  2:7835803/1‑93 (MQ=255)

                                              gTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAg                                           >  1:6930737/1‑93 (MQ=255)

                                                  tGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTa                                       >  1:7610516/1‑93 (MQ=255)

                                                   gATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTaa                                      >  2:5102816/1‑93 (MQ=255)

                                                   gATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTaa                                      >  1:2952275/1‑93 (MQ=255)

                                                    aTTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAAc                                     <  1:7871104/93‑1 (MQ=255)

                                                      tGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGc                                   >  1:1031303/1‑93 (MQ=255)

                                                      tGAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGc                                   <  1:3268047/93‑1 (MQ=255)

                                                       gAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGcc                                  >  2:1137775/1‑93 (MQ=255)

                                                       gAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGcc                                  <  1:5519130/93‑1 (MQ=255)

                                                       gAAGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGcc                                  >  1:4911110/1‑93 (MQ=255)

                                                         aGGCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGa                                <  1:1540793/93‑1 (MQ=255)

                                                          ggCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGAc                               >  2:480449/1‑93 (MQ=255)

                                                           gCAGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACg                              >  1:1898180/1‑93 (MQ=255)

                                                             aGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGc                            >  1:7432804/1‑93 (MQ=255)

                                                             aGTATTGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGc                            >  1:1284317/1‑93 (MQ=255)

                                                               tATTGATAACGATCGCGTTACGCTGACCAACCCGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCac                          >  1:6450971/1‑92 (MQ=255)

                                                                  tGATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTg                       >  1:4636897/1‑93 (MQ=255)

                                                                   gATAACGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTgg                      <  2:6522378/93‑1 (MQ=255)

                                                                      aaCGATCGCGTTACGCTGACCAAGCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCaa                   <  2:4347717/93‑1 (MQ=255)

                                                                      aaCGATCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCaa                   <  1:3349089/93‑1 (MQ=255)

                                                                          aTCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTgg               <  1:7930275/93‑1 (MQ=255)

                                                                          aTCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTgg               >  2:5352169/1‑93 (MQ=255)

                                                                           tCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGa              >  2:3637153/1‑93 (MQ=255)

                                                                           tCGCGTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGa              >  1:4061641/1‑93 (MQ=255)

                                                                            cgcgTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGaa             >  2:7279351/1‑93 (MQ=255)

                                                                            cgcgTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGaa             <  1:3131837/93‑1 (MQ=255)

                                                                             gcgTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGaaa            >  1:1209218/1‑93 (MQ=255)

                                                                               gTTACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAAc          >  2:4196766/1‑93 (MQ=255)

                                                                                ttACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAAcc         <  1:6650637/93‑1 (MQ=255)

                                                                                ttACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAAcc         <  2:757357/93‑1 (MQ=255)

                                                                                 tACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCt        <  1:3156917/93‑1 (MQ=255)

                                                                                 tACGCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCt        <  2:352745/93‑1 (MQ=255)

                                                                                    gCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAAccgcgc     >  2:7705506/1‑89 (MQ=255)

                                                                                    gCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcg     >  1:477723/1‑93 (MQ=255)

                                                                                    gCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcg     <  2:4911935/93‑1 (MQ=255)

                                                                                    gCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcg     <  2:2703037/93‑1 (MQ=255)

                                                                                    gCTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcg     >  1:1855587/1‑93 (MQ=255)

                                                                                     cTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGACACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTCGAAACCCTgcgc    >  2:6890439/1‑93 (MQ=255)

                                                                                     cTGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTgcgc    >  1:684572/1‑93 (MQ=255)

                                                                                      tGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCa   <  1:4842636/93‑1 (MQ=255)

                                                                                      tGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCa   <  2:3393074/93‑1 (MQ=255)

                                                                                      tGACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCa   <  1:2333567/93‑1 (MQ=255)

                                                                                       gACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCAt  <  1:2579229/93‑1 (MQ=255)

                                                                                       gACCAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCAt  <  1:5784434/93‑1 (MQ=255)

GTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCCCGCGGGACATTAACCGGAAACGCGCAGCGTAACGCCGACGGCAGCTGGCAAGTGGAAAACCTGCGCAT  >  NC_000913/3828584‑3828763

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