BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,928,950	C→T	Q218Q (CAG→CAA)	viaA ←	protein ViaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,928,950	0	C	T	99.7%	1314.0 / NA	306	Q218Q (CAG→CAA)	viaA	protein ViaA
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (129/176); total (130/176)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

TGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCGGCGCTCA  >  NC_000913/3928852‑3929047

tGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCGAgttg                                                                                                 <  1:7616400/101‑1 (MQ=255)

tGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgttg                                                                                                 >  1:7616998/1‑101 (MQ=255)

tGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAgttg                                                                                                 <  2:7616485/101‑1 (MQ=255)

   tCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGctc                                                                                              >  1:7617001/1‑101 (MQ=255)

    cAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGctct                                                                                             <  1:7616370/101‑1 (MQ=255)

       gAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGcc                                                                                          <  1:7616388/101‑1 (MQ=255)

        aaCCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCa                                                                                         >  2:7617007/1‑101 (MQ=255)

        aaCCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCa                                                                                         <  1:7616632/101‑1 (MQ=255)

           cGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGc                                                                                      <  2:7616366/101‑1 (MQ=255)

            gTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCg                                                                                     >  2:7617011/1‑101 (MQ=255)

             tCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGt                                                                                    <  1:7616293/101‑1 (MQ=255)

             tCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGt                                                                                    >  1:7617014/1‑101 (MQ=255)

             tCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGt                                                                                    <  2:7616317/101‑1 (MQ=255)

                ccGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTc                                                                                 <  1:7616713/101‑1 (MQ=255)

                 cGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCa                                                                                <  2:7616358/101‑1 (MQ=255)

                  ggTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAg                                                                               >  2:7617023/1‑101 (MQ=255)

                   gTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGt                                                                              >  1:7617025/1‑101 (MQ=255)

                    ttCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGtt                                                                             <  1:7616344/101‑1 (MQ=255)

                     tCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGtt                                                                                   >  2:7617029/1‑94 (MQ=255)

                     tCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTc                                                                            <  2:7616490/101‑1 (MQ=255)

                       gcgcACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCg                                                                          <  2:7616230/101‑1 (MQ=255)

                       gcgcACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCg                                                                          <  2:7616410/101‑1 (MQ=255)

                           aCCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTg                                                                      >  1:7617035/1‑101 (MQ=255)

                           aCCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTg                                                                      >  2:7617034/1‑101 (MQ=255)

                              aTGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTc                                                                   <  1:7616417/101‑1 (MQ=255)

                               tGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCg                                                                  >  1:7617042/1‑101 (MQ=255)

                               tGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCg                                                                  <  2:7616710/101‑1 (MQ=255)

                               tGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCg                                                                  >  2:7617041/1‑101 (MQ=255)

                                   gCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTaa                                                              <  2:7616554/101‑1 (MQ=255)

                                   gCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTaa                                                              <  1:7616552/101‑1 (MQ=255)

                                    cGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTaa                                                              <  1:7616537/100‑1 (MQ=255)

                                    cGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAg                                                             >  2:7617045/1‑101 (MQ=255)

                                    cGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAg                                                             <  1:7616565/101‑1 (MQ=255)

                                    cGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAg                                                             <  1:7616458/101‑1 (MQ=255)

                                     ggAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGt                                                                                    <  2:7616450/77‑1 (MQ=255)

                                     ggAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGa                                                            <  1:7616754/101‑1 (MQ=255)

                                        aGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAAt                                                         >  1:7617051/1‑101 (MQ=255)

                                         ggTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAAtt                                                        >  1:7617054/1‑101 (MQ=255)

                                             tCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCAcc                                                    <  1:7616484/101‑1 (MQ=255)

                                              ccATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCg                                                   >  1:7617059/1‑101 (MQ=255)

                                               cATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGt                                                  <  1:7616646/101‑1 (MQ=255)

                                                aTCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTa                                                 <  2:7616460/101‑1 (MQ=255)

                                                aTCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTa                                                 <  1:7616828/101‑1 (MQ=255)

                                                 tCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACCCCCCAGTTGCTCTGCCAg                                                                                        >  1:7617064/1‑61 (MQ=255)

                                                  cTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAtt                                               >  2:7617065/1‑101 (MQ=255)

                                                    gcgcATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTc                                             <  2:7616443/101‑1 (MQ=255)

                                                      gcATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTcac                                           <  2:7616648/101‑1 (MQ=255)

                                                       cATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTAttt                                              <  2:7616341/97‑1 (MQ=255)

                                                       cATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTcaca                                          >  2:7617069/1‑101 (MQ=255)

                                                        aTCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACaa                                         >  1:7617072/1‑101 (MQ=255)

                                                        aTCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACaa                                         <  2:7616280/101‑1 (MQ=255)

                                                         tCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAAt                                        <  1:7616084/101‑1 (MQ=255)

                                                         tCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAAt                                        <  2:7616574/101‑1 (MQ=255)

                                                          cGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATc                                       >  2:7617077/1‑101 (MQ=255)

                                                          cGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATc                                       >  1:7617075/1‑101 (MQ=255)

                                                           gTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCa                                      <  1:7616909/101‑1 (MQ=255)

                                                           gTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCa                                      <  1:7616475/101‑1 (MQ=255)

                                                           gTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCa                                      >  1:7617080/1‑101 (MQ=255)

                                                            ttGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCaa                                     >  1:7617083/1‑101 (MQ=255)

                                                            ttGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCaa                                     <  1:7616335/101‑1 (MQ=255)

                                                               cgcgGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTg                                  <  2:7616159/101‑1 (MQ=255)

                                                               cgcgGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTg                                  >  1:7617087/1‑101 (MQ=255)

                                                               cgcgGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTg                                  <  2:7616421/101‑1 (MQ=255)

                                                                gcgGTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGa                                 >  2:7617091/1‑101 (MQ=255)

                                                                  ggTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATa                               >  2:7617097/1‑101 (MQ=255)

                                                                  ggTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATa                               >  2:7617094/1‑101 (MQ=255)

                                                                   gTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAg                              <  2:7616633/101‑1 (MQ=255)

                                                                   gTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAg                              <  2:7616739/101‑1 (MQ=255)

                                                                   gTATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAg                              <  2:7616887/101‑1 (MQ=255)

                                                                    tATTGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCATTTCCGGCTg                                                                      >  1:7617099/1‑60 (MQ=255)

                                                                      ttGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCg                           <  2:7616390/101‑1 (MQ=255)

                                                                      ttGATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCg                           >  2:7617102/1‑101 (MQ=255)

                                                                        gATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGcc                         <  2:7616512/101‑1 (MQ=255)

                                                                        gATTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGcc                         <  2:7616859/101‑1 (MQ=255)

                                                                         aTTTGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCa                        >  1:7617106/1‑101 (MQ=255)

                                                                          tttGGCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCAc                       <  1:7616584/101‑1 (MQ=255)

                                                                             ggCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGtt                    <  2:7616592/101‑1 (MQ=255)

                                                                             ggCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGtt                    <  2:7616064/101‑1 (MQ=255)

                                                                              gCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTc                            >  2:7617114/1‑92 (MQ=255)

                                                                              gCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGttt                   >  2:7617111/1‑101 (MQ=255)

                                                                              gCTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGttt                   <  2:7617055/101‑1 (MQ=255)

                                                                               cTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTa                  >  2:7617117/1‑101 (MQ=255)

                                                                               cTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTa                  >  1:7617118/1‑101 (MQ=255)

                                                                               cTTCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTa                  <  2:7616814/101‑1 (MQ=255)

                                                                                ttCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTaa                 >  1:7617121/1‑101 (MQ=255)

                                                                                 tCCCGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAg                <  1:7616433/101‑1 (MQ=255)

                                                                                  cccGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGc               >  2:7617125/1‑101 (MQ=255)

                                                                                  cccGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGc               <  2:7616526/101‑1 (MQ=255)

                                                                                  cccGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGc               <  1:7616506/101‑1 (MQ=255)

                                                                                    cGAGAACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTg             <  1:7616243/101‑1 (MQ=255)

                                                                                      agaACGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGAc           >  1:7617128/1‑101 (MQ=255)

                                                                                        aaCGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCg         >  2:7617130/1‑101 (MQ=255)

                                                                                        aaCGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCg         <  1:7616114/101‑1 (MQ=255)

                                                                                        aaCGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCg         <  1:7616608/101‑1 (MQ=255)

                                                                                          cGCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCGgc       <  2:7616954/101‑1 (MQ=255)

                                                                                           gCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCGgcg      >  1:7617133/1‑101 (MQ=255)

                                                                                           gCCCCAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCGgcg      >  2:7617136/1‑101 (MQ=255)

                                                                                             cccAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCGGCGCt    >  1:7617140/1‑101 (MQ=255)

                                                                                              ccAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCGGCGCTc   <  2:7616694/101‑1 (MQ=255)

                                                                                               cAGTTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCGGCGCTCa  <  1:7616545/101‑1 (MQ=255)

TGCTCAGGAACCGTCGCCGGTTCGCGCACCATGGTGCGGAAGGTTTCCATCTGCGCATCGTTGCGCGGTATTGATTTGGCTTCCCGAGAACGCCCCAGCTGCTCTGCCAGGCGTTTCAGTTCCGGCTGTTCGTTAAGAAATTCACCGTATTTCACAATCAACTGATAGTCGCCACGTTTAAGCTGACCGGCGCTCA  >  NC_000913/3928852‑3929047

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