BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,905,202	T→C	D2G (GAC→GGC)	phoU ←	PhoU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,905,202	0	T	C	98.9%	1006.6 / ‑2.1	283	D2G (GAC→GGC)	phoU	PhoU
Reads supporting (aligned to +/- strand): ref base T (3/0); new base C (143/135); total (148/135)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 293 total aligned reads displayed.

ACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGT  >  NC_000913/3905112‑3905291

accaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGccc                                                                                         >  1:4175263/1‑93 (MQ=255)

 ccaTGCCGACCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCa                                                                                        >  1:1629658/1‑93 (MQ=255)

  caTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCAt                                                                                       <  2:5162707/93‑1 (MQ=255)

   aTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATa                                                                                      <  2:4857973/93‑1 (MQ=255)

    tGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATaa                                                                                     <  1:2420896/93‑1 (MQ=255)

    tGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATaa                                                                                     <  1:4681278/93‑1 (MQ=255)

    tGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATaa                                                                                     >  1:710718/1‑93 (MQ=255)

    tGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATaa                                                                                     <  2:3252519/93‑1 (MQ=255)

    tGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATaa                                                                                     <  2:4832775/93‑1 (MQ=255)

        gccCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCa                                                                                 <  2:552984/93‑1 (MQ=255)

          ccATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACt                                                                               >  1:5391134/1‑93 (MQ=255)

           cATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTc                                                                              >  1:4270313/1‑93 (MQ=255)

           cATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTc                                                                              >  2:5400236/1‑93 (MQ=255)

            aTGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTcc                                                                             >  1:2963986/1‑93 (MQ=255)

             tGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCt                                                                            <  2:5222137/93‑1 (MQ=255)

             tGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCt                                                                            >  1:2528249/1‑93 (MQ=255)

             tGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCt                                                                            >  1:2858906/1‑93 (MQ=255)

              ggTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTg                                                                           >  2:4598725/1‑93 (MQ=255)

               gTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGa                                                                          <  1:2766926/93‑1 (MQ=255)

                tcatcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGaa                                                                         >  2:699971/1‑93 (MQ=255)

                tcatcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGaa                                                                         <  2:1830952/93‑1 (MQ=255)

                  atcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATc                                                                       <  1:4812625/93‑1 (MQ=255)

                  atcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATc                                                                       >  2:2736111/1‑93 (MQ=255)

                  atcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATc                                                                       >  1:2099234/1‑93 (MQ=255)

                    caCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCaa                                                                     <  1:3224694/93‑1 (MQ=255)

                    caCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCaa                                                                     >  2:5331607/1‑93 (MQ=255)

                     aCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAAc                                                                    >  2:2515481/1‑93 (MQ=255)

                       cTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCg                                                                  <  2:875594/93‑1 (MQ=255)

                          cgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTaa                                                               <  2:2393855/93‑1 (MQ=255)

                          cgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTaa                                                               <  1:1105533/93‑1 (MQ=255)

                          cgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTaa                                                               <  2:4393197/93‑1 (MQ=255)

                            tgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACg                                                             >  1:1057593/1‑93 (MQ=255)

                            tgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACg                                                             >  2:2602037/1‑93 (MQ=255)

                               ggATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGAcc                                                          <  1:358791/93‑1 (MQ=255)

                                gATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCg                                                         >  1:2266742/1‑93 (MQ=255)

                                gATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCg                                                         <  2:1023970/93‑1 (MQ=255)

                                gATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCg                                                         >  1:513573/1‑93 (MQ=255)

                                  tACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGt                                                       <  2:1966588/93‑1 (MQ=255)

                                   aCTTTCCAGTTAGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTg                                                      >  1:1336291/1‑93 (MQ=255)

                                     tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGAt                                                    >  1:588979/1‑93 (MQ=255)

                                     tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGAt                                                    <  1:5833603/93‑1 (MQ=255)

                                      ttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATg                                                   >  1:1802391/1‑93 (MQ=255)

                                      gtCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATg                                                   <  1:5460314/92‑1 (MQ=255)

                                       tCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGt                                                  <  1:6144092/93‑1 (MQ=255)

                                        ccAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTa                                                 >  1:3778625/1‑93 (MQ=255)

                                        ccAGGTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTa                                                 <  1:3433128/93‑1 (MQ=255)

                                         cAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAg                                                >  2:4647466/1‑93 (MQ=255)

                                          aGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGt                                               >  2:1412885/1‑93 (MQ=255)

                                           gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtc                                              <  2:78556/93‑1 (MQ=255)

                                           gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtc                                              >  2:1574386/1‑93 (MQ=255)

                                           gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtc                                              >  2:5798127/1‑93 (MQ=255)

                                           gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtc                                              >  2:3692860/1‑93 (MQ=255)

                                             tCGGCGTTGAACTTACCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtctt                                            >  2:2825904/1‑93 (MQ=255)

                                              cGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtcttc                                           <  2:3048128/93‑1 (MQ=255)

                                              cGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtcttc                                           <  1:3983901/93‑1 (MQ=255)

                                              cGGCGTTGAACTGACCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtcttc                                           <  1:59931/93‑1 (MQ=255)

                                               ggCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtcttct                                          <  2:6092480/93‑1 (MQ=255)

                                               ggCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtcttct                                          >  1:5518572/1‑93 (MQ=255)

                                                gCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtg                                         >  2:4007282/1‑93 (MQ=255)

                                                 cGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgt                                        >  1:3171860/1‑93 (MQ=255)

                                                  gTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtt                                       <  1:2927078/93‑1 (MQ=255)

                                                  gTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtt                                       >  1:385106/1‑93 (MQ=255)

                                                   ttGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgttt                                      >  1:1871901/1‑93 (MQ=255)

                                                   ttGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgttt                                      >  1:3560099/1‑93 (MQ=255)

                                                    tGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttg                                     <  2:2211515/93‑1 (MQ=255)

                                                       aCTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTgtttgttt                                  <  2:2151829/93‑1 (MQ=255)

                                                       aCTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTgtttgttt                                  <  2:2075917/93‑1 (MQ=255)

                                                         tGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttct                                <  1:333351/93‑1 (MQ=255)

                                                          ggCCGGAAATATGTTTATTAAGATTGAGGCTGCCCATAACGCGCTCCTGAATCAACCGTAACGGCCGGTGATGTAGTCTTCTGTTTGTttctt                               <  2:3372805/93‑1 (MQ=255)

                                                            ccGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCg                             <  1:4437675/93‑1 (MQ=255)

                                                            ccGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCg                             <  2:4910753/93‑1 (MQ=255)

                                                             ctgAAATATGTTTATTAAGATTGAGACTGCACATAACGCACTCCTGAATCAGCCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGc                            <  1:1234892/91‑1 (MQ=255)

                                                             cGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGc                            >  2:2965217/1‑93 (MQ=255)

                                                             cGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGc                            >  1:2195211/1‑93 (MQ=255)

                                                              ggAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCt                           >  1:1151648/1‑93 (MQ=255)

                                                               gAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTg                          <  2:2583229/93‑1 (MQ=255)

                                                               gAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTg                          >  1:1529498/1‑93 (MQ=255)

                                                                aaaTATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTgg                         >  1:4768613/1‑93 (MQ=255)

                                                                  atatGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCt                       >  2:4174007/1‑93 (MQ=255)

                                                                   tatGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCtt                      <  1:10392/93‑1 (MQ=255)

                                                                     tGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTgg                    >  1:4123143/1‑93 (MQ=255)

                                                                     tGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTgg                    >  1:4945451/1‑93 (MQ=255)

                                                                     tGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGCAGTCTTCTGTTTGTTTCTTCTGTGGCTTgg                    >  2:126388/1‑93 (MQ=255)

                                                                       tttATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTg                  >  2:4719647/1‑93 (MQ=255)

                                                                       tttATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTg                  >  2:1349040/1‑93 (MQ=255)

                                                                       tttATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCAGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTg                  >  2:3867996/1‑93 (MQ=255)

                                                                         tattaAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGaa                <  2:1119209/93‑1 (MQ=255)

                                                                         tattaAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGaa                <  1:3529135/93‑1 (MQ=255)

                                                                          attaAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAAc               <  2:5863812/93‑1 (MQ=255)

                                                                            taAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAg             >  2:6145552/1‑93 (MQ=255)

                                                                             aaGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGa            >  1:4562267/1‑93 (MQ=255)

                                                                              aGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGAt           <  2:3487160/93‑1 (MQ=255)

                                                                              aGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGAt           <  2:5947563/93‑1 (MQ=255)

                                                                               aaTTGAGACGGCCCATAACGCACGCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATc          <  1:4050728/92‑1 (MQ=255)

                                                                                aTTGAGACTGCCCATAGCGCACTCCTGAATCAACCGTAACGACCGGTCATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCg         <  2:6015085/93‑1 (MQ=255)

                                                                                 ttGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGt        >  1:2031662/1‑93 (MQ=255)

                                                                                  tGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTTATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTc       <  2:1447324/93‑1 (MQ=255)

                                                                                    agaCTGACCATAACGCAATCCTGAATCAACCGTAACGAACGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCg     >  1:3396046/1‑93 (MQ=255)

                                                                                     gaCTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCgt    <  1:325497/93‑1 (MQ=255)

                                                                                     gaCTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCgt    >  2:5002807/1‑93 (MQ=255)

ACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGT  >  NC_000913/3905112‑3905291

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