BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,974,008	C→T	A40V (GCG→GTG)	wzxE →	putative cytochrome

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,974,008	0	C	T	99.4%	550.1 / NA	162	A40V (GCG→GTG)	wzxE	putative cytochrome
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (72/89); minor base G (1/0); total (73/89)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 166 total aligned reads displayed.

CGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACC  >  NC_000913/3973921‑3974097

cggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCaa                                                                                      <  2:4443257/93‑1 (MQ=255)

cggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCaa                                                                                      <  1:2342409/93‑1 (MQ=255)

 ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCaaa                                                                                     >  1:2357570/1‑93 (MQ=255)

 ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCaaa                                                                                     >  2:3355187/1‑93 (MQ=255)

 ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCaaa                                                                                     >  1:1526551/1‑93 (MQ=255)

  gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAAt                                                                                    >  1:3454373/1‑93 (MQ=255)

   ccAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAAtt                                                                                   <  2:1133748/93‑1 (MQ=255)

     aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGGGTCGTTTGGTCCCGCGGGGCTTGGGCTGGTGGGAAAtctc                                                                                 >  2:3579750/1‑93 (MQ=255)

      gTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTcg                                                                                >  1:3877021/1‑92 (MQ=255)

      gTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTcc                                                                                >  2:109375/1‑93 (MQ=255)

      gTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTcc                                                                                >  2:1639237/1‑93 (MQ=255)

         cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGcc                                                                             >  2:3491472/1‑93 (MQ=255)

         cacTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATGTCCGcc                                                                             >  2:1470380/1‑93 (MQ=255)

          acTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCa                                                                            <  2:3946475/93‑1 (MQ=255)

           cTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTCCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAg                                                                           <  1:897984/93‑1 (MQ=255)

            tGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGt                                                                          <  1:2986711/93‑1 (MQ=255)

            tGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGt                                                                          >  1:2030925/1‑93 (MQ=255)

            tGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGt                                                                          >  2:1156498/1‑93 (MQ=255)

             ggTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGtt                                                                         <  1:3226660/93‑1 (MQ=255)

                cAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGAt                                                                      <  1:267539/93‑1 (MQ=255)

                cAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGAt                                                                      <  1:3074604/93‑1 (MQ=255)

                cAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGAt                                                                      <  2:3584395/93‑1 (MQ=255)

                 aaGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGAtt                                                                     <  1:201122/93‑1 (MQ=255)

                 aaGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGAtt                                                                     <  2:865951/93‑1 (MQ=255)

                  aGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTa                                                                    <  1:2977743/93‑1 (MQ=255)

                  aGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTa                                                                    >  2:3503065/1‑93 (MQ=255)

                  aGATTGGTGCCGGGGTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCAAGTTCATTa                                                                    <  2:1679991/93‑1 (MQ=255)

                   gATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTAc                                                                   >  1:2927902/1‑93 (MQ=255)

                        gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGc                                                              <  1:987844/93‑1 (MQ=255)

                        gTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGc                                                              <  2:2657224/93‑1 (MQ=255)

                         ggCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCt                                                             <  2:99605/92‑1 (MQ=255)

                             gggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCggg                                                         >  2:844539/1‑92 (MQ=255)

                             gggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGc                                                         >  1:1194566/1‑93 (MQ=255)

                              ggTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCg                                                        >  2:29658/1‑93 (MQ=255)

                               gTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGt                                                       <  1:3194128/93‑1 (MQ=255)

                                ttACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCCAATTTTCGCCAGTTGATTACCGTGCTCGGCGTg                                                      >  2:1033926/1‑93 (MQ=255)

                                ttACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTg                                                      >  2:3011892/1‑93 (MQ=255)

                                 tACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGc                                                     >  1:1990496/1‑93 (MQ=255)

                                  aCTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCt                                                    <  2:93540/93‑1 (MQ=255)

                                  aCTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCt                                                    <  2:1715469/93‑1 (MQ=255)

                                  aCTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCt                                                    <  2:515682/93‑1 (MQ=255)

                                   cTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                   >  1:124249/1‑93 (MQ=255)

                                   cTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                   <  1:3171701/93‑1 (MQ=255)

                                      gTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGcc                                                <  1:2268305/93‑1 (MQ=255)

                                      gTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGcc                                                >  2:438318/1‑93 (MQ=255)

                                      gTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGcc                                                <  2:76695/93‑1 (MQ=255)

                                      gTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGcc                                                <  1:47712/93‑1 (MQ=255)

                                      gTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGcc                                                >  1:370717/1‑93 (MQ=255)

                                       tCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCg                                               <  1:4250076/93‑1 (MQ=255)

                                        cGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                              <  1:1310371/93‑1 (MQ=255)

                                        cGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                              <  1:1821509/93‑1 (MQ=255)

                                        cGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgg                                              <  2:3260947/93‑1 (MQ=255)

                                         ggTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                             <  2:1054238/93‑1 (MQ=255)

                                         ggTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                             >  2:4382529/1‑93 (MQ=255)

                                          gTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                            >  2:3656853/1‑93 (MQ=255)

                                           tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           <  1:336369/93‑1 (MQ=255)

                                           tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           <  1:1535728/93‑1 (MQ=255)

                                           tAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           >  2:1868215/1‑93 (MQ=255)

                                           gaaGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                           <  1:883012/92‑1 (MQ=255)

                                               ttGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGc                                       <  2:1595554/93‑1 (MQ=255)

                                                tGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCa                                      <  2:2941394/93‑1 (MQ=255)

                                                tGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCa                                      <  1:1639237/93‑1 (MQ=255)

                                                  cTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATc                                    <  1:1883007/93‑1 (MQ=255)

                                                   tGGCGGTGTCATTTGGGCCGGCGGGGCTTGGGCTGGTGGGAAATTTCCGCCAGCTGATTACCCTGTTCAGCGTGCTTGGCGGGGCTGGCATCt                                   >  2:3975756/1‑93 (MQ=255)

                                                    ggCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                  >  2:1379932/1‑93 (MQ=255)

                                                    ggCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                  <  2:3256159/93‑1 (MQ=255)

                                                     gCGGTGTCATTTTGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                 <  2:3554781/93‑1 (MQ=255)

                                                     gCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                 >  2:575019/1‑93 (MQ=255)

                                                     gCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                 >  1:3936369/1‑93 (MQ=255)

                                                       ggTGTCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTaa                               <  2:150674/93‑1 (MQ=255)

                                                        gtgtCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                              >  1:4380494/1‑93 (MQ=255)

                                                         tgtCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACg                             >  1:4676596/1‑93 (MQ=255)

                                                         tgtCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACg                             >  2:222599/1‑93 (MQ=255)

                                                          gtCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                            >  1:4151726/1‑93 (MQ=255)

                                                           tCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                           >  1:2857389/1‑93 (MQ=255)

                                                           tCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                           <  2:2353509/93‑1 (MQ=255)

                                                           tCATTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                           <  1:204757/93‑1 (MQ=255)

                                                             aTTTGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtgt                         >  2:4241564/1‑93 (MQ=255)

                                                              tttGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgtgga                        >  2:270749/1‑91 (MQ=255)

                                                              tttGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                        >  1:1206872/1‑93 (MQ=255)

                                                              tttGGTCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                        >  2:661476/1‑93 (MQ=255)

                                                                 ggTCCGTCGGGACTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAcc                     <  2:307696/93‑1 (MQ=255)

                                                                   tCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                   <  1:690372/93‑1 (MQ=255)

                                                                   tCCGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                   <  1:2209858/93‑1 (MQ=255)

                                                                    ccGGCGGGGCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaaa                  <  2:2189790/93‑1 (MQ=255)

                                                                          gggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt            >  1:3393883/1‑93 (MQ=255)

                                                                          gggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt            <  1:709033/93‑1 (MQ=255)

                                                                          gggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt            <  1:4497196/93‑1 (MQ=255)

                                                                           ggCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTg           >  1:989356/1‑93 (MQ=255)

                                                                            gCTTGGGCTGGTGGCAAATTTCCGCCAGTTTATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc          <  1:2631571/93‑1 (MQ=255)

                                                                            gCTTGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc          <  1:4244517/93‑1 (MQ=255)

                                                                              ttGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGccc        >  2:722536/1‑93 (MQ=255)

                                                                               tGGGCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCa       <  2:1833772/93‑1 (MQ=255)

                                                                                gggCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg      >  1:2772407/1‑93 (MQ=255)

                                                                                  gCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa    <  2:2115917/93‑1 (MQ=255)

                                                                                  gCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa    >  1:2546461/1‑93 (MQ=255)

                                                                                  gCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa    <  1:3112297/93‑1 (MQ=255)

                                                                                  gCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa    >  1:822197/1‑93 (MQ=255)

                                                                                  gCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa    >  1:1014901/1‑93 (MQ=255)

                                                                                  gCTGGTGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTAGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa    <  1:942841/93‑1 (MQ=255)

CGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACC  >  NC_000913/3973921‑3974097

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