BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,724,581	T→C	Q137R (CAG→CGG)	yiaA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,724,581	0	T	C	99.5%	702.9 / ‑4.8	217	Q137R (CAG→CGG)	yiaA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (104/112); total (104/113)
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.62e-01

Only 100 of 229 total aligned reads displayed.

AGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCA  >  NC_000913/3724491‑3724670

 gAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCcgtg                                                                                        >  1:2235105/1‑93 (MQ=255)

  aaaCAGATTAGTGAGTTGAGTAAAAGCCCGGTCACATTGGACTAACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGAAACCCgtgt                                                                                       <  2:4220896/93‑1 (MQ=255)

  aaaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCgtgt                                                                                       <  1:145871/93‑1 (MQ=255)

  aaaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCgtgt                                                                                       >  2:1861563/1‑93 (MQ=255)

     cAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTc                                                                                    <  2:2069180/93‑1 (MQ=255)

      aGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCt                                                                                   >  1:865834/1‑93 (MQ=255)

       gATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtt                                                                                  <  2:2246490/93‑1 (MQ=255)

       gATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtt                                                                                  >  2:1210976/1‑93 (MQ=255)

        aatAGAGAGTTAAGTAAAAGCCCGATCACATTGGACTGACAAAGCTTACGTGAGTAATTCGCTGTATTCTTCCTGGGTAACCCGTGAATCttt                                                                                 <  2:4334871/91‑1 (MQ=255)

         ttAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtttt                                                                                >  1:3973918/1‑93 (MQ=255)

         ttAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtttt                                                                                <  2:1847819/93‑1 (MQ=255)

          tAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtg                                                                               <  1:1887061/93‑1 (MQ=255)

          tAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtg                                                                               <  2:104256/93‑1 (MQ=255)

           aGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtgg                                                                              >  1:4298527/1‑93 (MQ=255)

            gTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggt                                                                             <  1:2166280/93‑1 (MQ=255)

            gTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggt                                                                             <  2:4511131/93‑1 (MQ=255)

            gTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggt                                                                             >  2:839703/1‑93 (MQ=255)

              gAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtgg                                                                           <  2:1689080/93‑1 (MQ=255)

              gAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtgg                                                                           <  2:2195744/93‑1 (MQ=255)

               aGTTAAGTAAAAGCCCGGTCACATTGGACTGACAGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTggg                                                                          >  2:153655/1‑93 (MQ=255)

                gTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGt                                                                         >  2:1259852/1‑93 (MQ=255)

                gTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGt                                                                         <  1:1162161/93‑1 (MQ=255)

                 ttAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGtt                                                                        <  1:3449948/93‑1 (MQ=255)

                 ttAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGtt                                                                        <  2:3657759/93‑1 (MQ=255)

                  tAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGttt                                                                       <  1:4388200/93‑1 (MQ=255)

                   aaGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTa                                                                      <  2:2780599/93‑1 (MQ=255)

                   aaGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTa                                                                      <  1:3002014/93‑1 (MQ=255)

                    aGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTAt                                                                     >  1:4616770/1‑93 (MQ=255)

                       aaaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTcc                                                                  <  2:2169589/93‑1 (MQ=255)

                       aaaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTcc                                                                  >  2:3477787/1‑93 (MQ=255)

                        aaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCg                                                                 >  1:1077574/1‑93 (MQ=255)

                         aaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCgg                                                                <  1:1328659/93‑1 (MQ=255)

                         aaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCgg                                                                >  2:41431/1‑93 (MQ=255)

                          aGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGc                                                               >  2:3536211/1‑93 (MQ=255)

                           gCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCa                                                              >  2:873703/1‑93 (MQ=255)

                           gCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCa                                                              <  1:2786698/93‑1 (MQ=255)

                           gCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCa                                                              <  1:3100575/93‑1 (MQ=255)

                           gCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCa                                                              >  2:3694601/1‑93 (MQ=255)

                            cccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCAt                                                             <  1:1521938/93‑1 (MQ=255)

                            cccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCAt                                                             <  2:1766877/93‑1 (MQ=255)

                             ccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATc                                                            <  2:2710201/93‑1 (MQ=255)

                             ccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATc                                                            >  1:263586/1‑93 (MQ=255)

                              cGGTCACATTGAACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCa                                                           >  2:557803/1‑93 (MQ=255)

                               ggTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCAc                                                          <  2:2229926/93‑1 (MQ=255)

                                gTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACg                                                         >  1:1774716/1‑93 (MQ=255)

                                 tCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGa                                                        <  2:1984439/93‑1 (MQ=255)

                                    caTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAata                                                     >  1:424459/1‑93 (MQ=255)

                                     aTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAatat                                                    >  2:3979341/1‑93 (MQ=255)

                                     aTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAatat                                                    >  1:2713351/1‑93 (MQ=255)

                                     aTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAatat                                                    >  1:238447/1‑93 (MQ=255)

                                      ttGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAtt                                                   >  1:1993828/1‑93 (MQ=255)

                                         gACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAttttt                                                >  2:1420416/1‑93 (MQ=255)

                                          aCTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTc                                               >  2:4099774/1‑93 (MQ=255)

                                           cTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCt                                              <  1:4648707/93‑1 (MQ=255)

                                           cTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCt                                              <  1:476958/93‑1 (MQ=255)

                                           cTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCt                                              >  2:1506456/1‑93 (MQ=255)

                                            tGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTg                                             <  2:537666/93‑1 (MQ=255)

                                                cGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCAcc                                         <  2:2962899/93‑1 (MQ=255)

                                                cGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCAcc                                         >  2:64467/1‑93 (MQ=255)

                                                 gggCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCg                                        <  1:80230/93‑1 (MQ=255)

                                                  ggCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGc                                       <  2:730371/93‑1 (MQ=255)

                                                    cTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTa                                     >  2:488420/1‑93 (MQ=255)

                                                     ttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTAc                                    <  2:4587532/93‑1 (MQ=255)

                                                     ttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTAc                                    <  2:2139180/93‑1 (MQ=255)

                                                     ttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTAc                                    >  1:937569/1‑93 (MQ=255)

                                                     ttACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTAc                                    <  1:4110671/93‑1 (MQ=255)

                                                        cGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGc                                 >  1:3582609/1‑93 (MQ=255)

                                                           gAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCAcc                              <  2:3039220/93‑1 (MQ=255)

                                                             gTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaa                            >  2:4677060/1‑93 (MQ=255)

                                                             gTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaa                            <  1:779357/93‑1 (MQ=255)

                                                              ttATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaaa                           <  1:3119787/93‑1 (MQ=255)

                                                               tATTCCGTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCCCGAATATTTTTCTGCCCCGCTACTGCACCAAAc                          <  2:150458/93‑1 (MQ=255)

                                                                aTTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTTGGTTTATTCCGGCATTACGAATATCTTTCTTCATCGCTACTGCTCCATACa                         >  1:4207580/1‑93 (MQ=255)

                                                                 ttCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACaa                        >  2:1081432/1‑93 (MQ=255)

                                                                 ttCGCTGTATTCTTCCTGCGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCCCGAATATTTTTCTGCACCGCTATTGCACCAAACaa                        >  1:3310374/1‑93 (MQ=255)

                                                                   cGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGc                      >  2:629350/1‑93 (MQ=255)

                                                                     gtGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCtt                    <  1:3430699/92‑1 (MQ=255)

                                                                     cTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCtt                    <  2:3560185/93‑1 (MQ=255)

                                                                      tGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTa                   >  1:3930742/1‑93 (MQ=255)

                                                                      tGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTa                   >  2:2894478/1‑93 (MQ=255)

                                                                       gTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaa                  >  2:322555/1‑93 (MQ=255)

                                                                       gTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaa                  <  1:991753/93‑1 (MQ=255)

                                                                       gTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaa                  <  1:2665182/93‑1 (MQ=255)

                                                                            tttcCTGGGTATCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAagaaga             <  2:1954676/92‑1 (MQ=255)

                                                                             ttcCTGGGTAACCCGTGTTTCTTTTGGTTGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAgaagaa            <  1:642885/93‑1 (MQ=255)

                                                                             ttcCTGGGTAACCCGTGTTTCTGTTGGTGGGTTTATTCCGGCATCACTAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAgaagaa            <  1:642883/93‑1 (MQ=255)

                                                                              tcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGaaa           <  2:3907735/93‑1 (MQ=255)

                                                                              tcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGaaa           <  1:3533568/93‑1 (MQ=255)

                                                                               ccTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAg          >  1:1635765/1‑93 (MQ=255)

                                                                                cTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGc         >  1:3814121/1‑93 (MQ=255)

                                                                                cTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGc         >  1:2938282/1‑93 (MQ=255)

                                                                                  gggTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCa       <  1:3616996/93‑1 (MQ=255)

                                                                                  gggTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCa       <  1:1710339/93‑1 (MQ=255)

                                                                                  gggTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCa       <  1:4488957/93‑1 (MQ=255)

                                                                                     tAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTc    <  1:1223288/93‑1 (MQ=255)

                                                                                     tAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTc    <  2:1152336/93‑1 (MQ=255)

                                                                                      aaCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTcc   <  2:4452231/93‑1 (MQ=255)

                                                                                      aaCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTcc   <  1:2805904/93‑1 (MQ=255)

                                                                                       aCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCa  <  1:3891460/93‑1 (MQ=255)

                                                                                       aCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCa  >  1:3946503/1‑93 (MQ=255)

AGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCA  >  NC_000913/3724491‑3724670

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