BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)_6→7	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	96.9%	348.3 / 4.5	98	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (2/1); new base C (46/49); total (48/50)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 116 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                   <  1:56069/100‑1 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c                                                                                                >  1:2703975/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑cc                                                                                               >  1:744844/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑cc                                                                                               >  2:803938/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑cc                                                                                               <  2:292757/100‑1 (MQ=255)

    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑cccn                                                                                             <  2:588653/100‑2 (MQ=255)

    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑cccc                                                                                             <  1:514947/100‑1 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCACATTCTTCCCACCTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑cccccc                                                                                           <  2:1132249/100‑1 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑caccc                                                                                            >  2:2400565/1‑99 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccccc                                                                                          <  1:1639033/100‑2 (MQ=255)

        tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccccct                                                                                         >  1:1317331/1‑98 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccccctgt                                                                                       >  1:977597/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccccctgt                                                                                       <  1:426094/100‑5 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccccctgtt                                                                                      >  2:2198660/1‑95 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑ccccccctgt                                                                                       >  1:1129448/1‑95 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTc                                                                                     >  2:1960750/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTc                                                                                     <  1:2978775/100‑1 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTc                                                                                     <  1:2548411/100‑1 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTc                                                                                     >  1:246757/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTc                                                                                     <  2:1562773/100‑1 (MQ=255)

             ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCt                                                                                    <  1:1235570/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCt                                                                                    <  2:1191119/99‑1 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACAGAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGa                                                                                 <  1:3020859/100‑1 (MQ=255)

                   gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGaaaa                                                                              >  1:1663205/1‑100 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATg                                                                            <  2:458108/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGAc                                                                          <  1:2034896/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGAc                                                                          <  1:1319008/100‑1 (MQ=255)

                        ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCCC‑TGTTCTGGAAAATGACTa                                                                        >  1:1297834/1‑100 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTa                                                                        <  2:426058/100‑1 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTa                                                                        >  1:139943/1‑100 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAg                                                                       <  2:744844/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGcc                                                                     >  1:2505197/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGcc                                                                     >  2:3116296/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGcc                                                                     <  2:3339219/100‑1 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:2661287/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:3234406/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTACTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:2061976/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:1983505/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:1916248/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  2:237387/1‑100 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             <  2:2665609/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:216567/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  2:2258239/1‑100 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       <  1:1996183/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  2:1919921/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  2:175507/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1351188/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:2871213/100‑1 (MQ=255)

                                                 caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  1:2416055/100‑1 (MQ=255)

                                                  aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGATAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  2:394674/99‑1 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              >  1:1212586/1‑100 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             >  1:2999897/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:656751/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  1:3047446/1‑100 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:3261211/100‑1 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  2:3288501/1‑100 (MQ=255)

                                                         tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        <  2:3234393/100‑1 (MQ=255)

                                                          cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  1:803938/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  2:2358213/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:2167445/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:20607/100‑1 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGGAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  1:2403399/1‑100 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTCGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  2:1150253/1‑100 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  2:782627/1‑100 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    <  2:1136545/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:175507/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  1:3255662/1‑100 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  <  1:3337715/100‑1 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  2:1916244/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:1989198/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:1432018/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:246756/100‑1 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:3110774/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:1598854/1‑100 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           >  1:1412989/1‑100 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          >  2:467931/1‑100 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         >  1:2445940/1‑100 (MQ=255)

                                                                         ttGCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        >  1:1451684/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  2:2696684/1‑100 (MQ=255)

                                                                            naaaGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGTTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:1496515/2‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  1:1228158/1‑100 (MQ=255)

                                                                               aGGAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  <  1:1598857/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:3036301/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  1:1998665/100‑1 (MQ=255)

                                                                                 gagaAGACTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                <  2:248049/100‑1 (MQ=255)

                                                                                   gaagaCTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  2:3255649/100‑1 (MQ=255)

                                                                                   gaagaCTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  2:886160/1‑100 (MQ=255)

                                                                                   gaagaCTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  2:1870341/100‑1 (MQ=255)

                                                                                     agaCTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  1:3168875/100‑1 (MQ=255)

                                                                                     agaCTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATACTGATTCGGCAGGTGATTg            <  2:1081345/100‑1 (MQ=255)

                                                                                      gaCTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  1:3265354/1‑100 (MQ=255)

                                                                                        cTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:1655472/99‑1 (MQ=255)

                                                                                        cTATGTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  1:357328/100‑1 (MQ=255)

                                                                                          aTGTTGAAAT‑‑CCCCCCATGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       >  2:1394340/1‑100 (MQ=255)

                                                                                            gTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:3082367/1‑100 (MQ=255)

                                                                                            gTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:2736768/1‑100 (MQ=255)

                                                                                            gTTGAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:1665731/1‑100 (MQ=255)

                                                                                               gAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:1960754/100‑1 (MQ=255)

                                                                                               gAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  1:1676306/1‑100 (MQ=255)

                                                                                               gAAATC‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:1488936/100‑1 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑CCCCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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