BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,354,373	G→A	G254G (GGC→GGT)	sapA ←	SapA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,354,373	0	G	A	98.8%	292.6 / NA	85	G254G (GGC→GGT)	sapA	SapA
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (34/50); minor base C (0/1); total (34/51)
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

TCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCGGC  >  NC_000913/1354284‑1354462

tcgtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGccaccg                                                                                        >  1:4367911/1‑93 (MQ=255)

tcgtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGccaccg                                                                                        >  1:1415720/1‑93 (MQ=255)

 cgtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCgg                                                                                       >  1:3609873/1‑93 (MQ=255)

 cgtcACGCAAAATGGATAGCTGGATGGCAGCAGGCCAGGCCGGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGCTGCCACCgg                                                                                       <  2:758452/93‑1 (MQ=255)

  gtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGa                                                                                      >  1:3186375/1‑93 (MQ=255)

  gtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGa                                                                                      >  1:394909/1‑93 (MQ=255)

  gtcACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGa                                                                                      >  2:1704554/1‑93 (MQ=255)

    cACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGc                                                                                    <  1:1322560/93‑1 (MQ=255)

    cACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGc                                                                                    >  1:61007/1‑93 (MQ=255)

    cACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGc                                                                                    <  1:1157877/93‑1 (MQ=255)

     aCGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGCGcc                                                                                   >  1:1944743/1‑93 (MQ=255)

     aCGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGcc                                                                                   <  2:3479734/93‑1 (MQ=255)

         aaaaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTaaa                                                                               <  1:249624/93‑1 (MQ=255)

         aaaaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTaaa                                                                               <  1:1544863/93‑1 (MQ=255)

          aaaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAAt                                                                              <  1:1333566/93‑1 (MQ=255)

           aaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGCCGTCCGGTCGCACCGGAGCCTAAATc                                                                             >  1:4380954/1‑93 (MQ=255)

           aaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATc                                                                             <  2:1652412/93‑1 (MQ=255)

           aaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATc                                                                             <  2:2235821/93‑1 (MQ=255)

           aaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATc                                                                             >  1:4435076/1‑93 (MQ=255)

           aaTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATc                                                                             >  1:2130758/1‑93 (MQ=255)

           aaTGGATAGCTGGCTGCCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATc                                                                             <  2:1440327/93‑1 (MQ=255)

            aTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCCCCGGAGCCTAAATcc                                                                            <  2:243182/93‑1 (MQ=255)

            aTGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATcc                                                                            >  1:1905891/1‑93 (MQ=255)

              ggATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATccac                                                                          >  1:883103/1‑93 (MQ=255)

               gATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATccacc                                                                         <  2:3370327/93‑1 (MQ=255)

                 tAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCcaccac                                                                       <  1:3714183/93‑1 (MQ=255)

                 tAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCcaccac                                                                       <  2:3487076/93‑1 (MQ=255)

                 tAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCcaccac                                                                       <  1:4550679/93‑1 (MQ=255)

                    ctggctggCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTAc                                                                    <  1:647055/93‑1 (MQ=255)

                    ctggctggCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTAc                                                                    >  1:4086514/1‑93 (MQ=255)

                     tggctggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTAcc                                                                   >  1:4493908/1‑93 (MQ=255)

                      ggctggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCt                                                                  >  2:2383169/1‑93 (MQ=255)

                      ggctggcAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCt                                                                  >  1:4411851/1‑93 (MQ=255)

                             agcagGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCAt                                                           >  1:4371740/1‑93 (MQ=255)

                              gcagGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTGCCTGCGGCAtt                                                          <  2:3922893/93‑1 (MQ=255)

                              gcagGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCAtt                                                          >  1:160928/1‑93 (MQ=255)

                              gcagGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCAtt                                                          >  2:3739670/1‑93 (MQ=255)

                                aggccaggccAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTaa                                                        >  2:2271048/1‑93 (MQ=255)

                                aggccaggccAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTaa                                                        <  2:3806267/93‑1 (MQ=255)

                                 ggccaggccaGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAAc                                                       <  2:4219104/93‑1 (MQ=255)

                                 ggccaggccaCAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAAc                                                       <  2:42326/93‑1 (MQ=255)

                                  gccagtccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACg                                                      >  1:3515044/1‑93 (MQ=255)

                                  gccaggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACg                                                      <  2:3690252/93‑1 (MQ=255)

                                  gccaggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACg                                                      <  1:180607/93‑1 (MQ=255)

                                      ggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGAGGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGttt                                                  >  1:399/1‑93 (MQ=255)

                                      ggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGttt                                                  <  1:2956341/93‑1 (MQ=255)

                                      ggccagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGttt                                                  >  2:1905209/1‑93 (MQ=255)

                                      ggccagAACGTCGCATTCCCAGGACAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCATAAATCCACCACTACCAGCGGCATTAACAGttt                                                  >  2:4566160/1‑93 (MQ=255)

                                       gccagAACTCCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTa                                                 <  2:4217943/93‑1 (MQ=255)

                                         cagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAcc                                               <  2:2625535/93‑1 (MQ=255)

                                         cagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAcc                                               >  2:44544/1‑93 (MQ=255)

                                         cagAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAcc                                               >  2:254287/1‑93 (MQ=255)

                                         cagAACGTCGCATCCCCCGGGCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAcc                                               <  1:1401905/93‑1 (MQ=255)

                                         cagAACGCCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAcc                                               <  1:1033286/93‑1 (MQ=255)

                                          agAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTAccc                                              >  1:2003467/1‑92 (MQ=255)

                                           gAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACcgc                                             >  2:4250166/1‑93 (MQ=255)

                                           gAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACcgc                                             >  1:3382917/1‑93 (MQ=255)

                                             aCGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCgcgc                                           >  2:3847827/1‑93 (MQ=255)

                                              cGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGcc                                          <  1:2275808/93‑1 (MQ=255)

                                                 cGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGAATTAACGGTTTACCGCGCCAGa                                       <  2:2518023/93‑1 (MQ=255)

                                                  gCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGaa                                      <  1:4433597/93‑1 (MQ=255)

                                                     ttCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtc                                   >  2:167009/1‑93 (MQ=255)

                                                     ttCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtc                                   >  1:1634616/1‑93 (MQ=255)

                                                      tCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtca                                  <  1:1317966/93‑1 (MQ=255)

                                                         ccGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtcatca                               <  1:1184731/93‑1 (MQ=255)

                                                         ccGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGtcatca                               <  1:1102637/93‑1 (MQ=255)

                                                          cGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTcatcat                              <  2:3519137/93‑1 (MQ=255)

                                                          cGGTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTcatcat                              <  2:2547939/93‑1 (MQ=255)

                                                           ggTCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATg                             >  2:1628629/1‑93 (MQ=255)

                                                             tCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGAc                           >  2:579768/1‑93 (MQ=255)

                                                             tCAGGAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGAc                           <  1:2599590/93‑1 (MQ=255)

                                                                ggAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGtt                        <  1:1135181/93‑1 (MQ=255)

                                                                ggAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGtt                        <  1:2298927/93‑1 (MQ=255)

                                                                ggAGTTTCGACAGACGTCCAGTACCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGtt                        <  1:2735834/93‑1 (MQ=255)

                                                                 gAGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTg                       <  1:634431/93‑1 (MQ=255)

                                                                  aGTTTCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGt                      <  1:3649495/93‑1 (MQ=255)

                                                                   gTTTCGTCAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTa                     <  1:3620262/93‑1 (MQ=255)

                                                                      tCGACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGc                  <  2:2646249/93‑1 (MQ=255)

                                                                        gACAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGa                >  2:2165766/1‑93 (MQ=255)

                                                                         aCAGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGCTGTAGGCGaa               >  2:266037/1‑93 (MQ=255)

                                                                           aGACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATa             <  2:4113135/93‑1 (MQ=255)

                                                                            gACGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATaa            <  1:2174955/93‑1 (MQ=255)

                                                                              cGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACTGTTTACAGCCCCAGAAGTCATCATGACGTTGTAGGCGAATAAAt          >  2:1963597/1‑93 (MQ=255)

                                                                              cGTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAAt          <  1:1207928/93‑1 (MQ=255)

                                                                               gTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAAtt         <  2:488217/93‑1 (MQ=255)

                                                                               gTCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAAtt         <  1:2777126/93‑1 (MQ=255)

                                                                                tCCGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTg        <  2:4315051/93‑1 (MQ=255)

                                                                                 ccGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGc       <  1:1436623/93‑1 (MQ=255)

                                                                                  cGGTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGcc      >  2:623071/1‑93 (MQ=255)

                                                                                   ggTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGccc     <  1:1110132/93‑1 (MQ=255)

                                                                                   ggTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGccc     <  1:4420852/93‑1 (MQ=255)

                                                                                    gTGCCACCGGAGCCTAAATCCACCACTGCCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCg    <  2:1435311/93‑1 (MQ=255)

                                                                                    gTGCCACCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCg    <  2:1338621/93‑1 (MQ=255)

                                                                                      gccaccGGAGCCTAAATCCACCACTACCTGCGGCATTAATGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCGgc  <  2:1901235/93‑1 (MQ=255)

TCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCCGGAGCCTAAATCCACCACTACCTGCGGCATTAACGGTTTACCGCGCCAGAAGTCATCATGACGTTGTAGGCGAATAAATTGCCCGGC  >  NC_000913/1354284‑1354462

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