BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,717	+C	coding (86/903 nt)	rarD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,001,717	1	.	C	99.6%	816.2 / ‑4.2	257	coding (86/903 nt)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (105/151); total (106/151)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 296 total aligned reads displayed.

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑G‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGTAGTAAATCAACTTGAAGTAC‑g                                                                                         >  2:2218888/1‑93 (MQ=255)

 gCTCATCAGCACCACCATAAAGAAAATCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gc                                                                                        <  1:189591/93‑2 (MQ=255)

 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gc                                                                                        <  1:380729/93‑2 (MQ=255)

 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAAATTGAAGTAC‑gc                                                                                        >  1:1080541/1‑92 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gctg                                                                                      >  2:244278/1‑90 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAATAC‑gctg                                                                                      >  1:3574098/1‑90 (MQ=255)

    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGTAGTAAATCAAATTGAAGTAC‑gctgg                                                                                     >  1:5490662/1‑89 (MQ=255)

    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gctgg                                                                                     <  2:1100839/93‑5 (MQ=255)

    catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑gctgg                                                                                     <  1:5901963/93‑5 (MQ=255)

        aGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCt                                                                                 <  2:2006704/93‑1 (MQ=255)

         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCta                                                                                >  2:4215997/1‑93 (MQ=255)

         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCta                                                                                >  1:4698394/1‑93 (MQ=255)

         gCACAACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCta                                                                                <  1:4262242/93‑1 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtata                                                                              <  2:5213734/93‑1 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtata                                                                              >  2:5291860/1‑93 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCtata                                                                              >  2:4334076/1‑93 (MQ=255)

               ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCGCGTAGTAAATCAACTTGAATTAA‑GCTGGCGCTATAcccc                                                                          >  2:2720157/1‑93 (MQ=255)

               ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATAcccc                                                                          >  2:1825131/1‑93 (MQ=255)

               ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATAaccc                                                                          >  2:5405751/1‑93 (MQ=255)

                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCa                                                                         >  2:554035/1‑93 (MQ=255)

                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaa                                                                        >  1:4594734/1‑93 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaa                                                                       >  2:6470683/1‑93 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaa                                                                       <  1:5179347/93‑1 (MQ=255)

                   aaaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaat                                                                      <  2:3688050/93‑1 (MQ=255)

                      gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaa                                                                   <  1:2997419/93‑1 (MQ=255)

                       aaagaCGACCAGATCACGCGATGCGTAAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaaa                                                                  <  1:4967258/93‑1 (MQ=255)

                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTAC‑GCTGGCGCTATACCCCaaataaaaa                                                                 >  1:91054/1‑92 (MQ=255)

                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaaat                                                                 <  2:5813525/93‑1 (MQ=255)

                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaaat                                                                 <  2:5117873/93‑1 (MQ=255)

                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCaaataaaata                                                                >  2:5084/1‑93 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAaataaaataa                                                               <  1:6398902/93‑1 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  1:4727354/93‑1 (MQ=255)

                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCg                                                            <  1:3331271/93‑1 (MQ=255)

                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCACCTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCa                                                         <  2:6690069/93‑1 (MQ=255)

                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCa                                                         <  2:6182329/93‑1 (MQ=255)

                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCa                                                         <  1:5754951/93‑1 (MQ=255)

                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCa                                                         >  2:2826153/1‑93 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCCGcca                                                        >  1:5276487/1‑91 (MQ=255)

                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                                       >  1:4864782/1‑93 (MQ=255)

                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                                       <  1:1831254/93‑1 (MQ=255)

                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAag                                                       <  2:5001269/93‑1 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGAAT‑GCTTGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      >  2:6340733/1‑93 (MQ=255)

                                     cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGG‑TGGCGCTATACCCCCAATAAAATAAAAGGCAagagg                                                    >  2:4473898/1‑92 (MQ=255)

                                      aCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCg                                                   <  1:4935213/93‑1 (MQ=255)

                                       cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                  <  2:3285029/93‑1 (MQ=255)

                                       cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGa                                                  <  1:4438471/93‑1 (MQ=255)

                                         cgATGCGTCAGGGTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGGAAGAGCGAga                                                >  1:3126256/1‑92 (MQ=255)

                                         cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAga                                                >  1:4479489/1‑92 (MQ=255)

                                         cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGAGGCAAGAGCGAGt                                                >  1:2813501/1‑93 (MQ=255)

                                          gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGGAAGAGCGAgga                                               >  1:3267379/1‑91 (MQ=255)

                                          gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                                               <  2:4417862/93‑1 (MQ=255)

                                          gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                                               >  2:945202/1‑93 (MQ=255)

                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGCAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACg                                          >  2:4821269/1‑93 (MQ=255)

                                               gTCAGGATTTCATCGGCGGGCACGTAGGAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGAGGCAAGAACGAGTAATaac                                          >  2:3365563/1‑91 (MQ=255)

                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                         >  1:1539872/1‑93 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  2:1270586/93‑1 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  1:1208402/93‑1 (MQ=255)

                                                   gtaTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATTCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                      <  1:2574917/91‑1 (MQ=255)

                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGcccc                                      >  1:5662217/1‑92 (MQ=255)

                                                        tCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgg                                 >  1:2065935/1‑92 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACAC‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              >  2:6547568/1‑93 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              <  2:4560389/93‑1 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              <  1:1687386/93‑1 (MQ=255)

                                                            cggcggGCACGTAGTAAATCATCTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                             <  2:1729133/93‑1 (MQ=255)

                                                             ggcggGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                            <  2:750234/93‑1 (MQ=255)

                                                             ggcggGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                            <  2:410735/93‑1 (MQ=255)

                                                               cggGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                          <  1:6183514/93‑1 (MQ=255)

                                                               cggGCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                          >  1:6024582/1‑93 (MQ=255)

                                                               cggGCACGTAGTAAATCAACGTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                          <  2:1547322/93‑1 (MQ=255)

                                                                gggCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt                         <  1:5232078/93‑1 (MQ=255)

                                                                 ggCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                        <  2:2516809/93‑1 (MQ=255)

                                                                 ggCACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                        <  2:4770920/93‑1 (MQ=255)

                                                                   cACGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc                      <  1:1611789/93‑1 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     <  1:6112778/93‑1 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  2:6086306/1‑93 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  2:2979072/1‑93 (MQ=255)

                                                                     cGTAGTAAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat                    >  1:6266958/1‑93 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  <  1:3887974/93‑1 (MQ=255)

                                                                        agtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                 <  1:806948/93‑1 (MQ=255)

                                                                         gtaAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                >  1:4027879/1‑93 (MQ=255)

                                                                          taAATCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccatc               <  2:5955336/93‑1 (MQ=255)

                                                                           aaaTCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg              <  1:2869360/93‑1 (MQ=255)

                                                                            aaTCAACTTGAAGTAC‑GCTGGCGCTATCCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             <  2:136783/93‑1 (MQ=255)

                                                                            aaTCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             >  1:6589935/1‑93 (MQ=255)

                                                                            aaTCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             >  1:1447457/1‑93 (MQ=255)

                                                                             aTCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCggg            >  2:1025503/1‑93 (MQ=255)

                                                                             aTCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCggg            >  2:3919319/1‑93 (MQ=255)

                                                                              tCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa           <  1:3465172/93‑1 (MQ=255)

                                                                              tCAACTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa           <  2:2613192/93‑1 (MQ=255)

                                                                              tCAACCTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa           <  1:1026984/93‑1 (MQ=255)

                                                                                aaCTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         <  1:6647845/93‑1 (MQ=255)

                                                                                 aCTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa        >  2:5723780/1‑93 (MQ=255)

                                                                                  cTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       >  2:3232915/1‑93 (MQ=255)

                                                                                  cTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       >  2:5597742/1‑93 (MQ=255)

                                                                                  cTTGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       <  1:2203799/93‑1 (MQ=255)

                                                                                   ttGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACt      >  2:3968265/1‑93 (MQ=255)

                                                                                   ttGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACt      <  1:375549/93‑1 (MQ=255)

                                                                                    tGAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTc     >  1:2261157/1‑93 (MQ=255)

                                                                                     gAAGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCa    >  2:3083657/1‑93 (MQ=255)

                                                                                      aaGTAC‑GCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAt   <  2:4857708/93‑1 (MQ=255)

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAC‑G‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

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