BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,289,578	G→A	L285L (CTG→TTG)	ccmH ←	CcmH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,289,578	0	G	A	99.5%	647.9 / ‑5.4	187	L285L (CTG→TTG)	ccmH	CcmH
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (98/88); total (98/89)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 190 total aligned reads displayed.

CGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGT  >  NC_000913/2289488‑2289667

 ggCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTg                                                                                        <  1:4063712/93‑1 (MQ=255)

  gCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATTGTCCGTTCTCACCAACTGa                                                                                       >  2:2380030/1‑93 (MQ=255)

  gCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGa                                                                                       >  1:3107122/1‑93 (MQ=255)

     aCGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGt                                                                                    >  2:3082648/1‑93 (MQ=255)

      cGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTa                                                                                   <  1:1063985/93‑1 (MQ=255)

        gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGc                                                                                 >  1:3843650/1‑93 (MQ=255)

        gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGc                                                                                 >  2:140176/1‑93 (MQ=255)

         cTTCGCCAAATCGCTGCTGCTCAAATGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCa                                                                                <  2:331278/93‑1 (MQ=255)

         cTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCa                                                                                >  2:2996368/1‑93 (MQ=255)

           tCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGt                                                                              <  2:1858832/93‑1 (MQ=255)

             gCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTc                                                                            <  1:331001/93‑1 (MQ=255)

             gCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTc                                                                            >  1:2045139/1‑93 (MQ=255)

             gCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTc                                                                            <  2:3048762/93‑1 (MQ=255)

                aaaTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCAcc                                                                         >  1:1128506/1‑93 (MQ=255)

                aaaTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCAcc                                                                         <  1:3772809/93‑1 (MQ=255)

                aaaTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCAcc                                                                         <  1:1695251/93‑1 (MQ=255)

                aaaTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCAcc                                                                         >  2:3287303/1‑93 (MQ=255)

                    cGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCAccgccg                                                                     >  1:2243855/1‑93 (MQ=255)

                     gctgctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGa                                                                    <  1:3811300/93‑1 (MQ=255)

                     gctgctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGa                                                                    >  2:4193201/1‑93 (MQ=255)

                       tgctgctCAAAGGCATTAAACGCATACATGCTTAGCACATGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAgg                                                                  <  1:4069403/93‑1 (MQ=255)

                       tgctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGCgg                                                                  >  1:232188/1‑93 (MQ=255)

                       tgctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAgg                                                                  >  1:1189158/1‑93 (MQ=255)

                       tgctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAgg                                                                  >  2:1671350/1‑93 (MQ=255)

                       cgctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAgg                                                                  <  1:3157907/92‑1 (MQ=255)

                           gctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGt                                                              <  2:2737638/93‑1 (MQ=255)

                           gctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGt                                                              >  1:2742021/1‑93 (MQ=255)

                            ctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGtt                                                             <  2:3081215/93‑1 (MQ=255)

                              cAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGt                                                           >  1:304971/1‑93 (MQ=255)

                              cAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGt                                                           >  1:3892548/1‑93 (MQ=255)

                               aaaGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTc                                                          >  2:1560948/1‑93 (MQ=255)

                               aaaGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTc                                                          <  1:63642/93‑1 (MQ=255)

                                 aGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGt                                                        >  2:2796942/1‑93 (MQ=255)

                                 aGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGt                                                        <  2:3167099/93‑1 (MQ=255)

                                  ggCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGGCGtt                                                       >  2:3142662/1‑93 (MQ=255)

                                    cATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTgg                                                     <  2:811517/93‑1 (MQ=255)

                                    cATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTgg                                                     >  2:1183889/1‑93 (MQ=255)

                                    cATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTgg                                                     >  2:4740300/1‑93 (MQ=255)

                                     aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTggg                                                    >  2:4533771/1‑93 (MQ=255)

                                      ttAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGa                                                   >  2:3501528/1‑93 (MQ=255)

                                       tAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGAt                                                  <  1:4607102/93‑1 (MQ=255)

                                        aaaCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATc                                                 <  1:4208789/93‑1 (MQ=255)

                                        aaaCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATc                                                 <  1:538551/93‑1 (MQ=255)

                                         aaCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCa                                                <  2:4396608/93‑1 (MQ=255)

                                         aaCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCa                                                >  2:2087513/1‑93 (MQ=255)

                                          aCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAg                                               <  1:2687991/93‑1 (MQ=255)

                                           cGCATACATGCTTCGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGa                                              <  1:881296/93‑1 (MQ=255)

                                           cGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGa                                              >  1:254910/1‑93 (MQ=255)

                                           cGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGa                                              <  1:4688997/93‑1 (MQ=255)

                                            gCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGAt                                             >  1:2492600/1‑93 (MQ=255)

                                              aTACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGa                                           >  2:794759/1‑93 (MQ=255)

                                               tACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGaa                                          >  1:503249/1‑93 (MQ=255)

                                                 cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACg                                        >  2:2925706/1‑93 (MQ=255)

                                                 cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACg                                        >  2:1692107/1‑93 (MQ=255)

                                                 cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACg                                        <  2:1953152/93‑1 (MQ=255)

                                                   tGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAg                                      >  1:3340795/1‑93 (MQ=255)

                                                    gCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAga                                     >  2:1216573/1‑92 (MQ=255)

                                                    gCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGt                                     >  1:2341761/1‑93 (MQ=255)

                                                       tAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCaa                                  <  2:3299804/93‑1 (MQ=255)

                                                        aGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAAc                                 >  1:2976728/1‑93 (MQ=255)

                                                        aGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAAc                                 >  1:4296877/1‑93 (MQ=255)

                                                        aGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAAc                                 <  1:3251222/93‑1 (MQ=255)

                                                          cacaCGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGc                               <  2:4180632/93‑1 (MQ=255)

                                                          cacaCGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGc                               <  1:4334046/93‑1 (MQ=255)

                                                           acacGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCt                              >  2:2702402/1‑93 (MQ=255)

                                                            cacGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCtt                             >  1:4429711/1‑93 (MQ=255)

                                                            cacGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGATTGTCGTTGGGATCAGATGAACGAGTCAACGCtt                             >  2:3938384/1‑93 (MQ=255)

                                                             acGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTc                            <  1:1482507/93‑1 (MQ=255)

                                                              cGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCa                           <  1:226871/93‑1 (MQ=255)

                                                              cGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCa                           >  1:3486314/1‑93 (MQ=255)

                                                                gATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGc                         <  2:2869345/93‑1 (MQ=255)

                                                                  tatTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGt                       >  1:1755855/1‑93 (MQ=255)

                                                                    ttGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTAt                     >  1:3442601/1‑93 (MQ=255)

                                                                    ttGCTATGGTCCGTTCTCACCAACTGACGTAGCAGCTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTAt                     >  2:4676173/1‑93 (MQ=255)

                                                                      gCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATcc                   >  2:3613452/1‑93 (MQ=255)

                                                                       cTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATccc                  <  1:1275700/93‑1 (MQ=255)

                                                                       cTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATccc                  >  2:2529474/1‑93 (MQ=255)

                                                                       cTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATccc                  >  1:565311/1‑93 (MQ=255)

                                                                         ttGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAg                <  1:1622533/92‑1 (MQ=255)

                                                                         aTGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAg                <  1:2606755/93‑1 (MQ=255)

                                                                         aTGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAg                >  2:45561/1‑93 (MQ=255)

                                                                         aTGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAg                <  1:182205/93‑1 (MQ=255)

                                                                          ttgTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGt               <  1:3418963/91‑1 (MQ=255)

                                                                          tGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGt               >  1:1251056/1‑93 (MQ=255)

                                                                          tGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGt               >  2:2959979/1‑93 (MQ=255)

                                                                          tGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGt               >  1:4439744/1‑93 (MQ=255)

                                                                          tGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGt               <  1:4369523/93‑1 (MQ=255)

                                                                            gTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgc             >  2:629261/1‑93 (MQ=255)

                                                                              ccGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgcag           >  2:1285977/1‑93 (MQ=255)

                                                                              ccGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgcag           <  1:4652903/93‑1 (MQ=255)

                                                                               cGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgcagc          <  1:596748/93‑1 (MQ=255)

                                                                                gTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgcagca         >  2:2325034/1‑93 (MQ=255)

                                                                                  tctcACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATc       >  1:3423428/1‑93 (MQ=255)

                                                                                   ctcACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCa      <  2:4063061/93‑1 (MQ=255)

                                                                                    tcACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCAc     >  1:2862131/1‑93 (MQ=255)

                                                                                     cACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACt    <  1:2414832/93‑1 (MQ=255)

                                                                                      aCCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTg   >  2:3742395/1‑93 (MQ=255)

                                                                                      aCCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTg   <  1:4150123/93‑1 (MQ=255)

                                                                                       ccAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGt  <  1:165450/93‑1 (MQ=255)

                                                                                       ccAACTGACGTAGCAGTTCACCGACGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGt  <  2:2593881/93‑1 (MQ=255)

CGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGT  >  NC_000913/2289488‑2289667

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