BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,793,926	G→A	A107A (GCG→GCA)	rfaC →	heptosyl transferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,793,926	0	G	A	83.3%	635.9 / 109.8	312	A107A (GCG→GCA)	rfaC	heptosyl transferase I
Reads supporting (aligned to +/- strand): ref base G (26/26); new base A (138/122); total (164/148)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-08

Only 100 of 318 total aligned reads displayed.

GTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCA  >  NC_000913/3793836‑3794018

gTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGcaca                                                                                            >  2:2362445/1‑93 (MQ=255)

gTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGcaca                                                                                            >  1:6404918/1‑93 (MQ=255)

  ttCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGTGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATg                                                                                          <  1:6951896/93‑1 (MQ=255)

  ttCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATg                                                                                          >  2:2426463/1‑93 (MQ=255)

   tCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATgg                                                                                         <  1:3096551/93‑1 (MQ=255)

    cGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGc                                                                                        <  1:1371630/93‑1 (MQ=255)

       gAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTa                                                                                     >  1:6624042/1‑93 (MQ=255)

       gAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGAGCTGGTGACGCGTCTGGCGCATGGCGTa                                                                                     >  2:6194552/1‑93 (MQ=255)

        aaGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTaa                                                                                    <  1:5532038/93‑1 (MQ=255)

         aGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTaaa                                                                                   >  1:284928/1‑93 (MQ=255)

          gcgcTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAg                                                                                  >  1:5626180/1‑93 (MQ=255)

          gcgcTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAg                                                                                  >  2:1708557/1‑93 (MQ=255)

            gcTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgca                                                                                >  1:4405766/1‑93 (MQ=255)

             cTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcat                                                                               <  2:6401977/93‑1 (MQ=255)

             cTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcat                                                                               <  2:5391762/93‑1 (MQ=255)

              tACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcatg                                                                              >  1:5680266/1‑93 (MQ=255)

              tACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcatg                                                                              <  1:4161647/93‑1 (MQ=255)

              tACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcatg                                                                              >  1:4012685/1‑93 (MQ=255)

               aCAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcatgg                                                                             >  2:6534334/1‑93 (MQ=255)

                cAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcatggc                                                                            <  2:7272731/93‑1 (MQ=255)

                 aaGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcatggca                                                                           <  2:1093921/93‑1 (MQ=255)

                  aGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcatggcat                                                                          >  2:5789393/1‑93 (MQ=255)

                   gCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAgcatggcatg                                                                         >  1:2119590/1‑93 (MQ=255)

                   gCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAgcatggcatg                                                                         >  1:184385/1‑93 (MQ=255)

                    cAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGCAAAGcatggcatgg                                                                        <  2:5185639/93‑1 (MQ=255)

                      gagaACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGAc                                                                      <  1:4261437/93‑1 (MQ=255)

                       agaACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACt                                                                     <  1:3334158/93‑1 (MQ=255)

                       agaACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACt                                                                     >  2:5504674/1‑93 (MQ=255)

                       agaACTATGACGCAGTTATCGACGCTCAGGGGCTGGAAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACt                                                                     >  1:6331010/1‑93 (MQ=255)

                        gaACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTg                                                                    >  2:251624/1‑93 (MQ=255)

                          aCTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGc                                                                  <  2:4955127/93‑1 (MQ=255)

                           cTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCa                                                                 >  1:6244392/1‑93 (MQ=255)

                            tATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCaa                                                                >  2:4618317/1‑93 (MQ=255)

                            tATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCaa                                                                >  1:7237633/1‑93 (MQ=255)

                            tATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCaa                                                                >  2:3423553/1‑93 (MQ=255)

                              tGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAAc                                                              >  1:2695361/1‑93 (MQ=255)

                               gACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAAcc                                                             >  1:2629908/1‑93 (MQ=255)

                               gACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAAcc                                                             >  2:7349229/1‑93 (MQ=255)

                                 cGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGc                                                           <  1:7413394/93‑1 (MQ=255)

                                 cGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGc                                                           <  2:4360797/93‑1 (MQ=255)

                                  gCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCt                                                          <  2:3678622/93‑1 (MQ=255)

                                  gCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCt                                                          <  2:3903630/93‑1 (MQ=255)

                                   cAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTc                                                         <  1:4760603/93‑1 (MQ=255)

                                    aGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTcg                                                        >  1:3547351/1‑93 (MQ=255)

                                     gTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTcgc                                                       >  2:1479560/1‑93 (MQ=255)

                                      ttATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTcgcg                                                      >  1:6110483/1‑93 (MQ=255)

                                       tATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGa                                                     >  2:4885777/1‑93 (MQ=255)

                                        aTCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGaa                                                    <  2:6626346/93‑1 (MQ=255)

                                          cGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAAcc                                                  <  2:5235076/93‑1 (MQ=255)

                                           gACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCt                                                 <  2:3364301/93‑1 (MQ=255)

                                             gcttCAGGGGCTGGAAAAAACCCCGGCCCGGTTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCttt                                               <  1:1209693/90‑1 (MQ=255)

                                               cTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTag                                             <  2:6438487/93‑1 (MQ=255)

                                                tCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTagc                                            >  1:5360987/1‑93 (MQ=255)

                                                  aGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTagcca                                          >  1:2722462/1‑93 (MQ=255)

                                                   ggggCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTagccag                                         >  2:3594507/1‑93 (MQ=255)

                                                     ggCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTagccagcc                                       <  1:6686649/93‑1 (MQ=255)

                                                      gCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCt                                      >  2:139352/1‑93 (MQ=255)

                                                       cTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTg                                     >  2:3049618/1‑93 (MQ=255)

                                                        tGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGt                                    >  1:5095/1‑93 (MQ=255)

                                                         ggTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGtt                                   >  1:6061201/1‑93 (MQ=255)

                                                          gTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGttt                                  <  1:98332/93‑1 (MQ=255)

                                                           tAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGtttt                                 >  2:4849745/1‑93 (MQ=255)

                                                            aaaaaGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTa                                <  2:599414/93‑1 (MQ=255)

                                                             aaaaGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTAc                               <  1:719219/93‑1 (MQ=255)

                                                              aaaGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACa                              >  2:6681323/1‑93 (MQ=255)

                                                              aaaGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACa                              >  1:484716/1‑93 (MQ=255)

                                                              aaaGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACa                              >  1:1240526/1‑93 (MQ=255)

                                                                aGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAAt                            <  2:7051674/93‑1 (MQ=255)

                                                                aGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAAt                            >  1:1894768/1‑93 (MQ=255)

                                                                aGCGCGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAAt                            >  1:3424756/1‑93 (MQ=255)

                                                                 gcgcGGCGCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATc                           <  1:41459/93‑1 (MQ=255)

                                                                   gcggcgCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGt                         <  1:880036/93‑1 (MQ=255)

                                                                    cggcgCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTa                        >  1:2020177/1‑93 (MQ=255)

                                                                    cggcgCTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTa                        >  2:3803416/1‑93 (MQ=255)

                                                                     ggcgCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTaa                       <  1:588399/93‑1 (MQ=255)

                                                                      gcgcTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAg                      >  1:3226941/1‑93 (MQ=255)

                                                                      gcgcTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAg                      <  1:1030292/93‑1 (MQ=255)

                                                                       cgcTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGc                     <  2:2079680/93‑1 (MQ=255)

                                                                         cTGGTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGcat                   >  1:6503252/1‑93 (MQ=255)

                                                                           ggTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGcatca                 <  2:7529448/93‑1 (MQ=255)

                                                                           ggTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGcatca                 <  1:3900240/93‑1 (MQ=255)

                                                                           ggTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGcatca                 <  1:5234471/93‑1 (MQ=255)

                                                                           ggTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGcatca                 >  1:854874/1‑93 (MQ=255)

                                                                            gTGACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGcatcat                <  2:3215783/93‑1 (MQ=255)

                                                                              gACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCatat              <  2:1582601/93‑1 (MQ=255)

                                                                              gACGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCatat              >  1:3042210/1‑93 (MQ=255)

                                                                               aCGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATAtt             >  2:1227735/1‑93 (MQ=255)

                                                                               aCGCGTCTGGCACATGGCGTAAAGCATGGCATGGGCTGGCAAACCCCTCGCGAACCTTTAGCCAGCCCGTTTTACAATCGTAAACATCATAtt             >  1:3665594/1‑93 (MQ=255)

                                                                               aCGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATAtt             <  1:812422/93‑1 (MQ=255)

                                                                               aCGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATAtt             <  2:2758092/93‑1 (MQ=255)

                                                                               aCGCGTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATAtt             <  2:4088396/93‑1 (MQ=255)

                                                                                cgcgTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTg            >  1:2328045/1‑93 (MQ=255)

                                                                                  cgTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCa          <  2:1663817/93‑1 (MQ=255)

                                                                                  cgTCTGGCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCa          <  1:494285/93‑1 (MQ=255)

                                                                                       ggCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca     >  2:6297146/1‑93 (MQ=255)

                                                                                       ggCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca     <  2:5607927/93‑1 (MQ=255)

                                                                                       ggCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca     >  2:7194618/1‑93 (MQ=255)

                                                                                       ggCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca     <  1:1648195/93‑1 (MQ=255)

                                                                                       ggCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca     <  1:1148558/93‑1 (MQ=255)

                                                                                       ggCACATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAAca     >  1:1056773/1‑93 (MQ=255)

GTTTCGTGAAGCGCTACAAGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGGACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCA  >  NC_000913/3793836‑3794018

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