BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	200,658	T→C	R59R (CGT→CGC)	hlpA →	histone‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	200,658	0	T	C	97.5%	598.0 / 3.0	197	R59R (CGT→CGC)	hlpA	histone‑like protein
Reads supporting (aligned to +/- strand): ref base T (3/2); new base C (84/108); total (87/110)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 201 total aligned reads displayed.

GGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGC  >  NC_000913/200568‑200747

gggCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCAt                                                                                         <  1:885542/93‑1 (MQ=255)

 ggCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATg                                                                                        <  1:4417898/93‑1 (MQ=255)

 ggCAGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATg                                                                                        <  1:837934/93‑1 (MQ=255)

    aGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGaa                                                                                     >  1:1010699/1‑93 (MQ=255)

    aGCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGaa                                                                                     >  1:2475917/1‑93 (MQ=255)

     gCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGaaa                                                                                    >  2:1073516/1‑93 (MQ=255)

     gCCTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGaaa                                                                                    >  1:1090543/1‑93 (MQ=255)

     accTGTTCCAGCGGGTAGCGCAGAAACCCGGTGTTTCTAACACGCTGAAAAATGAGTTAAAAGCCCGTGCAAGCGAACTGCAGCGCATGGaaa                                                                                    <  2:1232160/92‑1 (MQ=255)

       cTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAAcc                                                                                  >  1:1851294/1‑93 (MQ=255)

       cTGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAAcc                                                                                  >  1:4655865/1‑93 (MQ=255)

        tGTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCg                                                                                 >  2:1438383/1‑93 (MQ=255)

         gTTCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGa                                                                                <  1:1530110/93‑1 (MQ=255)

          ttCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGAt                                                                               >  2:3994962/1‑93 (MQ=255)

           tCCAGCGGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATc                                                                              <  2:4325088/93‑1 (MQ=255)

              agcggGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGc                                                                           <  1:3570006/93‑1 (MQ=255)

               gcgggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCa                                                                          >  1:3238947/4‑93 (MQ=255)

               gcgggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCa                                                                          <  2:3430063/90‑1 (MQ=255)

                cgggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAg                                                                         <  1:269174/91‑1 (MQ=255)

                cgggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAg                                                                         <  1:2600603/91‑1 (MQ=255)

                cgggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAg                                                                         <  2:2035069/91‑1 (MQ=255)

                cgggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAg                                                                         <  2:1617699/91‑1 (MQ=255)

                cgggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAg                                                                         <  1:1289022/91‑1 (MQ=255)

                 gggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAgg                                                                        <  1:2187683/92‑1 (MQ=255)

                  ggTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGc                                                                       <  2:226254/93‑1 (MQ=255)

                    tAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTa                                                                     >  2:4214731/1‑93 (MQ=255)

                      gcgcAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTaaa                                                                   >  1:2819191/1‑93 (MQ=255)

                       cgcAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTaaaa                                                                  <  2:2259416/93‑1 (MQ=255)

                       cgcAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTaaaa                                                                  <  2:2318823/93‑1 (MQ=255)

                         cAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATg                                                                >  1:3488602/1‑93 (MQ=255)

                         cAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATg                                                                >  2:2780740/1‑93 (MQ=255)

                          aGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGa                                                               >  1:638345/1‑93 (MQ=255)

                          aGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGa                                                               >  2:1277347/1‑93 (MQ=255)

                          aGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGa                                                               <  2:4267679/93‑1 (MQ=255)

                          aGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGa                                                               <  1:600803/93‑1 (MQ=255)

                            aaaaCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGaaa                                                             >  1:3831484/1‑93 (MQ=255)

                            aaaaCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGaaa                                                             >  2:772536/1‑93 (MQ=255)

                            aaaaCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGaaa                                                             >  1:327975/1‑93 (MQ=255)

                             aaaCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGaaaa                                                            <  2:4559905/93‑1 (MQ=255)

                             aaaCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGaaaa                                                            >  2:1337813/1‑93 (MQ=255)

                              aaCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGaaaaa                                                           <  1:3175102/93‑1 (MQ=255)

                              aaCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGaaaaa                                                           <  1:4185824/93‑1 (MQ=255)

                              aaCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGaaaaa                                                           <  2:1039234/93‑1 (MQ=255)

                               aCCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAg                                                          >  2:3667888/1‑93 (MQ=255)

                                ccGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGc                                                         <  2:1923497/93‑1 (MQ=255)

                                   gtgtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGc                                                      >  1:1193163/1‑93 (MQ=255)

                                   gtgtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGc                                                      <  2:285230/93‑1 (MQ=255)

                                   gtgtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGc                                                      >  2:1992281/1‑93 (MQ=255)

                                   gtgtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGc                                                      <  2:3391529/93‑1 (MQ=255)

                                   gtgtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGc                                                      <  1:3676018/93‑1 (MQ=255)

                                    tgtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCa                                                     >  2:3241877/1‑93 (MQ=255)

                                    tgtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCa                                                     >  2:263857/1‑93 (MQ=255)

                                    tgtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCACATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCa                                                     <  2:257408/93‑1 (MQ=255)

                                     gtTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAg                                                    >  2:2825715/1‑93 (MQ=255)

                                      tttCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGt                                                   >  1:4210135/1‑93 (MQ=255)

                                      tttCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGt                                                   >  2:3282002/1‑93 (MQ=255)

                                         cTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCa                                                >  1:3403282/1‑93 (MQ=255)

                                          tAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCAt                                               <  1:3981568/93‑1 (MQ=255)

                                           aaCACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATg                                              >  2:71940/1‑93 (MQ=255)

                                           aaCACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATg                                              >  1:1696203/1‑93 (MQ=255)

                                            acacGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAGACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGa                                             <  2:3829041/93‑1 (MQ=255)

                                            acacGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGa                                             >  2:2676632/1‑93 (MQ=255)

                                              acGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGaaa                                           <  2:653962/93‑1 (MQ=255)

                                                gCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGc                                         >  1:2553846/1‑93 (MQ=255)

                                                  tGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCgg                                       <  2:3056542/93‑1 (MQ=255)

                                                     aaaaTGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCa                                    <  2:898825/93‑1 (MQ=255)

                                                     aaaaTGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCa                                    <  1:2985809/93‑1 (MQ=255)

                                                      aaaTGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAg                                   >  2:1859877/1‑93 (MQ=255)

                                                       aaTGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGc                                  >  2:4669908/1‑93 (MQ=255)

                                                       aaTGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGc                                  <  1:1255569/93‑1 (MQ=255)

                                                       aaTGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGc                                  <  1:716419/93‑1 (MQ=255)

                                                        aTGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCg                                 >  2:3853076/1‑93 (MQ=255)

                                                         tGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGa                                <  1:1626993/93‑1 (MQ=255)

                                                         tGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGa                                <  1:3886690/93‑1 (MQ=255)

                                                          gAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGGGAt                               >  2:1363374/1‑93 (MQ=255)

                                                           aGTTCAAAGGCCGTGCCAGCGAACTGCAGGGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATc                              <  2:4602007/93‑1 (MQ=255)

                                                           aGTTCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATc                              <  1:387260/93‑1 (MQ=255)

                                                              tCAAAGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATCAAAAAGCTGCAGTCAATGAAAGAGGGCAGAGATCgaa                           >  2:2487881/1‑91 (MQ=255)

                                                                aaaGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACt                         >  1:4040706/1‑93 (MQ=255)

                                                                aaaGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACt                         >  2:3315890/1‑93 (MQ=255)

                                                                  aGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTaa                       >  1:2849634/1‑93 (MQ=255)

                                                                  aGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTaa                       >  2:1513129/1‑93 (MQ=255)

                                                                  aGGCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTaa                       >  1:2277850/1‑93 (MQ=255)

                                                                    gCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGc                     <  1:1057121/93‑1 (MQ=255)

                                                                    gCCGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGc                     <  2:709215/93‑1 (MQ=255)

                                                                     ccGTGCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCt                    <  1:425657/93‑1 (MQ=255)

                                                                       gTGCCAGCTAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTgg                  <  1:2101798/93‑1 (MQ=255)

                                                                         gCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGaa                >  1:3227561/1‑93 (MQ=255)

                                                                         gCCAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGaa                <  1:763566/93‑1 (MQ=255)

                                                                           cAGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGaaaa              >  2:1752123/1‑93 (MQ=255)

                                                                            aGCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGaaaaa             >  1:502708/1‑93 (MQ=255)

                                                                             gCGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAg            >  1:3317634/1‑93 (MQ=255)

                                                                              cGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGa           <  1:2359228/93‑1 (MQ=255)

                                                                              cGAACTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGa           <  1:2076391/93‑1 (MQ=255)

                                                                                aaCTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACg         >  2:2983016/1‑93 (MQ=255)

                                                                                 aCTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGt        >  2:4175663/1‑93 (MQ=255)

                                                                                 aCTGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGt        <  2:2092781/93‑1 (MQ=255)

                                                                                   tGCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGa      <  2:48934/93‑1 (MQ=255)

                                                                                    gCAGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGAt     <  1:4506445/93‑1 (MQ=255)

                                                                                      aGCGCATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATgg   >  2:2211239/1‑93 (MQ=255)

                                                                                       gcgcATGGAATCCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGc  <  2:4480938/93‑1 (MQ=255)

GGGCAGCCTGTTCCAGCAGGTAGCGCAGAAAACCGGTGTTTCTAACACGCTGGAAAATGAGTTCAAAGGCCGTGCCAGCGAACTGCAGCGTATGGAAACCGATCTGCAGGCTAAAATGAAAAAGCTGCAGTCCATGAAAGCGGGCAGCGATCGCACTAAGCTGGAAAAAGACGTGATGGC  >  NC_000913/200568‑200747

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