BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,211,571	A→G	R44R (CGT→CGC)	yjaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,211,571	0	A	G	97.2%	534.4 / 2.6	176	R44R (CGT→CGC)	yjaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/5); new base G (105/66); total (105/71)
Fisher's exact test for biased strand distribution p-value = 9.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 207 total aligned reads displayed.

TCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTC  >  NC_000913/4211482‑4211660

  cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaac                                                                                      >  2:3741613/1‑93 (MQ=255)

    tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTa                                                                                    <  2:3247492/93‑1 (MQ=255)

    tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTa                                                                                    >  2:1908642/1‑93 (MQ=255)

    tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTa                                                                                    >  2:795508/1‑93 (MQ=255)

    tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAAATa                                                                                    >  2:1908644/1‑93 (MQ=255)

    tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTa                                                                                    <  2:1580220/93‑1 (MQ=255)

      tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGTGCTTCCGGCAGGAAGGAGCGAACTAgg                                                                                  <  2:3103208/93‑1 (MQ=255)

      tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAgg                                                                                  <  2:2526246/93‑1 (MQ=255)

      tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAgg                                                                                  >  2:4651182/1‑93 (MQ=255)

       tctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGt                                                                                 <  1:2336070/92‑1 (MQ=255)

         tgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcc                                                                               <  1:1508386/93‑1 (MQ=255)

           cccACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcctc                                                                             >  2:363253/1‑93 (MQ=255)

           cccACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcctc                                                                             >  1:3466644/1‑93 (MQ=255)

            ccACTTAGCAACATAACTCCAACGGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcctcc                                                                            <  2:4499568/93‑1 (MQ=255)

            ccACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcctcc                                                                            <  1:1617592/93‑1 (MQ=255)

            ccACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcctcc                                                                            >  1:2137918/1‑93 (MQ=255)

             cACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCa                                                                           >  1:3591334/1‑93 (MQ=255)

             cACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCa                                                                           <  1:1479671/93‑1 (MQ=255)

               cTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGc                                                                         >  2:264142/1‑93 (MQ=255)

               acgcgCAACATACATCTAACCGGCTGATCCCGCTCATTAACCTCCATCCACAACGGAGCGTCGGCCAGGAAGGAGCGAATTAGGTCCTCCAGc                                                                         <  2:2519986/89‑1 (MQ=255)

                ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCt                                                                        >  1:3359427/1‑93 (MQ=255)

                ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCt                                                                        <  1:4198593/93‑1 (MQ=255)

                ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCt                                                                        <  1:1660241/93‑1 (MQ=255)

                ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCt                                                                        >  1:2020594/1‑93 (MQ=255)

                 tAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTc                                                                       <  2:4299072/93‑1 (MQ=255)

                  aGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCg                                                                      >  1:4484411/1‑93 (MQ=255)

                   gCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCgg                                                                     >  2:1747676/1‑93 (MQ=255)

                   gCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCgg                                                                     >  1:1333032/1‑93 (MQ=255)

                     aaCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTc                                                                   <  2:4251334/93‑1 (MQ=255)

                      aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTcc                                                                  >  2:1518310/1‑93 (MQ=255)

                       cATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCg                                                                 >  1:260503/1‑93 (MQ=255)

                       cATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCg                                                                 >  1:4264581/1‑93 (MQ=255)

                          aaaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATa                                                              <  1:586580/93‑1 (MQ=255)

                          aaaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATa                                                              >  1:1247168/1‑93 (MQ=255)

                           aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATAc                                                             >  1:3165073/1‑93 (MQ=255)

                            aTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATAct                                                            <  2:2233092/93‑1 (MQ=255)

                            aTCCAACCCGCTGATCCCGCTCATTAACCGCGACCCACAACGCCGCTTCCGGCAGAAAGAAACGAATTAGGTCCTCCAGCTCGGTCCGATAct                                                            <  2:1980299/93‑1 (MQ=255)

                             tCCAACCGGCTGATCCCGCTTATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATActc                                                           >  2:1431006/1‑93 (MQ=255)

                             tCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATActc                                                           >  1:3009593/1‑93 (MQ=255)

                              ccACCCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATActct                                                          <  2:2395033/93‑1 (MQ=255)

                              ccAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATActct                                                          >  1:2919267/1‑93 (MQ=255)

                                aaCCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctgc                                                        >  1:1891494/1‑93 (MQ=255)

                                 aCCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctgct                                                       >  1:2072527/1‑93 (MQ=255)

                                 aCCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctgct                                                       <  1:4024904/93‑1 (MQ=255)

                                    ggCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGCCCGATACCCTGCTGAt                                                    >  2:2200622/1‑93 (MQ=255)

                                      cTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAg                                                  >  1:1695107/1‑93 (MQ=255)

                                      cTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAg                                                  >  2:652865/1‑93 (MQ=255)

                                      cTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGGCCGATACTCTGCTGATAg                                                  >  1:1177289/1‑93 (MQ=255)

                                       tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCGGGAAGGAGTGAACTAGGTCCTCGAGCTCGGTCAGATACTCTGATGATAGa                                                 >  2:1315529/1‑93 (MQ=255)

                                       tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGACCTCCAGCTCGGTCCGATACTCTGCTGATAGa                                                 >  1:3291463/1‑93 (MQ=255)

                                        gATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaa                                                >  2:3417197/1‑93 (MQ=255)

                                        gATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaa                                                >  2:1334174/1‑93 (MQ=255)

                                        gATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaa                                                >  1:245364/1‑93 (MQ=255)

                                          tCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaaaa                                              >  1:1048094/1‑93 (MQ=255)

                                           cccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAAt                                             <  2:826686/93‑1 (MQ=255)

                                           cccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAAt                                             <  1:3565792/93‑1 (MQ=255)

                                            ccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATc                                            >  1:3857306/1‑93 (MQ=255)

                                            ccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATc                                            >  1:12257/1‑93 (MQ=255)

                                            ccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATc                                            <  2:4667711/93‑1 (MQ=255)

                                             cGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCg                                           >  2:517330/1‑93 (MQ=255)

                                             cGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCg                                           <  2:1171357/93‑1 (MQ=255)

                                             cGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCg                                           >  2:2159096/1‑93 (MQ=255)

                                              gCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGt                                          >  1:2488517/1‑93 (MQ=255)

                                              cgtCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGt                                          <  2:1893130/91‑1 (MQ=255)

                                               cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTg                                         >  2:2868162/1‑93 (MQ=255)

                                                tCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGa                                        <  2:613059/93‑1 (MQ=255)

                                                tCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGa                                        >  2:1042409/1‑93 (MQ=255)

                                                 cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       <  2:3558999/93‑1 (MQ=255)

                                                 cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       <  1:3920588/93‑1 (MQ=255)

                                                  aTTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGt                                      >  1:677879/1‑93 (MQ=255)

                                                  aTTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGt                                      <  2:1248436/93‑1 (MQ=255)

                                                   ttAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTg                                     >  1:2183605/1‑93 (MQ=255)

                                                    tAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTgg                                    >  1:2517230/1‑93 (MQ=255)

                                                      aCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCa                                  >  2:3893749/1‑93 (MQ=255)

                                                         gcgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCg                               >  2:397399/1‑93 (MQ=255)

                                                          cgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGa                              <  2:1835469/93‑1 (MQ=255)

                                                           gACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGAc                             >  1:2562377/1‑93 (MQ=255)

                                                               cacaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGaa                         <  1:779617/93‑1 (MQ=255)

                                                                acaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAAc                        >  2:3313737/1‑93 (MQ=255)

                                                                acaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAAc                        >  1:1580855/1‑93 (MQ=255)

                                                                 caACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTTCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACg                       >  2:4131506/1‑93 (MQ=255)

                                                                   aCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGac                     <  1:3341329/93‑1 (MQ=255)

                                                                   aCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGac                     <  1:172502/93‑1 (MQ=255)

                                                                         cTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCaaa               >  1:3689746/1‑93 (MQ=255)

                                                                          ttCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAAt              >  2:3188478/1‑93 (MQ=255)

                                                                          ttCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGGGAGTGGCATCGACAGAACGACACCAAAt              >  1:3485811/1‑93 (MQ=255)

                                                                             cGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGc           >  1:4603404/1‑93 (MQ=255)

                                                                              ggCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCa          <  1:1963879/93‑1 (MQ=255)

                                                                               gCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa         <  2:588232/93‑1 (MQ=255)

                                                                               gCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa         <  2:2109122/93‑1 (MQ=255)

                                                                               gCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa         <  1:1068967/93‑1 (MQ=255)

                                                                                gaggaaggaGCGAAATAGGTCCTCCAGCTCGGTCCGATTCTTTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAAc        <  1:2448209/92‑1 (MQ=255)

                                                                                cAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAAc        >  1:4634800/1‑93 (MQ=255)

                                                                                  ggaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGa      >  1:2289702/1‑93 (MQ=255)

                                                                                  ggaaggagCGAACTAGGTCATCCAGCTCGGTCCGACACTCTGCTGATAGAAAATCGTGAGAGGCATAGACAGAAGGACACCAAATCGACACGa      >  1:710141/1‑93 (MQ=255)

                                                                                   gaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAg     >  1:3047844/1‑93 (MQ=255)

                                                                                   gaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAg     <  1:3246297/93‑1 (MQ=255)

                                                                                    aaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGt    <  2:4164549/93‑1 (MQ=255)

                                                                                      ggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTTTCCTGTTAGAAAATCGTGGGGGGCTTGGACAGAACGAAACCAAATCGCAACGAGTtc  >  2:269385/1‑93 (MQ=255)

                                                                                      ggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTtc  >  1:1951127/1‑93 (MQ=255)

TCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTC  >  NC_000913/4211482‑4211660

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