BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,061,961	Δ1 bp	coding (1386/1458 nt)	b3876 ←	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,061,959	0	C	.	100.0%	1192.4 / NA	260	coding (1388/1458 nt)	b3876	putative permease
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (146/114); total (146/114)

Only 100 of 276 total aligned reads displayed.

TTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTA  >  NC_000913/4061867‑4062047

tttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGc                                                                                          >  1:2073003/1‑93 (MQ=255)

tttttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGc                                                                                          >  1:1981634/1‑93 (MQ=255)

  tttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccg                                                                                        >  2:224145/1‑92 (MQ=255)

  tttCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGccg                                                                                        <  2:2278293/93‑2 (MQ=255)

      ataataAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGtttt                                                                                   >  1:1679956/1‑93 (MQ=255)

       taataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGttttt                                                                                  <  1:3296921/93‑1 (MQ=255)

        aataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCCATTGACATTTATAACGTTACGGAAG‑CCGtttttt                                                                                 >  1:2000090/1‑93 (MQ=255)

         ataaGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttc                                                                                >  1:615668/1‑93 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttt                                                                             >  2:437653/1‑92 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttc                                                                             <  1:3173492/93‑1 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttc                                                                             <  1:3507538/93‑1 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttc                                                                             >  2:885377/1‑93 (MQ=255)

            aGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTttcttc                                                                             <  1:2460425/93‑1 (MQ=255)

             gTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                            <  2:1600046/93‑1 (MQ=255)

             gTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                            >  1:2316064/1‑93 (MQ=255)

             gTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                            >  1:3137118/1‑93 (MQ=255)

             gTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTtcttct                                                                            >  1:2647998/1‑93 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                           >  2:768814/1‑93 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTg                                                                           >  2:4260255/1‑93 (MQ=255)

              tCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGGTTTTtcttctt                                                                           >  1:2480903/1‑92 (MQ=255)

               ccATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGt                                                                          <  2:1845983/93‑1 (MQ=255)

                  ttttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGc                                                                       >  2:3269865/1‑93 (MQ=255)

                    ttttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCa                                                                     >  1:426980/1‑93 (MQ=255)

                     tttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAt                                                                    <  1:298983/93‑1 (MQ=255)

                      ttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtc                                                                   >  1:2399013/1‑93 (MQ=255)

                      ttAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGATTTTTCTTCTGTCGCCAtc                                                                   <  2:4559486/93‑1 (MQ=255)

                       tAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtct                                                                  <  2:479404/93‑1 (MQ=255)

                        aaaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtctc                                                                 <  1:922095/93‑1 (MQ=255)

                        aaaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtctc                                                                 <  1:1854245/93‑1 (MQ=255)

                        aaaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtctc                                                                 <  1:1286057/93‑1 (MQ=255)

                        aaaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCAtctc                                                                 <  1:3695758/93‑1 (MQ=255)

                         aaCGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATctct                                                                >  2:3226413/1‑93 (MQ=255)

                           cGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCttt                                                              <  1:2250650/93‑1 (MQ=255)

                            gTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTa                                                             <  2:2405150/93‑1 (MQ=255)

                             tAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCCTCTCTTTAt                                                            >  2:2014605/1‑93 (MQ=255)

                              aGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtt                                                           >  2:734733/1‑93 (MQ=255)

                              aGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAtt                                                           <  1:2574120/93‑1 (MQ=255)

                               gggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttt                                                          <  2:2136149/93‑1 (MQ=255)

                               gggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttt                                                          >  1:549625/1‑93 (MQ=255)

                               gggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttt                                                          >  2:4027059/1‑93 (MQ=255)

                               gggAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttt                                                          >  1:4552999/1‑93 (MQ=255)

                                 gAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTAttttt                                                        <  2:1546656/93‑1 (MQ=255)

                                   ggCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGc                                                      >  1:2879923/1‑93 (MQ=255)

                                      cgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgtt                                                   >  2:4529692/1‑93 (MQ=255)

                                      cgccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgtt                                                   <  2:2481135/93‑1 (MQ=255)

                                       gccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttg                                                  >  2:1626436/1‑93 (MQ=255)

                                       gccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttg                                                  >  1:3840263/1‑93 (MQ=255)

                                       gccACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttg                                                  >  2:1074206/1‑93 (MQ=255)

                                         cacccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggt                                                <  2:2683124/93‑1 (MQ=255)

                                         cacccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggt                                                <  1:314411/93‑1 (MQ=255)

                                         cacccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggt                                                <  2:2847317/93‑1 (MQ=255)

                                          acccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggtt                                               <  1:3870805/93‑1 (MQ=255)

                                          acccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggtt                                               >  1:2157828/1‑93 (MQ=255)

                                           cccaccAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCgttggttg                                              <  2:3373435/93‑1 (MQ=255)

                                                 aGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTcc                                        >  2:3499802/1‑93 (MQ=255)

                                                  ggATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCt                                       >  2:4456690/1‑93 (MQ=255)

                                                  ggATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCt                                       >  2:633431/1‑93 (MQ=255)

                                                   gATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTc                                      <  2:3576169/93‑1 (MQ=255)

                                                   gATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTc                                      >  1:4601222/1‑93 (MQ=255)

                                                     tGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaa                                    <  2:926238/93‑1 (MQ=255)

                                                       gAAATAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaata                                  >  1:4121741/1‑93 (MQ=255)

                                                        aaaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataa                                 >  1:3950360/1‑93 (MQ=255)

                                                         aaTAGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCaataat                                <  2:3688403/93‑1 (MQ=255)

                                                            aGCGGGTATTTACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCaa                             <  2:1956251/93‑1 (MQ=255)

                                                            aGCGGCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCaa                             <  1:1906689/93‑1 (MQ=255)

                                                               ggCTATTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGc                          >  1:4440331/1‑93 (MQ=255)

                                                                   aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAc                      <  2:2834663/93‑1 (MQ=255)

                                                                   aTTGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAc                      >  2:1200367/1‑93 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     <  1:1617960/93‑1 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     >  2:1738873/1‑93 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     <  1:450057/93‑1 (MQ=255)

                                                                    ttGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAAcc                     >  2:4108919/1‑93 (MQ=255)

                                                                     tGACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCg                    >  1:2848859/1‑93 (MQ=255)

                                                                      gACATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGt                   >  1:37271/1‑93 (MQ=255)

                                                                       aCATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtt                  >  2:3614311/1‑93 (MQ=255)

                                                                        cATTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttt                 >  1:1217929/1‑93 (MQ=255)

                                                                         aTTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                >  2:2379160/1‑93 (MQ=255)

                                                                         aTTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                <  1:1536299/93‑1 (MQ=255)

                                                                         aTTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                >  2:1381913/1‑93 (MQ=255)

                                                                         aTTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                <  2:1002959/93‑1 (MQ=255)

                                                                         aTTTATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                >  2:3927246/1‑93 (MQ=255)

                                                                         aTTTATAACGTGACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGtttt                <  2:2953047/93‑1 (MQ=255)

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttttt               >  2:2115300/1‑93 (MQ=255)

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttttt               >  1:4526665/1‑93 (MQ=255)

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttttt               <  2:1902080/93‑1 (MQ=255)

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGttttt               >  2:1474413/1‑93 (MQ=255)

                                                                          tttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGgtttt               >  1:4235741/1‑93 (MQ=255)

                                                                           ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc              <  2:3046954/93‑1 (MQ=255)

                                                                           ttATAACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTc              <  1:364424/93‑1 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           <  2:811136/93‑1 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           <  2:1874104/93‑1 (MQ=255)

                                                                              taACGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCAtt           <  1:2496593/93‑1 (MQ=255)

                                                                               aaCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTg          <  1:1011743/93‑1 (MQ=255)

                                                                                aCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGa         <  1:1094770/93‑1 (MQ=255)

                                                                                aCGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGa         <  2:4180379/93‑1 (MQ=255)

                                                                                 cGTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAg        <  2:4161269/93‑1 (MQ=255)

                                                                                  gTTACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgt       <  2:3764769/93‑1 (MQ=255)

                                                                                   ttACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGATTTTCATTGAgtg      >  1:684131/1‑93 (MQ=255)

                                                                                    tACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgt     >  1:2750638/1‑93 (MQ=255)

                                                                                    tACGGAAG‑CCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAgtgt     <  2:4288416/93‑1 (MQ=255)

TTTTTCATAATAAGTCCATTTTTTAAACGTAGGGAGGCCGCCACCCACCAGGATGGAAATAGCGGCTATTGACATTTATAACGTTACGGAAGCCCGTTTTTTCTTCTGTCGCCATCTCTTTATTTTTGCGTTGGTTGATTTCCTCAATAATCAATGCAAACCGTTTTTCATTGAGTGTGTA  >  NC_000913/4061867‑4062047

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: `‑`