BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,815,038	T→C	Q41R (CAG→CGG)	ydjC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,815,038	0	T	C	97.8%	618.7 / ‑0.4	181	Q41R (CAG→CGG)	ydjC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/3); new base C (86/91); total (87/94)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 195 total aligned reads displayed.

GTCAGTGGCTTACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCTGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCTTAAGCCAAAA  >  NC_000913/1814948‑1815127

 tCAGTGGCTTACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGcc                                                                                        <  1:5223738/93‑1 (MQ=255)

   aGTGGCTTACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCa                                                                                      <  2:1825997/93‑1 (MQ=255)

    gTGGCTTACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGAACCGCATGGTCAATAGCCCGCCCAt                                                                                     >  2:4298941/1‑93 (MQ=255)

      ggCTTACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTc                                                                                   <  2:2209327/93‑1 (MQ=255)

         ttACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCAcc                                                                                <  1:1001380/93‑1 (MQ=255)

         ttACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCAcc                                                                                >  1:1341936/1‑93 (MQ=255)

         ttACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCAcc                                                                                <  2:5303386/93‑1 (MQ=255)

          tACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCa                                                                               <  2:3546584/93‑1 (MQ=255)

           aCCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAg                                                                              <  1:4407246/93‑1 (MQ=255)

            cccATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAgc                                                                             >  1:4674683/1‑93 (MQ=255)

            cccATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCACGCCCATTCACAAgc                                                                             >  2:1976642/1‑93 (MQ=255)

             ccATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAgcg                                                                            <  1:1265769/93‑1 (MQ=255)

              cATAGTAAGGACAAAGTGCATCCCTATGGCCAGAATCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCACATTCACCAgcgc                                                                           <  1:4491974/93‑1 (MQ=255)

                tAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTGACCAGCGCcg                                                                         >  1:3948583/1‑93 (MQ=255)

                tAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCcg                                                                         >  2:3277540/1‑93 (MQ=255)

                tAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCcg                                                                         >  2:2183992/1‑93 (MQ=255)

                 aGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCcgt                                                                        <  2:1454671/93‑1 (MQ=255)

                  gTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCcgtc                                                                       <  2:1431894/93‑1 (MQ=255)

                   tAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTTATCACGACTCAACTGCACCGCATGGTCAACAGCCCGCCCATTCACCAGCGCcggcg                                                                      >  1:209694/1‑93 (MQ=255)

                   tAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCcgtcg                                                                      >  2:4234038/1‑93 (MQ=255)

                    aaGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCcgtcgt                                                                     <  2:3178390/93‑1 (MQ=255)

                    aaGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCcgtcgc                                                                     >  2:5830672/1‑92 (MQ=255)

                         cAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACg                                                                >  1:2747525/1‑93 (MQ=255)

                         cAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACg                                                                >  2:1048197/1‑93 (MQ=255)

                          aaaGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCTGCCCATTCACCAGCGCCGTCGTCGACGt                                                               <  2:4128632/93‑1 (MQ=255)

                             gTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGAc                                                            >  2:3065110/1‑93 (MQ=255)

                               gCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACaa                                                          <  2:5755661/93‑1 (MQ=255)

                               gCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACaa                                                          <  2:3817221/93‑1 (MQ=255)

                               gCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACaa                                                          >  2:4643427/1‑93 (MQ=255)

                               gCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACaa                                                          <  2:1806463/93‑1 (MQ=255)

                                 aTCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATc                                                        <  2:1740566/93‑1 (MQ=255)

                                  tCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATcc                                                       >  2:644349/1‑93 (MQ=255)

                                  tCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATcc                                                       >  1:2955246/1‑93 (MQ=255)

                                  tCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATcc                                                       <  2:4688390/93‑1 (MQ=255)

                                  tCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATcc                                                       >  2:3347423/1‑93 (MQ=255)

                                    ccTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCNCNGTCGTCGACGTGACAATCCCa                                                     >  2:1940917/1‑93 (MQ=255)

                                    ccTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCa                                                     >  1:5423712/1‑93 (MQ=255)

                                    ccTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCa                                                     >  2:4793867/1‑93 (MQ=255)

                                     cTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCAt                                                    <  1:5539003/93‑1 (MQ=255)

                                     cTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCAt                                                    <  2:3989604/93‑1 (MQ=255)

                                     cTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCAt                                                    <  1:1914930/93‑1 (MQ=255)

                                      tATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCAtt                                                   >  2:729777/1‑93 (MQ=255)

                                      tATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCAtt                                                   >  2:1636141/1‑93 (MQ=255)

                                       aTGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTg                                                  <  1:4345266/93‑1 (MQ=255)

                                        tGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGc                                                 <  2:1142806/93‑1 (MQ=255)

                                        tGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGc                                                 >  1:1657169/1‑93 (MQ=255)

                                        tGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGc                                                 <  1:4966503/93‑1 (MQ=255)

                                        tGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGc                                                 >  2:3682948/1‑93 (MQ=255)

                                         ggCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTTGACGTGACAATCCCATTGCg                                                >  1:6027565/1‑93 (MQ=255)

                                            cAGACTCGGTTCATCACGACTCAACTGCAACGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACa                                             <  2:2212065/93‑1 (MQ=255)

                                              gACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAgg                                           <  1:2383997/93‑1 (MQ=255)

                                               aCTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGc                                          >  2:4386511/1‑93 (MQ=255)

                                               aCTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGc                                          >  1:2606054/1‑93 (MQ=255)

                                                cTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACNTNACAATCCCATNGCGACAggac                                         >  2:5713127/1‑91 (MQ=255)

                                                 tCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGcat                                        >  2:1603130/1‑91 (MQ=255)

                                                 tCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCt                                        <  2:1220472/93‑1 (MQ=255)

                                                  cGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTc                                       >  1:101676/1‑93 (MQ=255)

                                                   ggTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCg                                      >  1:2054986/1‑93 (MQ=255)

                                                    gTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGa                                     <  1:1197279/93‑1 (MQ=255)

                                                    gTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGa                                     >  2:5250889/1‑93 (MQ=255)

                                                     ttCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGAt                                    <  1:4292078/93‑1 (MQ=255)

                                                     ttCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGAt                                    <  2:3728834/93‑1 (MQ=255)

                                                     ttCATCACGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGAt                                    <  1:586212/93‑1 (MQ=255)

                                                        atcaCGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAAt                                 <  2:203138/93‑1 (MQ=255)

                                                        atcaCGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAAt                                 <  2:4180268/93‑1 (MQ=255)

                                                        atcaCGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAAt                                 <  2:3317449/93‑1 (MQ=255)

                                                          caCGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGc                               >  2:5954840/1‑93 (MQ=255)

                                                            cGACTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCg                             <  1:2668257/93‑1 (MQ=255)

                                                              gcTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTa                           <  2:2076150/92‑1 (MQ=255)

                                                               cTCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAg                          <  1:74983/93‑1 (MQ=255)

                                                                tCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGt                         <  2:3239102/93‑1 (MQ=255)

                                                                tCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGt                         >  2:1173768/1‑93 (MQ=255)

                                                                tCAACTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGt                         <  1:3275742/93‑1 (MQ=255)

                                                                  aaCTGCACCGCATGGTCAATAGCCTGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTc                       <  2:6012109/93‑1 (MQ=255)

                                                                  aaCTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTc                       <  2:2380793/93‑1 (MQ=255)

                                                                  aaCCGCTCCGAATGGACAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTc                       <  2:1921055/93‑1 (MQ=255)

                                                                   aCTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCt                      <  1:930024/93‑1 (MQ=255)

                                                                   aCTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCt                      >  2:2814906/1‑93 (MQ=255)

                                                                    cTGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTg                     >  1:3131754/1‑93 (MQ=255)

                                                                     tGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGa                    >  1:890610/1‑93 (MQ=255)

                                                                     tGCACCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGa                    >  1:792907/1‑93 (MQ=255)

                                                                        aCCGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCt                 >  2:6144853/1‑93 (MQ=255)

                                                                         ccGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCtt                <  2:2646445/93‑1 (MQ=255)

                                                                          cGCATGGTGAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCttt               <  2:5102807/93‑1 (MQ=255)

                                                                          cGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCttt               >  1:1298023/1‑93 (MQ=255)

                                                                          cGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCttt               <  1:2532165/93‑1 (MQ=255)

                                                                          cGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCttt               <  2:4966848/93‑1 (MQ=255)

                                                                          cGCATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCttt               >  2:2270457/1‑93 (MQ=255)

                                                                            cATGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGc             <  1:5404009/93‑1 (MQ=255)

                                                                             aTGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCt            >  1:5464624/1‑93 (MQ=255)

                                                                             aTGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCt            <  2:4071526/93‑1 (MQ=255)

                                                                              tGGTCAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAAGCCTCGATAATGCCGTAGTTCTGACCTTTGCtt           >  2:5931005/1‑93 (MQ=255)

                                                                              tGGACAATAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCtt           <  1:4748873/93‑1 (MQ=255)

                                                                                gTCAATAGCCCGCCCATTCTCCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCTTaa         <  2:839895/93‑1 (MQ=255)

                                                                                   aaTAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCTTAAGcc      >  1:2192094/1‑93 (MQ=255)

                                                                                   aaTAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCTTAAGcc      <  1:3547987/93‑1 (MQ=255)

                                                                                    aTAGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCTTAAGCCa     <  1:3477230/93‑1 (MQ=255)

                                                                                      aGCCCGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGCTAATGCCGTAGTTCTGACCTTTTCTTAAGCCaaa   >  1:2184957/1‑93 (MQ=255)

                                                                                       gcccgcccATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCTTAAGCCaaaa  <  1:1037369/93‑1 (MQ=255)

                                                                                       gCCTGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTCACCTTTGCTTAAGCCaaaa  >  1:6081160/1‑93 (MQ=255)

GTCAGTGGCTTACCCATAGTAAGGACAAAGTGCATCCCTATGGCCAGACTCGGTTCATCACGACTCAACTGCACCGCATGGTCAATAGCCTGCCCATTCACCAGCGCCGTCGTCGACGTGACAATCCCATTGCGACAGGCCTCGATAATGCCGTAGTTCTGACCTTTGCTTAAGCCAAAA  >  NC_000913/1814948‑1815127

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