BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,896,435	G→A	A196V (GCG→GTG)	b3718 ←	putative isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,896,435	0	G	A	99.4%	1238.4 / ‑5.8	350	A196V (GCG→GTG)	b3718	putative isomerase
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (205/143); total (206/144)
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 380 total aligned reads displayed.

ATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCG  >  NC_000913/3896345‑3896524

aTCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCcacg                                                                                         >  2:4269967/1‑93 (MQ=255)

aTCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCCTGCAGCACATTTTTCAGCGCCTGCGCTTTACCcacg                                                                                         >  1:4009906/1‑93 (MQ=255)

 tCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCacgc                                                                                        >  2:6606672/1‑93 (MQ=255)

 tCAGCGAAGGATGCAGCGGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCacgc                                                                                        <  1:5448748/93‑1 (MQ=255)

  cAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCAcgcc                                                                                       <  2:8239172/93‑1 (MQ=255)

  cAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACACTTTTCAGTGCCTGCGCTTTACCCAcgcc                                                                                       >  2:7802900/1‑93 (MQ=255)

    gCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCAcgccgc                                                                                     >  2:4209592/1‑93 (MQ=255)

    gCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCAcgccgc                                                                                     >  1:1803066/1‑93 (MQ=255)

     cGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCt                                                                                    >  2:1874883/1‑93 (MQ=255)

      gAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCGCGCCGCTg                                                                                   >  1:8096164/1‑93 (MQ=255)

       aaGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGa                                                                                  >  2:8792333/1‑93 (MQ=255)

       aaGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGa                                                                                  <  2:5800784/93‑1 (MQ=255)

        aGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGAc                                                                                 >  2:3136091/1‑93 (MQ=255)

         ggATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACa                                                                                >  2:263690/1‑93 (MQ=255)

            tGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAAtt                                                                             >  1:2855292/1‑93 (MQ=255)

             gCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACGTTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTa                                                                            <  2:1462817/93‑1 (MQ=255)

             gCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTa                                                                            <  1:2970957/93‑1 (MQ=255)

             gCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTa                                                                            >  1:791319/1‑93 (MQ=255)

               aGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATc                                                                          >  1:2433528/1‑93 (MQ=255)

                gCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCa                                                                         >  1:7355004/1‑93 (MQ=255)

                 cTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAg                                                                        >  1:6747523/1‑93 (MQ=255)

                  tGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAga                                                                       >  1:2782453/1‑92 (MQ=255)

                  tGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAga                                                                       >  1:6275136/1‑92 (MQ=255)

                  tGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGc                                                                       >  2:1038704/1‑93 (MQ=255)

                   gCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCa                                                                      >  1:8191299/1‑93 (MQ=255)

                   gCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCa                                                                      >  2:5632209/1‑93 (MQ=255)

                   gCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCa                                                                      >  1:3420003/1‑93 (MQ=255)

                    cAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCaa                                                                     <  1:2597360/93‑1 (MQ=255)

                     aaCACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCGTTACCCACGCCGCTGCCAATTATCAGCAAg                                                                    >  1:4263058/1‑93 (MQ=255)

                     aaCACGGAAACTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAg                                                                    >  1:7257418/1‑93 (MQ=255)

                      acacGGAAGCTGGCACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGt                                                                   >  2:3473200/1‑93 (MQ=255)

                      acacGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGt                                                                   >  1:265752/1‑93 (MQ=255)

                      acacGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGt                                                                   >  2:2342908/1‑93 (MQ=255)

                       cacGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGtt                                                                  >  2:5459761/1‑93 (MQ=255)

                          ggAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGttttt                                                               >  2:8682033/1‑93 (MQ=255)

                            aaGCTGGAACATCTTCGGTAACCGGACCTTGCAGAATATTTTTCAGTGCCTGCGCTATACCCACTCCGTTGACAATCATCAGCAATTTTTTcg                                                             >  2:8074619/1‑93 (MQ=255)

                             aGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTcgc                                                            >  1:924804/1‑93 (MQ=255)

                               cTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTcgccg                                                          <  1:7015989/93‑1 (MQ=255)

                                tGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTcgccgc                                                         >  1:6817638/1‑93 (MQ=255)

                                 ggAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCgccgcc                                                        >  2:374538/1‑93 (MQ=255)

                                  gAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCgccgcct                                                       >  1:7599287/1‑92 (MQ=255)

                                  gAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCa                                                       <  1:3235526/93‑1 (MQ=255)

                                  gAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCa                                                       <  2:6967159/93‑1 (MQ=255)

                                    aCATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCatg                                                     <  1:6121872/93‑1 (MQ=255)

                                     cATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGTCGCTGACAATTATCAGCAAGTTTTTCGCCGCCatga                                                    <  1:3811827/93‑1 (MQ=255)

                                     cATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCatga                                                    <  1:4516683/93‑1 (MQ=255)

                                       tCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCatgatg                                                  >  1:7858836/1‑93 (MQ=255)

                                       tCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCatgatg                                                  >  1:1220282/1‑93 (MQ=255)

                                         ttCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCt                                                >  2:6745142/1‑93 (MQ=255)

                                         ttCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCt                                                <  1:4296158/93‑1 (MQ=255)

                                         ttCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCt                                                <  2:4560298/93‑1 (MQ=255)

                                          tCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCtt                                               <  2:2396084/93‑1 (MQ=255)

                                           cGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCttt                                              <  2:6462407/93‑1 (MQ=255)

                                              tAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCg                                           <  1:1949193/93‑1 (MQ=255)

                                               aaCCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCgg                                          <  2:6226913/93‑1 (MQ=255)

                                               aaCCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCgg                                          >  1:8121702/1‑93 (MQ=255)

                                               aaCCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCgg                                          >  2:8495179/1‑93 (MQ=255)

                                                aCCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGc                                         >  2:6889803/1‑93 (MQ=255)

                                                 ccGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGcc                                        <  1:8275735/93‑1 (MQ=255)

                                                  cGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGccc                                       <  1:3588932/93‑1 (MQ=255)

                                                   cgACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGcccc                                      <  1:6477530/92‑1 (MQ=255)

                                                    gACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCa                                     >  2:1363952/1‑93 (MQ=255)

                                                    gACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCa                                     >  1:5035462/1‑93 (MQ=255)

                                                     aCCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCAt                                    <  1:8431601/93‑1 (MQ=255)

                                                     aCCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCAt                                    <  1:1727955/93‑1 (MQ=255)

                                                     aCCTTGCAGAACATTTTTCAGTGACTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCAt                                    <  1:1458311/93‑1 (MQ=255)

                                                      ccTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATc                                   >  1:6046941/1‑93 (MQ=255)

                                                       cTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCg                                  >  2:5071957/1‑93 (MQ=255)

                                                           cAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCGGCCATGATGCTTTTCGGCCCCATCGTCAc                              >  2:2237349/1‑93 (MQ=255)

                                                           cAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCAc                              <  1:10656/93‑1 (MQ=255)

                                                           cAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCAc                              >  2:6366554/1‑93 (MQ=255)

                                                           cAGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCAc                              >  1:7120626/1‑93 (MQ=255)

                                                            aGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCTGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACg                             >  2:1719445/1‑93 (MQ=255)

                                                            aGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACg                             <  1:8709398/93‑1 (MQ=255)

                                                            aGAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACg                             >  2:1994980/1‑93 (MQ=255)

                                                             gAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGt                            <  2:6492511/93‑1 (MQ=255)

                                                             gAACATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGt                            >  2:1684519/1‑93 (MQ=255)

                                                               aCATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAg                          >  1:6921163/1‑93 (MQ=255)

                                                               aCATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAg                          >  1:8652327/1‑93 (MQ=255)

                                                               aCATTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAg                          >  1:5928616/1‑93 (MQ=255)

                                                                 aTTTTTCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCt                        >  2:8310198/1‑93 (MQ=255)

                                                                   ttttCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGt                      >  2:7106762/1‑93 (MQ=255)

                                                                    tttCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTc                     <  1:8541007/93‑1 (MQ=255)

                                                                      tCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTCGCTGTCAg                   >  1:1113986/1‑93 (MQ=255)

                                                                      tCAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAg                   <  1:4834501/93‑1 (MQ=255)

                                                                       cAGTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAgg                  <  2:3321824/93‑1 (MQ=255)

                                                                         gTGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGcc                >  1:8766617/1‑92 (MQ=255)

                                                                          tGCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCAc               <  2:1171648/93‑1 (MQ=255)

                                                                           gCCTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACt              >  1:6646066/1‑93 (MQ=255)

                                                                             cTGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTaa            >  1:4083512/1‑93 (MQ=255)

                                                                              tGCGCTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAAc           >  2:4796649/1‑93 (MQ=255)

                                                                               gcgcTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACg          >  1:2243536/1‑93 (MQ=255)

                                                                                cgcTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGa         >  1:8551450/1‑93 (MQ=255)

                                                                                 gcTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGaa        >  2:2996210/1‑93 (MQ=255)

                                                                                  cTTTACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGaaa       <  1:5240233/93‑1 (MQ=255)

                                                                                    ttACCCACGCCGCTGACAATTGTCAGCAAGTTTTTCGCCGCCATGCTGCTTTTCGGCCCCATCGTCACGTAGCTGTCAAGCACTAACGAAAAg     >  2:4683681/1‑93 (MQ=255)

                                                                                    ttACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCCTCGTCACGTAGCTGTCAGGCACTAACGAAAAg     <  2:7324963/93‑1 (MQ=255)

                                                                                    ttACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAg     <  1:639516/93‑1 (MQ=255)

                                                                                    ttACCCACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAg     >  1:5363376/1‑93 (MQ=255)

                                                                                       cccACGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGGcg  >  2:2536286/1‑93 (MQ=255)

ATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCAGTGCCTGCGCTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCG  >  NC_000913/3896345‑3896524

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