BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,560,887	T→C	V115A (GTG→GCG)	b2447 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,560,887	0	T	C	99.6%	818.3 / ‑5.3	252	V115A (GTG→GCG)	b2447	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (119/132); total (119/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 256 total aligned reads displayed.

ATCCAGAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGTGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGACATGGAACGGC  >  NC_000913/2560797‑2560976

 tCCAGAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCt                                                                                        <  1:1118074/93‑1 (MQ=255)

 tCCAGAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCt                                                                                        <  2:765074/93‑1 (MQ=255)

 tCCAGAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCt                                                                                        <  1:6712625/93‑1 (MQ=255)

 tCCAGAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCt                                                                                        >  1:4409789/1‑93 (MQ=255)

    aGAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGc                                                                                     >  1:3395786/1‑93 (MQ=255)

     gAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAACCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGcc                                                                                    <  1:5244158/93‑1 (MQ=255)

     gAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGcc                                                                                    >  1:5320557/1‑93 (MQ=255)

     gAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGcc                                                                                    >  1:1567341/1‑93 (MQ=255)

     gAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGcc                                                                                    <  1:3757741/93‑1 (MQ=255)

      aaTATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAACCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCt                                                                                   <  1:2791633/93‑1 (MQ=255)

      aaTATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCt                                                                                   <  1:1865945/93‑1 (MQ=255)

      aaTATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCt                                                                                   <  2:5931133/93‑1 (MQ=255)

       atatGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTg                                                                                  >  2:4624834/1‑93 (MQ=255)

       atatGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTg                                                                                  <  1:324825/93‑1 (MQ=255)

       atatGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTg                                                                                  <  1:3263097/93‑1 (MQ=255)

       atatGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTg                                                                                  <  2:4451805/93‑1 (MQ=255)

       atatGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTg                                                                                  <  2:4402342/93‑1 (MQ=255)

        tatGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTgg                                                                                 >  1:766259/1‑93 (MQ=255)

          tGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTgggg                                                                               >  1:739836/1‑92 (MQ=255)

            cGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGa                                                                             <  2:4147598/93‑1 (MQ=255)

              aaGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGc                                                                           >  1:4008587/1‑93 (MQ=255)

               aGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCt                                                                          >  1:2135781/1‑93 (MQ=255)

               aGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCt                                                                          >  2:4238996/1‑93 (MQ=255)

                gCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCtt                                                                         >  1:5350033/1‑93 (MQ=255)

                 cACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTa                                                                        <  1:1450976/93‑1 (MQ=255)

                  aCTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTAt                                                                       <  1:421529/93‑1 (MQ=255)

                   cTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATc                                                                      >  2:3545305/1‑93 (MQ=255)

                   cTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATc                                                                      <  2:3707217/93‑1 (MQ=255)

                   cTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATc                                                                      >  2:3945814/1‑93 (MQ=255)

                    tGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCg                                                                     >  1:5840758/1‑93 (MQ=255)

                     ggAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGc                                                                    >  2:5837877/1‑93 (MQ=255)

                     ggAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGc                                                                    >  2:2474374/1‑93 (MQ=255)

                       aaCATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGcac                                                                  <  2:6593025/93‑1 (MQ=255)

                        aCATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGcaca                                                                 >  2:4219357/1‑93 (MQ=255)

                        aCATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGcaca                                                                 >  1:6601918/1‑93 (MQ=255)

                            gAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGa                                                             >  1:4369814/1‑93 (MQ=255)

                             aaaGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGac                                                            <  2:6382485/93‑1 (MQ=255)

                             aaaGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGac                                                            <  2:3396918/93‑1 (MQ=255)

                             aaaGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGac                                                            >  2:1368847/1‑93 (MQ=255)

                              aaGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTTGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGaca                                                           >  1:5381895/1‑93 (MQ=255)

                              aaGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGaca                                                           <  1:2545518/93‑1 (MQ=255)

                              aaGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGaca                                                           >  1:2433016/1‑93 (MQ=255)

                                 gCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCt                                                        <  2:2417343/93‑1 (MQ=255)

                                 gCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCt                                                        <  2:608169/93‑1 (MQ=255)

                                  cGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCtt                                                       <  2:3293309/93‑1 (MQ=255)

                                   gTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTc                                                      <  1:1794387/93‑1 (MQ=255)

                                   gTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTc                                                      >  1:2470272/1‑93 (MQ=255)

                                    tCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCa                                                     >  2:2302271/1‑93 (MQ=255)

                                    tCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCa                                                     >  2:835563/1‑93 (MQ=255)

                                     cATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCAc                                                    >  1:1030047/1‑93 (MQ=255)

                                     cATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCAc                                                    <  2:3581327/93‑1 (MQ=255)

                                       tgaatgaCACACCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACgg                                                  <  2:1320376/93‑1 (MQ=255)

                                         aatgaCAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGt                                                >  2:485869/1‑93 (MQ=255)

                                         aatgaCAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGt                                                >  2:2601546/1‑93 (MQ=255)

                                          atgaCAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTg                                               >  2:2077489/1‑93 (MQ=255)

                                          atgaCAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTg                                               >  1:3290527/1‑93 (MQ=255)

                                           tgaCAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGc                                              >  1:3799777/1‑93 (MQ=255)

                                             aCAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGa                                            <  2:588388/93‑1 (MQ=255)

                                               aaaaCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGaaa                                          <  1:4785555/93‑1 (MQ=255)

                                               aaaaCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGaaa                                          <  2:6219396/93‑1 (MQ=255)

                                               aaaaCCCGCTTCCACGCCGACAAGTGGGGAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGaac                                          >  1:5458632/1‑92 (MQ=255)

                                                aaaCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAAc                                         >  2:2709342/1‑93 (MQ=255)

                                                aaaCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAAc                                         >  2:5192838/1‑93 (MQ=255)

                                                aaaCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAAc                                         <  1:6234434/93‑1 (MQ=255)

                                                aaaCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAAc                                         >  2:714989/1‑93 (MQ=255)

                                                 aaCCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACg                                        >  1:4500610/1‑93 (MQ=255)

                                                   cccGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGt                                      >  1:1932189/1‑93 (MQ=255)

                                                    ccGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGtt                                     >  2:5612832/1‑93 (MQ=255)

                                                     cGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGttt                                    <  2:4800840/93‑1 (MQ=255)

                                                      gCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTc                                   <  1:5040273/93‑1 (MQ=255)

                                                        ttCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCaa                                 <  1:4208185/93‑1 (MQ=255)

                                                         tCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAAt                                >  1:6432920/1‑93 (MQ=255)

                                                         tCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAAt                                <  1:4891642/93‑1 (MQ=255)

                                                         tCCACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAAt                                >  2:5396628/1‑93 (MQ=255)

                                                          ccACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATa                               <  1:4560633/93‑1 (MQ=255)

                                                          ccACGCCGACAAGTGGGTAGGGTACAGCAAGGCGC‑TGCCTGGGTGAGCTTATCGCACAAAGAAAACTTCACGGGGGCGAAACGGTTTCCATAc                              >  2:4323404/1‑93 (MQ=255)

                                                           cACGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATAc                              <  2:3836315/93‑1 (MQ=255)

                                                            aCGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACt                             <  2:1093668/93‑1 (MQ=255)

                                                            aCGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACt                             <  1:3017808/93‑1 (MQ=255)

                                                            aCGCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACt                             <  2:1686725/93‑1 (MQ=255)

                                                              gCCGACAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTc                           <  1:6343904/93‑1 (MQ=255)

                                                              gCCGACAAGTGGGTAGGGGACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTc                           <  2:641445/93‑1 (MQ=255)

                                                                   cAAGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAAc                      <  2:6542712/93‑1 (MQ=255)

                                                                   cAAGTGGCTAGAGTACAGCAAGGCGCTTGCCTGGGTGATCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAAc                      <  1:4656530/93‑1 (MQ=255)

                                                                    aaGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACg                     <  1:1501569/93‑1 (MQ=255)

                                                                     aGTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGa                    >  2:4570777/1‑93 (MQ=255)

                                                                      gTGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGAt                   <  2:2798312/93‑1 (MQ=255)

                                                                       tGGGTAGGGTACAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATg                  <  2:6305063/93‑1 (MQ=255)

                                                                        gggtagggtaCAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGt                 <  2:1599861/93‑1 (MQ=255)

                                                                        gggtagggtaCAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGt                 <  2:2939926/93‑1 (MQ=255)

                                                                            agggtaCAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGa             >  1:2243820/1‑93 (MQ=255)

                                                                            agggtaCAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGGTTCAATACTTCCAAACGATGTACGa             >  2:1488961/1‑93 (MQ=255)

                                                                              ggtaCAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAATGATGTACGACa           >  1:2896347/1‑93 (MQ=255)

                                                                                 aCAGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGACATgg        <  1:2600814/93‑1 (MQ=255)

                                                                                 aCAGCAAGGCGCTTGCCTGGGTGACCTTATCGCACAAAGACGCCTCCCCGGGTCCGAAACGGTTTCAATACTTCCAAGCGATGTACGACATgg        <  1:596613/93‑1 (MQ=255)

                                                                                   aGCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGACATGGaa      >  2:4734642/1‑93 (MQ=255)

                                                                                    gCAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGACATGGAAc     <  1:1393384/93‑1 (MQ=255)

                                                                                     cAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGACATGGAACg    >  2:551270/1‑93 (MQ=255)

                                                                                     cAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGACATGGAACg    <  2:3142004/93‑1 (MQ=255)

                                                                                     cAAGGCGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGACATGGAACg    >  2:4948242/1‑93 (MQ=255)

ATCCAGAATATGCGAAGCACTGGAACATGAAAGGCGTCATGAATGACAAAACCCGCTTCCACGCCGACAAGTGGGTAGGGTACAGCAAGGTGCTTGCCTGGGTGAGCTTATCGCACAAAGACACCTTCACGGGTGCGAAACGGTTTCAATACTTCCAAACGATGTACGACATGGAACGGC  >  NC_000913/2560797‑2560976

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: `‑`