BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,739,740	A→G	M141T (ATG→ACG)	yiaJ ←	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,739,740	0	A	G	99.8%	1385.0 / NA	441	M141T (ATG→ACG)	yiaJ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (222/218); minor base C (0/1); total (222/219)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 459 total aligned reads displayed.

TAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTC  >  NC_000913/3739650‑3739829

   tCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCg                                                                                      <  2:5810852/93‑1 (MQ=255)

   tCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCg                                                                                      >  2:4494078/1‑93 (MQ=255)

    ccGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCgg                                                                                     <  2:4968770/93‑1 (MQ=255)

      gggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGtt                                                                                   >  1:7527373/1‑93 (MQ=255)

      gggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGtt                                                                                   <  2:5688545/93‑1 (MQ=255)

        gTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgt                                                                                 >  1:4923031/1‑93 (MQ=255)

         tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                >  1:1044255/1‑93 (MQ=255)

         tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                >  2:436035/1‑93 (MQ=255)

         tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                <  2:4664942/93‑1 (MQ=255)

         tGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgtg                                                                                <  2:476288/93‑1 (MQ=255)

          gACCAAACGCGATGTAGATCTTGCCCATGGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTgg                                                                               <  2:229793/93‑1 (MQ=255)

          gACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTgg                                                                               <  1:4413428/93‑1 (MQ=255)

             cAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGtt                                                                            >  2:2911865/1‑93 (MQ=255)

              aaaCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTc                                                                           <  2:3102774/93‑1 (MQ=255)

               aaCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTcc                                                                          <  2:2093980/93‑1 (MQ=255)

                aCGCCATGTAGATCTTGCCCATTGCGGAGCAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCa                                                                         <  1:2042794/93‑1 (MQ=255)

                aCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCa                                                                         <  2:5591906/93‑1 (MQ=255)

                aCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCa                                                                         <  2:6616378/93‑1 (MQ=255)

                aCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCa                                                                         <  2:1179328/93‑1 (MQ=255)

                  gCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGc                                                                       >  2:331107/1‑93 (MQ=255)

                    cATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCtt                                                                     >  1:2903779/1‑93 (MQ=255)

                     aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    >  1:56597/1‑93 (MQ=255)

                     aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    >  1:7444078/1‑93 (MQ=255)

                     aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGAACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    >  2:5255362/1‑93 (MQ=255)

                      tGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtat                                                                   <  2:2387195/93‑1 (MQ=255)

                      tGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtat                                                                   <  1:7856380/93‑1 (MQ=255)

                      tGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtat                                                                   >  1:961953/1‑93 (MQ=255)

                      tGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTaat                                                                   >  1:7761425/1‑93 (MQ=255)

                       gTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtata                                                                  >  1:4641106/1‑93 (MQ=255)

                        tAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATGTAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaa                                                                 >  1:860882/1‑93 (MQ=255)

                        tAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaa                                                                 >  2:6813647/1‑93 (MQ=255)

                        tAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATaa                                                                 >  1:6106311/1‑93 (MQ=255)

                          gATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAAt                                                               >  1:6423257/1‑93 (MQ=255)

                           aTCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATc                                                              >  2:3023862/1‑93 (MQ=255)

                           aTCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATc                                                              <  2:2658464/93‑1 (MQ=255)

                             cTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCGCGGGTTCGCAGCGTCCGGGGTGTGGGTTCCAGCTTATAAATCaa                                                            <  2:7238723/93‑1 (MQ=255)

                              ttGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaaa                                                           <  1:3283338/93‑1 (MQ=255)

                              ttGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaaa                                                           >  1:3062907/1‑93 (MQ=255)

                                gCCCATGGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAAt                                                         <  2:4155539/93‑1 (MQ=255)

                                 cccATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATa                                                        <  1:2759119/93‑1 (MQ=255)

                                   cATTGCGGAACAGTAGAGCGGNATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      <  1:1129909/93‑1 (MQ=255)

                                   cATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      >  2:499348/1‑93 (MQ=255)

                                   cATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      >  1:4730151/1‑93 (MQ=255)

                                   cATTGCGGAACAGTAGAGCGGCATATGCTCGCCAATATAGGCATGGGTTCGCAGCCTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      <  1:5222325/93‑1 (MQ=255)

                                    aTTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCg                                                     <  1:5777908/93‑1 (MQ=255)

                                     ttGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGt                                                    <  2:1928470/93‑1 (MQ=255)

                                       gCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGa                                                  <  2:6002121/93‑1 (MQ=255)

                                        acgAACAGTAGAGGCGGATATTCGGGCCAACATTGGCACGGGTTCGCACCGCCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGAt                                                 <  1:2514643/91‑1 (MQ=255)

                                         ggAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATc                                                >  2:619361/1‑93 (MQ=255)

                                         ggAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATc                                                <  2:6713249/93‑1 (MQ=255)

                                         ggAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATc                                                >  1:4792685/1‑93 (MQ=255)

                                         ggAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATc                                                >  1:3565106/1‑93 (MQ=255)

                                          gAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCg                                               <  1:1588121/93‑1 (MQ=255)

                                          gAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCg                                               >  2:1631031/1‑93 (MQ=255)

                                          gAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCg                                               >  1:3151819/1‑93 (MQ=255)

                                             cAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCt                                            <  1:1781544/93‑1 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          <  1:349207/93‑1 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          <  2:6761152/93‑1 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          >  1:6899564/1‑93 (MQ=255)

                                                tAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcg                                         <  1:3924928/93‑1 (MQ=255)

                                                 agagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcgc                                        <  2:7629041/93‑1 (MQ=255)

                                                 agagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcgc                                        >  1:6221694/1‑93 (MQ=255)

                                                  gagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTtgcg                                       >  2:5378784/1‑93 (MQ=255)

                                                  gagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcgcg                                       >  2:1447508/1‑93 (MQ=255)

                                                    gCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGc                                     >  2:6325182/1‑93 (MQ=255)

                                                      ggCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTg                                   <  2:1352412/93‑1 (MQ=255)

                                                       gCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTgg                                  >  2:7404592/1‑93 (MQ=255)

                                                        cATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGga                                 >  2:6560054/1‑93 (MQ=255)

                                                        cATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGga                                 <  2:3711759/93‑1 (MQ=255)

                                                         atatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgag                                >  1:124806/1‑93 (MQ=255)

                                                          tatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgaga                               >  1:2414760/1‑93 (MQ=255)

                                                          tatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgaga                               >  1:4335121/1‑93 (MQ=255)

                                                          tatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgaga                               >  2:3575158/1‑93 (MQ=255)

                                                                ggCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTaa                         <  2:1042897/93‑1 (MQ=255)

                                                                gcccATTATTGGCACGGTTTCGCGGCGTCCCGTTTGTGGGTCCCGGC‑TATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        <  2:642929/91‑1 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        >  2:5142618/1‑93 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        <  1:6553933/93‑1 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        <  1:3717900/93‑1 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        <  2:713557/93‑1 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        >  1:1005327/1‑93 (MQ=255)

                                                                  ccAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATg                       <  1:7416083/93‑1 (MQ=255)

                                                                    aaTATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGt                     <  1:5430777/93‑1 (MQ=255)

                                                                    aaTATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGt                     >  1:6471139/1‑93 (MQ=255)

                                                                      tataGGCACGGGTTCGCAGCGTCCCGGTTGCGGGGTCCAGCTTATAAATCAAAATAGCGCGATCGTCTTCACGCGTGGGGTAGTTAATGGttt                   >  1:5943256/1‑93 (MQ=255)

                                                                        taGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCa                 <  1:2243897/93‑1 (MQ=255)

                                                                             aCGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCatt            <  2:4234789/93‑3 (MQ=255)

                                                                             aCGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGt            <  2:961650/93‑1 (MQ=255)

                                                                             aCGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGt            <  1:5155562/93‑1 (MQ=255)

                                                                              cGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTg           >  2:1557764/1‑93 (MQ=255)

                                                                              cGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTg           >  2:846662/1‑93 (MQ=255)

                                                                               gggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGGTAATGGTTTCACCAGTgg          >  2:2059761/1‑93 (MQ=255)

                                                                                ggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCCGTggg         >  2:6978479/1‑92 (MQ=255)

                                                                                ggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGc         <  2:7862983/93‑1 (MQ=255)

                                                                                ggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGGTTCACCAGTGGc         >  1:1463681/1‑93 (MQ=255)

                                                                                ggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGGTAATTGTTTCACCAggggg         >  1:4171424/1‑92 (MQ=255)

                                                                                 gTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCg        >  1:6799665/1‑93 (MQ=255)

                                                                                  ttCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGa       <  2:3382503/93‑1 (MQ=255)

                                                                                   tCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGAt      <  2:2853800/93‑1 (MQ=255)

                                                                                     gcagcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGt    >  1:7155271/1‑93 (MQ=255)

                                                                                       agcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTc  >  1:53630/1‑93 (MQ=255)

TAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTC  >  NC_000913/3739650‑3739829

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