BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	641,051	T→C	M14V (ATG→GTG)	ybdQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	641,051	0	T	C	97.3%	686.5 / 4.7	219	M14V (ATG→GTG)	ybdQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/3); new base C (91/122); total (94/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

Only 100 of 225 total aligned reads displayed.

CTGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTATAAT  >  NC_000913/640961‑641140

ctgACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCActt                                                                                         <  1:974445/93‑1 (MQ=255)

    cccGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCaaa                                                                                     >  1:1593952/1‑93 (MQ=255)

    cccGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCaaa                                                                                     >  2:323913/1‑93 (MQ=255)

     ccGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAAt                                                                                    >  1:3693670/1‑93 (MQ=255)

       gggTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATAc                                                                                  <  2:4528409/93‑1 (MQ=255)

         gtagtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACAt                                                                                >  1:1024946/1‑93 (MQ=255)

         gtagtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACAt                                                                                >  1:1401118/1‑93 (MQ=255)

            gtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAGCAGCTTTGTCGCTCAATTCCACTTCAAATACATCaa                                                                             >  2:940898/1‑93 (MQ=255)

            gtaCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCaa                                                                             >  2:809740/1‑93 (MQ=255)

             taCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAAc                                                                            >  1:4246392/1‑93 (MQ=255)

             taCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAAc                                                                            >  1:605441/1‑93 (MQ=255)

              aCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACt                                                                           <  1:1157200/93‑1 (MQ=255)

              aCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACt                                                                           <  1:410404/93‑1 (MQ=255)

              aCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACt                                                                           >  2:503651/1‑93 (MQ=255)

              aCGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGAGAACAGATCTGTCGCTCAATTCCACTTCAAATACATCAAcg                                                                           >  1:576637/1‑92 (MQ=255)

               cGTGAAGTAGATGAATAACTCCGTCATCCTGGGCTAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTg                                                                          <  1:2904110/93‑1 (MQ=255)

                gTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTgg                                                                         >  1:3341472/1‑93 (MQ=255)

                 tGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGc                                                                        >  2:768054/1‑93 (MQ=255)

                  gAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCa                                                                       <  2:3090342/93‑1 (MQ=255)

                   aaGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCAt                                                                      <  2:3819720/93‑1 (MQ=255)

                    aGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATa                                                                     >  2:710119/1‑93 (MQ=255)

                     gTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATaa                                                                    <  1:4713793/93‑1 (MQ=255)

                     gTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATaa                                                                    <  2:4616838/93‑1 (MQ=255)

                       aGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATg                                                                  >  2:1321630/1‑93 (MQ=255)

                        gATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGa                                                                 <  1:1980773/93‑1 (MQ=255)

                         aTGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAt                                                                <  2:1168376/93‑1 (MQ=255)

                         aTGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAt                                                                <  2:2158772/93‑1 (MQ=255)

                          tGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGAtt                                                               >  1:2025207/1‑93 (MQ=255)

                            aaTAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACATCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGt                                                             <  2:1848263/93‑1 (MQ=255)

                            aaTAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGt                                                             <  2:3160378/93‑1 (MQ=255)

                             aTAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTc                                                            <  1:3042215/93‑1 (MQ=255)

                             aTAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTc                                                            >  1:2667730/1‑93 (MQ=255)

                              tAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCt                                                           <  2:2424716/93‑1 (MQ=255)

                              tAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCt                                                           <  1:3772222/93‑1 (MQ=255)

                              tAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGATCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCt                                                           <  1:2191271/93‑1 (MQ=255)

                               aaCTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCtt                                                          <  1:4334924/93‑1 (MQ=255)

                                  tctgTCATCGGGGGCTAGGAATTGCGCGTGGGGAACAGCTGTGTCTCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTtata                                                       <  2:2600017/90‑1 (MQ=255)

                                  tCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTtata                                                       <  2:4313956/93‑1 (MQ=255)

                                   ccGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATAc                                                      <  1:3619606/93‑1 (MQ=255)

                                   ccGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATAc                                                      <  1:1240929/93‑1 (MQ=255)

                                    cGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACa                                                     >  2:3563904/1‑93 (MQ=255)

                                    cGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACa                                                     <  2:1724677/93‑1 (MQ=255)

                                     gTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACAt                                                    >  1:1266843/1‑93 (MQ=255)

                                         tCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAAcc                                                >  2:2503065/1‑93 (MQ=255)

                                         tCCTGGGAGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAAGACATCAACTGGCAGAATGATTGTCTTATACATAAcc                                                <  1:1786169/93‑1 (MQ=255)

                                          ccTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAAccc                                               >  2:3462355/1‑93 (MQ=255)

                                           cTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCt                                              <  1:3992288/93‑1 (MQ=255)

                                           cTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCt                                              <  1:160567/93‑1 (MQ=255)

                                              ggCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTtc                                           >  1:2308890/1‑93 (MQ=255)

                                              ggCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCGCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTtc                                           >  1:852625/1‑93 (MQ=255)

                                               gCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTtct                                          >  2:2154557/1‑93 (MQ=255)

                                                  aGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTccc                                       <  2:261842/93‑1 (MQ=255)

                                                  aGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTccc                                       <  1:1373161/93‑1 (MQ=255)

                                                   ggAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCt                                      <  1:2629266/93‑1 (MQ=255)

                                                     aaTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGt                                    <  1:3293297/93‑1 (MQ=255)

                                                      aTTCAGCGTGGCGCACAGCTTTGTCGCTCCATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGtt                                   <  1:1462962/93‑1 (MQ=255)

                                                      aTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGtt                                   <  2:3784189/93‑1 (MQ=255)

                                                      aTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGtt                                   <  1:4374473/93‑1 (MQ=255)

                                                      aTTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGtt                                   <  1:2961916/93‑1 (MQ=255)

                                                       ttCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTa                                  <  1:4428573/93‑1 (MQ=255)

                                                           gCGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTCTCTCCCTGTTAATct                              >  2:91625/1‑92 (MQ=255)

                                                            cGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCAt                             <  1:1748034/93‑1 (MQ=255)

                                                            cGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCAt                             <  1:249065/93‑1 (MQ=255)

                                                            cGTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCAt                             <  2:3534185/93‑1 (MQ=255)

                                                             gTGGCGATCACCTTTGTCGCTCAAGTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTACTCCCTGTTAATCATg                            <  2:2693137/93‑1 (MQ=255)

                                                             gTGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATg                            >  2:4594325/1‑93 (MQ=255)

                                                              tGGCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGa                           >  1:2790493/1‑93 (MQ=255)

                                                               ggCGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACCTAACCCTTTCTCCCTGTTCATCATGaa                          >  2:336385/1‑93 (MQ=255)

                                                                 cGAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACa                        >  2:4203393/1‑93 (MQ=255)

                                                                  gAACAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACaa                       >  2:4000708/1‑93 (MQ=255)

                                                                   aaCAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACaaa                      <  2:4328252/93‑1 (MQ=255)

                                                                   aaCAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACaaa                      >  2:565175/1‑93 (MQ=255)

                                                                   aaCAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACaaa                      >  2:2025067/1‑93 (MQ=255)

                                                                    aCAGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAAt                     >  1:1859506/1‑93 (MQ=255)

                                                                      aGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCa                   >  2:1608643/1‑93 (MQ=255)

                                                                      aGCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCa                   <  2:3723235/93‑1 (MQ=255)

                                                                       gCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCAt                  >  2:1645957/1‑93 (MQ=255)

                                                                       gCTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCAt                  <  2:3339184/93‑1 (MQ=255)

                                                                        cTTTGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCAtt                 <  1:3084514/93‑1 (MQ=255)

                                                                        cTTTGTCACTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCAtt                 <  1:3814406/93‑1 (MQ=255)

                                                                         tttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTc                <  2:1831504/93‑1 (MQ=255)

                                                                         tttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTc                <  1:1323082/93‑1 (MQ=255)

                                                                         tttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTc                <  2:3513086/93‑1 (MQ=255)

                                                                         tttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTc                <  1:2641837/93‑1 (MQ=255)

                                                                         tttGTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTc                >  1:3784180/1‑93 (MQ=255)

                                                                            gTCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGcc             >  2:597421/1‑93 (MQ=255)

                                                                             tCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCa            <  2:303474/93‑1 (MQ=255)

                                                                             tCGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCa            <  2:2850692/93‑1 (MQ=255)

                                                                              cGCTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCAt           <  2:851382/93‑1 (MQ=255)

                                                                               gCTCAATTCCACTTCAAATACATCAACTGGCATAATAATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATg          >  1:1075877/1‑93 (MQ=255)

                                                                                cTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGa         >  2:1912944/1‑93 (MQ=255)

                                                                                cTCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGa         <  1:1815407/93‑1 (MQ=255)

                                                                                 tCAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGAt        >  1:3149720/1‑93 (MQ=255)

                                                                                  cAATTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGAtt       >  1:789644/1‑93 (MQ=255)

                                                                                   aaTTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATta      <  2:2894600/93‑1 (MQ=255)

                                                                                    aTTCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATtat     <  2:1144802/93‑1 (MQ=255)

                                                                                     ttCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATtata    <  2:3929590/93‑1 (MQ=255)

                                                                                     ttCCACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATtata    <  1:4519857/93‑1 (MQ=255)

                                                                                       ccACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTataat  <  1:3879460/93‑1 (MQ=255)

                                                                                       ccACTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTataat  <  1:1485182/93‑1 (MQ=255)

CTGACCCGGGTAGTACGTGAAGTAGATGAATAACTCCGTCATCCTGGGCGAGGAATTCAGCGTGGCGAACAGCTTTGTCGCTCAATTCCATTTCAAATACATCAACTGGCATAATGATTGTCTTATACATAACCCTTTCTCCCTGTTAATCATGAACAAATCATTCGCCATGATTATAAT  >  NC_000913/640961‑641140

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