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.6%	795.3 / NA	244	L396P (CTC→CCC)	ycaQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (113/130); minor base G (0/1); total (113/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

CGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTG  >  NC_000913/969708‑969887

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

  cTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCttt                                                                                       >  1:6422073/1‑93 (MQ=255)

  cTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCttt                                                                                       >  2:3829110/1‑93 (MQ=255)

   tGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTa                                                                                      <  1:3240209/93‑1 (MQ=255)

    gCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTAc                                                                                     <  1:1134499/93‑1 (MQ=255)

      aaaaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTg                                                                                   >  1:3097310/1‑93 (MQ=255)

       aaaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCACGGCCCCTTTACTGa                                                                                  >  1:871578/1‑93 (MQ=255)

       aaaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGa                                                                                  <  1:3055278/93‑1 (MQ=255)

         aaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGAtt                                                                                <  2:3412119/93‑1 (MQ=255)

         aaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCACCTTTACTGAtt                                                                                >  2:6291642/1‑93 (MQ=255)

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

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

             gCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCg                                                                            <  2:3207629/93‑1 (MQ=255)

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

                aCGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCAc                                                                         <  1:613437/93‑1 (MQ=255)

                aCGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGACCCTTTACTGATTGCCGCAc                                                                         >  1:5127313/1‑93 (MQ=255)

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

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

                    cAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGc                                                                     <  2:551469/93‑1 (MQ=255)

                    cAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGc                                                                     <  2:1622332/93‑1 (MQ=255)

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

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

                        cGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTggg                                                                 >  1:2467955/1‑93 (MQ=255)

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

                          aTTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGaa                                                               >  1:6211347/1‑93 (MQ=255)

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

                           ttACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGaaa                                                              >  2:5921766/1‑93 (MQ=255)

                           ttACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGaaa                                                              <  1:831287/93‑1 (MQ=255)

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

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

                              cTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATaa                                                           >  1:4595000/1‑92 (MQ=255)

                              cTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAg                                                           <  2:23143/93‑1 (MQ=255)

                               tGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGa                                                          >  2:1229743/1‑93 (MQ=255)

                               tGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGa                                                          >  1:5559318/1‑93 (MQ=255)

                               tGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGa                                                          >  2:4293123/1‑93 (MQ=255)

                                gATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAc                                                         <  1:762155/93‑1 (MQ=255)

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

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

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

                                  tttCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGAccc                                                       >  1:5711223/1‑93 (MQ=255)

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

                                      gCTAGCTGACAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCGGCCCCGGCTGGGAAATAGAACCCGTc                                                   <  1:3828359/93‑1 (MQ=255)

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

                                         aaCTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCa                                                >  1:4802939/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                  cagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGAATGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAtt                                       <  2:4696583/93‑1 (MQ=255)

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

                                                    ggCAACGCGGGTGACATTAGGNTGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTAt                                     <  2:2727666/93‑1 (MQ=255)

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

                                                      cAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCGCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGa                                   <  2:130923/93‑1 (MQ=255)

                                                       aaCGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGAt                                  <  1:2627479/93‑1 (MQ=255)

                                                       aaCGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGAt                                  >  2:5592458/1‑93 (MQ=255)

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

                                                        aCGCGGGTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATa                                 <  2:4484968/93‑1 (MQ=255)

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

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

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

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

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

                                                            gggTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAt                             <  1:388590/93‑1 (MQ=255)

                                                            gggTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAt                             >  1:2282508/1‑93 (MQ=255)

                                                            gggTGACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAt                             >  2:2900559/1‑93 (MQ=255)

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

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

                                                                gACATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTAt                         <  2:700750/93‑1 (MQ=255)

                                                                 aCATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATg                        <  2:1096411/93‑1 (MQ=255)

                                                                 aCATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATg                        <  2:5993635/93‑1 (MQ=255)

                                                                  cATTAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCTGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGc                       <  2:2104206/93‑1 (MQ=255)

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

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

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

                                                                     tAGGATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCAtt                    <  2:2343996/93‑1 (MQ=255)

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

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

                                                                       ggATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCa                  <  1:3608722/93‑1 (MQ=255)

                                                                       ggATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCa                  <  2:4789415/93‑1 (MQ=255)

                                                                       ggATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCa                  <  2:4917051/93‑1 (MQ=255)

                                                                       ggATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCa                  >  2:5102492/1‑93 (MQ=255)

                                                                       ggATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCa                  <  1:1026694/93‑1 (MQ=255)

                                                                        gATGCTGCCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCAt                 <  1:2893435/93‑1 (MQ=255)

                                                                          tgctgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATaa               <  2:2935477/93‑1 (MQ=255)

                                                                          tgctgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATaa               <  2:3565163/93‑1 (MQ=255)

                                                                           gctgcCCGCAAGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAAc              <  2:2484169/93‑1 (MQ=255)

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

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

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

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

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

                                                                                    aaGGCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCAt     >  1:5900803/1‑93 (MQ=255)

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

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

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

                                                                                      ggCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCt   <  1:6653247/93‑1 (MQ=255)

                                                                                      ggCCCCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCt   >  1:2273591/1‑93 (MQ=255)

CGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTG  >  NC_000913/969708‑969887

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