BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,691,941	T→C	Y116C (TAC→TGC)	uidB ←	glucuronide permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,691,941	0	T	C	90.4%	626.1 / 48.5	229	Y116C (TAC→TGC)	uidB	glucuronide permease
Reads supporting (aligned to +/- strand): ref base T (14/8); new base C (107/100); total (121/108)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 230 total aligned reads displayed.

GACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAG  >  NC_000913/1691849‑1692031

  cGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAg                                                                                          <  2:5690392/93‑1 (MQ=255)

  cGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAg                                                                                          >  2:7161161/1‑93 (MQ=255)

   gggCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAgg                                                                                         >  1:2307473/1‑93 (MQ=255)

    ggCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGt                                                                                        >  2:2847517/1‑93 (MQ=255)

    ggCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGt                                                                                        >  2:5481528/1‑93 (MQ=255)

     gcgcGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTc                                                                                       <  1:4715091/93‑1 (MQ=255)

      cgcgGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCa                                                                                      <  1:973573/93‑1 (MQ=255)

        cgGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCaaa                                                                                    <  2:7066402/93‑1 (MQ=255)

          ggATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAata                                                                                  >  1:366917/1‑93 (MQ=255)

           gATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCGCCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAatat                                                                                 <  2:6552188/93‑1 (MQ=255)

           gATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAatat                                                                                 <  2:956854/93‑1 (MQ=255)

           gATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAatat                                                                                 <  1:1552679/93‑1 (MQ=255)

            aTTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATg                                                                                <  1:5884120/93‑1 (MQ=255)

             ttGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGc                                                                               >  2:195033/1‑93 (MQ=255)

             ttGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGc                                                                               >  1:1379658/1‑93 (MQ=255)

               gtgGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCAt                                                                             <  1:2126913/93‑1 (MQ=255)

               gtgGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCAt                                                                             <  2:6184537/93‑1 (MQ=255)

                tgGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATa                                                                            >  1:3186037/1‑93 (MQ=255)

                tgGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATa                                                                            <  2:384550/93‑1 (MQ=255)

                  gttgttGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATaca                                                                          <  2:6443028/93‑1 (MQ=255)

                    tgttgGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATAcacc                                                                        >  1:6796406/1‑93 (MQ=255)

                     gttgGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATAcacca                                                                       >  2:6687489/1‑93 (MQ=255)

                      ttgGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATAcaccac                                                                      <  2:5400998/93‑1 (MQ=255)

                        gggTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACtt                                                                    <  1:469145/93‑1 (MQ=255)

                           tCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGc                                                                 <  1:177892/93‑1 (MQ=255)

                           tCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGc                                                                 <  1:5586566/93‑1 (MQ=255)

                           tCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGc                                                                 >  2:1642800/1‑93 (MQ=255)

                           tCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGc                                                                 <  1:5178191/93‑1 (MQ=255)

                           tCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGc                                                                 <  2:4160596/93‑1 (MQ=255)

                             aTCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTa                                                               >  1:879747/1‑93 (MQ=255)

                             aTCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTa                                                               >  1:7329638/1‑93 (MQ=255)

                             aTCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTa                                                               <  1:165408/93‑1 (MQ=255)

                              tCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTAc                                                              <  2:465214/93‑1 (MQ=255)

                              tCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTAc                                                              <  2:1401467/93‑1 (MQ=255)

                              tCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCACACACCACTTTGCTAc                                                              >  1:1057061/1‑93 (MQ=255)

                                gcgGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCa                                                            <  1:239393/93‑1 (MQ=255)

                                gcgGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCa                                                            <  1:5416599/93‑1 (MQ=255)

                                  ggTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATg                                                          >  1:1420412/1‑93 (MQ=255)

                                   gTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATgg                                                         <  1:596561/93‑1 (MQ=255)

                                   gTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATgg                                                         >  2:2021458/1‑93 (MQ=255)

                                    tAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGc                                                        >  1:3399430/1‑93 (MQ=255)

                                       cAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTcc                                                     <  1:4468293/93‑1 (MQ=255)

                                       cAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTcc                                                     <  1:4021435/93‑1 (MQ=255)

                                       cAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTcc                                                     <  2:5249022/93‑1 (MQ=255)

                                        aaGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCa                                                    >  2:3341981/1‑93 (MQ=255)

                                            gAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCg                                                >  2:4097283/1‑93 (MQ=255)

                                            gAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCg                                                <  2:7423356/93‑1 (MQ=255)

                                               ccATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAACTATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTc                                             <  1:6210399/93‑1 (MQ=255)

                                               ccATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTc                                             >  1:1183045/1‑93 (MQ=255)

                                               ccATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTc                                             <  1:3790118/93‑1 (MQ=255)

                                               ccATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTc                                             >  1:6442352/1‑93 (MQ=255)

                                                cATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCa                                            >  1:5022855/1‑93 (MQ=255)

                                                cATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCa                                            >  1:2523697/1‑93 (MQ=255)

                                                cATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCa                                            >  2:1586542/1‑93 (MQ=255)

                                                cATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCa                                            >  2:4776198/1‑93 (MQ=255)

                                                 aTAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAg                                           >  2:713757/1‑93 (MQ=255)

                                                  tAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGc                                          <  1:2950612/93‑1 (MQ=255)

                                                  tAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGc                                          >  1:4844948/1‑93 (MQ=255)

                                                   aaGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCa                                         >  2:6337064/1‑93 (MQ=255)

                                                     ggAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCAcc                                       >  1:4129353/1‑93 (MQ=255)

                                                      gAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCAccc                                      >  2:5722837/1‑93 (MQ=255)

                                                       aaTATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCa                                     >  2:4438581/1‑93 (MQ=255)

                                                         tatTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGa                                   >  1:612844/1‑93 (MQ=255)

                                                         tatTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGa                                   >  1:475157/1‑93 (MQ=255)

                                                         tatTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGa                                   >  1:4764346/1‑93 (MQ=255)

                                                           ttCACCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAAt                                 <  2:1102074/93‑1 (MQ=255)

                                                              aCCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATAcc                              <  2:304077/93‑1 (MQ=255)

                                                              aCCAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATAcc                              >  2:2548666/1‑93 (MQ=255)

                                                               ccAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCa                             <  2:5064304/93‑1 (MQ=255)

                                                               ccAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCa                             <  2:2909485/93‑1 (MQ=255)

                                                               ccAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCa                             >  1:6232197/1‑93 (MQ=255)

                                                                cAGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAg                            >  1:4432648/1‑93 (MQ=255)

                                                                 aGGCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGc                           >  1:4812990/1‑93 (MQ=255)

                                                                  ggCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCa                          >  1:5316244/1‑93 (MQ=255)

                                                                  ggCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCa                          >  2:5022586/1‑93 (MQ=255)

                                                                  ggCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCa                          <  1:2076984/93‑1 (MQ=255)

                                                                  ggCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCa                          >  2:6819084/1‑93 (MQ=255)

                                                                   gCTGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCAc                         <  1:5558649/93‑1 (MQ=255)

                                                                     tGTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGc                       <  2:5877904/93‑1 (MQ=255)

                                                                      gTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCt                      >  2:6599495/1‑93 (MQ=255)

                                                                      gTAGCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCt                      <  2:3608701/93‑1 (MQ=255)

                                                                       tAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTg                     >  2:4510163/1‑93 (MQ=255)

                                                                         gCAAAGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGaa                   >  1:7122519/1‑93 (MQ=255)

                                                                           aaaGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGa                 >  1:2023442/1‑93 (MQ=255)

                                                                           aaaGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGa                 >  1:4260553/1‑93 (MQ=255)

                                                                             aGCCCGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATc               >  1:5084607/1‑93 (MQ=255)

                                                                                 cGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTa           >  2:5736962/1‑93 (MQ=255)

                                                                                 cGAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTa           <  1:6942951/93‑1 (MQ=255)

                                                                                  gAGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTaa          >  2:5332640/1‑93 (MQ=255)

                                                                                   aGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAAc         <  2:6429413/93‑1 (MQ=255)

                                                                                   aGGCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAAc         <  2:1226915/93‑1 (MQ=255)

                                                                                    ggCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACg        >  2:825974/1‑93 (MQ=255)

                                                                                     gCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACgg       <  2:2280095/93‑1 (MQ=255)

                                                                                     gCCCATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACgg       <  2:6798672/93‑1 (MQ=255)

                                                                                      cccATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGgc      <  2:1807081/93‑1 (MQ=255)

                                                                                        cATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGgcgc    >  2:550917/1‑93 (MQ=255)

                                                                                        cATGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGgcgc    >  2:7371234/1‑93 (MQ=255)

                                                                                         aTGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCa   >  1:2472897/1‑93 (MQ=255)

                                                                                         aTGCAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCa   <  1:7460578/93‑1 (MQ=255)

                                                                                          tGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAg  >  2:1147184/1‑93 (MQ=255)

GACGGGCGCGGGATTGTGGTTGTTGGGTCATCGCGGTAGCAAGTGAACCATAAGGAATATTCACCAGGCTGTAGCAAAGCCCGAGGCCCATGTAGGTCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAG  >  NC_000913/1691849‑1692031

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