BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	969,798	T→C	L396P (CTC→CCC)	ycaQ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	969,798	0	T	C	99.3%	451.5 / ‑4.9	142	L396P (CTC→CCC)	ycaQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (78/63); total (78/64)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 144 total aligned reads displayed.

GCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCT  >  NC_000913/969709‑969886

gCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCtt                                                                                       >  2:672160/1‑93 (MQ=255)

 cTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCttt                                                                                      <  2:1690464/93‑1 (MQ=255)

 cTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCttt                                                                                      <  2:4065617/93‑1 (MQ=255)

   gCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTAc                                                                                    >  1:1722592/1‑93 (MQ=255)

       aaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGAt                                                                                >  1:83477/1‑93 (MQ=255)

       aaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGAt                                                                                >  1:3260850/1‑93 (MQ=255)

         aGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTg                                                                              <  2:4495031/93‑1 (MQ=255)

         aGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTg                                                                              >  2:1234616/1‑93 (MQ=255)

         aGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGAATg                                                                              >  2:3086846/1‑93 (MQ=255)

          gggCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGc                                                                             >  1:856355/1‑93 (MQ=255)

          gggCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGc                                                                             <  1:92214/93‑1 (MQ=255)

           ggCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGcc                                                                            >  1:2816251/1‑93 (MQ=255)

             cTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGc                                                                          >  1:2339058/1‑93 (MQ=255)

              tACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCa                                                                         <  2:1089383/93‑1 (MQ=255)

              tACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCa                                                                         >  1:2955499/1‑93 (MQ=255)

               aCGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCAc                                                                        >  1:2650471/1‑93 (MQ=255)

                cGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCAcc                                                                       <  1:1879166/93‑1 (MQ=255)

                cGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCAcc                                                                       <  2:2039879/93‑1 (MQ=255)

                  tCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCgg                                                                     >  1:2028536/1‑93 (MQ=255)

                  tCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCgg                                                                     >  2:4013414/1‑93 (MQ=255)

                   cAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGc                                                                    >  1:1394092/1‑93 (MQ=255)

                    aGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCt                                                                   >  1:1888114/1‑93 (MQ=255)

                    aGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCt                                                                   <  1:2438496/93‑1 (MQ=255)

                    aGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCt                                                                   >  2:3759830/1‑93 (MQ=255)

                    aGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCt                                                                   <  2:3363593/93‑1 (MQ=255)

                     ggCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTg                                                                  <  1:340574/93‑1 (MQ=255)

                     ggCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTg                                                                  >  2:831663/1‑93 (MQ=255)

                        gATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGa                                                               >  2:1440222/1‑93 (MQ=255)

                        gATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGAGGCTGCCCGCAAGGCCCCTTTACAGATTGCCGCACCGGCGGGGa                                                               >  1:3595321/1‑93 (MQ=255)

                         aTTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGaa                                                              <  2:1807980/93‑1 (MQ=255)

                           tACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAAt                                                            >  2:500730/1‑93 (MQ=255)

                           tACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAAt                                                            <  1:4411968/93‑1 (MQ=255)

                           tACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAAt                                                            >  1:44556/1‑93 (MQ=255)

                           tACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAAt                                                            <  1:2498206/93‑1 (MQ=255)

                            aCTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATa                                                           <  2:1393092/93‑1 (MQ=255)

                             cTGATTTTGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAg                                                          >  2:4197107/1‑93 (MQ=255)

                             cTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAg                                                          >  2:1513081/1‑93 (MQ=255)

                             cTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAAAAg                                                          >  1:1378262/1‑93 (MQ=255)

                              tGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGa                                                         <  1:3930074/93‑1 (MQ=255)

                              tGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGa                                                         <  2:4026016/93‑1 (MQ=255)

                               gATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAc                                                        >  1:3587685/1‑93 (MQ=255)

                                aTTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAcc                                                       >  2:4185300/1‑93 (MQ=255)

                                aTTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAcc                                                       >  2:4185299/1‑93 (MQ=255)

                                   tCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCg                                                    <  2:761432/93‑1 (MQ=255)

                                    ggCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGt                                                   <  1:1084286/92‑1 (MQ=255)

                                     gCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTc                                                  <  1:1068507/93‑1 (MQ=255)

                                       tAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGc                                                >  2:2771470/1‑93 (MQ=255)

                                          cTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCAGCTGGGAAATAGACCCCGTCGCATa                                             <  1:1976972/93‑1 (MQ=255)

                                            ggCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATaaa                                           <  2:4591943/93‑1 (MQ=255)

                                             gcagcaGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAAt                                          <  1:2852199/93‑1 (MQ=255)

                                             gcagcaGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAAt                                          <  1:529920/93‑1 (MQ=255)

                                               agcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGa                                        >  2:3529038/1‑93 (MQ=255)

                                               agcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGa                                        >  1:599702/1‑93 (MQ=255)

                                               agcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGa                                        <  2:1486969/93‑1 (MQ=255)

                                                gcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAt                                       >  1:4248091/1‑93 (MQ=255)

                                                gcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAt                                       <  1:4543055/93‑1 (MQ=255)

                                                gcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTCCCGCACCGCCTGGGAAATAGACCCCGTC‑CCTAAAAGAtt                                      >  1:332174/1‑93 (MQ=255)

                                                 cagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAtt                                      >  1:4340405/1‑93 (MQ=255)

                                                 cagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAtt                                      <  1:4646589/93‑1 (MQ=255)

                                                 cagGCAAAGCGGGTGACAGTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAtt                                      <  1:3189934/93‑1 (MQ=255)

                                                  agGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTa                                     >  1:2564673/1‑93 (MQ=255)

                                                  agGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTa                                     >  2:1411514/1‑93 (MQ=255)

                                                   ggCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTAt                                    >  2:3757769/1‑93 (MQ=255)

                                                     caccgcgGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCGCCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGa                                  <  1:3792035/90‑1 (MQ=255)

                                                     cAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGa                                  <  2:3644480/93‑1 (MQ=255)

                                                       aCGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATa                                >  2:3101520/1‑93 (MQ=255)

                                                       aCGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATa                                >  2:4106992/1‑93 (MQ=255)

                                                       aCGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATa                                >  2:183485/1‑93 (MQ=255)

                                                         gcgGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAg                              <  1:1859261/93‑1 (MQ=255)

                                                          cgGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGa                             >  2:152347/1‑93 (MQ=255)

                                                          cgGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGa                             <  1:2648870/93‑1 (MQ=255)

                                                          cgGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGa                             <  2:2209191/93‑1 (MQ=255)

                                                            ggTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAtt                           >  1:3015984/1‑93 (MQ=255)

                                                              tGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTa                         <  1:473493/93‑1 (MQ=255)

                                                                 cATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATTATTATGATAAGATTTATGc                      >  1:2766159/1‑93 (MQ=255)

                                                                 cATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGc                      >  1:3504936/1‑93 (MQ=255)

                                                                  aTTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCa                     <  2:2513370/93‑1 (MQ=255)

                                                                  aTTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCa                     >  1:3527484/1‑93 (MQ=255)

                                                                  aTTAGGACGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCa                     <  2:1596990/93‑1 (MQ=255)

                                                                   ttAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCAt                    <  1:1201542/93‑1 (MQ=255)

                                                                   ttAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCAt                    >  1:3297498/1‑93 (MQ=255)

                                                                     aGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTc                  >  2:3392461/1‑93 (MQ=255)

                                                                        aTGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATa               <  1:1648665/93‑1 (MQ=255)

                                                                         tgctgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATaa              >  1:4732614/1‑93 (MQ=255)

                                                                         tgctgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATaa              >  2:1229890/1‑93 (MQ=255)

                                                                           ctgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAAcc            <  2:1414638/93‑1 (MQ=255)

                                                                            tgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAAccc           >  1:2460201/1‑93 (MQ=255)

                                                                            tgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAAccc           <  2:1218543/93‑1 (MQ=255)

                                                                              cccGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGa         <  1:2530380/93‑1 (MQ=255)

                                                                              cccGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGa         <  1:199265/93‑1 (MQ=255)

                                                                               ccGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAAAAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGAt        >  1:1824880/1‑93 (MQ=255)

                                                                                 gCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATcc      >  1:2080934/1‑93 (MQ=255)

                                                                                 gCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATcc      <  1:4111679/93‑1 (MQ=255)

                                                                                  cAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATcct     >  1:604612/1‑92 (MQ=255)

                                                                                  cAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCa     <  1:4104436/93‑1 (MQ=255)

                                                                                    aGGCCCCTTTACTGATTGCCGCACCGTCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATc   <  1:1647365/93‑1 (MQ=255)

                                                                                    aGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAGGATTTATGCATTCATAAACCGATCCATc   >  2:3461719/1‑93 (MQ=255)

                                                                                    aGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATc   <  2:2816251/93‑1 (MQ=255)

                                                                                    aGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATc   >  1:4050238/1‑93 (MQ=255)

                                                                                     ggCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAAACCTCTCCATCt  >  1:2554165/1‑93 (MQ=255)

GCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCT  >  NC_000913/969709‑969886

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