BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,630,942	T→C	K383K (AAA→AAG)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,630,942	0	T	C	99.3%	457.6 / ‑4.9	144	K383K (AAA→AAG)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (78/65); total (78/66)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 148 total aligned reads displayed.

CACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCC  >  NC_000913/2630852‑2631031

cACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGActt                                                                                         <  1:519297/93‑1 (MQ=255)

 aCCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTg                                                                                        >  2:3640908/1‑93 (MQ=255)

 aCCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTg                                                                                        >  1:4242271/1‑93 (MQ=255)

  ccAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGt                                                                                       >  1:670917/1‑93 (MQ=255)

  acAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGt                                                                                       <  1:1968576/92‑1 (MQ=255)

     gTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAg                                                                                    >  2:2739021/1‑93 (MQ=255)

      tttGTCGTCTGCGTTATCTGTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCGCGGTAAGACTTGTAAGa                                                                                   <  2:2148561/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGATCCTTAGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  1:4423150/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAgca                                                                                  >  2:1775569/1‑91 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  1:1570796/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  2:757077/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  2:1557422/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  1:3575736/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  2:4087632/93‑1 (MQ=255)

        tGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAAc                                                                                 <  2:550515/93‑1 (MQ=255)

         gTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACg                                                                                <  1:1091951/93‑1 (MQ=255)

         gTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACCGTAAGACTTGTAAGAACg                                                                                >  1:2200448/1‑93 (MQ=255)

           cGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGc                                                                              >  1:1382478/1‑93 (MQ=255)

            ggCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGcc                                                                             <  2:4131452/93‑1 (MQ=255)

             gCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGcccc                                                                            >  2:158024/1‑93 (MQ=255)

             gCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGccc                                                                            <  2:1907744/93‑1 (MQ=255)

               aGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTg                                                                          >  1:1207997/1‑93 (MQ=255)

               aGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTg                                                                          >  2:468541/1‑93 (MQ=255)

               aGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTg                                                                          <  1:1940773/93‑1 (MQ=255)

                gCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTgg                                                                         >  1:503169/1‑93 (MQ=255)

                gCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTgg                                                                         <  1:4617484/93‑1 (MQ=255)

                gCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTgg                                                                         >  2:3139755/1‑93 (MQ=255)

                 cGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGt                                                                        >  1:3659763/1‑93 (MQ=255)

                  gTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAGCCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTa                                                                       <  1:3190229/93‑1 (MQ=255)

                  gTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTa                                                                       <  2:484712/93‑1 (MQ=255)

                   ttATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTag                                                                      >  1:3994254/1‑93 (MQ=255)

                    tATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGTCCCTGGGagg                                                                     >  2:2333429/1‑92 (MQ=255)

                    tATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTaga                                                                     <  1:2252467/93‑1 (MQ=255)

                      tCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGt                                                                   >  2:4438163/1‑93 (MQ=255)

                      tCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGt                                                                   >  2:277976/1‑93 (MQ=255)

                      tCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGt                                                                   >  2:36087/1‑93 (MQ=255)

                        gCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTc                                                                 >  1:2453561/1‑93 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               >  2:3192973/1‑93 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               <  1:36008/93‑1 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               <  2:602995/93‑1 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               <  2:1599846/93‑1 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               <  1:762241/93‑1 (MQ=255)

                           ctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAt                                                              <  1:1260727/93‑1 (MQ=255)

                             ggAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAttt                                                            >  1:3794366/1‑93 (MQ=255)

                               aaaTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCa                                                          >  1:2380387/1‑93 (MQ=255)

                               aaaTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCa                                                          <  2:4602136/93‑1 (MQ=255)

                                 aTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAcc                                                        >  2:3910779/1‑93 (MQ=255)

                                  tAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAccc                                                       <  2:956552/93‑1 (MQ=255)

                                  tAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAccc                                                       <  1:3103630/93‑1 (MQ=255)

                                  tAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAccc                                                       >  1:1699958/1‑93 (MQ=255)

                                     cGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGga                                                    <  1:81956/93‑1 (MQ=255)

                                     cGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGga                                                    >  1:3732829/1‑93 (MQ=255)

                                        tCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAt                                                 >  2:1449084/1‑93 (MQ=255)

                                         cAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAtt                                                <  2:4322484/93‑1 (MQ=255)

                                           gagGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttct                                              >  2:3501953/1‑93 (MQ=255)

                                           gagGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttct                                              <  1:1039383/93‑1 (MQ=255)

                                           gagGAACCTTTGGACATCGCGCCCAGGGAAACCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttct                                              >  2:142096/1‑93 (MQ=255)

                                            agGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttctt                                             <  1:2711824/93‑1 (MQ=255)

                                            agGAACCTTTGGACATCGCCCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttctt                                             <  1:834804/93‑1 (MQ=255)

                                             ggAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttcttc                                            >  1:3040657/1‑93 (MQ=255)

                                             ggAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttcttc                                            <  2:3512932/93‑1 (MQ=255)

                                              gAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGGTTCTTCa                                           >  1:4215593/1‑93 (MQ=255)

                                                   tttGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTAcc                                      <  1:752796/93‑1 (MQ=255)

                                                    ttGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTAcca                                     <  1:797895/93‑2 (MQ=255)

                                                       gACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACGCGGAGATTCTTCAGTACACGcc                                  >  1:1518994/1‑93 (MQ=255)

                                                        aCATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCa                                 <  2:66895/93‑1 (MQ=255)

                                                         cATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTTAGAACGGCCCTGGTAGAGTTCGATTTCAACCGGAGATTCTTCAGTACCCGCCAg                                >  2:4293449/1‑93 (MQ=255)

                                                         cATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAg                                >  1:1036518/1‑93 (MQ=255)

                                                          aTCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGc                               >  1:2812926/1‑93 (MQ=255)

                                                          aTCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGc                               >  2:1873091/1‑93 (MQ=255)

                                                            cgcgCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCAt                             <  2:4237623/93‑1 (MQ=255)

                                                            cgcgCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCAt                             <  1:4036159/93‑1 (MQ=255)

                                                               gcCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGa                          <  2:1117672/93‑1 (MQ=255)

                                                               gcCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGa                          <  2:2870712/93‑1 (MQ=255)

                                                                cccAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGaa                         <  1:948193/93‑1 (MQ=255)

                                                                cccAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGaa                         >  1:3294690/1‑93 (MQ=255)

                                                                 ccAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAAc                        <  1:3228615/93‑1 (MQ=255)

                                                                  cAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAAcc                       >  2:3913222/1‑93 (MQ=255)

                                                                  cAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAAcc                       >  2:759908/1‑93 (MQ=255)

                                                                   aGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCt                      >  2:3375802/1‑93 (MQ=255)

                                                                    gggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTa                     >  1:604999/1‑93 (MQ=255)

                                                                    gggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTa                     >  1:4516636/1‑93 (MQ=255)

                                                                      gAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTAcc                   >  1:1702340/1‑93 (MQ=255)

                                                                       aaCCCATCCCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCa                  <  2:4069353/93‑1 (MQ=255)

                                                                        aCCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCAt                 >  2:1123699/1‑93 (MQ=255)

                                                                          ccATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCa               >  2:880132/1‑93 (MQ=255)

                                                                           cATCCCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCGGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAc              <  2:240806/93‑1 (MQ=255)

                                                                           cATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAc              >  2:1483678/1‑93 (MQ=255)

                                                                           cATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAc              >  1:1381888/1‑93 (MQ=255)

                                                                            aTACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAcc             >  1:2187681/1‑93 (MQ=255)

                                                                            aTACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAcc             >  1:3880183/1‑93 (MQ=255)

                                                                              aCCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACcgc           >  2:4150252/1‑93 (MQ=255)

                                                                                cACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCgcgc         >  1:3485761/1‑93 (MQ=255)

                                                                                cACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCgcgc         >  1:385446/1‑93 (MQ=255)

                                                                                 aCGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCt        >  1:3462660/1‑93 (MQ=255)

                                                                                 aCGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCt        >  2:1307210/1‑93 (MQ=255)

                                                                                    gTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGACCCTACCATCACCGCGCTTgc     <  1:4479654/93‑1 (MQ=255)

                                                                                     tAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTgcg    >  2:2277208/1‑93 (MQ=255)

                                                                                     tAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTgcg    >  1:895700/1‑93 (MQ=255)

                                                                                       aGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGcc  <  2:2953146/93‑1 (MQ=255)

CACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCC  >  NC_000913/2630852‑2631031

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`