BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	408,017	T→C	L42P (CTG→CCG)	b0392 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	408,017	0	T	C	99.6%	1502.6 / NA	470	L42P (CTG→CCG)	b0392	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); major base C (234/231); minor base A (1/1); total (237/233)
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 500 total aligned reads displayed.

CGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAG  >  NC_000913/407927‑408109

 ggCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAgccggc                                                                                           <  2:427789/93‑1 (MQ=255)

    aTTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGc                                                                                        >  2:4484067/1‑93 (MQ=255)

     ttGCGCAACGCACATACCCCGTCGGTGTTGTTAACGCCATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCAGAGCCGGCTGcc                                                                                       >  2:2640999/1‑93 (MQ=255)

     ttGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGcc                                                                                       >  1:311519/1‑93 (MQ=255)

     ttGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGcc                                                                                       >  1:5088089/1‑93 (MQ=255)

     ttGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGcc                                                                                       >  1:6931438/1‑93 (MQ=255)

     ttGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGcc                                                                                       >  1:3164145/1‑93 (MQ=255)

     ttGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGcc                                                                                       >  2:3630953/1‑93 (MQ=255)

     ttGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGATGcc                                                                                       >  2:4945255/1‑93 (MQ=255)

       gcgcAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGa                                                                                     <  2:6643697/93‑1 (MQ=255)

       gcgcAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGa                                                                                     <  2:8211351/93‑1 (MQ=255)

       gcgcAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGa                                                                                     <  1:1383023/93‑1 (MQ=255)

         gcAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGAtg                                                                                   >  1:6784384/1‑93 (MQ=255)

          cAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGAtgt                                                                                  <  1:8449661/93‑1 (MQ=255)

          cAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGAtgt                                                                                  <  1:3414269/93‑1 (MQ=255)

            aCGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGGGa                                                                                >  1:9238754/1‑93 (MQ=255)

             cGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGaa                                                                               <  1:943784/93‑1 (MQ=255)

               cacaTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGaaaa                                                                             <  1:8080877/93‑1 (MQ=255)

               cacaTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGaaaa                                                                             <  2:1292993/93‑1 (MQ=255)

                acaTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAAc                                                                            >  1:8358155/1‑93 (MQ=255)

                acaTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAAc                                                                            <  2:1661165/93‑1 (MQ=255)

                acaTACCCCAGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAAc                                                                            >  1:8192524/1‑93 (MQ=255)

                 caTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACg                                                                           >  1:1987667/1‑93 (MQ=255)

                 caTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACg                                                                           >  2:10009492/1‑93 (MQ=255)

                 caTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACg                                                                           >  2:890276/1‑93 (MQ=255)

                  aTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGCAAACgg                                                                          >  1:1019926/1‑93 (MQ=255)

                  aTACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACgg                                                                          <  1:4373108/93‑1 (MQ=255)

                   tACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGc                                                                         >  1:7062611/1‑93 (MQ=255)

                   tACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGc                                                                         <  1:4260914/93‑1 (MQ=255)

                     ccccGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGt                                                                       <  2:9380132/93‑1 (MQ=255)

                      cccGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTg                                                                      >  1:4812487/1‑93 (MQ=255)

                        cGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGt                                                                    >  1:6191019/1‑93 (MQ=255)

                          gCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTa                                                                  >  2:6017448/1‑93 (MQ=255)

                          gCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTa                                                                  >  2:6142321/1‑93 (MQ=255)

                             ttgttgttAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTaaaa                                                               >  1:3250733/1‑93 (MQ=255)

                             ttgttgttAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTaaaa                                                               >  2:7225504/1‑93 (MQ=255)

                              tgttgttACCAGCATTGACTGAGCCTCAGTATCGGTCGCTGGCTATTAATAGTCCACAGCGGGCTGCCGATGTGAAAACGGCGTGGTTAAAag                                                              <  1:8852508/93‑1 (MQ=255)

                               gttgttAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAaga                                                             <  2:5318742/93‑1 (MQ=255)

                                ttgttAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAagag                                                            >  1:1561540/1‑93 (MQ=255)

                                      aCGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCca                                                      <  1:4919780/93‑1 (MQ=255)

                                      aCGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCca                                                      <  2:1908747/93‑1 (MQ=255)

                                        gACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCcatc                                                    <  2:1989171/93‑1 (MQ=255)

                                         aCATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCcatca                                                   <  1:1120658/93‑1 (MQ=255)

                                         aCATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCcatca                                                   <  2:2876059/93‑1 (MQ=255)

                                         aCATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCcatca                                                   <  1:2241826/93‑1 (MQ=255)

                                         aCATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCcatca                                                   <  2:9531181/93‑1 (MQ=255)

                                         aCATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCcatca                                                   >  1:6451219/1‑93 (MQ=255)

                                         aCATTGACTGAGCCTCAGAATCGGCCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCcatca                                                   <  2:3733889/93‑1 (MQ=255)

                                             tGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTAc                                               >  1:5880984/1‑93 (MQ=255)

                                              gACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTAct                                              >  2:9648438/1‑93 (MQ=255)

                                              gACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTAct                                              >  2:7836805/1‑93 (MQ=255)

                                                cTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTActct                                            <  2:693767/93‑1 (MQ=255)

                                                cTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTActct                                            <  2:1424812/93‑1 (MQ=255)

                                                cTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTActct                                            >  2:8284473/1‑93 (MQ=255)

                                                 tGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCtt                                           >  1:3502634/1‑93 (MQ=255)

                                                 tGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCtt                                           >  1:1918636/1‑93 (MQ=255)

                                                  gAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTa                                          <  1:747222/93‑1 (MQ=255)

                                                   ggcttCAGGATCGGTCGCTGCCTATTAATAAACCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTAt                                         <  2:9048291/89‑1 (MQ=255)

                                                    gCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTAtt                                        >  2:5546612/1‑93 (MQ=255)

                                                     ccTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTAttt                                       >  1:6387017/1‑93 (MQ=255)

                                                      cTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTg                                      >  2:9835517/1‑93 (MQ=255)

                                                      cTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTg                                      >  2:316217/1‑93 (MQ=255)

                                                      cTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTg                                      <  1:6031022/93‑1 (MQ=255)

                                                       tCAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGt                                     >  1:5035206/1‑93 (MQ=255)

                                                        cAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTc                                    >  2:2510825/1‑93 (MQ=255)

                                                        cAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTc                                    <  1:4049768/93‑1 (MQ=255)

                                                        cAGGATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTc                                    <  1:9788224/93‑1 (MQ=255)

                                                          ggATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAACAGAGGAGCCATCATTGCTTTTATTTGTCGa                                  >  1:3065043/1‑93 (MQ=255)

                                                          ggATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGa                                  >  1:2608432/1‑93 (MQ=255)

                                                          ggATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGa                                  >  1:5493669/1‑93 (MQ=255)

                                                          ggATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGa                                  >  2:9105259/1‑93 (MQ=255)

                                                          ggATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGa                                  >  2:4741949/1‑93 (MQ=255)

                                                          ggATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGa                                  >  2:4004747/1‑93 (MQ=255)

                                                           gATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGaa                                 <  1:7902349/93‑1 (MQ=255)

                                                           gATCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGaa                                 >  1:9076101/1‑93 (MQ=255)

                                                             tCGGTCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGaaca                               >  2:8596241/1‑93 (MQ=255)

                                                                 tCGCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAAccc                           >  2:7548654/1‑93 (MQ=255)

                                                                   gCTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGa                         >  1:5721081/1‑93 (MQ=255)

                                                                   gCTGACTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGa                         <  2:6390158/93‑1 (MQ=255)

                                                                    cTGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGAt                        >  1:269999/1‑93 (MQ=255)

                                                                     tGGCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGAtt                       >  1:2101954/1‑92 (MQ=255)

                                                                      ggCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCt                      >  2:2414653/1‑93 (MQ=255)

                                                                      ggCTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCt                      <  1:762201/93‑1 (MQ=255)

                                                                        cTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCttt                    <  2:5760541/93‑1 (MQ=255)

                                                                        cTATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCttt                    <  2:8437843/93‑1 (MQ=255)

                                                                         tATTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTc                   <  2:3341084/93‑1 (MQ=255)

                                                                          aTTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCg                  >  2:2169450/1‑93 (MQ=255)

                                                                          aTTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCg                  >  1:3722016/1‑93 (MQ=255)

                                                                          aTTAATAATCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCg                  >  1:7714952/1‑93 (MQ=255)

                                                                            taataaTCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCt                <  1:1673128/93‑1 (MQ=255)

                                                                            taataaTCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCAATACTCTTATTTGTCGAACAACACGATCTTTCGCt                >  2:5931655/1‑93 (MQ=255)

                                                                              ataatCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGt              >  1:1001698/1‑93 (MQ=255)

                                                                                aatCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTa            >  2:1117072/1‑93 (MQ=255)

                                                                                aatCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTa            <  1:9561726/93‑1 (MQ=255)

                                                                                 atCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTAc           <  2:806169/93‑1 (MQ=255)

                                                                                  tCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACg          >  2:3476120/1‑93 (MQ=255)

                                                                                  tCCACAGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACg          <  2:772395/93‑1 (MQ=255)

                                                                                    cacaGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTg        <  1:8652424/93‑1 (MQ=255)

                                                                                      caGCCGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGAt      >  1:4544251/1‑93 (MQ=255)

                                                                                          tggctcGCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTATCAAACAACCCGATCTTTCGCTGTTACGTGATTTAg  <  1:3837949/93‑1 (MQ=255)

CGGCATTGCGCAACGCACATACCCCGGCGTTGTTGTTAACGACATTGACTGAGCCTCAGGATCGGTCGCTGGCTATTAATAATCCACAGCTGGCTGCCGATGTGAAAACGGCGTGGTTAAAAGAGGATCCATCATTACTCTTATTTGTCGAACAACCCGATCTTTCGCTGTTACGTGATTTAG  >  NC_000913/407927‑408109

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