BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,421,120	Δ1 bp	noncoding (60/121 nt)	rrfF ←	5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,421,119	0	G	.	94.9%	314.6 / 17.6	137	noncoding (61/121 nt)	rrfF	5S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (4/3); new base . (51/79); total (55/82)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

Only 100 of 155 total aligned reads displayed.

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207

cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                          >  1:5638065/1‑93 (MQ=255)

 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         <  1:1835278/93‑2 (MQ=39)

 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  1:6537269/1‑92 (MQ=255)

 aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                         >  1:545218/1‑92 (MQ=39)

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:3060183/1‑91 (MQ=255)

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:2603891/1‑91 (MQ=255)

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:1903451/1‑91 (MQ=255)

  gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                        >  2:6144596/1‑91 (MQ=255)

   tCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                                                       >  1:1594612/1‑93 (MQ=255)

    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      <  2:7071437/93‑2 (MQ=25)

    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      >  1:4271369/1‑92 (MQ=255)

    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      >  1:5284591/1‑92 (MQ=255)

    cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      >  1:6072536/1‑92 (MQ=255)

    cTTCCGACTGAGCCTTTTATTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                      <  1:1229960/93‑2 (MQ=255)

     ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCa                                                                                    >  1:3564915/1‑93 (MQ=255)

        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACtt                                                                                 <  2:3785996/93‑1 (MQ=9)

        cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACACCACACTACCATCGGCGCTACGGC‑GTTTCACtt                                                                                 >  1:1148956/1‑93 (MQ=255)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                <  2:3524350/93‑1 (MQ=9)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                <  2:1276953/93‑1 (MQ=9)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                <  1:4299138/93‑1 (MQ=9)

         gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                >  1:3208226/1‑93 (MQ=16)

          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCt                                                                               >  1:41072/1‑93 (MQ=16)

          aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCt                                                                               <  1:3816825/93‑1 (MQ=9)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              >  2:4443340/1‑93 (MQ=16)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              >  1:720346/1‑93 (MQ=16)

           cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                              <  1:6585166/93‑1 (MQ=9)

            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGa                                                                             >  1:3700891/1‑93 (MQ=12)

            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGa                                                                             >  2:2031192/1‑93 (MQ=16)

            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGa                                                                             >  2:6290061/1‑93 (MQ=16)

            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGa                                                                             >  1:1118079/1‑93 (MQ=16)

            tGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGa                                                                             <  1:4650330/93‑1 (MQ=9)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGTTCCTACGGC‑GTTTCACTTCTGat                                                                            >  1:6359698/1‑92 (MQ=11)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            <  2:1689289/93‑1 (MQ=9)

             gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                            <  2:6838004/93‑1 (MQ=9)

              aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGt                                                                           <  1:4769971/93‑1 (MQ=9)

                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTGCTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         <  1:7601460/93‑1 (MQ=9)

                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         <  2:7054386/93‑1 (MQ=9)

                ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                         >  1:7290025/1‑93 (MQ=16)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        >  1:4959554/1‑93 (MQ=9)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        <  1:2678187/93‑1 (MQ=9)

                 cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                        <  1:4707285/93‑1 (MQ=9)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:511520/1‑93 (MQ=12)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       <  2:2708218/93‑1 (MQ=9)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:2234794/1‑93 (MQ=9)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  2:3758677/1‑93 (MQ=12)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  2:6906401/1‑93 (MQ=12)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  2:7651633/1‑93 (MQ=12)

                  ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                       >  1:6066963/1‑93 (MQ=12)

                   tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                      >  2:649082/1‑93 (MQ=12)

                   tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                      >  1:2374832/1‑93 (MQ=12)

                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                     <  2:7546743/93‑1 (MQ=0)

                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                     >  1:4890105/1‑93 (MQ=9)

                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                     <  2:1472452/93‑1 (MQ=0)

                    ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                     <  2:5490071/93‑1 (MQ=0)

                                                                   cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAACTTCttt                      >  1:6745650/1‑93 (MQ=0)

                                                                   cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  1:877168/93‑1 (MQ=0)

                                                                   cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  2:4489916/93‑1 (MQ=0)

                                                                   cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      <  2:2359887/93‑1 (MQ=0)

                                                                   cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                      >  1:6745648/1‑93 (MQ=0)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:3255797/93‑1 (MQ=12)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:5136916/93‑1 (MQ=12)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:3229797/93‑1 (MQ=16)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:6757431/93‑1 (MQ=16)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:1861520/93‑1 (MQ=16)

                                                                     cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  2:1114237/93‑1 (MQ=16)

                                                                     cacaCTACCACCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                    <  1:3925836/93‑1 (MQ=255)

                                                                      atacTACCATCGGAGCTAAGG‑AGTTTCATTTCTGAGTGCGGCATGGGGTCAGGTTGGACCACTGCGCTGCTGTCGCGAGACGAATTCtttttt                   >  2:4189099/3‑91 (MQ=0)

                                                                      acacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                   >  1:868308/1‑93 (MQ=16)

                                                                       cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  <  1:7523001/93‑1 (MQ=19)

                                                                       cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                  <  1:1181776/93‑1 (MQ=255)

                                                                        acTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  1:689021/93‑1 (MQ=255)

                                                                        acTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 <  1:1699196/93‑1 (MQ=255)

                                                                        acTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                 >  1:4471679/1‑93 (MQ=255)

                                                                         cTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                <  1:6911186/93‑1 (MQ=255)

                                                                         cTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                <  1:3639614/93‑1 (MQ=255)

                                                                         cTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                <  1:4620426/93‑1 (MQ=255)

                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTATTc               >  2:1919392/1‑93 (MQ=255)

                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               <  1:7375419/93‑1 (MQ=255)

                                                                          tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc               <  2:351625/93‑1 (MQ=255)

                                                                           aCCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCt              <  1:6879265/93‑1 (MQ=255)

                                                                           aCCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCt              <  2:1008141/93‑1 (MQ=255)

                                                                            ccATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTg             <  2:7465520/93‑1 (MQ=255)

                                                                            ccATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCAAATCTg             <  1:2621945/93‑1 (MQ=255)

                                                                              aTCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCCGGTGGGACCCCCCCCCTACTGCCGCCAGACAACTTCTCTTCTAATCTccc           >  1:926885/1‑93 (MQ=11)

                                                                              aTCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGGCAAATTCTTTTCTAATCTGcc           >  2:4345897/1‑93 (MQ=255)

                                                                              aTCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGcc           >  1:384053/1‑93 (MQ=255)

                                                                               tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg          >  1:3693978/1‑93 (MQ=255)

                                                                                cGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGa         <  2:5232487/93‑1 (MQ=255)

                                                                                cGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGa         <  2:4374655/93‑1 (MQ=255)

                                                                                 ggCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGaa        <  1:5835301/93‑1 (MQ=255)

                                                                                 ggCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGaa        <  1:4113232/93‑1 (MQ=255)

                                                                                    gcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt     <  2:382283/93‑1 (MQ=255)

                                                                                     cTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACttt    <  1:549539/93‑1 (MQ=255)

                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   >  2:3731204/1‑93 (MQ=255)

                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   <  2:7736241/93‑1 (MQ=255)

                                                                                      tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa   >  2:3219084/1‑93 (MQ=255)

                                                                                       aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  <  2:5507590/93‑1 (MQ=255)

                                                                                       aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  <  2:1425096/93‑1 (MQ=255)

                                                                                       aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  <  2:2318620/93‑1 (MQ=255)

                                                                                       aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa  <  2:994388/93‑1 (MQ=255)

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAA  >  NC_000913/3421027‑3421207

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