BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,012,620	T→C	G132G (GGA→GGG)	b2876 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,012,620	0	T	C	97.8%	537.6 / 0.4	183	G132G (GGA→GGG)	b2876	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (86/91); total (90/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

Only 100 of 195 total aligned reads displayed.

CAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCC  >  NC_000913/3012530‑3012709

cAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGccc                                                                                         <  2:2079211/93‑1 (MQ=255)

  aGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGccctc                                                                                       >  1:3428708/1‑93 (MQ=255)

  aGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGccctc                                                                                       >  1:2317521/1‑93 (MQ=255)

   ggTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGccctcc                                                                                      <  1:1122063/93‑1 (MQ=255)

   ggTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGccctcc                                                                                      >  2:80934/1‑93 (MQ=255)

    gTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTccc                                                                                     <  1:1237916/93‑1 (MQ=255)

    gTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTccc                                                                                     <  2:1838827/93‑1 (MQ=255)

    gTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTccc                                                                                     >  2:1649014/1‑93 (MQ=255)

      tGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCAt                                                                                   >  1:562626/1‑93 (MQ=255)

      tGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCAt                                                                                   >  1:3918906/1‑93 (MQ=255)

       gCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATc                                                                                  <  2:920437/93‑1 (MQ=255)

       gCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATc                                                                                  <  1:1756126/93‑1 (MQ=255)

         aTCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATcac                                                                                >  2:968128/1‑93 (MQ=255)

         aTCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATcac                                                                                >  1:1351129/1‑93 (MQ=255)

          tCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATcaca                                                                               >  2:4074544/1‑93 (MQ=255)

          tCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATcaca                                                                               >  2:2838871/1‑93 (MQ=255)

          tCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATcaca                                                                               >  1:719461/1‑93 (MQ=255)

             ggTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCa                                                                            >  1:1035530/1‑93 (MQ=255)

              gtgtTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCaa                                                                           >  1:1851709/1‑93 (MQ=255)

               tgtTACCCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAAt                                                                          <  2:3321753/93‑1 (MQ=255)

               tgtTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAAt                                                                          >  1:3842037/1‑93 (MQ=255)

                 ttAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGCCATTTTCTGTGCTCACTTTCGTGCCCAAAGTATGCACTCCCATCACAGCAATCa                                                                        >  2:393441/1‑93 (MQ=255)

                  tAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCAc                                                                       <  1:4259285/93‑1 (MQ=255)

                    aCCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACgc                                                                     >  2:2440472/1‑93 (MQ=255)

                     cccGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACgca                                                                    >  2:3493489/1‑93 (MQ=255)

                      ccGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACgcag                                                                   >  1:1798280/1‑93 (MQ=255)

                      ccGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACgcag                                                                   >  2:25589/1‑93 (MQ=255)

                       cGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACgcagc                                                                  >  2:2626596/1‑93 (MQ=255)

                       cGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACgcagc                                                                  <  2:2111944/93‑1 (MQ=255)

                       cGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCCCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACgcagc                                                                  >  1:452743/1‑93 (MQ=255)

                        ggTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACgcagca                                                                 >  2:470281/1‑93 (MQ=255)

                        ggTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGACCTCCCATCACAGCAATCACgcagca                                                                 >  1:1681811/1‑93 (MQ=255)

                         gTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCaa                                                                >  1:4449132/1‑93 (MQ=255)

                           aaTATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACt                                                              <  1:2688742/93‑1 (MQ=255)

                            atatCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTg                                                             <  2:3432288/93‑1 (MQ=255)

                             tatCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCCACTGc                                                            >  1:1899581/1‑93 (MQ=255)

                             tatCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGc                                                            >  2:1252001/1‑93 (MQ=255)

                             tatCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGc                                                            <  1:2283018/93‑1 (MQ=255)

                               tCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCtt                                                          >  2:317809/1‑93 (MQ=255)

                                  gCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTg                                                       >  1:1207077/1‑93 (MQ=255)

                                  gCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTg                                                       >  2:187603/1‑93 (MQ=255)

                                  gCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTg                                                       <  2:1885609/93‑1 (MQ=255)

                                    aaaCTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTggg                                                     <  1:2805452/93‑1 (MQ=255)

                                      aCTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTa                                                   <  2:1598229/93‑1 (MQ=255)

                                      aCTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTa                                                   <  2:179021/93‑1 (MQ=255)

                                        tGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATg                                                 <  2:2104463/93‑1 (MQ=255)

                                         gcgcCCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGc                                                <  2:3669464/93‑1 (MQ=255)

                                         gCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGc                                                >  1:813416/1‑93 (MQ=255)

                                         gCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGc                                                <  1:1803516/93‑1 (MQ=255)

                                         gCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGc                                                >  2:4023330/1‑93 (MQ=255)

                                         gCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGc                                                <  1:4672808/93‑1 (MQ=255)

                                           gACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCaa                                              <  1:2998969/93‑1 (MQ=255)

                                           gACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCaa                                              <  1:4621205/93‑1 (MQ=255)

                                               aGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGtt                                          <  1:2475467/93‑1 (MQ=255)

                                               aGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGtt                                          >  2:3142124/1‑93 (MQ=255)

                                               aGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGtt                                          <  2:1804744/93‑1 (MQ=255)

                                                 cGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCg                                        <  2:4180915/93‑1 (MQ=255)

                                                  gATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGt                                       <  2:330263/93‑1 (MQ=255)

                                                  gATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGt                                       >  1:3938249/1‑93 (MQ=255)

                                                   aTGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTg                                      >  2:844791/1‑93 (MQ=255)

                                                    tGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGc                                     <  2:2653061/93‑1 (MQ=255)

                                                    tGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGc                                     >  1:3010746/1‑93 (MQ=255)

                                                    tGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGc                                     <  1:597791/93‑1 (MQ=255)

                                                     tgACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCt                                    <  2:544893/92‑1 (MQ=255)

                                                        cATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCtcat                                 >  2:2250454/1‑93 (MQ=255)

                                                         aTTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCtcatc                                >  2:1098823/1‑93 (MQ=255)

                                                         aTTTTCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCtcatc                                <  2:3894664/93‑1 (MQ=255)

                                                          ttttCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCtcatca                               >  1:910671/1‑93 (MQ=255)

                                                          ttttCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCtcatca                               >  1:1744108/1‑93 (MQ=255)

                                                           tttCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAg                              <  1:779008/93‑1 (MQ=255)

                                                           tttCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAg                              <  1:3085308/93‑1 (MQ=255)

                                                            ttCTGTGCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAgg                             <  1:2936838/93‑1 (MQ=255)

                                                                gtgCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCt                         >  1:1493228/1‑93 (MQ=255)

                                                                gtgCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCt                         <  2:1195735/93‑1 (MQ=255)

                                                                  gCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCttt                       <  2:2425823/93‑1 (MQ=255)

                                                                  gCTCACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCttt                       <  1:2056660/93‑1 (MQ=255)

                                                                     cACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTaa                    <  1:1453344/93‑1 (MQ=255)

                                                                     cACTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTaa                    <  1:2178291/93‑1 (MQ=255)

                                                                      aCTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAAc                   <  2:329249/93‑1 (MQ=255)

                                                                       cTTTCGCGCGCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACg                  <  2:4482460/93‑1 (MQ=255)

                                                                       cTTTCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACg                  <  1:3401029/93‑1 (MQ=255)

                                                                        tttCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACgg                 >  2:4399793/1‑93 (MQ=255)

                                                                          tCGCGCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCACCTGCTTTTGGGTATGCAAGGTTCGTGCTCGTCAGGCGCTTTTAACGGCa               <  2:4569919/93‑1 (MQ=255)

                                                                           cgcgCCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCAt              <  2:1919936/93‑1 (MQ=255)

                                                                             cgcCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTc            >  1:3283431/1‑93 (MQ=255)

                                                                             cgcCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTc            <  1:4143391/93‑1 (MQ=255)

                                                                              gcCCAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTcc           >  1:2849273/1‑93 (MQ=255)

                                                                               cccAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTGTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCa          <  1:3800083/93‑1 (MQ=255)

                                                                                ccAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAc         <  2:2330788/93‑1 (MQ=255)

                                                                                ccAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAc         >  2:1713525/1‑93 (MQ=255)

                                                                                ccAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAc         >  2:1608214/1‑93 (MQ=255)

                                                                                ccAAAGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAc         <  2:3004524/93‑1 (MQ=255)

                                                                                  aaaGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcgc       <  1:1215521/93‑1 (MQ=255)

                                                                                   aaGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcgcg      >  2:144209/1‑93 (MQ=255)

                                                                                   aaGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcgcg      >  1:4191923/1‑93 (MQ=255)

                                                                                   aaGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCAcgcg      >  1:2086189/1‑93 (MQ=255)

                                                                                    aGTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGa     <  1:3678570/93‑1 (MQ=255)

                                                                                     gTATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAg    >  1:1338839/1‑93 (MQ=255)

                                                                                      tATGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGc   <  1:2574149/93‑1 (MQ=255)

                                                                                       aTGCCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGcc  <  1:3768619/93‑1 (MQ=255)

CAAGGTTGCATCAGGTGTTAACCCGGTAATATCAGCAAACTGCGACCAGCGATGGACATTTTCTGTGCTCACTTTCGCGCCCAAAGTATGTCCTCCCATCACAGCAATCACGCAGCAACTGCTTTTGGGTATGCAAGGTTCGTGCTCATCAGGCGCTTTTAACGGCATTCCACGCGAGCC  >  NC_000913/3012530‑3012709

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