BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,456,098	A→G	Q41Q (CAA→CAG)	yheG →	YheG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,456,098	0	A	G	99.7%	1235.6 / ‑5.0	379	Q41Q (CAA→CAG)	yheG	YheG
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (197/181); total (197/182)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

Only 100 of 392 total aligned reads displayed.

GGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGG  >  NC_000913/3456008‑3456187

ggCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCg                                                                                         >  1:4934936/1‑93 (MQ=255)

   gcTGGCACAGCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCa                                                                                      <  1:3066735/93‑1 (MQ=255)

   gcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCa                                                                                      <  2:1441240/93‑1 (MQ=255)

    cTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAg                                                                                     >  2:9425387/1‑93 (MQ=255)

     tGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAgg                                                                                    <  1:4266175/93‑1 (MQ=255)

       gCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGtt                                                                                  <  2:935461/93‑1 (MQ=255)

       gCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGtt                                                                                  >  1:2726145/1‑93 (MQ=255)

        cacaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCCTTCGAGCGGTTACTCGATATTCAGCCTCAGGTTg                                                                                 >  2:9221668/1‑93 (MQ=255)

        cacaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTg                                                                                 <  1:2168596/93‑1 (MQ=255)

         acaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTgg                                                                                >  2:6223849/1‑93 (MQ=255)

         acaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTgg                                                                                >  1:3236044/1‑93 (MQ=255)

            aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCa                                                                             >  2:4309996/1‑93 (MQ=255)

            aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCa                                                                             <  1:4533923/93‑1 (MQ=255)

            aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCa                                                                             <  2:293206/93‑1 (MQ=255)

            aCGAATTGGTTATCTATATTGAGAGGGAGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCa                                                                             <  1:1722118/93‑1 (MQ=255)

              gAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGt                                                                           >  1:8632133/1‑93 (MQ=255)

                aTTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGc                                                                         >  1:5091877/1‑93 (MQ=255)

                 ttGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTggc                                                                        >  1:6408158/1‑91 (MQ=255)

                 ttGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGcc                                                                        >  1:33595/1‑93 (MQ=255)

                  tGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCa                                                                       <  1:260644/93‑1 (MQ=255)

                    gTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGa                                                                     >  1:9710021/1‑93 (MQ=255)

                    gTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGa                                                                     >  1:3783057/1‑93 (MQ=255)

                       aTCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCa                                                                  <  2:1100191/93‑1 (MQ=255)

                        tctcTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCAt                                                                 <  1:4221287/93‑1 (MQ=255)

                         ctctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATg                                                                <  2:8957774/93‑1 (MQ=255)

                          tctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGa                                                               >  2:4184233/1‑93 (MQ=255)

                          tctATTCAGAGGGCGACATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGa                                                               <  2:8305163/93‑1 (MQ=255)

                            tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                             <  2:8873882/93‑1 (MQ=255)

                            tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                             >  1:734355/1‑93 (MQ=255)

                            tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                             <  1:5959058/93‑1 (MQ=255)

                             aTTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                            >  2:2454550/1‑93 (MQ=255)

                             aTTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                            >  1:9549838/1‑93 (MQ=255)

                             aTTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                            <  1:2073005/93‑1 (MQ=255)

                              ttCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAg                                                           <  2:4624906/93‑1 (MQ=255)

                               tCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGc                                                          >  1:15029/1‑93 (MQ=255)

                               tCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGc                                                          >  2:7153925/1‑93 (MQ=255)

                               tCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGc                                                          >  1:6583660/1‑93 (MQ=255)

                               tCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGc                                                          <  1:8478828/93‑1 (MQ=255)

                                cAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCt                                                         >  2:5241320/1‑93 (MQ=255)

                                 agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                        >  2:7454086/1‑93 (MQ=255)

                                      gCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcac                                                   >  2:3496577/1‑93 (MQ=255)

                                       cGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcaca                                                  >  1:2317118/1‑93 (MQ=255)

                                       cGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcaca                                                  <  1:10025896/93‑1 (MQ=255)

                                        gagaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAg                                                 <  1:8179086/93‑1 (MQ=255)

                                        gagaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAg                                                 >  1:1907619/1‑93 (MQ=255)

                                         agaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAgg                                                >  1:9192926/1‑93 (MQ=255)

                                          gaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGc                                               >  2:5588139/1‑93 (MQ=255)

                                          gaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGc                                               <  1:9933518/93‑1 (MQ=255)

                                           attattTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCt                                              >  2:8648654/1‑93 (MQ=255)

                                           attattTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCt                                              >  2:9810907/1‑93 (MQ=255)

                                            ttatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTg                                             >  2:6004842/1‑93 (MQ=255)

                                            ttatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTg                                             >  2:3990696/1‑93 (MQ=255)

                                             tatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGa                                            >  2:1639159/1‑93 (MQ=255)

                                              atttatCTCGCGGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAt                                           <  2:6769/93‑1 (MQ=255)

                                              atttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAt                                           >  1:9850572/1‑93 (MQ=255)

                                               tttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtt                                          <  2:2549766/93‑1 (MQ=255)

                                                ttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAttt                                         >  1:5661807/1‑93 (MQ=255)

                                                 tatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtttt                                        >  1:1096950/1‑93 (MQ=255)

                                                   tctcGCCGATACGCCATTCTAGCTGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGa                                      <  1:7019156/93‑1 (MQ=255)

                                                   tctcGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGa                                      >  1:4072624/1‑93 (MQ=255)

                                                   tatcGCCGATACGCCATTCAAGCGGTTACGCGATATTCAGCGGCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGa                                      <  2:8122175/91‑1 (MQ=255)

                                                      cGCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGc                                   <  1:1256847/93‑1 (MQ=255)

                                                       gCCGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCATGTCACAGGCTGATTTTGAAGCt                                  >  2:7247786/1‑93 (MQ=255)

                                                        ccGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTc                                 <  2:3704747/93‑1 (MQ=255)

                                                        ccGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGCTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTc                                 <  1:4420851/93‑1 (MQ=255)

                                                         cGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCg                                <  2:2427386/93‑1 (MQ=255)

                                                         cGATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCg                                <  1:5149932/93‑1 (MQ=255)

                                                          gATACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                               >  1:3958150/1‑93 (MQ=255)

                                                          gATACGCCATTCGAGCGGTTACTCGATATTCAGCGGCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                               >  1:1849210/1‑93 (MQ=255)

                                                           aTACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGc                              <  1:6327370/93‑1 (MQ=255)

                                                            tACGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCt                             >  1:7465132/1‑93 (MQ=255)

                                                             aCGCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTg                            >  1:5529361/1‑93 (MQ=255)

                                                              cGCCATTCGCGCGGTTACTCGATATTCAGCGTCATGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTgg                           <  1:1206540/93‑1 (MQ=255)

                                                               gCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGa                          >  2:5016661/1‑93 (MQ=255)

                                                               gCCATTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTATCACAGGCTGATTTTGAAGCTTGGCTGGa                          <  2:2985108/93‑1 (MQ=255)

                                                                  aTTCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGc                       >  1:7779639/1‑93 (MQ=255)

                                                                   ttCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCg                      >  2:1920915/1‑93 (MQ=255)

                                                                   ttCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCg                      >  2:6937762/1‑93 (MQ=255)

                                                                    tCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCgg                     <  2:4812230/93‑1 (MQ=255)

                                                                    tCGAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCgg                     >  2:7798727/1‑93 (MQ=255)

                                                                      gAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                   <  2:9565867/93‑1 (MQ=255)

                                                                      gAGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                   <  1:8590260/93‑1 (MQ=255)

                                                                       aGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTAt                  <  1:3478387/93‑1 (MQ=255)

                                                                       aGCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTAt                  <  2:3390237/93‑1 (MQ=255)

                                                                        gCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTAtt                 >  2:7590780/1‑93 (MQ=255)

                                                                        gCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTAtt                 <  2:5091854/93‑1 (MQ=255)

                                                                        gCGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGGAtt                 >  2:2217654/1‑93 (MQ=255)

                                                                         cGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTc                <  1:7957864/93‑1 (MQ=255)

                                                                         cGGTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTc                <  2:3548554/93‑1 (MQ=255)

                                                                          ggTTACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc               >  1:2957158/1‑93 (MQ=255)

                                                                            ttACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCAt             >  1:4757835/1‑93 (MQ=255)

                                                                            ttACTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCAt             >  2:785464/1‑93 (MQ=255)

                                                                              aCTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCa           >  1:8868020/1‑93 (MQ=255)

                                                                              aCTCGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCa           <  2:8468751/93‑1 (MQ=255)

                                                                                 cGATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGaa        <  1:4636259/93‑1 (MQ=255)

                                                                                  gATATTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAAt       >  1:2488637/1‑93 (MQ=255)

                                                                                    tatTCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAc     >  2:6560356/1‑93 (MQ=255)

                                                                                      ttCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCg   >  2:1831328/1‑93 (MQ=255)

                                                                                      ttCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCg   >  2:3094049/1‑93 (MQ=255)

                                                                                       tCAGCGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCgg  >  1:9414509/1‑93 (MQ=255)

GGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGG  >  NC_000913/3456008‑3456187

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