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.4%	1136.9 / NA	315	L285L (CTG→TTG)	ccmH	CcmH
Reads supporting (aligned to +/- strand): ref base G (1/0); major base A (152/160); minor base C (1/1); total (154/161)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

CGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGT  >  NC_000913/2289488‑2289667

cgGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACt                                                                                         <  2:4123862/93‑1 (MQ=255)

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

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

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

    aaCGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACg                                                                                     <  1:259879/93‑1 (MQ=255)

    aaCGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACg                                                                                     <  1:2711373/93‑1 (MQ=255)

     aCGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCACCTGACGt                                                                                    >  1:4047764/1‑93 (MQ=255)

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

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

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

       ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAg                                                                                  <  2:2012807/93‑1 (MQ=255)

       ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATAGTCCGTTCTCACCAACTGACGTAg                                                                                  >  1:640020/1‑93 (MQ=255)

        gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGc                                                                                 <  2:3998631/93‑1 (MQ=255)

        gCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGc                                                                                 <  1:7332857/93‑1 (MQ=255)

         cTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCa                                                                                <  1:3623699/93‑1 (MQ=255)

         cTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCa                                                                                <  1:9361045/93‑1 (MQ=255)

          ttCTCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAg                                                                               <  2:2955618/93‑1 (MQ=255)

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

           tCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGt                                                                              >  1:8063703/1‑93 (MQ=255)

              ccAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCa                                                                           >  1:2341871/1‑93 (MQ=255)

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

                 aaTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCAccg                                                                        <  1:3077488/93‑1 (MQ=255)

                  aTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCAccgc                                                                       <  1:373065/93‑1 (MQ=255)

                  aTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCAccgc                                                                       >  1:5579227/1‑93 (MQ=255)

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

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

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

                        gctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGc                                                                 >  2:3481559/1‑93 (MQ=255)

                           gctCAAAGGCATTAAACGCATACATGCTTAGCACTCGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGt                                                              <  2:8403351/93‑1 (MQ=255)

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

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

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

                                aaGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCg                                                         <  2:3115687/93‑1 (MQ=255)

                                aaGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCg                                                         >  1:6295190/1‑93 (MQ=255)

                                 aGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGCTGTCGt                                                        <  2:9760691/93‑1 (MQ=255)

                                  ggttttAAACGCATACATGCGTTGTACACGGATATTGCTTTGGTCTGTACGCACTAATTGACTTATCAGTTCACTGCCGAGGCTGTTGTCGtt                                                       >  1:8948275/5‑93 (MQ=255)

                                  ggCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGtt                                                       <  2:8623607/93‑1 (MQ=255)

                                   gCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTg                                                      <  1:9791528/93‑1 (MQ=255)

                                   gCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTg                                                      <  1:6064927/93‑1 (MQ=255)

                                   gCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTg                                                      <  2:5600656/93‑1 (MQ=255)

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

                                     aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTggg                                                    <  1:6646297/93‑1 (MQ=255)

                                     aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGATGTCGTTggg                                                    >  2:1120601/1‑93 (MQ=255)

                                      ttAAGCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGa                                                   >  2:8033980/1‑93 (MQ=255)

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

                                        aaaCGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATc                                                 >  2:890547/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                               tACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGaa                                          <  1:5650463/93‑1 (MQ=255)

                                               tACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGaa                                          <  2:9005903/93‑1 (MQ=255)

                                                 cATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGGTGTCGTTGGGGTCCGATGAACg                                        >  2:7508952/1‑93 (MQ=255)

                                                  aTGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGa                                       <  1:4969611/93‑1 (MQ=255)

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

                                                      ttAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCa                                   <  2:4896419/93‑1 (MQ=255)

                                                      ttAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCa                                   >  2:7856698/1‑93 (MQ=255)

                                                      ttAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCa                                   <  2:6945116/93‑1 (MQ=255)

                                                      ttAGCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCa                                   <  1:9468215/93‑1 (MQ=255)

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

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

                                                         gCACACGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACg                                <  1:8701861/93‑1 (MQ=255)

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

                                                           acacGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCt                              >  1:6970152/1‑93 (MQ=255)

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

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

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

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

                                                             acGGATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTc                            >  2:9570587/1‑93 (MQ=255)

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

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

                                                               ggATATTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAg                          >  1:6713198/1‑93 (MQ=255)

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

                                                                 atatTTGTATGGTCCGTTCTCACCAACTGACGTAGCAGGTCCCCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCg                        <  1:1123498/93‑1 (MQ=255)

                                                                 atatTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCg                        <  2:3012362/93‑1 (MQ=255)

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

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

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

                                                                   atTGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTa                      >  2:9165415/1‑93 (MQ=255)

                                                                    ttGCTATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTAt                     <  2:1544622/93‑1 (MQ=255)

                                                                        tATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCa                 <  2:728156/93‑1 (MQ=255)

                                                                        tATGGTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCa                 <  1:4844338/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                           ggTCCGTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTg              <  1:2115940/93‑1 (MQ=255)

                                                                             cccGTTCGCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTGGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgca            <  2:7074058/92‑1 (MQ=255)

                                                                                gtcctcACCACCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgcagca         <  2:10037587/90‑1 (MQ=255)

                                                                                gTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgcagca         <  1:380470/93‑1 (MQ=255)

                                                                                gTTCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTgcagca         <  2:8297994/93‑1 (MQ=255)

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

                                                                                 ttCTCACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCAt        <  1:8324214/93‑1 (MQ=255)

                                                                                  tctcACCAACTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATc       <  2:3709494/93‑1 (MQ=255)

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

CGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGTTGTCGTTGGGATCAGATGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGT  >  NC_000913/2289488‑2289667

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